BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSPhysicsCherenkov.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 "BDSPhysicsCherenkov.hh"
20
21#include "globals.hh" // geant4 types / globals
22#include "G4AutoDelete.hh"
23#include "G4Cerenkov.hh"
24#include "G4Gamma.hh"
25#include "G4LeptonConstructor.hh"
26#include "G4OpticalPhoton.hh"
27#include "G4ParticleDefinition.hh"
28#include "G4PhysicsListHelper.hh"
29#include "G4Version.hh"
30
31BDSPhysicsCherenkov::BDSPhysicsCherenkov(G4double maxPhotonsPerStepIn,
32 G4double maxBetaChangePerStepIn):
33 G4VPhysicsConstructor("BDSPhysicsCherenkov"),
34 maxPhotonsPerStep(maxPhotonsPerStepIn),
35 maxBetaChangePerStep(maxBetaChangePerStepIn)
36{;}
37
38BDSPhysicsCherenkov::~BDSPhysicsCherenkov()
39{;}
40
41void BDSPhysicsCherenkov::ConstructParticle()
42{
43 G4LeptonConstructor leptons;
44 leptons.ConstructParticle();
45
46 G4Gamma::Gamma();
47 G4OpticalPhoton::OpticalPhotonDefinition();
48}
49
50void BDSPhysicsCherenkov::ConstructProcess()
51{
52 if (Activated())
53 {return;}
54
55 G4Cerenkov* cherenkov = new G4Cerenkov();
56 G4AutoDelete::Register(cherenkov);
57
58 // reduce memory profile
59 cherenkov->SetTrackSecondariesFirst(true);
60
61 // common settings (similar to optical physics)
62 if (maxPhotonsPerStep > 0)
63 {cherenkov->SetMaxNumPhotonsPerStep(maxPhotonsPerStep);}
64 cherenkov->SetMaxBetaChangePerStep(maxBetaChangePerStep);
65
66 G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
67
68#if G4VERSION_NUMBER > 1029
69 auto aParticleIterator = GetParticleIterator();
70#endif
71 aParticleIterator->reset();
72 while( (*aParticleIterator)() )
73 {
74 G4ParticleDefinition* particle = aParticleIterator->value();
75
76 if (cherenkov->IsApplicable(*particle))
77 {ph->RegisterProcess(cherenkov, particle);}
78 }
79
81}
G4double maxPhotonsPerStep
Parameter from global constants to control cherenkov physics.
G4double maxBetaChangePerStep
Parameter from global constants to control cherenkov physics.
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