G4EmPenelopePhysics Class Reference

#include <G4EmPenelopePhysics.hh>

Inheritance diagram for G4EmPenelopePhysics:

Public Member Functions
	G4EmPenelopePhysics (G4int ver=1, const G4String &name="")
	~G4EmPenelopePhysics () override
void	ConstructParticle () override
void	ConstructProcess () override
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
virtual void	ConstructParticle ()=0
virtual void	ConstructProcess ()=0
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Types inherited from G4VPhysicsConstructor
using	PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
PhysicsBuilder_V	GetBuilders () const
void	AddBuilder (G4PhysicsBuilderInterface *bld)
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 35 of file G4EmPenelopePhysics.hh.

Constructor & Destructor Documentation

G4EmPenelopePhysics()

G4EmPenelopePhysics::G4EmPenelopePhysics ( G4int  ver = 1, const G4String &  name = ""  )

explicit

Definition at line 107 of file G4EmPenelopePhysics.cc.

  : G4VPhysicsConstructor("G4EmPenelope"), verbose(ver)
{
  G4EmParameters* param = G4EmParameters::Instance();
  param->SetDefaults();
  param->SetVerbose(verbose);
  param->SetMinEnergy(100*CLHEP::eV);
  param->SetLowestElectronEnergy(100*CLHEP::eV);
  param->SetNumberOfBinsPerDecade(20);
  param->ActivateAngularGeneratorForIonisation(true);
  param->SetStepFunction(0.2, 10*CLHEP::um);
  param->SetStepFunctionMuHad(0.1, 50*CLHEP::um);
  param->SetStepFunctionLightIons(0.1, 20*CLHEP::um);
  param->SetStepFunctionIons(0.1, 1*CLHEP::um);
  param->SetUseMottCorrection(true);
  param->SetMscStepLimitType(fUseSafetyPlus);
  param->SetMscSkin(3);
  param->SetMscRangeFactor(0.08);
  param->SetMuHadLateralDisplacement(true);
  param->SetFluo(true);
  param->SetMaxNIELEnergy(1*CLHEP::MeV);
  param->SetPIXEElectronCrossSectionModel("Penelope");
  SetPhysicsType(bElectromagnetic);
}

~G4EmPenelopePhysics()

G4EmPenelopePhysics::~G4EmPenelopePhysics ( )

override

Definition at line 134 of file G4EmPenelopePhysics.cc.

{}

Member Function Documentation

ConstructParticle()

void G4EmPenelopePhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 139 of file G4EmPenelopePhysics.cc.

{
  // minimal set of particles for EM physics
  G4EmBuilder::ConstructMinimalEmSet();
}

ConstructProcess()

void G4EmPenelopePhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 147 of file G4EmPenelopePhysics.cc.

