G4EmStandardPhysicsGS Class Reference

#include <G4EmStandardPhysicsGS.hh>

Inheritance diagram for G4EmStandardPhysicsGS:

Public Member Functions
	G4EmStandardPhysicsGS (G4int ver=1, const G4String &name="")
	~G4EmStandardPhysicsGS () 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 ()
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const

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 = 0
G4String	namePhysics = ""
G4int	typePhysics = 0
G4ParticleTable *	theParticleTable = nullptr
G4int	g4vpcInstanceID = 0
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 51 of file G4EmStandardPhysicsGS.hh.

Constructor & Destructor Documentation

G4EmStandardPhysicsGS()

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

explicit

Definition at line 108 of file G4EmStandardPhysicsGS.cc.

  : G4VPhysicsConstructor("G4EmStandardGS")
{
  SetVerboseLevel(ver);
  G4EmParameters* param = G4EmParameters::Instance();
  param->SetDefaults();
  param->SetVerbose(ver);
  // use a denser discrete kinetic energy grid for more accurate interpolation
  param->SetNumberOfBinsPerDecade(16);
  // set the continuous step limit to: 0.2*Range that goes to Range below 10 um
  param->SetStepFunction(0.2, 10*CLHEP::um);
  // set the GS MSC model for e-/e+ to be used below 1.0 GeV with its (Mott,
  // screening, scattering power) corrections activated and with the accurate
  // stepping and boundary crossing algorithms (no other options since 11.4)
  // with a skin of 3 elastic MFP near boundary
  param->SetMscEnergyLimit(1.0*CLHEP::GeV);
  param->SetUseMottCorrection(true);
  param->SetMscStepLimitType(fUseSafetyPlus);
  param->SetMscSkin(3);
  param->SetMscRangeFactor(0.08);
  // activate fluoresence, i.e. emission of characteristic X-ray
  param->SetFluo(true);
  // set the energy loss fluctuation type
  param->SetFluctuationType(fUrbanFluctuation);
  SetPhysicsType(bElectromagnetic);
}

~G4EmStandardPhysicsGS()

G4EmStandardPhysicsGS::~G4EmStandardPhysicsGS ( )

override

Definition at line 137 of file G4EmStandardPhysicsGS.cc.

{}

Member Function Documentation

ConstructParticle()

void G4EmStandardPhysicsGS::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 142 of file G4EmStandardPhysicsGS.cc.

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

ConstructProcess()

void G4EmStandardPhysicsGS::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 150 of file G4EmStandardPhysicsGS.cc.

