#include <ExGflashPhysicsList.hh>

Inheritance diagram for ExGflashPhysicsList:

Public Member Functions
	ExGflashPhysicsList ()

virtual	~ExGflashPhysicsList ()

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 ()

virtual void	ConstructProcess ()

virtual void	SetCuts ()

virtual void	ConstructBosons ()

virtual void	ConstructLeptons ()

virtual void	ConstructMesons ()

virtual void	ConstructBaryons ()

virtual void	ConstructIons ()

void	AddParameterisation ()

virtual void	ConstructGeneral ()

virtual void	ConstructEM ()

virtual void	AddTransportation ()

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 38 of file ExGflashPhysicsList.hh.

Constructor & Destructor Documentation

ExGflashPhysicsList::ExGflashPhysicsList ( )

Definition at line 51 of file ExGflashPhysicsList.cc.

References G4VUserPhysicsList::SetVerboseLevel().

                                         :  G4VUserPhysicsList()
 {
   SetVerboseLevel(0);
 }

ExGflashPhysicsList::~ExGflashPhysicsList ( )

virtual

Definition at line 56 of file ExGflashPhysicsList.cc.

57 {

58 }

Member Function Documentation

void ExGflashPhysicsList::AddParameterisation ( )

protected

Definition at line 303 of file ExGflashPhysicsList.cc.

References G4ProcessManager::AddDiscreteProcess(), G4cout, G4endl, G4ParticleDefinition::GetProcessManager(), and theParticleIterator.

Referenced by ConstructProcess().

 {
   G4FastSimulationManagerProcess* 
     theFastSimulationManagerProcess = 
       new G4FastSimulationManagerProcess();
   G4cout << "FastSimulationManagerProcess" <<G4endl;
   theParticleIterator->reset();
   //std::cout<<"---"<<std::endl;
   while( (*theParticleIterator)() ){
     //std::cout<<"+++"<<std::endl;
   
     G4ParticleDefinition* particle = theParticleIterator->value();
     // std::cout<<"--- particle "<<particle->GetParticleName()<<std::endl;
     G4ProcessManager* pmanager = particle->GetProcessManager();
     // The fast simulation process becomes a discrete process only since 9.0:
     pmanager->AddDiscreteProcess(theFastSimulationManagerProcess);
   }
 }

void ExGflashPhysicsList::AddTransportation ( )

protectedvirtual

Definition at line 133 of file ExGflashPhysicsList.cc.

References G4VUserPhysicsList::AddTransportation().

Referenced by ConstructProcess().

 {
   G4VUserPhysicsList::AddTransportation();
 }

void ExGflashPhysicsList::ConstructBaryons ( )

protectedvirtual

Definition at line 106 of file ExGflashPhysicsList.cc.

References G4BaryonConstructor::ConstructParticle().

Referenced by ConstructParticle().

 {
   //  Construct all barions
   G4BaryonConstructor  pConstructor;
   pConstructor.ConstructParticle();  
 }

void ExGflashPhysicsList::ConstructBosons ( )

protectedvirtual

Definition at line 76 of file ExGflashPhysicsList.cc.

References G4ChargedGeantino::ChargedGeantinoDefinition(), G4Gamma::GammaDefinition(), G4Geantino::GeantinoDefinition(), and G4OpticalPhoton::OpticalPhotonDefinition().

Referenced by ConstructParticle().

 {
   // pseudo-particles
   G4Geantino::GeantinoDefinition();
   G4ChargedGeantino::ChargedGeantinoDefinition();
 
   // gamma
   G4Gamma::GammaDefinition();
 
   // optical photon
   G4OpticalPhoton::OpticalPhotonDefinition();
 }

void ExGflashPhysicsList::ConstructEM ( )

protectedvirtual

Definition at line 157 of file ExGflashPhysicsList.cc.

References G4ProcessManager::AddDiscreteProcess(), G4ProcessManager::AddProcess(), G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetProcessManager(), idxAlongStep, idxAtRest, idxPostStep, G4ParticleDefinition::IsShortLived(), SetCuts(), G4ProcessManager::SetProcessOrdering(), G4ProcessManager::SetProcessOrderingToFirst(), and theParticleIterator.

Referenced by ConstructProcess().

 {
 
    G4cout<<"Physics List constructor"<<G4endl;
    SetCuts();
   theParticleIterator->reset();
   while( (*theParticleIterator)() ){
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* pmanager = particle->GetProcessManager();
     G4String particleName = particle->GetParticleName();
      
     if (particleName == "gamma") {
     // gamma
       // Construct processes for gamma
       G4VProcess* theGammaConversion = new G4GammaConversion();
       G4VProcess* theComptonScattering = new G4ComptonScattering();
       G4VProcess* thePhotoElectricEffect = new G4PhotoElectricEffect();     
       // G4VProcess* thegammacut = new G4UserLimits();
       //      thegammacut->SetUserMinEkine(1.0*MeV);
 
       pmanager->AddDiscreteProcess(theGammaConversion);
       pmanager->AddDiscreteProcess(theComptonScattering);      
       pmanager->AddDiscreteProcess(thePhotoElectricEffect);
       //  pmanager->AddProcess(thegammacut);
       //   G4cout <<"theGammaConversion" << theGammaConversion <<endl;
       //G4cout <<"theComptonScattering" << theComptonScattering <<endl;
       //G4cout <<"thePhotoElectricEffect" << thePhotoElectricEffect <<endl;
 
     } else if (particleName == "e-") {
     //electron
       // Construct processes for electron
       G4VProcess* theeminusMultipleScattering = new G4eMultipleScattering();
       G4VProcess* theeminusIonisation = new G4eIonisation();
       G4VProcess* theeminusBremsstrahlung = new G4eBremsstrahlung();
       //      G4VProcess* theeminuscut = new G4UserLimits();
       //  theeminuscut->SetUserMinEkine(1.0*MeV);
       // add processes
       pmanager->AddProcess(theeminusMultipleScattering);
       pmanager->AddProcess(theeminusIonisation);
       pmanager->AddProcess(theeminusBremsstrahlung);      
      
       //  pmanager->AddProcess( theeminuscut);
 
 
   
  // set ordering for AlongStepDoIt
       pmanager->SetProcessOrdering(theeminusMultipleScattering, idxAlongStep,  1);
       pmanager->SetProcessOrdering(theeminusIonisation, idxAlongStep,  2);
       // set ordering for PostStepDoIt
       pmanager->SetProcessOrdering(theeminusMultipleScattering, idxPostStep, 1);
       pmanager->SetProcessOrdering(theeminusIonisation, idxPostStep, 2);
       pmanager->SetProcessOrdering(theeminusBremsstrahlung, idxPostStep, 3);
 
       //G4cout <<"theeminusMultipleScattering" << theeminusMultipleScattering <<endl;
       //G4cout <<"theeminusIonisation" << theeminusIonisation <<endl;
       //G4cout <<"theeminusBremsstrahlung" << theeminusBremsstrahlung <<endl;
 
     } else if (particleName == "e+") {
     //positron
       // Construct processes for positron
       G4VProcess* theeplusMultipleScattering = new G4eMultipleScattering();
       G4VProcess* theeplusIonisation = new G4eIonisation();
       G4VProcess* theeplusBremsstrahlung = new G4eBremsstrahlung();
       G4VProcess* theeplusAnnihilation = new G4eplusAnnihilation();
       // add processes
       pmanager->AddProcess(theeplusMultipleScattering);
       pmanager->AddProcess(theeplusIonisation);
       pmanager->AddProcess(theeplusBremsstrahlung);
       pmanager->AddProcess(theeplusAnnihilation);
       // set ordering for AtRestDoIt
       pmanager->SetProcessOrderingToFirst(theeplusAnnihilation, idxAtRest);
       // set ordering for AlongStepDoIt
       pmanager->SetProcessOrdering(theeplusMultipleScattering, idxAlongStep,  1);
       pmanager->SetProcessOrdering(theeplusIonisation, idxAlongStep,  2);
       // set ordering for PostStepDoIt
       pmanager->SetProcessOrdering(theeplusMultipleScattering, idxPostStep, 1);
       pmanager->SetProcessOrdering(theeplusIonisation, idxPostStep, 2);
       pmanager->SetProcessOrdering(theeplusBremsstrahlung, idxPostStep, 3);
       pmanager->SetProcessOrdering(theeplusAnnihilation, idxPostStep, 4);
   
       //G4cout <<"theeplusMultipleScattering" << theeplusMultipleScattering <<endl;
       //G4cout <<"theeplusIonisation" << theeplusIonisation <<endl;
       //G4cout <<"theeplusBremsstrahlung" << theeplusBremsstrahlung <<endl;
 
     } else if( particleName == "mu+" || 
                particleName == "mu-"    ) {
     //muon  
      // Construct processes for muon+
      G4VProcess* aMultipleScattering = new G4MuMultipleScattering();
      G4VProcess* aBremsstrahlung = new G4MuBremsstrahlung();
      G4VProcess* aPairProduction = new G4MuPairProduction();
      G4VProcess* anIonisation = new G4MuIonisation();
       // add processes
      pmanager->AddProcess(anIonisation);
      pmanager->AddProcess(aMultipleScattering);
      pmanager->AddProcess(aBremsstrahlung);
      pmanager->AddProcess(aPairProduction);
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering(aMultipleScattering, idxAlongStep,  1);
      pmanager->SetProcessOrdering(anIonisation, idxAlongStep,  2);
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering(aMultipleScattering, idxPostStep, 1);
      pmanager->SetProcessOrdering(anIonisation, idxPostStep, 2);
      pmanager->SetProcessOrdering(aBremsstrahlung, idxPostStep, 3);
      pmanager->SetProcessOrdering(aPairProduction, idxPostStep, 4);
      
    } else if ((!particle->IsShortLived()) &&
         (particle->GetPDGCharge() != 0.0) && 
         (particle->GetParticleName() != "chargedgeantino")) {
      // all others charged particles except geantino
      G4VProcess* aMultipleScattering = new G4hMultipleScattering();
      G4VProcess* anIonisation = new G4hIonisation();
      // add processes
      pmanager->AddProcess(anIonisation);
      pmanager->AddProcess(aMultipleScattering);
      // set ordering for AlongStepDoIt
      pmanager->SetProcessOrdering(aMultipleScattering, idxAlongStep,  1);
      pmanager->SetProcessOrdering(anIonisation, idxAlongStep,  2);
      // set ordering for PostStepDoIt
      pmanager->SetProcessOrdering(aMultipleScattering, idxPostStep, 1);
      pmanager->SetProcessOrdering(anIonisation, idxPostStep, 2);
     }
   }
 }

