#include <CustomPhysicsList.h>

Inheritance diagram for CustomPhysicsList:

Public Member Functions
virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

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

virtual	~CustomPhysicsList ()

Protected Member Functions
void	addCustomPhysics ()

Private Member Functions
void	setupRHadronPhycis (G4ParticleDefinition *particle)

void	setupSUSYPhycis (G4ParticleDefinition *particle)

Private Attributes
double	dfactor

edm::ParameterSet	myConfig

G4ProcessHelper *	myHelper

std::string	particleDefFilePath

std::string	processDefFilePath

Detailed Description

Definition at line 14 of file CustomPhysicsList.h.

Constructor & Destructor Documentation

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

Definition at line 28 of file CustomPhysicsList.cc.

References dfactor, 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");
   particleDefFilePath = fp.fullPath();
   edm::LogInfo("CustomPhysics")<<"Path for custom particle definition file: "
                    <<particleDefFilePath<< "\n" << "      dark_factor= " << dfactor;
   myHelper = 0;  
 }

CustomPhysicsList::~CustomPhysicsList ( )

virtual

Definition at line 40 of file CustomPhysicsList.cc.

References myHelper.

                                       {
   delete myHelper;
 }

Member Function Documentation

void CustomPhysicsList::addCustomPhysics ( )

protected

Definition at line 52 of file CustomPhysicsList.cc.

References CommonMethods::cp(), dfactor, CustomParticle::GetCloud(), CustomParticle::GetSpectator(), GeV, CustomParticleFactory::isCustomParticle(), LogDebug, myConfig, and myHelper.

Referenced by ConstructProcess().

                                         {
 
   edm::LogInfo("CustomPhysics") << " CustomPhysicsList: adding CustomPhysics processes "
                 << "for the list of particles: \n";
   aParticleIterator->reset();
 
   while((*aParticleIterator)()) {
     G4ParticleDefinition* particle = aParticleIterator->value();
     if(CustomParticleFactory::isCustomParticle(particle)) {
       CustomParticle* cp = dynamic_cast<CustomParticle*>(particle);
       edm::LogInfo("CustomPhysics") << particle->GetParticleName()
                     <<"  PDGcode= "<<particle->GetPDGEncoding()
                     << "  Mass= "
                     <<particle->GetPDGMass()/GeV  <<" GeV.";
       if(cp->GetCloud()!=0) {
     edm::LogInfo("CustomPhysics") << particle->GetParticleName()
                       <<" CloudMass= "
                       <<cp->GetCloud()->GetPDGMass()/GeV
                       <<" GeV; SpectatorMass= "
                       << cp->GetSpectator()->GetPDGMass()/GeV
                       <<" GeV.";
       }
       G4ProcessManager* pmanager = particle->GetProcessManager();
       if(pmanager) {
     if(particle->GetPDGCharge() != 0.0) {
       pmanager->AddProcess(new G4hMultipleScattering,-1, 1, 1);
       pmanager->AddProcess(new G4hIonisation,        -1, 2, 2);
     }
     if(cp != 0) {
       if(particle->GetParticleType()=="rhadron" || 
          particle->GetParticleType()=="mesonino" || 
          particle->GetParticleType() == "sbaryon"){
         if(!myHelper) myHelper = new G4ProcessHelper(myConfig);
         pmanager->AddDiscreteProcess(new FullModelHadronicProcess(myHelper));
       }
           if(particle->GetParticleType()=="darkpho"){
             CMSDarkPairProductionProcess * darkGamma = new CMSDarkPairProductionProcess(dfactor);
             pmanager->AddDiscreteProcess(darkGamma);
           }   
     }
       }
       else      LogDebug("CustomPhysics") << "   No pmanager";
     }
   }
 }

void CustomPhysicsList::ConstructParticle ( )

virtual

Definition at line 44 of file CustomPhysicsList.cc.

References CustomParticleFactory::loadCustomParticles(), and particleDefFilePath.

                                          {
   CustomParticleFactory::loadCustomParticles(particleDefFilePath);     
 }

void CustomPhysicsList::ConstructProcess ( )

virtual

Definition at line 48 of file CustomPhysicsList.cc.

References addCustomPhysics().

                                          {
   addCustomPhysics();
 }

void CustomPhysicsList::setupRHadronPhycis ( G4ParticleDefinition * particle )

private

Definition at line 99 of file CustomPhysicsList.cc.

References CommonMethods::cp(), CustomParticle::GetCloud(), GeV, LogDebug, myConfig, and myHelper.

 {
   //    LogDebug("CustomPhysics")<<"Configuring rHadron: "
   //    <<cp->
 
   CustomParticle* cp = dynamic_cast<CustomParticle*>(particle);
   if(cp->GetCloud()!=0) {
     LogDebug("CustomPhysics")
       <<"Cloud mass is "
       <<cp->GetCloud()->GetPDGMass()/GeV
       <<" GeV. Spectator mass is "
       <<static_cast<CustomParticle*>(particle)->GetSpectator()->GetPDGMass()/GeV
       <<" GeV.";
   }
   G4ProcessManager* pmanager = particle->GetProcessManager();
   if(pmanager){
     if(!myHelper) myHelper = new G4ProcessHelper(myConfig);
     if(particle->GetPDGCharge()/eplus != 0){
       pmanager->AddProcess(new G4hMultipleScattering,-1, 1, 1);
       pmanager->AddProcess(new G4hIonisation,        -1, 2, 2);
     }
     pmanager->AddDiscreteProcess(new FullModelHadronicProcess(myHelper)); //GHEISHA
   }
   else      LogDebug("CustomPhysics") << "   No pmanager";
 }

void CustomPhysicsList::setupSUSYPhycis ( G4ParticleDefinition * particle )

private

Definition at line 125 of file CustomPhysicsList.cc.

References LogDebug.

 {
   G4ProcessManager* pmanager = particle->GetProcessManager();
   if(pmanager){
     if(particle->GetPDGCharge()/eplus != 0){
       pmanager->AddProcess(new G4hMultipleScattering,-1, 1, 1);
       pmanager->AddProcess(new G4hIonisation,        -1, 2, 2);
     }
     pmanager->AddProcess(new G4Decay, 1, -1, 3);
   }
   else      LogDebug("CustomPhysics") << "   No pmanager";
 }