#include <DicomPhysicsList.hh>

Inheritance diagram for DicomPhysicsList:

Public Member Functions
	DicomPhysicsList ()

	~DicomPhysicsList ()

Public Member Functions inherited from G4VUserPhysicsList
	G4VUserPhysicsList ()

virtual	~G4VUserPhysicsList ()

	G4VUserPhysicsList (const G4VUserPhysicsList &)

G4VUserPhysicsList &	operator= (const G4VUserPhysicsList &)

void	Construct ()

void	UseCoupledTransportation (G4bool vl=true)

void	SetDefaultCutValue (G4double newCutValue)

G4double	GetDefaultCutValue () const

void	BuildPhysicsTable ()

void	PreparePhysicsTable (G4ParticleDefinition *)

void	BuildPhysicsTable (G4ParticleDefinition *)

G4bool	StorePhysicsTable (const G4String &directory=".")

G4bool	IsPhysicsTableRetrieved () const

G4bool	IsStoredInAscii () const

const G4String &	GetPhysicsTableDirectory () const

void	SetPhysicsTableRetrieved (const G4String &directory="")

void	SetStoredInAscii ()

void	ResetPhysicsTableRetrieved ()

void	ResetStoredInAscii ()

void	DumpList () const

void	DumpCutValuesTable (G4int flag=1)

void	DumpCutValuesTableIfRequested ()

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

void	SetCutsWithDefault ()

void	SetCutValue (G4double aCut, const G4String &pname)

G4double	GetCutValue (const G4String &pname) const

void	SetCutValue (G4double aCut, const G4String &pname, const G4String &rname)

void	SetParticleCuts (G4double cut, G4ParticleDefinition particle, G4Region region=0)

void	SetParticleCuts (G4double cut, const G4String &particleName, G4Region *region=0)

void	SetCutsForRegion (G4double aCut, const G4String &rname)

void	ResetCuts ()
	obsolete methods More...

void	SetApplyCuts (G4bool value, const G4String &name)

G4bool	GetApplyCuts (const G4String &name) const

void	RemoveProcessManager ()

void	AddProcessManager (G4ParticleDefinition newParticle, G4ProcessManager newManager=0)

void	CheckParticleList ()

void	DisableCheckParticleList ()

G4int	GetInstanceID () const

void	InitializeWorker ()

Protected Member Functions
virtual void	ConstructParticle ()

void	ConstructBosons ()

void	ConstructLeptons ()

void	ConstructBaryons ()

virtual void	ConstructProcess ()

void	ConstructEM ()

void	ConstructHad ()

void	ConstructGeneral ()

virtual void	SetCuts ()

Protected Member Functions inherited from G4VUserPhysicsList
void	AddTransportation ()

G4bool	RegisterProcess (G4VProcess process, G4ParticleDefinition particle)

void	BuildIntegralPhysicsTable (G4VProcess , G4ParticleDefinition )

virtual void	RetrievePhysicsTable (G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)

void	InitializeProcessManager ()

Additional Inherited Members
Static Public Member Functions inherited from G4VUserPhysicsList
static const G4VUPLManager &	GetSubInstanceManager ()

Protected Attributes inherited from G4VUserPhysicsList
G4ParticleTable *	theParticleTable

G4int	verboseLevel

G4double	defaultCutValue

G4bool	isSetDefaultCutValue

G4ProductionCutsTable *	fCutsTable

G4bool	fRetrievePhysicsTable

G4bool	fStoredInAscii

G4bool	fIsCheckedForRetrievePhysicsTable

G4bool	fIsRestoredCutValues

G4String	directoryPhysicsTable

G4bool	fDisableCheckParticleList

G4int	g4vuplInstanceID

Static Protected Attributes inherited from G4VUserPhysicsList
static G4RUN_DLL G4VUPLManager	subInstanceManager

Detailed Description

Definition at line 39 of file DicomPhysicsList.hh.

Constructor & Destructor Documentation

DicomPhysicsList::DicomPhysicsList ( )

Definition at line 44 of file DicomPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, python.hepunit::micrometer, and G4VUserPhysicsList::SetVerboseLevel().

                                   :  G4VUserPhysicsList()
 {
   defaultCutValue = 0.01*micrometer;
   fCutForGamma     = defaultCutValue;
   fCutForElectron  = defaultCutValue;
   fCutForPositron  = defaultCutValue;
   SetVerboseLevel(1);
 }

DicomPhysicsList::~DicomPhysicsList ( )

Definition at line 54 of file DicomPhysicsList.cc.

55 {}

Member Function Documentation

void DicomPhysicsList::ConstructBaryons ( )

protected

Definition at line 84 of file DicomPhysicsList.cc.

References G4BaryonConstructor::ConstructParticle(), G4MesonConstructor::ConstructParticle(), and G4IonConstructor::ConstructParticle().

