#include <eRositaPhysicsList.hh>

Inheritance diagram for eRositaPhysicsList:

Public Member Functions
	eRositaPhysicsList ()

	~eRositaPhysicsList ()

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

void	ConstructProcess ()

void	SetCuts ()

void	ConstructBosons ()

void	ConstructLeptons ()

void	ConstructMesons ()

void	ConstructBaryons ()

void	ConstructGeneral ()

void	ConstructEM ()

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 36 of file eRositaPhysicsList.hh.

Constructor & Destructor Documentation

eRositaPhysicsList::eRositaPhysicsList ( )

Definition at line 59 of file eRositaPhysicsList.cc.

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

                                       :  G4VUserPhysicsList()
 {
   defaultCutValue = 0.001*mm;
    SetVerboseLevel(1);
 
    std::cout << "==============================================================================="
          << std::endl
          << "Geant4 eRosita example - based on a simplified version of eROSITA simulation"
          << std::endl
          << "Further details can be found in:"
          << std::endl
          << "M.G. Pia et al., 'PIXE Simulation With Geant4', "
          << "IEEE Trans. Nucl. Sci., vol. 56, no. 6, pp. 3614-3649, 2009"
          << std::endl
          << "N. Meidinger et al., 'Development of the focal plane PNCCD camera system for the X-ray space telescope eROSITA', " 
          << std::endl
          <<"NIM A 624, 321-329, 2010"
          << std::endl
          << "==============================================================================="
          << std::endl;
 
    std::cout<< std::endl;
    
    std::cout << "==============================================================================="
          << std::endl
          << " The use of G4LowEnergyIonisation, G4LowEnergyBremsstrahlung, "
          << std::endl
          << "G4LowEnergyPhotoElectric, G4LowEnergyCompton, G4LowEnergyGammaConversion"
          << std::endl
          << "in this example is intentional. These classes will be replaced by other classes"
          << std::endl
          << "appropriate to the problem domain in a forthcoming Geant4 version"
          << std::endl
          << "==============================================================================="
          << std::endl;
 }

eRositaPhysicsList::~eRositaPhysicsList ( )

Definition at line 97 of file eRositaPhysicsList.cc.

98 {}

Member Function Documentation

void eRositaPhysicsList::ConstructBaryons ( )

protected

Definition at line 157 of file eRositaPhysicsList.cc.

References G4AntiProton::AntiProtonDefinition(), and G4Proton::ProtonDefinition().

Referenced by ConstructParticle().

 {
   //  barions
   G4Proton::ProtonDefinition();
   G4AntiProton::AntiProtonDefinition();
 
   //G4Neutron::NeutronDefinition();
   //G4AntiNeutron::AntiNeutronDefinition();
 }

void eRositaPhysicsList::ConstructBosons ( )

protected

Definition at line 110 of file eRositaPhysicsList.cc.

References G4Gamma::GammaDefinition().

Referenced by ConstructParticle().

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

void eRositaPhysicsList::ConstructEM ( )

protected

Definition at line 178 of file eRositaPhysicsList.cc.

References G4LowEnergyPhotoElectric::ActivateAuger(), G4ProcessManager::AddProcess(), python.hepunit::eV, G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetProcessManager(), G4ParticleDefinition::IsShortLived(), python.hepunit::keV, python.hepunit::MeV, G4hImpactIonisation::SetCutForAugerElectrons(), G4LowEnergyPhotoElectric::SetCutForLowEnSecElectrons(), G4LowEnergyPhotoElectric::SetCutForLowEnSecPhotons(), G4hImpactIonisation::SetCutForSecondaryPhotons(), G4hImpactIonisation::SetPixeCrossSectionK(), G4hImpactIonisation::SetPixeCrossSectionL(), G4hImpactIonisation::SetPixeCrossSectionM(), G4hImpactIonisation::SetPixeProjectileMaxEnergy(), G4hImpactIonisation::SetPixeProjectileMinEnergy(), and theParticleIterator.

