ParametrisedPhysics.cc

Go to the documentation of this file.

#include "SimG4Core/GFlash/interface/ParametrisedPhysics.h"

#include "G4Electron.hh"
#include "G4FastSimulationManagerProcess.hh"
#include "G4ProcessManager.hh"

#include "G4LeptonConstructor.hh"
#include "G4MesonConstructor.hh"
#include "G4BaryonConstructor.hh"
#include "G4ShortLivedConstructor.hh"
#include "G4IonConstructor.hh"
#include "G4RegionStore.hh"

using namespace CLHEP;


ParametrisedPhysics::ParametrisedPhysics(std::string name, const edm::ParameterSet & p) :
  G4VPhysicsConstructor(name), theParSet(p) 
{
  theEMShowerModel = 0;
  theHadShowerModel = 0;
  theHadronShowerModel = 0;
}

ParametrisedPhysics::~ParametrisedPhysics() {
  delete theEMShowerModel;
  delete theHadShowerModel;
  delete theHadronShowerModel;
}

void ParametrisedPhysics::ConstructParticle() 
{
    G4LeptonConstructor pLeptonConstructor;
    pLeptonConstructor.ConstructParticle();

    G4MesonConstructor pMesonConstructor;
    pMesonConstructor.ConstructParticle();

    G4BaryonConstructor pBaryonConstructor;
    pBaryonConstructor.ConstructParticle();

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

void ParametrisedPhysics::ConstructProcess() {

  bool gem  = theParSet.getParameter<bool>("GflashEcal");
  bool ghad = theParSet.getParameter<bool>("GflashHcal");
  std::cout << "GFlash Construct: " << gem << "  " << ghad << std::endl;

  if(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){
    std::cout << "EcalRegion is not defined !!!" << std::endl;
    std::cout << "This means that GFlash will not be turned on." << std::endl;
    
      } 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) {
    std::cout << "HcalRegion is not defined !!!" << std::endl;
    std::cout << "This means that GFlash will not be turned on." << std::endl;
    
      } else {

    //Electromagnetic Shower Model for HCAL
    theHadShowerModel = 
      new GflashEMShowerModel("GflashHadShowerModel",aRegion,theParSet);
    std::cout << "GFlash is defined for HcalRegion" << std::endl;    
      }
    }
  }
}