Referenced by ConstructParticle().

 {
   //  baryons
   G4BaryonConstructor bConstructor;
   bConstructor.ConstructParticle();
 
   G4IonConstructor iConstructor;
   iConstructor.ConstructParticle();
 
   G4MesonConstructor mConstructor;
   mConstructor.ConstructParticle();
 }

void DicomPhysicsList::ConstructBosons ( )

protected

Definition at line 66 of file DicomPhysicsList.cc.

References G4Gamma::GammaDefinition(), and G4OpticalPhoton::OpticalPhotonDefinition().

Referenced by ConstructParticle().

 { 
   // gamma
   G4Gamma::GammaDefinition();
 
   // optical photon
   G4OpticalPhoton::OpticalPhotonDefinition();
 }

void DicomPhysicsList::ConstructEM ( )

protected

Definition at line 162 of file DicomPhysicsList.cc.

References G4ProcessManager::AddDiscreteProcess(), G4VEmProcess::AddEmModel(), G4VEnergyLossProcess::AddEmModel(), G4ProcessManager::AddProcess(), fUseDistanceToBoundary, G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetProcessManager(), G4VEnergyLossProcess::SetEmModel(), G4VEnergyLossProcess::SetStepFunction(), G4VMultipleScattering::SetStepLimitType(), and theParticleIterator.

Referenced by ConstructProcess().

 {
   theParticleIterator->reset();
 
   while( (*theParticleIterator)() ){
 
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* pmanager = particle->GetProcessManager();
     G4String particleName = particle->GetParticleName();
 
     if (particleName == "gamma") {
 
 
       G4PhotoElectricEffect* thePhotoElectricEffect =
                                                    new G4PhotoElectricEffect();
       G4LivermorePhotoElectricModel* theLivermorePhotoElectricModel
                                          = new G4LivermorePhotoElectricModel();
       thePhotoElectricEffect->AddEmModel(0, theLivermorePhotoElectricModel);
       pmanager->AddDiscreteProcess(thePhotoElectricEffect);
 
       G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
       G4LivermoreComptonModel* theLivermoreComptonModel =
                                                  new G4LivermoreComptonModel();
       theComptonScattering->AddEmModel(0, theLivermoreComptonModel);
       pmanager->AddDiscreteProcess(theComptonScattering);
 
       G4GammaConversion* theGammaConversion = new G4GammaConversion();
       G4LivermoreGammaConversionModel* theLivermoreGammaConversionModel
                                        = new G4LivermoreGammaConversionModel();
       theGammaConversion->AddEmModel(0, theLivermoreGammaConversionModel);
       pmanager->AddDiscreteProcess(theGammaConversion);
 
       G4RayleighScattering* theRayleigh = new G4RayleighScattering();
       G4LivermoreRayleighModel* theRayleighModel
                                               = new G4LivermoreRayleighModel();
       theRayleigh->AddEmModel(0, theRayleighModel);
       pmanager->AddDiscreteProcess(theRayleigh);
       
       pmanager->AddProcess(new G4StepLimiter(), -1, -1, 5);
       
     } else if (particleName == "e-") {
      
       G4eMultipleScattering* msc = new G4eMultipleScattering();
       msc->SetStepLimitType(fUseDistanceToBoundary);
       pmanager->AddProcess(msc,                   -1, 1, 1);
       
       // Ionisation
       G4eIonisation* eIoni = new G4eIonisation();
       eIoni->AddEmModel(0, new G4LivermoreIonisationModel(),
                            new G4UniversalFluctuation() );
       eIoni->SetStepFunction(0.2, 100*um); //     
       pmanager->AddProcess(eIoni,                 -1, 2, 2);
       
       // Bremsstrahlung
       G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
       eBrem->AddEmModel(0, new G4LivermoreBremsstrahlungModel());
       pmanager->AddProcess(eBrem,                   -1,-3, 3);
 
       pmanager->AddProcess(new G4StepLimiter(), -1, -1, 4);
       
     } else if (particleName == "e+") {
 
       // Identical to G4EmStandardPhysics_option3
     
       G4eMultipleScattering* msc = new G4eMultipleScattering();
       msc->SetStepLimitType(fUseDistanceToBoundary);
       pmanager->AddProcess(msc,                   -1, 1, 1);
 
       G4eIonisation* eIoni = new G4eIonisation();
       eIoni->SetStepFunction(0.2, 100*um);      
       pmanager->AddProcess(eIoni,                 -1, 2, 2);
       
       pmanager->AddProcess(new G4eBremsstrahlung, -1,-3, 3);
       
       pmanager->AddProcess(new G4eplusAnnihilation,0,-1, 4);
       
       pmanager->AddProcess(new G4StepLimiter(), -1, -1, 5);
 
     } else if (particleName == "GenericIon") {
 
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
 
       G4ionIonisation* ionIoni = new G4ionIonisation();
       ionIoni->SetEmModel(new G4IonParametrisedLossModel());
       ionIoni->SetStepFunction(0.1, 20*um);
       pmanager->AddProcess(ionIoni,                   -1, 2, 2);
 
       pmanager->AddProcess(new G4StepLimiter(), -1, -1, 3);
 
     } else if (particleName == "alpha" ||
                particleName == "He3" ) {
 
       // Identical to G4EmStandardPhysics_option3
       
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
 
       G4ionIonisation* ionIoni = new G4ionIonisation();
       ionIoni->SetStepFunction(0.1, 20*um);
       pmanager->AddProcess(ionIoni,                   -1, 2, 2);
 
       pmanager->AddProcess(new G4StepLimiter(), -1, -1, 3);
     }
 
    //
 
   }
 }

