#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
ElectronLimiter *	theElectronLimiter

GFlashEMShowerModel *	theEMShowerModel

GFlashEMShowerModel *	theHadShowerModel

edm::ParameterSet	theParSet

ElectronLimiter *	thePositronLimiter

Detailed Description

Definition at line 17 of file ParametrisedEMPhysics.h.

Constructor & Destructor Documentation

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

Definition at line 29 of file ParametrisedEMPhysics.cc.

References theEMShowerModel, and theHadShowerModel.

   : G4VPhysicsConstructor(name), theParSet(p) 
 {
   theEMShowerModel = 0;
   theHadShowerModel = 0;
 }

ParametrisedEMPhysics::~ParametrisedEMPhysics ( )

virtual

Definition at line 36 of file ParametrisedEMPhysics.cc.

References theEMShowerModel, and theHadShowerModel.

                                               {
   delete theEMShowerModel;
   delete theHadShowerModel;
 }

Member Function Documentation

void ParametrisedEMPhysics::ConstructParticle ( )

virtual

Definition at line 41 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 59 of file ParametrisedEMPhysics.cc.

References edm::ParameterSet::getParameter(), i, ElectronLimiter::SetFieldCheckFlag(), ElectronLimiter::SetRangeCheckFlag(), theElectronLimiter, theEMShowerModel, theHadShowerModel, theParSet, and thePositronLimiter.

                                              {
 
   // GFlash part 
   bool gem  = theParSet.getParameter<bool>("GflashEcal");
   bool ghad = theParSet.getParameter<bool>("GflashHcal");
 
   if(gem || ghad) {
     edm::LogInfo("SimG4CoreApplication") 
       << "ParametrisedEMPhysics: GFlash Construct: " 
       << gem << "  " << ghad;
     G4FastSimulationManagerProcess * theFastSimulationManagerProcess = 
       new G4FastSimulationManagerProcess();
     theParticleIterator->reset();
     while ((*theParticleIterator)()) {
       G4ParticleDefinition * particle = theParticleIterator->value();
       G4ProcessManager * pmanager = particle->GetProcessManager();
       G4String pname = particle->GetParticleName();
       if(pname == "e-" || pname == "e+") {
     pmanager->AddProcess(theFastSimulationManagerProcess, -1, -1, 1);
       }
     }
 
     if(gem) {
       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 {
 
     //Electromagnetic Shower Model for ECAL
     theEMShowerModel = 
       new GFlashEMShowerModel("GflashEMShowerModel",aRegion,theParSet);
     //std::cout << "GFlash is defined for EcalRegion" << std::endl;
       }    
     }
     if(ghad) {
       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 {
 
     //Electromagnetic Shower Model for HCAL
     theHadShowerModel = 
       new GFlashEMShowerModel("GflashHadShowerModel",aRegion,theParSet);
     //std::cout << "GFlash is defined for HcalRegion" << std::endl;
       }
     }
   }
   // Russian Roulette part 
   G4EmProcessOptions opt;
   const G4int NREG = 6; 
   const G4String rname[NREG] = {"EcalRegion", "HcalRegion", "MuonIron",
                 "PreshowerRegion","CastorRegion",
                 "DefaultRegionForTheWorld"};
   G4double rrfact[NREG] = { 1.0 };
 
   // Russian roulette for gamma
   G4double energyLim = theParSet.getParameter<double>("RusRoGammaEnergyLimit")*MeV;
   if(energyLim > 0.0) {
     rrfact[0] = theParSet.getParameter<double>("RusRoEcalGamma");
     rrfact[1] = theParSet.getParameter<double>("RusRoHcalGamma");
     rrfact[2] = theParSet.getParameter<double>("RusRoMuonIronGamma");
     rrfact[3] = theParSet.getParameter<double>("RusRoPreShowerGamma");
     rrfact[4] = theParSet.getParameter<double>("RusRoCastorGamma");
     rrfact[5] = theParSet.getParameter<double>("RusRoWorldGamma");
     for(int i=0; i<NREG; ++i) {
       if(rrfact[i] < 1.0) {
     opt.ActivateSecondaryBiasing("eBrem",rname[i],rrfact[i],energyLim);
     edm::LogInfo("SimG4CoreApplication") 
       << "ParametrisedEMPhysics: Russian Roulette"
       << " for gamma Prob= " << rrfact[i]  
       << " Elimit(MeV)= " << energyLim/CLHEP::MeV
       << " inside " << rname[i];
       }
     }
   }
   // Russian roulette for e-
   energyLim = theParSet.getParameter<double>("RusRoElectronEnergyLimit")*MeV;
   if(energyLim > 0.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) {
     opt.ActivateSecondaryBiasing("eIoni",rname[i],rrfact[i],energyLim);
     opt.ActivateSecondaryBiasing("hIoni",rname[i],rrfact[i],energyLim);
     //opt.ActivateSecondaryBiasing("muIoni",rname[i],rrfact[i],energyLim);
     //opt.ActivateSecondaryBiasing("ionIoni",rname[i],rrfact[i],energyLim);
     //opt.ActivateSecondaryBiasingForGamma("phot",rname[i],rrfact[i],energyLim);
     //opt.ActivateSecondaryBiasingForGamma("compt",rname[i],rrfact[i],energyLim);
     //opt.ActivateSecondaryBiasingForGamma("conv",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 ||  pLimiter) {
     G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
     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());
     }
   }
 }