ParametrisedEMPhysics Class Reference

#include <ParametrisedEMPhysics.h>

Inheritance diagram for ParametrisedEMPhysics:

Classes
struct	TLSmod

Public Member Functions
void	ConstructParticle () override
void	ConstructProcess () override
	ParametrisedEMPhysics (const std::string &name, const edm::ParameterSet &p)
	~ParametrisedEMPhysics () override

Private Member Functions
void	ModifyTransportation (const G4ParticleDefinition *, int ntry, double th1, double th2)

Private Attributes
edm::ParameterSet	theParSet

Static Private Attributes
static G4ThreadLocal TLSmod *	m_tpmod = 0

Detailed Description

Definition at line 16 of file ParametrisedEMPhysics.h.

Constructor & Destructor Documentation

ParametrisedEMPhysics::ParametrisedEMPhysics	(	const std::string &	name,
		const edm::ParameterSet &	p
	)

Definition at line 54 of file ParametrisedEMPhysics.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GeV, mps_fire::i, MeV, nt, and theParSet.

  : G4VPhysicsConstructor(name), theParSet(p) 
{
  // bremsstrahlung threshold and EM verbosity
  G4EmParameters* param = G4EmParameters::Instance();
  G4int verb = theParSet.getUntrackedParameter<int>("Verbosity",0);
  param->SetVerbose(verb);

  G4double bremth = theParSet.getParameter<double>("G4BremsstrahlungThreshold")*GeV; 
  param->SetBremsstrahlungTh(bremth);

  bool fluo = theParSet.getParameter<bool>("FlagFluo");
  param->SetFluo(fluo);

  bool modifyT = theParSet.getParameter<bool>("ModifyTransportation");
  double th1 = theParSet.getUntrackedParameter<double>("ThresholdWarningEnergy")*MeV;
  double th2 = theParSet.getUntrackedParameter<double>("ThresholdImportantEnergy")*MeV;
  int nt = theParSet.getUntrackedParameter<int>("ThresholdTrials");

  edm::LogVerbatim("SimG4CoreApplication") 
    << "ParametrisedEMPhysics::ConstructProcess: bremsstrahlung threshold Eth= "
    << bremth/GeV << " GeV" 
    << "\n                                         verbosity= " << verb
    << "  fluoFlag: " << fluo << "  modifyTransport: " << modifyT 
    << "  Ntrials= " << nt 
    << "\n                                         ThWarning(MeV)= " << th1
    << "  ThException(MeV)= " << th2;

  // Russian roulette and tracking cut for e+-
  double energyLim = 
    theParSet.getParameter<double>("RusRoElectronEnergyLimit")*MeV;
  if(energyLim > 0.0) {
    const G4int NREG = 6; 
    const G4String rname[NREG] = {"EcalRegion", "HcalRegion", "MuonIron",
                                  "PreshowerRegion","CastorRegion",
                                  "DefaultRegionForTheWorld"};
    G4double rrfact[NREG] = { 1.0 };

    rrfact[0] = theParSet.getParameter<double>("RusRoEcalElectron");
    rrfact[1] = theParSet.getParameter<double>("RusRoHcalElectron");
    rrfact[2] = theParSet.getParameter<double>("RusRoMuonIronElectron");
    rrfact[3] = theParSet.getParameter<double>("RusRoPreShowerElectron");
    rrfact[4] = theParSet.getParameter<double>("RusRoCastorElectron");
    rrfact[5] = theParSet.getParameter<double>("RusRoWorldElectron");
    for(int i=0; i<NREG; ++i) {
      if(rrfact[i] < 1.0) {
        param->ActivateSecondaryBiasing("eIoni",rname[i],rrfact[i],energyLim);
        param->ActivateSecondaryBiasing("hIoni",rname[i],rrfact[i],energyLim);
        edm::LogVerbatim("SimG4CoreApplication") 
          << "ParametrisedEMPhysics: Russian Roulette"
          << " for e- Prob= " << rrfact[i]  
          << " Elimit(MeV)= " << energyLim/CLHEP::MeV
          << " inside " << rname[i];
      }
    }
  }
}

