ParametrisedEMPhysics Class Reference

#include <ParametrisedEMPhysics.h>

List of all members.

Public Member Functions
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
	ParametrisedEMPhysics (std::string name, const edm::ParameterSet &p)
virtual	~ParametrisedEMPhysics ()
Private Attributes
GFlashEMShowerModel *	theEMShowerModel
GFlashEMShowerModel *	theHadShowerModel
edm::ParameterSet	theParSet

---

Detailed Description

Definition at line 16 of file ParametrisedEMPhysics.h.

---

Constructor & Destructor Documentation

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

Definition at line 24 of file ParametrisedEMPhysics.cc.

References theEMShowerModel, and theHadShowerModel.

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

ParametrisedEMPhysics::~ParametrisedEMPhysics

(

)

[virtual]

Definition at line 30 of file ParametrisedEMPhysics.cc.

References theEMShowerModel, and theHadShowerModel.

                                              {
  delete theEMShowerModel;
  delete theHadShowerModel;
}

---

Member Function Documentation

void ParametrisedEMPhysics::ConstructParticle

(

)

[virtual]

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

References edm::ParameterSet::getParameter(), theEMShowerModel, theHadShowerModel, and theParSet.

                                             {

  // 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;
  double gamEcal = theParSet.getParameter<double>("RusRoEcalGamma");
  if(gamEcal < 1.0) {
    double gamEcalLim = theParSet.getParameter<double>("RusRoEcalGammaLimit")
      *CLHEP::MeV;
    if(gamEcalLim > 0.0) {
      opt.ActivateSecondaryBiasing("eBrem","EcalRegion",gamEcal,gamEcalLim);
      edm::LogInfo("SimG4CoreApplication") 
        << "ParametrisedEMPhysics: Russian Roulette "
        << " for gamma in ECAL Prob= " 
        << gamEcal << "  Elimit(MeV)= " << gamEcalLim/CLHEP::MeV;
    }
  }
  double gamHcal = theParSet.getParameter<double>("RusRoHcalGamma");
  if(gamHcal < 1.0) {
    double gamHcalLim = theParSet.getParameter<double>("RusRoHcalGammaLimit")
      *CLHEP::MeV;
    if(gamHcalLim > 0.0) {
      opt.ActivateSecondaryBiasing("eBrem","HcalRegion",gamHcal,gamHcalLim);
      edm::LogInfo("SimG4CoreApplication") 
        << "ParametrisedEMPhysics: Russian Roulette "
        << " for gamma in HCAL Prob= " 
        << gamHcal << "  Elimit(MeV)= " << gamHcalLim/CLHEP::MeV;
    }
  }
  double eEcal = theParSet.getParameter<double>("RusRoEcalElectron");
  if(eEcal < 1.0) {
    double eEcalLim = theParSet.getParameter<double>("RusRoEcalElectronLimit")
      *CLHEP::MeV;
    if(eEcalLim > 0.0) {
      opt.ActivateSecondaryBiasing("eIoni","EcalRegion",eEcal,eEcalLim);
      opt.ActivateSecondaryBiasing("hIoni","EcalRegion",eEcal,eEcalLim);
      opt.ActivateSecondaryBiasing("muIoni","EcalRegion",eEcal,eEcalLim);
      opt.ActivateSecondaryBiasing("ionIoni","EcalRegion",eEcal,eEcalLim);
      // opt.ActivateSecondaryBiasingForGamma("phot","EcalRegion",eEcal,eEcalLim);
      // opt.ActivateSecondaryBiasingForGamma("compt","EcalRegion",eEcal,eEcalLim);
      // opt.ActivateSecondaryBiasingForGamma("conv","EcalRegion",eEcal,eEcalLim);
      edm::LogInfo("SimG4CoreApplication") 
        << "ParametrisedEMPhysics: Russian Roulette "
        << " for electrons in ECAL Prob= " 
        << eEcal << "  Elimit(MeV)= " << eEcalLim/CLHEP::MeV;
    }
  }
  double eHcal = theParSet.getParameter<double>("RusRoHcalElectron");
  if(eHcal < 1.0) {
    double eHcalLim = theParSet.getParameter<double>("RusRoHcalElectronLimit")
      *CLHEP::MeV;
    if(eHcalLim > 0.0) {
      opt.ActivateSecondaryBiasing("eIoni","HcalRegion",eHcal,eHcalLim);
      opt.ActivateSecondaryBiasing("hIoni","HcalRegion",eHcal,eHcalLim);
      opt.ActivateSecondaryBiasing("muIoni","HcalRegion",eHcal,eHcalLim);
      opt.ActivateSecondaryBiasing("ionIoni","HcalRegion",eHcal,eHcalLim);
      //opt.ActivateSecondaryBiasingForGamma("phot","HcalRegion",eHcal,eHcalLim);
      // opt.ActivateSecondaryBiasingForGamma("compt","HcalRegion",eHcal,eHcalLim);
      // opt.ActivateSecondaryBiasingForGamma("conv","HcalRegion",eHcal,eHcalLim);
      edm::LogInfo("SimG4CoreApplication") 
        << "ParametrisedEMPhysics: Russian Roulette "
        << " for electrons in HCAL Prob= " 
        << eHcal << "  Elimit(MeV)= " << eHcalLim/CLHEP::MeV;
    }
  }
}

---

Member Data Documentation

GFlashEMShowerModel* ParametrisedEMPhysics::theEMShowerModel [private]

Definition at line 30 of file ParametrisedEMPhysics.h.

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

GFlashEMShowerModel* ParametrisedEMPhysics::theHadShowerModel [private]

Definition at line 31 of file ParametrisedEMPhysics.h.

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

edm::ParameterSet ParametrisedEMPhysics::theParSet [private]

Definition at line 28 of file ParametrisedEMPhysics.h.

Referenced by ConstructProcess().