eRositaPhysicsList Class Reference

#include <eRositaPhysicsList.hh>

Inheritance diagram for eRositaPhysicsList:

Collaboration 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 ()
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const

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::eRositaPhysicsList

(

)

Definition at line 59 of file eRositaPhysicsList.cc.

                                      :  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;
}

Here is the call graph for this function:

~eRositaPhysicsList()

eRositaPhysicsList::~eRositaPhysicsList

(

)

Definition at line 97 of file eRositaPhysicsList.cc.

{}

Member Function Documentation

ConstructBaryons()

void eRositaPhysicsList::ConstructBaryons ( )

protected

Definition at line 157 of file eRositaPhysicsList.cc.

{
  //  barions
  G4Proton::ProtonDefinition();
  G4AntiProton::AntiProtonDefinition();

  //G4Neutron::NeutronDefinition();
  //G4AntiNeutron::AntiNeutronDefinition();
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructBosons()

void eRositaPhysicsList::ConstructBosons ( )

protected

Definition at line 110 of file eRositaPhysicsList.cc.

{
  // pseudo-particles
  //G4Geantino::GeantinoDefinition();
  //G4ChargedGeantino::ChargedGeantinoDefinition();

  // gamma
  G4Gamma::GammaDefinition();
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructEM()

void eRositaPhysicsList::ConstructEM ( )

protected

Definition at line 178 of file eRositaPhysicsList.cc.

{
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->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);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructGeneral()

void eRositaPhysicsList::ConstructGeneral ( )

protected

Definition at line 272 of file eRositaPhysicsList.cc.

{
  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->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);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructLeptons()

void eRositaPhysicsList::ConstructLeptons ( )

protected

Definition at line 121 of file eRositaPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructMesons()

void eRositaPhysicsList::ConstructMesons ( )

protected

Definition at line 139 of file eRositaPhysicsList.cc.

{
  //  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();
}

Here is the caller graph for this function:

ConstructParticle()

void eRositaPhysicsList::ConstructParticle ( void  )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 101 of file eRositaPhysicsList.cc.

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

Here is the call graph for this function:

ConstructProcess()

void eRositaPhysicsList::ConstructProcess ( void  )

protectedvirtual

Implements G4VUserPhysicsList.

Definition at line 168 of file eRositaPhysicsList.cc.

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

Here is the call graph for this function:

SetCuts()

void eRositaPhysicsList::SetCuts ( )

protectedvirtual

Reimplemented from G4VUserPhysicsList.

Definition at line 316 of file eRositaPhysicsList.cc.

{
  //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();
}

Here is the call graph for this function:

---

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

• eRositaPhysicsList.hh
• eRositaPhysicsList.cc