Referenced by ConstructProcess().

 {
   theParticleIterator->reset();
   while( (*theParticleIterator)() ){
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* processManager = particle->GetProcessManager();
     G4String particleName = particle->GetParticleName();
      
     if (particleName == "gamma") {
 
       // photon   
 
       G4LowEnergyPhotoElectric* photoelectric = new G4LowEnergyPhotoElectric;
       photoelectric->ActivateAuger(true);
       photoelectric->SetCutForLowEnSecPhotons(0.250 * keV);
       photoelectric->SetCutForLowEnSecElectrons(0.250 * keV);
       G4LowEnergyCompton* compton = new G4LowEnergyCompton;
       G4LowEnergyGammaConversion* gammaConversion = new G4LowEnergyGammaConversion;
       G4LowEnergyRayleigh* rayleigh = new G4LowEnergyRayleigh;
 
       processManager -> AddDiscreteProcess(photoelectric);
       processManager -> AddDiscreteProcess(compton);
       processManager -> AddDiscreteProcess(gammaConversion);
       processManager -> AddDiscreteProcess(rayleigh);
     
     } else if (particleName == "e-") {
 
       // electron
 
       G4eMultipleScattering* eMultipleScattering = new G4eMultipleScattering();
       G4LowEnergyIonisation* eIonisation = new G4LowEnergyIonisation();
       G4LowEnergyBremsstrahlung* eBremsstrahlung = new G4LowEnergyBremsstrahlung();
 
       processManager -> AddProcess(eMultipleScattering, -1, 1, 1);
       processManager -> AddProcess(eIonisation, -1, 2, 2);
       processManager -> AddProcess(eBremsstrahlung, -1, -1, 3);   
 
     } else if (particleName == "e+") {
       // positron
       processManager->AddProcess(new G4eMultipleScattering, -1, 1, 1);
       processManager->AddProcess(new G4eIonisation,         -1, 2, 2);
       processManager->AddProcess(new G4eBremsstrahlung,     -1, 3, 3);
       processManager->AddProcess(new G4eplusAnnihilation,    0,-1, 4);
 
       //} else if( particleName == "mu+" || 
       //        particleName == "mu-"    ) {
       //muon  
       //processManager->AddProcess(new G4MuMultipleScattering, -1, 1, 1);
       //processManager->AddProcess(new G4MuIonisation,         -1, 2, 2);
       //processManager->AddProcess(new G4MuBremsstrahlung,     -1, 3, 3);
       //processManager->AddProcess(new G4MuPairProduction,     -1, 4, 4);       
              
     } else if( particleName == "proton" ||
                particleName == "pi-" ||
                particleName == "pi+"    ) {
       //proton  
 
       G4hImpactIonisation* hIonisation = new G4hImpactIonisation();
       hIonisation->SetPixeCrossSectionK("ecpssr");
       hIonisation->SetPixeCrossSectionL("ecpssr");
       hIonisation->SetPixeCrossSectionM("ecpssr");
       hIonisation->SetPixeProjectileMinEnergy(1.* keV);
       hIonisation->SetPixeProjectileMaxEnergy(200. * MeV);
       hIonisation->SetCutForSecondaryPhotons(250. * eV);
       hIonisation->SetCutForAugerElectrons(250. * eV);
 
       G4hMultipleScattering* hMultipleScattering = new G4hMultipleScattering();
 
       processManager -> AddProcess(hMultipleScattering, -1, 1, 1);   
       processManager -> AddProcess(hIonisation, -1, 2, 2);
      
     } else if( particleName == "alpha" || 
            particleName == "He3" ||
            particleName == "pi-" ||
                particleName == "pi+" ||
            particleName == "GenericIon" ) {
 
       // pions, alpha, ions (should never occur in the current example) 
       processManager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       processManager->AddProcess(new G4ionIonisation,       -1, 2, 2);
                       
     } else if ((!particle->IsShortLived()) &&
            (particle->GetPDGCharge() != 0.0) && 
            (particle->GetParticleName() != "chargedgeantino")) {
       //all others charged particles except geantino
       processManager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       processManager->AddProcess(new G4hIonisation,         -1, 2, 2);
     }
   }
 }

void eRositaPhysicsList::ConstructGeneral ( )

protected

Definition at line 271 of file eRositaPhysicsList.cc.

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

Referenced by ConstructProcess().

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

void eRositaPhysicsList::ConstructLeptons ( )

protected

Definition at line 121 of file eRositaPhysicsList.cc.

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

Referenced by ConstructParticle().

 {
   // leptons
   //  e+/-
   G4Electron::ElectronDefinition();
   G4Positron::PositronDefinition();
   // mu+/-
   //G4MuonPlus::MuonPlusDefinition();
   //G4MuonMinus::MuonMinusDefinition();
   // nu_e
   //G4NeutrinoE::NeutrinoEDefinition();
   //G4AntiNeutrinoE::AntiNeutrinoEDefinition();
   // nu_mu
   //G4NeutrinoMu::NeutrinoMuDefinition();
   //G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
 }

void eRositaPhysicsList::ConstructMesons ( )

protected

Definition at line 139 of file eRositaPhysicsList.cc.

Referenced by ConstructParticle().

 {
   //  mesons
   //    light mesons
   //G4PionPlus::PionPlusDefinition();
   //G4PionMinus::PionMinusDefinition();
   //G4PionZero::PionZeroDefinition();
   //G4Eta::EtaDefinition();
   //G4EtaPrime::EtaPrimeDefinition();
   //G4KaonPlus::KaonPlusDefinition();
   //G4KaonMinus::KaonMinusDefinition();
   //G4KaonZero::KaonZeroDefinition();
   //G4AntiKaonZero::AntiKaonZeroDefinition();
   //G4KaonZeroLong::KaonZeroLongDefinition();
   //G4KaonZeroShort::KaonZeroShortDefinition();
 }

void eRositaPhysicsList::ConstructParticle ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 101 of file eRositaPhysicsList.cc.

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

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

void eRositaPhysicsList::ConstructProcess ( void )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 168 of file eRositaPhysicsList.cc.

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

 {
   AddTransportation();
   ConstructEM();
   ConstructGeneral();
   //AddStepMax();
 }

void eRositaPhysicsList::SetCuts ( )

protectedvirtual

Reimplemented from G4VUserPhysicsList.

Definition at line 313 of file eRositaPhysicsList.cc.

References G4VUserPhysicsList::DumpCutValuesTable(), python.hepunit::eV, G4ProductionCutsTable::GetProductionCutsTable(), python.hepunit::GeV, G4VUserPhysicsList::SetCutsWithDefault(), G4ProductionCutsTable::SetEnergyRange(), and G4VUserPhysicsList::verboseLevel.

 {
   //G4VUserPhysicsList::SetCutsWithDefault method sets 
   //the default cut value for all particle types 
   //
   SetCutsWithDefault();
 
   // Set the secondary production cut lower than 990. eV
   // Very important for processes at low energies
  
   G4double lowLimit = 250. * eV;
   G4double highLimit = 100. * GeV;
   G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowLimit, highLimit);
      
   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