BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSPhysicsAnnihiToMuMu.cc
1/*
2Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway,
3University of London 2001 - 2023.
4
5This file is part of BDSIM.
6
7BDSIM is free software: you can redistribute it and/or modify
8it under the terms of the GNU General Public License as published
9by the Free Software Foundation version 3 of the License.
10
11BDSIM is distributed in the hope that it will be useful, but
12WITHOUT ANY WARRANTY; without even the implied warranty of
13MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14GNU General Public License for more details.
15
16You should have received a copy of the GNU General Public License
17along with BDSIM. If not, see <http://www.gnu.org/licenses/>.
18*/
19#include "BDSDebug.hh"
20#include "BDSPhysicsAnnihiToMuMu.hh"
21
22#include "G4AutoDelete.hh"
23#include "G4AnnihiToMuPair.hh"
24#include "G4LeptonConstructor.hh"
25#include "G4ParticleDefinition.hh"
26#include "G4PhysicsListHelper.hh"
27#include "G4String.hh"
28#include "G4Version.hh"
29
30
31BDSPhysicsAnnihiToMuMu::BDSPhysicsAnnihiToMuMu():
32 G4VPhysicsConstructor("BDSPhysicsAnnihiToMuMu")
33{;}
34
35BDSPhysicsAnnihiToMuMu::~BDSPhysicsAnnihiToMuMu()
36{;}
37
39{
40 G4LeptonConstructor leptons;
41 leptons.ConstructParticle();
42}
43
45{
46 if (Activated())
47 {return;}
48
49 G4AnnihiToMuPair* ePlusToMuPair = new G4AnnihiToMuPair();
50 G4AutoDelete::Register(ePlusToMuPair);
51
52 G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
53
54#if G4VERSION_NUMBER > 1029
55 auto aParticleIterator = GetParticleIterator();
56#endif
57 aParticleIterator->reset();
58
59 while( (*aParticleIterator)() )
60 {
61 G4ParticleDefinition* particle = aParticleIterator->value();
62 G4String particleName = particle->GetParticleName();
63
64 if (particleName == "e+")
65 {ph->RegisterProcess(ePlusToMuPair, particle);}
66 }
67
69}
virtual void ConstructParticle()
Construct all leptions, photons (inc optical), and pion +- just in case.
virtual void ConstructProcess()
Construct and attach the processes to the relevant particles.
G4bool Activated() const
Get whether this instance has been activated.
Definition: BDSSingleUse.hh:37
void SetActivated()
Flag this instance as activated for later querying.
Definition: BDSSingleUse.hh:40