void DicomPhysicsList::ConstructGeneral ( )

protected

Definition at line 329 of file DicomPhysicsList.cc.

Referenced by ConstructProcess().

330 { }

void DicomPhysicsList::ConstructHad ( )

protected

Definition at line 279 of file DicomPhysicsList.cc.

References G4HadronicProcess::AddDataSet(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetProcessManager(), python.hepunit::GeV, python.hepunit::MeV, G4HadronicProcess::RegisterMe(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), and theParticleIterator.

Referenced by ConstructProcess().

 {
 
   G4HadronElasticProcess * theElasticProcess = new G4HadronElasticProcess;
   theElasticProcess->RegisterMe( new G4HadronElastic() );
 
   theParticleIterator->reset();
   while( (*theParticleIterator)() )
   {
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* pManager = particle->GetProcessManager();
 
     if (particle->GetParticleName() == "alpha") { 
 
       // INELASTIC SCATTERING
       // Binary Cascade
       G4BinaryLightIonReaction* theBC = new G4BinaryLightIonReaction();
       // DHW - change lower limit from 80 MeV to zero because
       //       G4LEDeuteronInelastic is deprecated
       theBC->SetMinEnergy(0.*MeV);
       theBC->SetMaxEnergy(40.*GeV);
   
       // TRIPATHI CROSS SECTION
       // Implementation of formulas in analogy to NASA technical paper 3621 by 
       // Tripathi, et al. Cross-sections for ion ion scattering
       G4TripathiCrossSection* TripathiCrossSection = new G4TripathiCrossSection;
   
       // IONS SHEN CROSS SECTION
       // Implementation of formulas 
       // Shen et al. Nuc. Phys. A 491 130 (1989) 
       // Total Reaction Cross Section for Heavy-Ion Collisions
       G4IonsShenCrossSection* aShen = new G4IonsShenCrossSection;
   
       G4AlphaInelasticProcess* theIPalpha = new G4AlphaInelasticProcess;                  
       theIPalpha->AddDataSet(TripathiCrossSection);
       theIPalpha->AddDataSet(aShen);
 
       // Register the Binary Cascade Model
       theIPalpha->RegisterMe(theBC);
           
       // Activate the alpha inelastic scattering using the binary cascade model
       pManager -> AddDiscreteProcess(theIPalpha);
           
       // Activate the Hadron Elastic Process
       pManager -> AddDiscreteProcess(theElasticProcess);           
     }
   }
 }

void DicomPhysicsList::ConstructLeptons ( )

protected

Definition at line 76 of file DicomPhysicsList.cc.

References G4Electron::ElectronDefinition(), and G4Positron::PositronDefinition().

Referenced by ConstructParticle().

 {
   // leptons
   G4Electron::ElectronDefinition();
   G4Positron::PositronDefinition();
 }

void DicomPhysicsList::ConstructParticle ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 58 of file DicomPhysicsList.cc.

References ConstructBaryons(), ConstructBosons(), and ConstructLeptons().

 {
   ConstructBosons();
   ConstructLeptons();
   ConstructBaryons();
 }

void DicomPhysicsList::ConstructProcess ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 98 of file DicomPhysicsList.cc.

References G4VUserPhysicsList::AddTransportation(), ConstructEM(), ConstructGeneral(), and ConstructHad().

 {
   AddTransportation();
   ConstructEM();
   ConstructHad();
   ConstructGeneral();
 }

void DicomPhysicsList::SetCuts ( )

protectedvirtual

Reimplemented from G4VUserPhysicsList.

Definition at line 333 of file DicomPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4VUserPhysicsList::DumpCutValuesTable(), G4BestUnit, G4cout, G4endl, G4VUserPhysicsList::SetCutValue(), and G4VUserPhysicsList::verboseLevel.

 {
   if (verboseLevel >0){
     G4cout << "DicomPhysicsList::SetCuts:";
     G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
   }  
   
   // set cut values for gamma at first and for e- second and next for e+,
   // because some processes for e+/e- need cut values for gamma 
   SetCutValue(fCutForGamma, "gamma");
   SetCutValue(fCutForElectron, "e-");
   SetCutValue(fCutForPositron, "e+");
   
   if (verboseLevel>0) DumpCutValuesTable();
 
 }

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation