#include <CustomPhysicsListSS.h>

Inheritance diagram for CustomPhysicsListSS:

Public Member Functions
virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

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

virtual	~CustomPhysicsListSS ()

Private Attributes
bool	fHadronicInteraction

edm::ParameterSet	myConfig

std::string	particleDefFilePath

std::string	processDefFilePath

Static Private Attributes
static G4ThreadLocal G4Decay *	fDecayProcess = 0

static G4ThreadLocal bool	fInitialized = false

static G4ThreadLocal G4ProcessHelper *	myHelper = 0

Detailed Description

Definition at line 15 of file CustomPhysicsListSS.h.

Constructor & Destructor Documentation

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

Definition at line 26 of file CustomPhysicsListSS.cc.

References fHadronicInteraction, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), myConfig, myHelper, AlCaHLTBitMon_ParallelJobs::p, and particleDefFilePath.

   :  G4VPhysicsConstructor(name) 
 {  
   myConfig = p;
   edm::FileInPath fp = p.getParameter<edm::FileInPath>("particlesDef");
   fHadronicInteraction = p.getParameter<bool>("rhadronPhysics");
   particleDefFilePath = fp.fullPath();
   edm::LogInfo("SimG4CoreCustomPhysics")
     <<"CustomPhysicsListSS: Path for custom particle definition file: \n" 
     <<particleDefFilePath;
   myHelper = 0;  
 }

CustomPhysicsListSS::~CustomPhysicsListSS ( )

virtual

Definition at line 39 of file CustomPhysicsListSS.cc.

39 {

40 }

Member Function Documentation

void CustomPhysicsListSS::ConstructParticle ( )

virtual

Definition at line 42 of file CustomPhysicsListSS.cc.

References CustomParticleFactory::loadCustomParticles(), and particleDefFilePath.

                                            {
   CustomParticleFactory::loadCustomParticles(particleDefFilePath);     
 }

void CustomPhysicsListSS::ConstructProcess ( )

virtual

Definition at line 46 of file CustomPhysicsListSS.cc.

References CommonMethods::cp(), fDecayProcess, fHadronicInteraction, fInitialized, CustomParticle::GetCloud(), CustomParticle::GetSpectator(), GeV, CustomParticleFactory::isCustomParticle(), myConfig, myHelper, and CustomPDGParser::s_isRHadron().

                                            {
 
   if(fInitialized) { return; }
   fInitialized = true;
 
   edm::LogInfo("SimG4CoreCustomPhysics") 
     <<"CustomPhysicsListSS: adding CustomPhysics processes";
 
   fDecayProcess = new G4Decay();
 
   aParticleIterator->reset();
 
   while((*aParticleIterator)()) {
     G4ParticleDefinition* particle = aParticleIterator->value();
     if(CustomParticleFactory::isCustomParticle(particle)) {
       CustomParticle* cp = dynamic_cast<CustomParticle*>(particle);
       G4ProcessManager* pmanager = particle->GetProcessManager();
       edm::LogInfo("SimG4CoreCustomPhysics") 
     <<"CustomPhysicsListSS: " << particle->GetParticleName()
     <<" PDGcode= " << particle->GetPDGEncoding()
     << " Mass= " << particle->GetPDGMass()/GeV  <<" GeV.";
       if(cp && pmanager) {
     if(particle->GetPDGCharge()/eplus != 0) {
       pmanager->AddProcess(new G4CoulombScattering,  -1,-1, 1);
       pmanager->AddProcess(new G4hIonisation,        -1, 1, 2);
     }
     if(fDecayProcess->IsApplicable(*particle)) {
       pmanager->AddProcess(new G4Decay, 0, -1, 3);
     }
     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) myHelper = new G4ProcessHelper(myConfig);
       pmanager->AddDiscreteProcess(new FullModelHadronicProcess(myHelper));
     }
       }
     }
   }
 }