{
  if(verbose > 1) {
    G4cout << "### " << GetPhysicsName() << " Construct Processes " << G4endl;
  }
  G4EmBuilder::PrepareEMPhysics();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
  
  // processes used by several particles
  G4ePairProduction* ee = new G4ePairProduction();
  G4hMultipleScattering* hmsc = new G4hMultipleScattering("ionmsc");
 
  // high energy limit for e+- scattering models
  G4double highEnergyLimit = G4EmParameters::Instance()->MscEnergyLimit();
 
  // nuclear stopping
  G4double nielEnergyLimit = G4EmParameters::Instance()->MaxNIELEnergy();
  G4NuclearStopping* pnuc = new G4NuclearStopping();
  pnuc->SetMaxKinEnergy(nielEnergyLimit);
 
  //Applicability range for Penelope models
  //for higher energies, the Standard models are used   
  G4double PenelopeHighEnergyLimit = 1.0*CLHEP::GeV;
 
  // Add gamma EM Processes
  G4ParticleDefinition* particle = G4Gamma::Gamma();
 
  //Photo-electric effect
  G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
  thePhotoElectricEffect->SetEmModel(new G4PEEffectFluoModel());
  G4PenelopePhotoElectricModel* thePEPenelopeModel = new G4PenelopePhotoElectricModel();   
  thePEPenelopeModel->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  thePhotoElectricEffect->AddEmModel(0, thePEPenelopeModel);
 
  //Compton scattering
  G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
  G4PenelopeComptonModel* theComptonPenelopeModel = new G4PenelopeComptonModel();
  theComptonScattering->SetEmModel(new G4KleinNishinaModel());
  theComptonPenelopeModel->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  theComptonScattering->AddEmModel(0, theComptonPenelopeModel);
 
  //Gamma conversion
  G4GammaConversion* theGammaConversion = new G4GammaConversion();
  G4PenelopeGammaConversionModel* theGCPenelopeModel = new G4PenelopeGammaConversionModel();
  theGCPenelopeModel->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  theGammaConversion->AddEmModel(0, theGCPenelopeModel);
 
  //Rayleigh scattering
  G4RayleighScattering* theRayleigh = new G4RayleighScattering();
  G4PenelopeRayleighModel* theRayleighPenelopeModel = new G4PenelopeRayleighModel();
  theRayleigh->SetEmModel(theRayleighPenelopeModel);
 
  ph->RegisterProcess(thePhotoElectricEffect, particle);
  ph->RegisterProcess(theComptonScattering, particle);
  ph->RegisterProcess(theGammaConversion, particle);
  ph->RegisterProcess(theRayleigh, particle);
 
  // e-
  particle = G4Electron::Electron();
 
  // multiple scattering
  G4eMultipleScattering* msc = new G4eMultipleScattering;
  G4GoudsmitSaundersonMscModel* msc1 = new G4GoudsmitSaundersonMscModel();
  G4WentzelVIModel* msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);
  msc->SetEmModel(msc1);
  msc->SetEmModel(msc2);
 
  G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel(); 
  G4CoulombScattering* ss = new G4CoulombScattering();
  ss->SetEmModel(ssm); 
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);
      
  //Ionisation
  G4eIonisation* eioni = new G4eIonisation();
  G4PenelopeIonisationModel* theIoniPenelope = new G4PenelopeIonisationModel();
  theIoniPenelope->SetHighEnergyLimit(PenelopeHighEnergyLimit);     
  eioni->AddEmModel(0, theIoniPenelope, new G4UniversalFluctuation());
      
  //Bremsstrahlung
  G4eBremsstrahlung* brem = new G4eBremsstrahlung();
  G4PenelopeBremsstrahlungModel* theBremPenelope = new G4PenelopeBremsstrahlungModel();
  theBremPenelope->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  brem->SetEmModel(theBremPenelope);
 
  // register processes
  ph->RegisterProcess(msc, particle);
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(brem, particle);
  ph->RegisterProcess(ee, particle);
  ph->RegisterProcess(ss, particle);
 
  // e+
  particle = G4Positron::Positron();
    
  // multiple scattering
  msc = new G4eMultipleScattering;
  msc1 = new G4GoudsmitSaundersonMscModel();
  msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);
  msc->SetEmModel(msc1);
  msc->SetEmModel(msc2);
 
  ssm = new G4eCoulombScatteringModel(); 
  ss = new G4CoulombScattering();
  ss->SetEmModel(ssm); 
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);
      
  //Ionisation
  eioni = new G4eIonisation();
  theIoniPenelope = new G4PenelopeIonisationModel();
  theIoniPenelope->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  eioni->AddEmModel(0,theIoniPenelope, new G4UniversalFluctuation());
 
  //Bremsstrahlung
  brem = new G4eBremsstrahlung();
  theBremPenelope = new G4PenelopeBremsstrahlungModel();
  theBremPenelope->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  brem->SetEmModel(theBremPenelope);
      
  //Annihilation
  G4eplusAnnihilation* anni = new G4eplusAnnihilation();
  G4PenelopeAnnihilationModel* theAnnPenelope = new G4PenelopeAnnihilationModel();
  theAnnPenelope->SetHighEnergyLimit(PenelopeHighEnergyLimit);
  anni->AddEmModel(0, theAnnPenelope);
 
  // register processes
  ph->RegisterProcess(msc, particle);
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(brem, particle);
  ph->RegisterProcess(ee, particle);
  ph->RegisterProcess(anni, particle);
  ph->RegisterProcess(ss, particle);
 
  // generic ion
  particle = G4GenericIon::GenericIon();
  G4ionIonisation* ionIoni = new G4ionIonisation();
  ionIoni->SetEmModel(new G4IonParametrisedLossModel());
  ph->RegisterProcess(hmsc, particle);
  ph->RegisterProcess(ionIoni, particle);
  ph->RegisterProcess(pnuc, particle);
 
  // muons, hadrons, ions
  G4EmBuilder::ConstructCharged(hmsc, pnuc);
 
  // extra configuration
  G4EmModelActivator mact(GetPhysicsName());
}

---

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

• G4EmPenelopePhysics.hh
• G4EmPenelopePhysics.cc