20#include "BDSGlobalConstants.hh"
21#include "BDSPhysicsMuon.hh"
23#include "G4AutoDelete.hh"
24#include "G4AnnihiToMuPair.hh"
25#include "G4eeToHadrons.hh"
27#include "G4GammaConversionToMuons.hh"
28#include "G4LeptonConstructor.hh"
29#include "G4MesonConstructor.hh"
30#include "G4MuBremsstrahlung.hh"
31#include "G4MuIonisation.hh"
32#include "G4MuMultipleScattering.hh"
33#include "G4MuPairProduction.hh"
34#include "G4OpticalPhoton.hh"
35#include "G4ParticleDefinition.hh"
36#include "G4PhysicsListHelper.hh"
37#include "G4Version.hh"
40BDSPhysicsMuon::BDSPhysicsMuon():
44BDSPhysicsMuon::BDSPhysicsMuon(G4bool emWillBeUsedIn):
45 G4VPhysicsConstructor(
"BDSPhysicsMuon"),
46 emWillBeUsed(emWillBeUsedIn)
49BDSPhysicsMuon::~BDSPhysicsMuon()
55 G4LeptonConstructor leptons;
56 leptons.ConstructParticle();
59 G4MesonConstructor mConstructor;
60 mConstructor.ConstructParticle();
64 G4OpticalPhoton::OpticalPhoton();
73 G4GammaConversionToMuons* gammaToMuPair =
new G4GammaConversionToMuons();
74 G4AutoDelete::Register(gammaToMuPair);
77 G4AnnihiToMuPair* ePlusToMuPair =
new G4AnnihiToMuPair();
78 G4AutoDelete::Register(ePlusToMuPair);
79 G4eeToHadrons* eeToHadrons =
new G4eeToHadrons();
80 G4AutoDelete::Register(eeToHadrons);
83 G4MuMultipleScattering* mumsc =
nullptr;
84 G4MuIonisation* muion =
nullptr;
85 G4MuBremsstrahlung* mubrm =
nullptr;
88 mumsc =
new G4MuMultipleScattering();
89 G4AutoDelete::Register(mumsc);
90 muion =
new G4MuIonisation();
91 G4AutoDelete::Register(muion);
92 mubrm =
new G4MuBremsstrahlung();
93 G4AutoDelete::Register(mubrm);
95 G4MuPairProduction* mupar =
new G4MuPairProduction();
96 G4AutoDelete::Register(mupar);
98 G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
100#if G4VERSION_NUMBER > 1029
101 auto aParticleIterator = GetParticleIterator();
103 aParticleIterator->reset();
105 while( (*aParticleIterator)() )
107 G4ParticleDefinition* particle = aParticleIterator->value();
108 G4String particleName = particle->GetParticleName();
110 if(particleName ==
"gamma")
112 ph->RegisterProcess(gammaToMuPair, particle);
115 if(particleName ==
"e+")
117 ph->RegisterProcess(ePlusToMuPair, particle);
118 ph->RegisterProcess(eeToHadrons, particle);
121 if(particleName ==
"mu+" || particleName ==
"mu-")
125 ph->RegisterProcess(mumsc, particle);
126 ph->RegisterProcess(muion, particle);
127 ph->RegisterProcess(mubrm, particle);
128 ph->RegisterProcess(mupar, particle);
High energy muon processes.
virtual void ConstructProcess()
Construct and attach the processes to the relevant particles.
virtual void ConstructParticle()
Construct all leptons, photons (inc optical), and pion +- just in case.
G4bool Activated() const
Get whether this instance has been activated.
void SetActivated()
Flag this instance as activated for later querying.