#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, bool useuni=false)

	~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< G4ProcessHelper >	myHelper

Detailed Description

Definition at line 12 of file CustomPhysicsList.h.

Constructor & Destructor Documentation

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

Definition at line 30 of file CustomPhysicsList.cc.

References dfactor, 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());
   myHelper.reset(nullptr);
 
   edm::LogVerbatim("SimG4CoreCustomPhysics") << "CustomPhysicsList: Path for custom particle definition file: \n"
                                              << particleDefFilePath << "\n"
                                              << "      dark_factor= " << dfactor;
 }

CustomPhysicsList::~CustomPhysicsList ( )

override

Definition at line 51 of file CustomPhysicsList.cc.

51 {}

Member Function Documentation

void CustomPhysicsList::ConstructParticle ( )

override

Definition at line 53 of file CustomPhysicsList.cc.

References fParticleFactory, and particleDefFilePath.

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

void CustomPhysicsList::ConstructProcess ( )

override

Definition at line 58 of file CustomPhysicsList.cc.

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

                                          {
   edm::LogVerbatim("SimG4CoreCustomPhysics") << "CustomPhysicsList: adding CustomPhysics processes "
                                              << "for the list of particles";
 
   G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
   for (auto particle : fParticleFactory.get()->GetCustomParticles()) {
     CustomParticle* cp = dynamic_cast<CustomParticle*>(particle);
     if (cp) {
       G4ProcessManager* pmanager = particle->GetProcessManager();
       edm::LogVerbatim("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 (cp->GetCloud() && fHadronicInteraction &&
             (CustomPDGParser::s_isgluinoHadron(particle->GetPDGEncoding()) ||
              (CustomPDGParser::s_isstopHadron(particle->GetPDGEncoding())) ||
              (CustomPDGParser::s_issbottomHadron(particle->GetPDGEncoding())))) {
           edm::LogVerbatim("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);
         }
         if (particle->GetParticleName() == "anti_sexaq") {
           // here the different sexaquark interactions get defined
           G4HadronicProcess* sqInelPr = new G4HadronicProcess();
           CMSSQNeutronAnnih* sqModel = new CMSSQNeutronAnnih(particle->GetPDGMass() / GeV);
           sqInelPr->RegisterMe(sqModel);
           CMSSQInelasticCrossSection* sqInelXS = new CMSSQInelasticCrossSection(particle->GetPDGMass() / GeV);
           sqInelPr->AddDataSet(sqInelXS);
           pmanager->AddDiscreteProcess(sqInelPr);
           // add also the looping needed to simulate flat interaction probability
           CMSSQLoopProcess* sqLoopPr = new CMSSQLoopProcess();
           pmanager->AddContinuousProcess(sqLoopPr);
           CMSSQLoopProcessDiscr* sqLoopPrDiscr = new CMSSQLoopProcessDiscr(particle->GetPDGMass() / GeV);
           pmanager->AddDiscreteProcess(sqLoopPrDiscr);
         } else if (particle->GetParticleName() == "sexaq") {
           edm::LogVerbatim("CustomPhysics") << "   No pmanager implemented for sexaq, only for anti_sexaq";
         }
       }
     }
   }
 }