CMSEmStandardPhysics Class Reference

#include <CMSEmStandardPhysics.h>

Inheritance diagram for CMSEmStandardPhysics:

Public Member Functions
	CMSEmStandardPhysics (G4int ver, const edm::ParameterSet &p)
void	ConstructParticle () override
void	ConstructProcess () override
	~CMSEmStandardPhysics () override

Private Attributes
G4double	fGeomFactor
G4double	fLambdaLimit
G4double	fRangeFactor
G4double	fSafetyFactor
G4MscStepLimitType	fStepLimitType
G4int	verbose

Detailed Description

Definition at line 10 of file CMSEmStandardPhysics.h.

Constructor & Destructor Documentation

CMSEmStandardPhysics()

CMSEmStandardPhysics::CMSEmStandardPhysics	(	G4int	ver,
		const edm::ParameterSet &	p
	)

Definition at line 46 of file CMSEmStandardPhysics.cc.

    : G4VPhysicsConstructor("CMSEmStandard_emm"), verbose(ver) {
  G4EmParameters* param = G4EmParameters::Instance();
  param->SetDefaults();
  param->SetVerbose(verbose);
  param->SetApplyCuts(true);
  param->SetStepFunction(0.8, 1 * CLHEP::mm);
  param->SetMscRangeFactor(0.2);
  param->SetMscStepLimitType(fMinimal);
  SetPhysicsType(bElectromagnetic);
  fRangeFactor = p.getParameter<double>("G4MscRangeFactor");
  fGeomFactor = p.getParameter<double>("G4MscGeomFactor");
  fSafetyFactor = p.getParameter<double>("G4MscSafetyFactor");
  fLambdaLimit = p.getParameter<double>("G4MscLambdaLimit") * CLHEP::mm;
  std::string msc = p.getParameter<std::string>("G4MscStepLimit");
  fStepLimitType = fUseSafety;
  if (msc == "UseSafetyPlus") {
    fStepLimitType = fUseSafetyPlus;
  }
  if (msc == "Minimal") {
    fStepLimitType = fMinimal;
  }
}

References fGeomFactor, fLambdaLimit, fRangeFactor, fSafetyFactor, fStepLimitType, AlCaHLTBitMon_ParallelJobs::p, AlCaHLTBitMon_QueryRunRegistry::string, and verbose.

~CMSEmStandardPhysics()

CMSEmStandardPhysics::~CMSEmStandardPhysics ( )

override

Definition at line 70 of file CMSEmStandardPhysics.cc.

{}

Member Function Documentation

ConstructParticle()

void CMSEmStandardPhysics::ConstructParticle ( )

override

Definition at line 72 of file CMSEmStandardPhysics.cc.

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

ConstructProcess()

void CMSEmStandardPhysics::ConstructProcess ( )

override