void ExGflashPhysicsList::ConstructGeneral ( )

protectedvirtual

Definition at line 284 of file ExGflashPhysicsList.cc.

References G4ProcessManager::AddProcess(), G4ParticleDefinition::GetProcessManager(), idxAtRest, idxPostStep, G4Decay::IsApplicable(), G4ProcessManager::SetProcessOrdering(), and theParticleIterator.

Referenced by ConstructProcess().

 {
   // Add Decay Process
   G4Decay* theDecayProcess = new G4Decay();
   //G4cout << "decay" <<theDecayProcess<<endl;
   theParticleIterator->reset();
   while( (*theParticleIterator)() ){
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* pmanager = particle->GetProcessManager();
     if (theDecayProcess->IsApplicable(*particle)) { 
       pmanager ->AddProcess(theDecayProcess);
       // set ordering for PostStepDoIt and AtRestDoIt
       pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
       pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
     }
   }  
 }

void ExGflashPhysicsList::ConstructIons ( )

protectedvirtual

Definition at line 114 of file ExGflashPhysicsList.cc.

References G4IonConstructor::ConstructParticle().

Referenced by ConstructParticle().

 {
   //  Construct light ions
   G4IonConstructor pConstructor;
   pConstructor.ConstructParticle();  
 }

void ExGflashPhysicsList::ConstructLeptons ( )

protectedvirtual

Definition at line 90 of file ExGflashPhysicsList.cc.

References G4LeptonConstructor::ConstructParticle().

Referenced by ConstructParticle().

 {
   // Construct all leptons
   G4LeptonConstructor pConstructor;
   pConstructor.ConstructParticle();
 }

void ExGflashPhysicsList::ConstructMesons ( )

protectedvirtual

Definition at line 98 of file ExGflashPhysicsList.cc.

References G4MesonConstructor::ConstructParticle().

Referenced by ConstructParticle().

 {
   //  Construct all mesons
   G4MesonConstructor pConstructor;
   pConstructor.ConstructParticle();
 }

void ExGflashPhysicsList::ConstructParticle ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 60 of file ExGflashPhysicsList.cc.

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

 {
   // In this method, static member functions should be called
   // for all particles which you want to use.
   // This ensures that objects of these particle types will be
   // created in the program. 
 
   std::cout<<"start construct particle"<<std::endl;
   ConstructBosons();
   ConstructLeptons();
   ConstructMesons();
   ConstructBaryons();
   ConstructIons();
   std::cout<<"end construct particle"<<std::endl;
 }

void ExGflashPhysicsList::ConstructProcess ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 121 of file ExGflashPhysicsList.cc.

References AddParameterisation(), AddTransportation(), ConstructEM(), and ConstructGeneral().

 {
   AddTransportation();
   AddParameterisation();
   std::cout<<"AddParameterisation"<<std::endl;
  
   ConstructEM();
   std::cout<<"ConstructEM"<<std::endl;
   ConstructGeneral();
 }

void ExGflashPhysicsList::SetCuts ( )

protectedvirtual

Reimplemented from G4VUserPhysicsList.

Definition at line 322 of file ExGflashPhysicsList.cc.

References G4VUserPhysicsList::DumpCutValuesTable(), G4cout, G4RegionStore::GetInstance(), G4ProductionCutsTable::GetProductionCutsTable(), G4RegionStore::GetRegion(), G4VUserPhysicsList::SetCutsWithDefault(), G4Region::SetProductionCuts(), and G4VUserPhysicsList::verboseLevel.

Referenced by ConstructEM().

 {
   if (verboseLevel >1){
     G4cout << "ExGflashPhysicsList::SetCuts:";
   }  
   //  " G4VUserPhysicsList::SetCutsWithDefault" method sets 
   //   the default cut value for all particle types 
   DumpCutValuesTable();
    SetCutsWithDefault();
    //  SetCutValue(100*mm, "gamma");
 //   SetCutValue(0*mm, "e-");
 //   SetCutValue(0*mm, "e+");
 
 //   SetCutValue(62*mm, "gamma");
 //   SetCutValue(0.73*mm, "e-");
 //   SetCutValue(0.78*mm, "e+");
 
 
 
  
  
   DumpCutValuesTable();
 // set cuts for region crystals with default Cuts
   G4Region* region = G4RegionStore::GetInstance()->GetRegion("crystals");
   region->SetProductionCuts(
           G4ProductionCutsTable::GetProductionCutsTable()->GetDefaultProductionCuts());
 }

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