#include <CustomPhysicsList.h>

Inheritance diagram for CustomPhysicsList:

Public Member Functions
void	ConstructParticle () override

void	ConstructProcess () override

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

	~CustomPhysicsList () override

Private Attributes
double	dfactor

bool	fHadronicInteraction

std::unique_ptr< CustomParticleFactory >	fParticleFactory

edm::ParameterSet	myConfig

std::string	particleDefFilePath

std::string	processDefFilePath

Static Private Attributes
static G4ThreadLocal std::unique_ptr< G4Decay >	fDecayProcess

static G4ThreadLocal std::unique_ptr< G4ProcessHelper >	myHelper

Detailed Description

Definition at line 13 of file CustomPhysicsList.h.

Constructor & Destructor Documentation

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

Definition at line 25 of file CustomPhysicsList.cc.

References dfactor, fDecayProcess, fHadronicInteraction, fParticleFactory, edm::ParameterSet::getParameter(), myConfig, myHelper, AlCaHLTBitMon_ParallelJobs::p, and particleDefFilePath.

   :  G4VPhysicsConstructor(name) 
 {  
   myConfig = p;
   dfactor = p.getParameter<double>("dark_factor");
   edm::FileInPath fp = p.getParameter<edm::FileInPath>("particlesDef");
   fHadronicInteraction = p.getParameter<bool>("rhadronPhysics");
 
   particleDefFilePath = fp.fullPath();
 
   fParticleFactory.reset(new CustomParticleFactory());
   fDecayProcess.reset(nullptr);
   myHelper.reset(nullptr);
 
   edm::LogInfo("SimG4CoreCustomPhysics") 
     << "CustomPhysicsList: Path for custom particle definition file: \n"
     <<particleDefFilePath << "\n" << "      dark_factor= " << dfactor;
 }

CustomPhysicsList::~CustomPhysicsList ( )

override

Definition at line 44 of file CustomPhysicsList.cc.

44 {

45 }

Member Function Documentation

void CustomPhysicsList::ConstructParticle ( )

override

Definition at line 47 of file CustomPhysicsList.cc.

References fParticleFactory, and particleDefFilePath.

                                          {
   edm::LogInfo("SimG4CoreCustomPhysics") 
     << "===== CustomPhysicsList::ConstructParticle ";
   fParticleFactory.get()->loadCustomParticles(particleDefFilePath);
 }

void CustomPhysicsList::ConstructProcess ( )

override

Definition at line 53 of file CustomPhysicsList.cc.

References SimDataFormats::CaloAnalysis::cp, dfactor, fDecayProcess, fHadronicInteraction, fParticleFactory, CustomParticle::GetCloud(), CustomParticle::GetSpectator(), GeV, myConfig, myHelper, and CustomPDGParser::s_isRHadron().

                                          {
   
   edm::LogInfo("SimG4CoreCustomPhysics") 
     <<"CustomPhysicsList: adding CustomPhysics processes "
     << "for the list of particles";
 
   fDecayProcess.reset(new G4Decay());
   G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
   for(auto particle : fParticleFactory.get()->GetCustomParticles()) {
     CustomParticle* cp = dynamic_cast<CustomParticle*>(particle);
     if(cp) {
       G4ProcessManager* pmanager = particle->GetProcessManager();
       edm::LogInfo("SimG4CoreCustomPhysics") 
     <<"CustomPhysicsList: " << particle->GetParticleName()
     <<"  PDGcode= " << particle->GetPDGEncoding()
     << "  Mass= " << particle->GetPDGMass()/GeV  <<" GeV.";
       if(pmanager) {
     if(particle->GetPDGCharge() != 0.0) {
       ph->RegisterProcess(new G4hMultipleScattering, particle);
       ph->RegisterProcess(new G4hIonisation, particle);
     }
     if(fDecayProcess.get()->IsApplicable(*particle)) {
       ph->RegisterProcess(fDecayProcess.get(), particle);
     }
     if(cp->GetCloud() && fHadronicInteraction && 
        CustomPDGParser::s_isRHadron(particle->GetPDGEncoding())) {
       edm::LogInfo("SimG4CoreCustomPhysics") 
         <<"CustomPhysicsList: " << particle->GetParticleName()
         <<" CloudMass= " <<cp->GetCloud()->GetPDGMass()/GeV
         <<" GeV; SpectatorMass= " << cp->GetSpectator()->GetPDGMass()/GeV<<" GeV.";
        
           if(!myHelper.get()) { myHelper.reset(new G4ProcessHelper(myConfig, fParticleFactory.get())); }
       pmanager->AddDiscreteProcess(new FullModelHadronicProcess(myHelper.get()));
     }
         if(particle->GetParticleType()=="darkpho"){
           CMSDarkPairProductionProcess * darkGamma = new CMSDarkPairProductionProcess(dfactor);
           pmanager->AddDiscreteProcess(darkGamma);
         }
       }
     }
   }
 }