Definition at line 77 of file CMSEmStandardPhysics.cc.

                                            {
  if (verbose > 0) {
    edm::LogVerbatim("PhysicsList") << "### " << GetPhysicsName() << " Construct EM Processes";
  }
 
  // This EM builder takes default models of Geant4 10 EMV.
  // Multiple scattering by Urban for all particles
  // except e+e- below 100 MeV for which the Urban93 model is used
  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 highEnergyLimit = G4EmParameters::Instance()->MscEnergyLimit();
 
  const G4Region* aRegion = G4RegionStore::GetInstance()->GetRegion("HcalRegion", false);
  const G4Region* bRegion = G4RegionStore::GetInstance()->GetRegion("HGCalRegion", false);
 
  // Add gamma EM Processes
  G4ParticleDefinition* particle = G4Gamma::Gamma();
 
  G4PhotoElectricEffect* pee = new G4PhotoElectricEffect();
  pee->SetEmModel(new G4LivermorePhotoElectricModel());
 
  if (G4EmParameters::Instance()->GeneralProcessActive()) {
    G4GammaGeneralProcess* sp = new G4GammaGeneralProcess();
    sp->AddEmProcess(pee);
    sp->AddEmProcess(new G4ComptonScattering());
    sp->AddEmProcess(new G4GammaConversion());
    G4LossTableManager::Instance()->SetGammaGeneralProcess(sp);
    ph->RegisterProcess(sp, particle);
 
  } else {
    ph->RegisterProcess(pee, particle);
    ph->RegisterProcess(new G4ComptonScattering(), particle);
    ph->RegisterProcess(new G4GammaConversion(), particle);
  }
 
  // e-
  particle = G4Electron::Electron();
 
  G4eIonisation* eioni = new G4eIonisation();
 
  G4eMultipleScattering* msc = new G4eMultipleScattering;
  G4UrbanMscModel* msc1 = new G4UrbanMscModel();
  G4WentzelVIModel* msc2 = new G4WentzelVIModel();
  G4UrbanMscModel* msc3 = new G4UrbanMscModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);
  msc3->SetHighEnergyLimit(highEnergyLimit);
  msc3->SetRangeFactor(fRangeFactor);
  msc3->SetGeomFactor(fGeomFactor);
  msc3->SetSafetyFactor(fSafetyFactor);
  msc3->SetLambdaLimit(fLambdaLimit);
  msc3->SetStepLimitType(fStepLimitType);
  msc3->SetLocked(true);
  msc->SetEmModel(msc1);
  msc->SetEmModel(msc2);
  if (nullptr != aRegion) {
    msc->AddEmModel(-1, msc3, aRegion);
  }
  if (nullptr != bRegion) {
    msc->AddEmModel(-1, msc3, bRegion);
  }
 
  G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel();
  G4CoulombScattering* ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);
 
  ph->RegisterProcess(msc, particle);
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(new G4eBremsstrahlung(), particle);
  ph->RegisterProcess(ss, particle);
 
  // e+
  particle = G4Positron::Positron();
  eioni = new G4eIonisation();
 
  msc = new G4eMultipleScattering();
  msc1 = new G4UrbanMscModel();
  msc2 = new G4WentzelVIModel();
  msc3 = new G4UrbanMscModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);
  msc3->SetHighEnergyLimit(highEnergyLimit);
  msc3->SetRangeFactor(fRangeFactor);
  msc3->SetGeomFactor(fGeomFactor);
  msc3->SetSafetyFactor(fSafetyFactor);
  msc3->SetLambdaLimit(fLambdaLimit);
  msc3->SetStepLimitType(fStepLimitType);
  msc3->SetLocked(true);
  msc->SetEmModel(msc1);
  msc->SetEmModel(msc2);
  if (nullptr != aRegion) {
    msc->AddEmModel(-1, msc3, aRegion);
  }
  if (nullptr != bRegion) {
    msc->AddEmModel(-1, msc3, bRegion);
  }
 
  ssm = new G4eCoulombScatteringModel();
  ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);
 
  ph->RegisterProcess(msc, particle);
  ph->RegisterProcess(eioni, particle);
  ph->RegisterProcess(new G4eBremsstrahlung(), 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);
}

References nanoDQM_cff::Electron, fGeomFactor, fLambdaLimit, fRangeFactor, fSafetyFactor, fStepLimitType, contentValuesCheck::ss, and verbose.

Member Data Documentation

fGeomFactor

G4double CMSEmStandardPhysics::fGeomFactor

private

Definition at line 20 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().

fLambdaLimit

G4double CMSEmStandardPhysics::fLambdaLimit

private

Definition at line 22 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().

fRangeFactor

G4double CMSEmStandardPhysics::fRangeFactor

private

Definition at line 19 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().

fSafetyFactor

G4double CMSEmStandardPhysics::fSafetyFactor

private

Definition at line 21 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().

fStepLimitType

G4MscStepLimitType CMSEmStandardPhysics::fStepLimitType

private

Definition at line 23 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().

verbose

G4int CMSEmStandardPhysics::verbose

private

Definition at line 24 of file CMSEmStandardPhysics.h.

Referenced by CMSEmStandardPhysics(), and ConstructProcess().