CMSEmStandardPhysicsXS Class Reference

#include <CMSEmStandardPhysicsXS.h>

Inheritance diagram for CMSEmStandardPhysicsXS:

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

Private Attributes
bool	fG4HepEmActive
G4double	fGeomFactor
G4double	fLambdaLimit
G4double	fRangeFactor
G4double	fSafetyFactor
G4MscStepLimitType	fStepLimitType

Detailed Description

Definition at line 19 of file CMSEmStandardPhysicsXS.h.

Constructor & Destructor Documentation

CMSEmStandardPhysicsXS()

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

Definition at line 48 of file CMSEmStandardPhysicsXS.cc.

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

    : G4VPhysicsConstructor("CMSEmStandard_emn") {
  SetVerboseLevel(ver);
  // EM parameters specific for this EM physics configuration
  G4EmParameters* param = G4EmParameters::Instance();
  param->SetDefaults();
  param->SetVerbose(ver);
  param->SetApplyCuts(true);
  param->SetStepFunction(0.8, 1 * CLHEP::mm);
  param->SetMscRangeFactor(0.2);
  param->SetMscStepLimitType(fMinimal);
  param->SetFluo(false);
  param->SetUseMottCorrection(true);  // use Mott-correction for e-/e+ msc gs
  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;
  }
  double tcut = p.getParameter<double>("G4TrackingCut") * CLHEP::MeV;
  param->SetLowestElectronEnergy(tcut);
  param->SetLowestMuHadEnergy(tcut);
  fG4HepEmActive = p.getParameter<bool>("G4HepEmActive");
  if (fG4HepEmActive) {
    // At the moment, G4HepEm supports only one configuration of MSC, so use
    // the most generic parameters everywhere.
    param->SetMscRangeFactor(fRangeFactor);
    param->SetMscGeomFactor(fGeomFactor);
    param->SetMscSafetyFactor(fSafetyFactor);
    param->SetMscLambdaLimit(fLambdaLimit);
    param->SetMscStepLimitType(fStepLimitType);
  }
}

~CMSEmStandardPhysicsXS()

CMSEmStandardPhysicsXS::~CMSEmStandardPhysicsXS ( )

overridedefault

Member Function Documentation

ConstructParticle()

void CMSEmStandardPhysicsXS::ConstructParticle ( )

override

Definition at line 89 of file CMSEmStandardPhysicsXS.cc.

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

ConstructProcess()

void CMSEmStandardPhysicsXS::ConstructProcess ( )

override

Definition at line 94 of file CMSEmStandardPhysicsXS.cc.

