PhysListEmPenelope Class Reference

#include <PhysListEmPenelope.hh>

Inheritance diagram for PhysListEmPenelope:

Collaboration diagram for PhysListEmPenelope:

Public Member Functions
	PhysListEmPenelope (const G4String &name="penelope")
	~PhysListEmPenelope ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
	PhysListEmPenelope (const G4String &name="penelope")
	~PhysListEmPenelope ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
	PhysListEmPenelope (const G4String &name="Penelope")
	~PhysListEmPenelope ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4int	g4vpcInstanceID
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 42 of file PhysListEmPenelope.hh.

Constructor & Destructor Documentation

PhysListEmPenelope::PhysListEmPenelope

(

const G4String &

name = "penelope"

)

Definition at line 82 of file PhysListEmPenelope.cc.

  :  G4VPhysicsConstructor(name)
{ }

PhysListEmPenelope::~PhysListEmPenelope

(

)

Definition at line 88 of file PhysListEmPenelope.cc.

{ }

PhysListEmPenelope::PhysListEmPenelope

(

const G4String &

name = "penelope"

)

PhysListEmPenelope::~PhysListEmPenelope

(

)

PhysListEmPenelope::PhysListEmPenelope

(

const G4String &

name = "Penelope"

)

PhysListEmPenelope::~PhysListEmPenelope

(

)

Member Function Documentation

virtual void PhysListEmPenelope::ConstructParticle ( void  )

inlinevirtual

Implements G4VPhysicsConstructor.

Definition at line 46 of file PhysListEmPenelope.hh.

{ };

virtual void PhysListEmPenelope::ConstructParticle ( void  )

inlinevirtual

Implements G4VPhysicsConstructor.

Definition at line 50 of file PhysListEmPenelope.hh.

{};

virtual void PhysListEmPenelope::ConstructParticle ( void  )

inlinevirtual

Implements G4VPhysicsConstructor.

Definition at line 50 of file PhysListEmPenelope.hh.

{};

virtual void PhysListEmPenelope::ConstructProcess ( )

virtual

Implements G4VPhysicsConstructor.

void PhysListEmPenelope::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 93 of file PhysListEmPenelope.cc.

{
  // Add standard EM Processes

  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();

    //Applicability range for Penelope models
    //for higher energies, the Standard models are used   
    G4double highEnergyLimit = 1*GeV;
         
    if (particleName == "gamma") {
      // gamma         

      G4PhotoElectricEffect* phot = new G4PhotoElectricEffect();
      G4PenelopePhotoElectricModel* 
      photModel = new G4PenelopePhotoElectricModel();
      photModel->SetHighEnergyLimit(highEnergyLimit);
      phot->AddEmModel(0, photModel);
      pmanager->AddDiscreteProcess(phot);

      G4ComptonScattering* compt = new G4ComptonScattering();
      G4PenelopeComptonModel* 
      comptModel = new G4PenelopeComptonModel();
      comptModel->SetHighEnergyLimit(highEnergyLimit);
      compt->AddEmModel(0, comptModel);
      pmanager->AddDiscreteProcess(compt);

      G4GammaConversion* conv = new G4GammaConversion();
      G4PenelopeGammaConversionModel* 
      convModel = new G4PenelopeGammaConversionModel();
      convModel->SetHighEnergyLimit(highEnergyLimit);
      conv->AddEmModel(0, convModel);
      pmanager->AddDiscreteProcess(conv);

      G4RayleighScattering* rayl = new G4RayleighScattering();
      G4PenelopeRayleighModel* 
      raylModel = new G4PenelopeRayleighModel();
      raylModel->SetHighEnergyLimit(highEnergyLimit);
      rayl->AddEmModel(0, raylModel);
      pmanager->AddDiscreteProcess(rayl);
      
    } else if (particleName == "e-") {
      //electron

      G4eIonisation* eIoni = new G4eIonisation();
      G4PenelopeIonisationModel* 
      eIoniModel = new G4PenelopeIonisationModel();
      eIoniModel->SetHighEnergyLimit(highEnergyLimit); 
      eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
      pmanager->AddProcess(eIoni,                   -1,-1, 1);
      
      G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
      G4PenelopeBremsstrahlungModel* 
      eBremModel = new G4PenelopeBremsstrahlungModel();
      eBremModel->SetHighEnergyLimit(highEnergyLimit);
      eBrem->AddEmModel(0, eBremModel);
      pmanager->AddProcess(eBrem,                   -1,-1, 2);
                  
    } else if (particleName == "e+") {
      //positron
      G4eIonisation* eIoni = new G4eIonisation();
      G4PenelopeIonisationModel* 
      eIoniModel = new G4PenelopeIonisationModel();
      eIoniModel->SetHighEnergyLimit(highEnergyLimit); 
      eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
      pmanager->AddProcess(eIoni,                   -1,-1, 1);
      
      G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
      G4PenelopeBremsstrahlungModel* 
      eBremModel = new G4PenelopeBremsstrahlungModel();
      eBremModel->SetHighEnergyLimit(highEnergyLimit);
      eBrem->AddEmModel(0, eBremModel);
      pmanager->AddProcess(eBrem,                   -1,-1, 2);      

      G4eplusAnnihilation* eAnni = new G4eplusAnnihilation();
      G4PenelopeAnnihilationModel* 
      eAnniModel = new G4PenelopeAnnihilationModel();
      eAnniModel->SetHighEnergyLimit(highEnergyLimit); 
      eAnni->AddEmModel(0, eAnniModel);
      pmanager->AddProcess(eAnni,                    0,-1, 3);
            
    } else if( particleName == "mu+" || 
               particleName == "mu-"    ) {
      //muon  
      pmanager->AddProcess(new G4MuIonisation,      -1,-1, 1);
      pmanager->AddProcess(new G4MuBremsstrahlung,  -1,-1, 2);
      pmanager->AddProcess(new G4MuPairProduction,  -1,-1, 3);       
     
    } else if( particleName == "alpha" || particleName == "GenericIon" ) { 
      pmanager->AddProcess(new G4ionIonisation,     -1,-1, 1);

    } else if ((!particle->IsShortLived()) &&
               (particle->GetPDGCharge() != 0.0) && 
               (particle->GetParticleName() != "chargedgeantino")) {
      //all others charged particles except geantino
      pmanager->AddProcess(new G4hIonisation,       -1,-1, 1);
    }
  }
}

Here is the call graph for this function:

virtual void PhysListEmPenelope::ConstructProcess ( )

virtual

Implements G4VPhysicsConstructor.

---

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

• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm13/include/PhysListEmPenelope.hh
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm14/include/PhysListEmPenelope.hh
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm18/include/PhysListEmPenelope.hh
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm13/src/PhysListEmPenelope.cc
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm14/src/PhysListEmPenelope.cc
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm18/src/PhysListEmPenelope.cc