#include "MyG4DecayPhysics.hh"

Inheritance diagram for MyG4DecayPhysics:

Public Member Functions
	MyG4DecayPhysics (G4int ver=1)

	MyG4DecayPhysics (const G4String &name, G4int ver=1)

virtual	~MyG4DecayPhysics ()

virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

virtual G4Decay *	GetDecayProcess ()

Private Attributes
G4int	verbose

Static Private Attributes
static G4ThreadLocal G4Decay *	fDecayProcess = 0

static G4ThreadLocal G4bool	wasActivated = false

Detailed Description

Definition at line 50 of file MyG4DecayPhysics.hh.

Constructor & Destructor Documentation

MyG4DecayPhysics::MyG4DecayPhysics ( G4int ver = 1 )

Definition at line 65 of file MyG4DecayPhysics.cc.

   :  G4VPhysicsConstructor("Decay"), verbose(ver)
 {
 }

MyG4DecayPhysics::MyG4DecayPhysics	(	const G4String &	name,
		G4int	ver = `1`
	)

Definition at line 70 of file MyG4DecayPhysics.cc.

   :  G4VPhysicsConstructor(name), verbose(ver)
 {
 }

MyG4DecayPhysics::~MyG4DecayPhysics ( )

virtual

Definition at line 75 of file MyG4DecayPhysics.cc.

76 {

77 }

Member Function Documentation

void MyG4DecayPhysics::ConstructParticle ( )

virtual

Definition at line 79 of file MyG4DecayPhysics.cc.

 {
 
 // G4cout << "MyG4DecayPhysics::ConstructParticle" << G4endl;
   G4BosonConstructor  pBosonConstructor;
   pBosonConstructor.ConstructParticle();
 
   G4LeptonConstructor pLeptonConstructor;
   pLeptonConstructor.ConstructParticle();
 
   G4MesonConstructor pMesonConstructor;
   pMesonConstructor.ConstructParticle();
 
   G4BaryonConstructor pBaryonConstructor;
   pBaryonConstructor.ConstructParticle();
 
   G4IonConstructor pIonConstructor;
   pIonConstructor.ConstructParticle();
 
   G4ShortLivedConstructor pShortLivedConstructor;
   pShortLivedConstructor.ConstructParticle();  
 }

void MyG4DecayPhysics::ConstructProcess ( )

virtual

Definition at line 102 of file MyG4DecayPhysics.cc.

References fDecayProcess, verbose, and wasActivated.

 {
   if(wasActivated) { return; }
   wasActivated = true;
 
   G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
   // Add Decay Process
   fDecayProcess = new G4Decay();
   auto myParticleIterator=GetParticleIterator();
   myParticleIterator->reset();
   G4ParticleDefinition* particle=0;
 
   while( (*myParticleIterator)() )
   {
     particle = myParticleIterator->value();
 
     // -- check for mu+
     bool isMuPlus = particle->GetParticleName() == "mu+";
     if( fDecayProcess->IsApplicable(*particle) && !isMuPlus ) 
     { 
       if(verbose > 1) {
         G4cout << "### Decays for " << particle->GetParticleName() << G4endl;
       }
       ph->RegisterProcess(fDecayProcess, particle);
     }
     if (isMuPlus) {
       G4cout << "MyG4DecayPhysics::ConstructProcess() mu+ decays will be disabled\n" << G4endl;
       particle->DumpTable();
       particle->GetProcessManager()->DumpInfo();
     }
   }
 }