References nanoDQM_cfi::Electron, fG4HepEmActive, fGeomFactor, fLambdaLimit, fRangeFactor, fSafetyFactor, fStepLimitType, g4SimHits_cfi::G4GammaGeneralProcess, MaterialEffects_cfi::MultipleScattering, contentValuesCheck::ss, and HLT_2023v12_cff::verboseLevel.

                                              {
  if (verboseLevel > 0) {
    edm::LogVerbatim("PhysicsList") << "### " << GetPhysicsName() << " Construct 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 Urban 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
  auto param = G4EmParameters::Instance();
  G4double highEnergyLimit = param->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();

  if (param->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();

  // multiple scattering
  G4UrbanMscModel* msc1 = new G4UrbanMscModel();
  G4WentzelVIModel* msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);

  // msc for HCAL using the Urban model
  G4UrbanMscModel* msc4 = nullptr;
  if (nullptr != aRegion) {
    msc4 = new G4UrbanMscModel();
    msc4->SetHighEnergyLimit(highEnergyLimit);
    msc4->SetRangeFactor(fRangeFactor);
    msc4->SetGeomFactor(fGeomFactor);
    msc4->SetSafetyFactor(fSafetyFactor);
    msc4->SetLambdaLimit(fLambdaLimit);
    msc4->SetStepLimitType(fStepLimitType);
    msc4->SetLocked(true);
  }

  // msc GS with Mott-correction
  G4GoudsmitSaundersonMscModel* msc3 = nullptr;
  if (nullptr != bRegion) {
    msc3 = new G4GoudsmitSaundersonMscModel();
    msc3->SetHighEnergyLimit(highEnergyLimit);
    msc3->SetRangeFactor(0.08);
    msc3->SetSkin(3.);
    msc3->SetStepLimitType(fUseSafetyPlus);
    msc3->SetLocked(true);
  }

  G4TransportationWithMscType transportationWithMsc = G4EmParameters::Instance()->TransportationWithMsc();
  if (transportationWithMsc != G4TransportationWithMscType::fDisabled) {
    G4ProcessManager* procManager = particle->GetProcessManager();
    // Remove default G4Transportation and replace with G4TransportationWithMsc.
    G4VProcess* removed = procManager->RemoveProcess(0);
    if (removed->GetProcessName() != "Transportation") {
      G4Exception("CMSEmStandardPhysics::ConstructProcess",
                  "em0050",
                  FatalException,
                  "replaced process is not G4Transportation!");
    }
    G4TransportationWithMsc* transportWithMsc =
        new G4TransportationWithMsc(G4TransportationWithMsc::ScatteringType::MultipleScattering);
    if (transportationWithMsc == G4TransportationWithMscType::fMultipleSteps) {
      transportWithMsc->SetMultipleSteps(true);
    }
    transportWithMsc->AddMscModel(msc1);
    transportWithMsc->AddMscModel(msc2);
    if (msc4 != nullptr) {
      transportWithMsc->AddMscModel(msc4, -1, aRegion);
    }
    if (msc3 != nullptr) {
      transportWithMsc->AddMscModel(msc3, -1, bRegion);
    }
    procManager->AddProcess(transportWithMsc, -1, 0, 0);
  } else {
    // Register as a separate process.
    G4eMultipleScattering* msc = new G4eMultipleScattering;
    msc->SetEmModel(msc1);
    msc->SetEmModel(msc2);
    if (msc4 != nullptr) {
      msc->AddEmModel(-1, msc4, aRegion);
    }
    if (msc3 != nullptr) {
      msc->AddEmModel(-1, msc3, bRegion);
    }
    ph->RegisterProcess(msc, particle);
  }

  // single scattering
  G4eCoulombScatteringModel* ssm = new G4eCoulombScatteringModel();
  G4CoulombScattering* ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);

  // register processes
  ph->RegisterProcess(new G4eIonisation(), particle);
  ph->RegisterProcess(new G4eBremsstrahlung(), particle);
  ph->RegisterProcess(ss, particle);

  // e+
  particle = G4Positron::Positron();
  msc3 = nullptr;
  msc4 = nullptr;

  // multiple scattering
  msc1 = new G4UrbanMscModel();
  msc2 = new G4WentzelVIModel();
  msc1->SetHighEnergyLimit(highEnergyLimit);
  msc2->SetLowEnergyLimit(highEnergyLimit);

  // msc for HCAL using the Urban model
  if (nullptr != aRegion) {
    msc4 = new G4UrbanMscModel();
    msc4->SetHighEnergyLimit(highEnergyLimit);
    msc4->SetRangeFactor(fRangeFactor);
    msc4->SetGeomFactor(fGeomFactor);
    msc4->SetSafetyFactor(fSafetyFactor);
    msc4->SetLambdaLimit(fLambdaLimit);
    msc4->SetStepLimitType(fStepLimitType);
    msc4->SetLocked(true);
  }

  // msc GS with Mott-correction
  if (nullptr != bRegion) {
    msc3 = new G4GoudsmitSaundersonMscModel();
    msc3->SetHighEnergyLimit(highEnergyLimit);
    msc3->SetRangeFactor(0.08);
    msc3->SetSkin(3.);
    msc3->SetStepLimitType(fUseSafetyPlus);
    msc3->SetLocked(true);
  }

  if (transportationWithMsc != G4TransportationWithMscType::fDisabled) {
    G4ProcessManager* procManager = particle->GetProcessManager();
    // Remove default G4Transportation and replace with G4TransportationWithMsc.
    G4VProcess* removed = procManager->RemoveProcess(0);
    if (removed->GetProcessName() != "Transportation") {
      G4Exception("CMSEmStandardPhysics::ConstructProcess",
                  "em0050",
                  FatalException,
                  "replaced process is not G4Transportation!");
    }
    G4TransportationWithMsc* transportWithMsc =
        new G4TransportationWithMsc(G4TransportationWithMsc::ScatteringType::MultipleScattering);
    if (transportationWithMsc == G4TransportationWithMscType::fMultipleSteps) {
      transportWithMsc->SetMultipleSteps(true);
    }
    transportWithMsc->AddMscModel(msc1);
    transportWithMsc->AddMscModel(msc2);
    if (msc4 != nullptr) {
      transportWithMsc->AddMscModel(msc4, -1, aRegion);
    }
    if (msc3 != nullptr) {
      transportWithMsc->AddMscModel(msc3, -1, bRegion);
    }
    procManager->AddProcess(transportWithMsc, -1, 0, 0);
  } else {
    // Register as a separate process.
    G4eMultipleScattering* msc = new G4eMultipleScattering;
    msc->SetEmModel(msc1);
    msc->SetEmModel(msc2);
    if (msc4 != nullptr) {
      msc->AddEmModel(-1, msc4, aRegion);
    }
    if (msc3 != nullptr) {
      msc->AddEmModel(-1, msc3, bRegion);
    }
    ph->RegisterProcess(msc, particle);
  }

  // single scattering
  ssm = new G4eCoulombScatteringModel();
  ss = new G4CoulombScattering();
  ss->SetEmModel(ssm);
  ss->SetMinKinEnergy(highEnergyLimit);
  ssm->SetLowEnergyLimit(highEnergyLimit);
  ssm->SetActivationLowEnergyLimit(highEnergyLimit);

  // register processes
  ph->RegisterProcess(new G4eIonisation(), particle);
  ph->RegisterProcess(new G4eBremsstrahlung(), particle);
  ph->RegisterProcess(new G4eplusAnnihilation(), particle);
  ph->RegisterProcess(ss, particle);

  if (fG4HepEmActive) {
    auto* hepEmTM = new CMSHepEmTrackingManager(highEnergyLimit);
    G4Electron::Electron()->SetTrackingManager(hepEmTM);
    G4Positron::Positron()->SetTrackingManager(hepEmTM);
  }

  // generic ion
  particle = G4GenericIon::GenericIon();
  G4ionIonisation* ionIoni = new G4ionIonisation();
  ph->RegisterProcess(hmsc, particle);
  ph->RegisterProcess(ionIoni, particle);

  // muons, hadrons, ions
  G4EmBuilder::ConstructCharged(hmsc, pnuc);
}

Member Data Documentation

fG4HepEmActive

bool CMSEmStandardPhysicsXS::fG4HepEmActive

private

Definition at line 33 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().

fGeomFactor

G4double CMSEmStandardPhysicsXS::fGeomFactor

private

Definition at line 29 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().

fLambdaLimit

G4double CMSEmStandardPhysicsXS::fLambdaLimit

private

Definition at line 31 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().

fRangeFactor

G4double CMSEmStandardPhysicsXS::fRangeFactor

private

Definition at line 28 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().

fSafetyFactor

G4double CMSEmStandardPhysicsXS::fSafetyFactor

private

Definition at line 30 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().

fStepLimitType

G4MscStepLimitType CMSEmStandardPhysicsXS::fStepLimitType

private

Definition at line 32 of file CMSEmStandardPhysicsXS.h.

Referenced by CMSEmStandardPhysicsXS(), and ConstructProcess().