CustomPhysicsListSS Class Reference

#include <CustomPhysicsListSS.h>

Inheritance diagram for CustomPhysicsListSS:

Public Member Functions
void	ConstructParticle () override
void	ConstructProcess () override
	CustomPhysicsListSS (const std::string &name, const edm::ParameterSet &p, bool useuni=false)
	~CustomPhysicsListSS () 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 12 of file CustomPhysicsListSS.h.

Constructor & Destructor Documentation

CustomPhysicsListSS::CustomPhysicsListSS	(	const std::string &	name,
		const edm::ParameterSet &	p,
		bool	useuni = false
	)

Definition at line 26 of file CustomPhysicsListSS.cc.

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

  :  G4VPhysicsConstructor(name) 
{  
  myConfig = p;
  if(apinew) {
    dfactor = p.getParameter<double>("DarkMPFactor");
    fHadronicInteraction = p.getParameter<bool>("RhadronPhysics");
  } else {
    // this is left for backward compatibility
    dfactor = p.getParameter<double>("dark_factor");
    fHadronicInteraction = p.getParameter<bool>("rhadronPhysics");
  }
  edm::FileInPath fp = p.getParameter<edm::FileInPath>("particlesDef");
  particleDefFilePath = fp.fullPath();
  fParticleFactory.reset(new CustomParticleFactory());
  fDecayProcess.reset(nullptr);
  myHelper.reset(nullptr);

  edm::LogInfo("SimG4CoreCustomPhysics")
    <<"CustomPhysicsListSS: Path for custom particle definition file: \n" 
    <<particleDefFilePath;
}

CustomPhysicsListSS::~CustomPhysicsListSS ( )

override

Definition at line 50 of file CustomPhysicsListSS.cc.

                                          {
}

Member Function Documentation

void CustomPhysicsListSS::ConstructParticle ( )

override

Definition at line 53 of file CustomPhysicsListSS.cc.

References fParticleFactory, and particleDefFilePath.

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

void CustomPhysicsListSS::ConstructProcess ( )

override

Definition at line 59 of file CustomPhysicsListSS.cc.

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

                                           {

  edm::LogInfo("SimG4CoreCustomPhysicsSS") 
    <<"CustomPhysicsListSS: adding CustomPhysics processes";

  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") 
    <<"CustomPhysicsListSS: " << 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") 
        <<"CustomPhysicsListSS: " << 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);
        }
      }
    }
  }
}

Member Data Documentation

double CustomPhysicsListSS::dfactor

private

Definition at line 35 of file CustomPhysicsListSS.h.

Referenced by ConstructProcess(), and CustomPhysicsListSS().

G4ThreadLocal std::unique_ptr< G4Decay > CustomPhysicsListSS::fDecayProcess

staticprivate

Definition at line 24 of file CustomPhysicsListSS.h.

Referenced by ConstructProcess(), and CustomPhysicsListSS().

bool CustomPhysicsListSS::fHadronicInteraction

private

Definition at line 29 of file CustomPhysicsListSS.h.

Referenced by ConstructProcess(), and CustomPhysicsListSS().

std::unique_ptr<CustomParticleFactory> CustomPhysicsListSS::fParticleFactory

private

Definition at line 27 of file CustomPhysicsListSS.h.

Referenced by ConstructParticle(), ConstructProcess(), and CustomPhysicsListSS().

edm::ParameterSet CustomPhysicsListSS::myConfig

private

Definition at line 31 of file CustomPhysicsListSS.h.

Referenced by ConstructProcess(), and CustomPhysicsListSS().

G4ThreadLocal std::unique_ptr< G4ProcessHelper > CustomPhysicsListSS::myHelper

staticprivate

Definition at line 25 of file CustomPhysicsListSS.h.

Referenced by ConstructProcess(), and CustomPhysicsListSS().

std::string CustomPhysicsListSS::particleDefFilePath

private

Definition at line 33 of file CustomPhysicsListSS.h.

Referenced by ConstructParticle(), and CustomPhysicsListSS().

std::string CustomPhysicsListSS::processDefFilePath

private

Definition at line 34 of file CustomPhysicsListSS.h.