{
  if(verboseLevel > 1) {
    G4cout << "### " << GetPhysicsName() << " Construct Processes " << G4endl;
  }
  G4EmBuilder::PrepareEMPhysics();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  // processes used by several particles
  G4hMultipleScattering* hmsc = new G4hMultipleScattering("ionmsc");
  G4NuclearStopping* pnuc(nullptr);
 
  // high energy limit for e+- scattering models and bremsstrahlung
  G4double mscEnergyLimit = G4EmParameters::Instance()->MscEnergyLimit();
 
  // Add gamma EM processes
 
  // gamma
  G4ParticleDefinition* particle = G4Gamma::Gamma();
 
  G4PhotoElectricEffect* pe = new G4PhotoElectricEffect();
  pe->SetEmModel(new G4LivermorePhotoElectricModel());
 
  G4ComptonScattering* cs = new G4ComptonScattering;
  G4GammaConversion* gc = new G4GammaConversion;
  G4RayleighScattering* rs = new G4RayleighScattering;
 
  if (G4EmParameters::Instance()->GeneralProcessActive()) {
    G4GammaGeneralProcess* sp = new G4GammaGeneralProcess();
    sp->AddEmProcess(pe);
    sp->AddEmProcess(cs);
    sp->AddEmProcess(gc);
    sp->AddEmProcess(rs);
    G4LossTableManager::Instance()->SetGammaGeneralProcess(sp);
    ph->RegisterProcess(sp, particle);
  } else {
    ph->RegisterProcess(pe, particle);
    ph->RegisterProcess(cs, particle);
    ph->RegisterProcess(gc, particle);
    ph->RegisterProcess(rs, particle);
  }
 
  // e-
  particle = G4Electron::Electron();
 
  // msc: GS[:100 MeV] + WentzelVI[100 MeV:]
  G4GoudsmitSaundersonMscModel* msc1 = new G4GoudsmitSaundersonMscModel();
  G4WentzelVIModel* msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(mscEnergyLimit);
  msc2->SetLowEnergyLimit(mscEnergyLimit);
  G4EmBuilder::ConstructElectronMscProcess(msc1, msc2, particle);
  // (WVI is a mixed model, i.e. needs single scattering)
  G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel();
  G4CoulombScattering* ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(mscEnergyLimit);
  ssm->SetLowEnergyLimit(mscEnergyLimit);
  ssm->SetActivationLowEnergyLimit(mscEnergyLimit);
 
  // ionisation: Penelope[:1.0 GeV] + Moller[1.0 GeV:]
  G4eIonisation* eioni = new G4eIonisation();
  eioni->SetFluctModel(G4EmStandUtil::ModelOfFluctuations());
  G4VEmModel* theIoniMod = new G4PenelopeIonisationModel();
  theIoniMod->SetHighEnergyLimit(1.0*CLHEP::GeV);
  eioni->AddEmModel(0, theIoniMod);
 
  // bremsstrahlung: Seltzer-Berger[:1.0 GeV] + extended Bethe–Heitler[1.0 GeV:]
  G4eBremsstrahlung* brem = new G4eBremsstrahlung();
  G4SeltzerBergerModel* br1 = new G4SeltzerBergerModel();
  G4eBremsstrahlungRelModel* br2 = new G4eBremsstrahlungRelModel();
  br1->SetAngularDistribution(new G4Generator2BS());
  br2->SetAngularDistribution(new G4Generator2BS());
  brem->SetEmModel(br1);
  brem->SetEmModel(br2);
  br1->SetHighEnergyLimit(1.0*CLHEP::GeV);
 
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(brem, particle);
  ph->RegisterProcess(ss, particle);
 
  // e+
  particle = G4Positron::Positron();
 
  // msc: GS[:100 MeV] + WentzelVI[100 MeV:]
  msc1 = new G4GoudsmitSaundersonMscModel();
  msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(mscEnergyLimit);
  msc2->SetLowEnergyLimit(mscEnergyLimit);
  G4EmBuilder::ConstructElectronMscProcess(msc1, msc2, particle);
  // (WVI is a mixed model, i.e. needs single scattering)
  ssm = new G4eCoulombScatteringModel();
  ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(mscEnergyLimit);
  ssm->SetLowEnergyLimit(mscEnergyLimit);
  ssm->SetActivationLowEnergyLimit(mscEnergyLimit);
 
  // ionisation: Penelope[:1.0 GeV] + Bhabha[1.0 GeV:]
  eioni = new G4eIonisation();
  eioni->SetFluctModel(G4EmStandUtil::ModelOfFluctuations());
  G4VEmModel* pen = new G4PenelopeIonisationModel();
  pen->SetHighEnergyLimit(1.0*CLHEP::GeV);
  eioni->AddEmModel(0, pen);
 
  // bremsstrahlung: Seltzer-Berger[:1.0 GeV] + extended Bethe–Heitler[1.0 GeV:]
  brem = new G4eBremsstrahlung();
  br1 = new G4SeltzerBergerModel();
  br2 = new G4eBremsstrahlungRelModel();
  br1->SetAngularDistribution(new G4Generator2BS());
  br2->SetAngularDistribution(new G4Generator2BS());
  brem->SetEmModel(br1);
  brem->SetEmModel(br2);
  br1->SetHighEnergyLimit(1.0*CLHEP::GeV);
 
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(brem, particle);
  ph->RegisterProcess(new G4eplusAnnihilation(), particle);
  ph->RegisterProcess(ss, particle);
 
  // generic ion
  particle = G4GenericIon::GenericIon();
  G4ionIonisation* ionIoni = new G4ionIonisation();
  ph->RegisterProcess(hmsc, particle);
  ph->RegisterProcess(ionIoni, particle);
 
  // muons, hadrons, ions
  G4EmBuilder::ConstructCharged(hmsc, pnuc);
 
  // extra configuration
  G4EmModelActivator mact(GetPhysicsName());
}

---

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

• G4EmStandardPhysicsGS.hh
• G4EmStandardPhysicsGS.cc