ParametrisedEMPhysics::~ParametrisedEMPhysics ( )

override

Definition at line 113 of file ParametrisedEMPhysics.cc.

References m_tpmod.

                                              {
  if(m_tpmod) {
    delete m_tpmod;
    m_tpmod = nullptr;
  } 
}

Member Function Documentation

void ParametrisedEMPhysics::ConstructParticle ( )

override

Definition at line 120 of file ParametrisedEMPhysics.cc.

{
  G4LeptonConstructor pLeptonConstructor;
  pLeptonConstructor.ConstructParticle();

  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

  G4ShortLivedConstructor pShortLivedConstructor;
  pShortLivedConstructor.ConstructParticle();  
    
  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

void ParametrisedEMPhysics::ConstructProcess ( )

override

Definition at line 138 of file ParametrisedEMPhysics.cc.

References nanoDQM_cff::Electron, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_tpmod, MeV, ModifyTransportation(), nt, ParametrisedEMPhysics::TLSmod::theEcalEMShowerModel, ParametrisedEMPhysics::TLSmod::theEcalHadShowerModel, ParametrisedEMPhysics::TLSmod::theElectronLimiter, ParametrisedEMPhysics::TLSmod::theFastSimulationManagerProcess, ParametrisedEMPhysics::TLSmod::theHcalEMShowerModel, ParametrisedEMPhysics::TLSmod::theHcalHadShowerModel, theParSet, and ParametrisedEMPhysics::TLSmod::thePositronLimiter.

                                             {

  // GFlash part 
  bool gem  = theParSet.getParameter<bool>("GflashEcal");
  bool ghad = theParSet.getParameter<bool>("GflashHcal");
  bool gemHad  = theParSet.getParameter<bool>("GflashEcalHad");
  bool ghadHad = theParSet.getParameter<bool>("GflashHcalHad");

  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
  if(gem || ghad || gemHad || ghadHad) {
    if(!m_tpmod) { m_tpmod = new TLSmod; }
    edm::LogVerbatim("SimG4CoreApplication") 
      << "ParametrisedEMPhysics: GFlash Construct for e+-: " 
      << gem << "  " << ghad << " for hadrons: " << gemHad << "  " << ghadHad;

    m_tpmod->theFastSimulationManagerProcess.reset(new G4FastSimulationManagerProcess());

    if(gem || ghad) {
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4Electron::Electron());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4Positron::Positron());
    }
    if(gemHad || ghadHad) {
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4Proton::Proton());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4AntiProton::AntiProton());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4PionPlus::PionPlus());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4PionMinus::PionMinus());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4KaonPlus::KaonPlus());
      ph->RegisterProcess(m_tpmod->theFastSimulationManagerProcess.get(), G4KaonMinus::KaonMinus());
    }

    if(gem || gemHad) {
      G4Region* aRegion = 
        G4RegionStore::GetInstance()->GetRegion("EcalRegion");
      
      if(!aRegion){
        edm::LogVerbatim("SimG4CoreApplication") 
          << "ParametrisedEMPhysics::ConstructProcess: " 
          << "EcalRegion is not defined, GFlash will not be enabled for ECAL!";
        
      } else {
        if(gem) {

          //Electromagnetic Shower Model for ECAL
          m_tpmod->theEcalEMShowerModel.reset(new GFlashEMShowerModel("GflashEcalEMShowerModel",
                                              aRegion,theParSet));
        }
        if(gemHad) {

          //Electromagnetic Shower Model for ECAL
          m_tpmod->theEcalHadShowerModel.reset(new GFlashHadronShowerModel("GflashEcalHadShowerModel",
                                               aRegion,theParSet));
        }    
      }
    }
    if(ghad || ghadHad) {
      G4Region* aRegion = 
        G4RegionStore::GetInstance()->GetRegion("HcalRegion");
      if(!aRegion) {
        edm::LogVerbatim("SimG4CoreApplication") 
          << "ParametrisedEMPhysics::ConstructProcess: " 
          << "HcalRegion is not defined, GFlash will not be enabled for HCAL!";
        
      } else {
        if(ghad) {

          //Electromagnetic Shower Model for HCAL
          m_tpmod->theHcalEMShowerModel.reset(new GFlashEMShowerModel("GflashHcalEMShowerModel",
                                              aRegion,theParSet));
        }
        if(ghadHad) {

          //Electromagnetic Shower Model for ECAL
          m_tpmod->theHcalHadShowerModel.reset(new GFlashHadronShowerModel("GflashHcalHadShowerModel",
                                               aRegion,theParSet));
        }    
      }
    }
  }

  // Step limiters for e+-
  bool eLimiter = theParSet.getParameter<bool>("ElectronStepLimit");
  bool rLimiter = theParSet.getParameter<bool>("ElectronRangeTest");
  bool pLimiter = theParSet.getParameter<bool>("PositronStepLimit");

  if(eLimiter || rLimiter) {
    if(!m_tpmod) { m_tpmod = new TLSmod; }
    m_tpmod->theElectronLimiter.reset(new ElectronLimiter(theParSet));
    m_tpmod->theElectronLimiter.get()->SetRangeCheckFlag(rLimiter);
    m_tpmod->theElectronLimiter.get()->SetFieldCheckFlag(eLimiter);
    ph->RegisterProcess(m_tpmod->theElectronLimiter.get(), G4Electron::Electron());
  }
  
  if(pLimiter){
    if(!m_tpmod) { m_tpmod = new TLSmod; }
    m_tpmod->thePositronLimiter.reset(new ElectronLimiter(theParSet));
    m_tpmod->thePositronLimiter.get()->SetFieldCheckFlag(pLimiter);
    ph->RegisterProcess(m_tpmod->theElectronLimiter.get(), G4Positron::Positron());
  }
  // enable fluorescence
  bool fluo = theParSet.getParameter<bool>("FlagFluo");
  if(fluo && !G4LossTableManager::Instance()->AtomDeexcitation()) {
    G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
    G4LossTableManager::Instance()->SetAtomDeexcitation(de);
  }
  // change parameters of transportation
  bool modifyT = theParSet.getParameter<bool>("ModifyTransportation");
  if(modifyT) {
    double th1 = theParSet.getUntrackedParameter<double>("ThresholdWarningEnergy")*MeV;
    double th2 = theParSet.getUntrackedParameter<double>("ThresholdImportantEnergy")*MeV;
    int nt = theParSet.getUntrackedParameter<int>("ThresholdTrials");
    ModifyTransportation(G4Electron::Electron(), nt, th1, th2);
    ModifyTransportation(G4Positron::Positron(), nt, th1, th2);
    ModifyTransportation(G4Proton::Proton(), nt, th1, th2);
  }
}

void ParametrisedEMPhysics::ModifyTransportation ( const G4ParticleDefinition *  part, int  ntry, double  th1, double  th2  )

private

Definition at line 254 of file ParametrisedEMPhysics.cc.

Referenced by ConstructProcess().

{
  G4ProcessManager* man = part->GetProcessManager();
  G4Transportation* trans = (G4Transportation*)((*(man->GetProcessList()))[0]);
  if(trans) {
    trans->SetThresholdWarningEnergy(th1); 
    trans->SetThresholdImportantEnergy(th2); 
    trans->SetThresholdTrials(ntry); 
    edm::LogVerbatim("SimG4CoreApplication") 
      << "ParametrisedEMPhysics: printout level changed for " << part->GetParticleName();
 
  }
}

Member Data Documentation

G4ThreadLocal ParametrisedEMPhysics::TLSmod * ParametrisedEMPhysics::m_tpmod = 0

staticprivate

Definition at line 33 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess(), and ~ParametrisedEMPhysics().

edm::ParameterSet ParametrisedEMPhysics::theParSet

private

Definition at line 31 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess(), and ParametrisedEMPhysics().