ParametrisedEMPhysics Class Reference

#include <ParametrisedEMPhysics.h>

Inheritance diagram for ParametrisedEMPhysics:

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

Private Attributes
GFlashEMShowerModel *	theEcalEMShowerModel
GFlashHadronShowerModel *	theEcalHadShowerModel
ElectronLimiter *	theElectronLimiter
GFlashEMShowerModel *	theHcalEMShowerModel
GFlashHadronShowerModel *	theHcalHadShowerModel
edm::ParameterSet	theParSet
ElectronLimiter *	thePositronLimiter

Detailed Description

Definition at line 18 of file ParametrisedEMPhysics.h.

Constructor & Destructor Documentation

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

Definition at line 41 of file ParametrisedEMPhysics.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GeV, theEcalEMShowerModel, theEcalHadShowerModel, theHcalEMShowerModel, theHcalHadShowerModel, and theParSet.

  : G4VPhysicsConstructor(name), theParSet(p) 
{
  theEcalEMShowerModel = nullptr;
  theEcalHadShowerModel = nullptr;
  theHcalEMShowerModel = nullptr;
  theHcalHadShowerModel = nullptr;

  // 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);

  edm::LogInfo("SimG4CoreApplication") 
    << "ParametrisedEMPhysics::ConstructProcess: bremsstrahlung threshold Eth= "
    << bremth/GeV << " GeV" 
    << "\n                                         verbosity= " << verb
    << "  fluoFlag: " << fluo; 
}

ParametrisedEMPhysics::~ParametrisedEMPhysics ( )

virtual

Definition at line 68 of file ParametrisedEMPhysics.cc.

References theEcalEMShowerModel, theEcalHadShowerModel, theHcalEMShowerModel, and theHcalHadShowerModel.

                                              {
  delete theEcalEMShowerModel;
  delete theEcalHadShowerModel;
  delete theHcalEMShowerModel;
  delete theHcalHadShowerModel;
}

Member Function Documentation

void ParametrisedEMPhysics::ConstructParticle ( )

virtual

Definition at line 75 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 ( )

virtual

Definition at line 93 of file ParametrisedEMPhysics.cc.

References configurableAnalysis::Electron, edm::ParameterSet::getParameter(), mps_fire::i, MeV, ElectronLimiter::SetFieldCheckFlag(), ElectronLimiter::SetRangeCheckFlag(), theEcalEMShowerModel, theEcalHadShowerModel, theElectronLimiter, theHcalEMShowerModel, theHcalHadShowerModel, theParSet, and 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) {
    edm::LogInfo("SimG4CoreApplication") 
      << "ParametrisedEMPhysics: GFlash Construct for e+-: " 
      << gem << "  " << ghad << " for hadrons: " << gemHad << "  " << ghadHad;

    G4FastSimulationManagerProcess * theFastSimulationManagerProcess = 
      new G4FastSimulationManagerProcess();

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

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

      //Electromagnetic Shower Model for ECAL
      theEcalEMShowerModel = 
        new GFlashEMShowerModel("GflashEcalEMShowerModel",aRegion,theParSet);
      //std::cout << "GFlash is defined for EcalRegion" << std::endl;
    }
    if(gemHad) {

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

      //Electromagnetic Shower Model for HCAL
      theHcalEMShowerModel = 
        new GFlashEMShowerModel("GflashHcalEMShowerModel",aRegion,theParSet);
      //std::cout << "GFlash is defined for HcalRegion" << std::endl;
    }
    if(ghadHad) {

      //Electromagnetic Shower Model for ECAL
      theHcalHadShowerModel = 
        new GFlashHadronShowerModel("GflashHcalHadShowerModel",aRegion,theParSet);
      //std::cout << "GFlash is defined for EcalRegion" << std::endl;
    }    
      }
    }
  }
  // Russian roulette and tracking cut for e+-
  const G4int NREG = 6; 
  const G4String rname[NREG] = {"EcalRegion", "HcalRegion", "MuonIron",
                "PreshowerRegion","CastorRegion",
                "DefaultRegionForTheWorld"};
  G4double rrfact[NREG] = { 1.0 };

  double energyLim = 
    theParSet.getParameter<double>("RusRoElectronEnergyLimit")*MeV;
  if(energyLim > 0.0) {
    G4EmProcessOptions opt;
    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) {
    opt.ActivateSecondaryBiasing("eIoni",rname[i],rrfact[i],energyLim);
    opt.ActivateSecondaryBiasing("hIoni",rname[i],rrfact[i],energyLim);
    edm::LogInfo("SimG4CoreApplication") 
      << "ParametrisedEMPhysics: Russian Roulette"
      << " for e- Prob= " << rrfact[i]  
      << " Elimit(MeV)= " << energyLim/CLHEP::MeV
      << " inside " << rname[i];
      }
    }
  }

  // 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) {
    theElectronLimiter = new ElectronLimiter(theParSet);
    theElectronLimiter->SetRangeCheckFlag(rLimiter);
    theElectronLimiter->SetFieldCheckFlag(eLimiter);
    ph->RegisterProcess(theElectronLimiter, G4Electron::Electron());
  }
  
  if(pLimiter){
    thePositronLimiter = new ElectronLimiter(theParSet);
    thePositronLimiter->SetFieldCheckFlag(pLimiter);
    ph->RegisterProcess(theElectronLimiter, G4Positron::Positron());
  }
  // enable fluorescence
  bool fluo = theParSet.getParameter<bool>("FlagFluo");
  if(fluo) {
    G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
    G4LossTableManager::Instance()->SetAtomDeexcitation(de);
    de->SetFluo(true);
  }
}

Member Data Documentation

GFlashEMShowerModel* ParametrisedEMPhysics::theEcalEMShowerModel

private

Definition at line 32 of file ParametrisedEMPhysics.h.

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

GFlashHadronShowerModel* ParametrisedEMPhysics::theEcalHadShowerModel

private

Definition at line 34 of file ParametrisedEMPhysics.h.

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

ElectronLimiter* ParametrisedEMPhysics::theElectronLimiter

private

Definition at line 37 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess().

GFlashEMShowerModel* ParametrisedEMPhysics::theHcalEMShowerModel

private

Definition at line 33 of file ParametrisedEMPhysics.h.

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

GFlashHadronShowerModel* ParametrisedEMPhysics::theHcalHadShowerModel

private

Definition at line 35 of file ParametrisedEMPhysics.h.

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

edm::ParameterSet ParametrisedEMPhysics::theParSet

private

Definition at line 30 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess(), and ParametrisedEMPhysics().

ElectronLimiter* ParametrisedEMPhysics::thePositronLimiter

private

Definition at line 38 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess().