ParametrisedPhysics.cc

Go to the documentation of this file.

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

#include "G4Electron.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;

G4ThreadLocal ParametrisedPhysics::ThreadPrivate* ParametrisedPhysics::tpdata = nullptr;

ParametrisedPhysics::ParametrisedPhysics(std::string name, const edm::ParameterSet & p) :
  G4VPhysicsConstructor(name), theParSet(p) 
{}

ParametrisedPhysics::~ParametrisedPhysics() {
  if(nullptr != tpdata) {
    delete tpdata->theEMShowerModel;
    delete tpdata->theHadShowerModel;
    delete tpdata->theHadronShowerModel;
    delete tpdata->theFastSimulationManagerProcess;
    tpdata = nullptr;
  }
}

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() {

  tpdata = new ThreadPrivate;
  tpdata->theEMShowerModel = nullptr;
  tpdata->theHadShowerModel = nullptr;
  tpdata->theHadronShowerModel = nullptr;
  tpdata->theFastSimulationManagerProcess = nullptr;

  bool gem  = theParSet.getParameter<bool>("GflashEcal");
  bool ghad = theParSet.getParameter<bool>("GflashHcal");
  G4cout << "GFlash Construct: " << gem << "  " << ghad << G4endl;

  if(gem || ghad) {
    tpdata->theFastSimulationManagerProcess =
      new G4FastSimulationManagerProcess();

    G4ParticleTable* table = G4ParticleTable::GetParticleTable();
    EmParticleList emList;
    for(const auto& particleName : emList.PartNames()) {
      G4ParticleDefinition* particle = table->FindParticle(particleName);
      G4ProcessManager * pmanager = particle->GetProcessManager();
      const G4String& pname = particle->GetParticleName();
      if(pname == "e-" || pname == "e+") {
    pmanager->AddDiscreteProcess(tpdata->theFastSimulationManagerProcess);
      }
    }

    if(gem) {
      G4Region* aRegion = G4RegionStore::GetInstance()->GetRegion("EcalRegion");

      if(!aRegion){
    G4cout << "EcalRegion is not defined !!!" << G4endl;
    G4cout << "This means that GFlash will not be turned on." << G4endl;
    
      } else {

    //Electromagnetic Shower Model for ECAL
    tpdata->theEMShowerModel = 
      new GflashEMShowerModel("GflashEMShowerModel",aRegion,theParSet);
    G4cout << "GFlash is defined for EcalRegion" << G4endl;
      }    
    }
    if(ghad) {
      G4Region* aRegion = G4RegionStore::GetInstance()->GetRegion("HcalRegion");
      if(!aRegion) {
    G4cout << "HcalRegion is not defined !!!" << G4endl;
    G4cout << "This means that GFlash will not be turned on." << G4endl;
    
      } else {

    //Electromagnetic Shower Model for HCAL
    tpdata->theHadShowerModel = 
      new GflashEMShowerModel("GflashHadShowerModel",aRegion,theParSet);
    G4cout << "GFlash is defined for HcalRegion" << G4endl;    
      }
    }
  }
}