CustomPhysicsList Class Reference

#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::CustomPhysicsList	(	const std::string &	name,
		const edm::ParameterSet &	p,
		bool	useuni = false
	)

Definition at line 33 of file CustomPhysicsList.cc.

References dfactor, fHadronicInteraction, personalPlayback::fp, fParticleFactory, 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 = std::make_unique<CustomParticleFactory>();
  myHelper.reset(nullptr);

  edm::LogVerbatim("SimG4CoreCustomPhysics") << "CustomPhysicsList: Path for custom particle definition file: \n"
                                             << particleDefFilePath << "\n"
                                             << "      dark_factor= " << dfactor;
}

~CustomPhysicsList()

CustomPhysicsList::~CustomPhysicsList ( )

override

Definition at line 54 of file CustomPhysicsList.cc.

{}

Member Function Documentation

ConstructParticle()

void CustomPhysicsList::ConstructParticle ( )

override

Definition at line 56 of file CustomPhysicsList.cc.

References fParticleFactory, and particleDefFilePath.

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

ConstructProcess()

void CustomPhysicsList::ConstructProcess ( )

override

Definition at line 61 of file CustomPhysicsList.cc.

References CMSQGSPSIMPBuilder::Build(), dfactor, fHadronicInteraction, fParticleFactory, 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()) {
    if (particle->GetParticleType() == "simp") {
      G4ProcessManager* pmanager = particle->GetProcessManager();
      if (pmanager) {
        CMSSIMPInelasticProcess* simpInelPr = new CMSSIMPInelasticProcess();
        CMSQGSPSIMPBuilder* theQGSPSIMPB = new CMSQGSPSIMPBuilder();
        theQGSPSIMPB->Build(simpInelPr);
        pmanager->AddDiscreteProcess(simpInelPr);
      } else
        edm::LogVerbatim("CustomPhysics") << "   No pmanager";
    }

    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 = std::make_unique<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";
        }
      }
    }
  }
}

Member Data Documentation

dfactor

double CustomPhysicsList::dfactor

private

Definition at line 30 of file CustomPhysicsList.h.

Referenced by ConstructProcess(), and CustomPhysicsList().

fHadronicInteraction

bool CustomPhysicsList::fHadronicInteraction

private

Definition at line 24 of file CustomPhysicsList.h.

Referenced by ConstructProcess(), and CustomPhysicsList().

fParticleFactory

std::unique_ptr<CustomParticleFactory> CustomPhysicsList::fParticleFactory

private

Definition at line 22 of file CustomPhysicsList.h.

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

myConfig

edm::ParameterSet CustomPhysicsList::myConfig

private

Definition at line 26 of file CustomPhysicsList.h.

Referenced by ConstructProcess(), and CustomPhysicsList().

myHelper

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

staticprivate

Definition at line 21 of file CustomPhysicsList.h.

Referenced by ConstructProcess(), and CustomPhysicsList().

particleDefFilePath

std::string CustomPhysicsList::particleDefFilePath

private

Definition at line 28 of file CustomPhysicsList.h.

Referenced by ConstructParticle(), and CustomPhysicsList().

processDefFilePath

std::string CustomPhysicsList::processDefFilePath

private

Definition at line 29 of file CustomPhysicsList.h.