| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // -------------------------------------------
28 //
29 // GEANT4 Class header file
30 //
31 //
32 // File name: G4RiGeMuPairProductionModel
33 //
34 // Authors: Girardo Depaola & Ricardo Pa
35 //
36 // Creation date: 29.10.2024
37 //
38 //
39 // -------------------------------------------
40 //
41
42 #ifndef G4RiGeMuPairProductionModel_h
43 #define G4RiGeMuPairProductionModel_h 1
44
45 #include "G4VEmModel.hh"
46 #include "G4NistManager.hh"
47 #include "G4ElementData.hh"
48 #include "G4Physics2DVector.hh"
49 #include "G4VEmAngularDistribution.hh"
50 #include <vector>
51
52 class G4Element;
53 class G4ParticleChangeForLoss;
54 class G4RiGeAngularGenerator;
55
56 class G4RiGeMuPairProductionModel : public G4V
57 {
58 public:
59
60 explicit G4RiGeMuPairProductionModel(const G
61
62 ~G4RiGeMuPairProductionModel() override = de
63
64 void Initialise(const G4ParticleDefinition*,
65
66 void InitialiseLocal(const G4ParticleDefinit
67 G4VEmModel* masterModel
68
69 G4double ComputeCrossSectionPerAtom(const G4
70 G4double kineticEnergy,
71 G4double Z, G4double A,
72 G4double cutEnergy,
73 G4double maxEnergy) override;
74
75 G4double ComputeDEDXPerVolume(const G4Materi
76 const G4Partic
77 G4double kinet
78 G4double cutEn
79
80 void SampleSecondaries(std::vector<G4Dynamic
81 const G4MaterialCutsC
82 const G4DynamicPartic
83 G4double tmin,
84 G4double maxEnergy) o
85
86 G4double MinPrimaryEnergy(const G4Material*,
87 const G4ParticleDe
88 G4double) override
89
90 G4double
91 ComputeDMicroscopicCrossSection(G4double tki
92 G4double pairEnergy);
93
94 inline void SetLowestKineticEnergy(G4double
95
96 inline void SetParticle(const G4ParticleDefi
97
98 // hide assignment operator and copy constru
99 G4RiGeMuPairProductionModel& operator=
100 (const G4RiGeMuPairProductionModel& right) =
101 G4RiGeMuPairProductionModel(const G4RiGeMuPa
102
103 protected:
104
105 G4double ComputMuPairLoss(G4double Z, G4doub
106 G4double tmax);
107
108 G4double ComputeMicroscopicCrossSection(G4do
109 G4do
110 G4do
111
112 G4double FindScaledEnergy(G4int Z, G4double
113 G4double yymin, G4double yymax);
114
115 inline G4double MaxSecondaryEnergyForElement
116 G4double Z);
117
118 void MakeSamplingTables();
119
120 void StoreTables() const;
121
122 G4bool RetrieveTables();
123
124 virtual void DataCorrupted(G4int Z, G4double
125
126 G4ParticleChangeForLoss* fParticleChange = n
127 const G4ParticleDefinition* particle = nullp
128 G4NistManager* nist = nullptr;
129
130 G4double factorForCross;
131 G4double sqrte;
132 G4double particleMass = 0.0;
133 G4double z13 = 0.0;
134 G4double z23 = 0.0;
135 G4double lnZ = 0.0;
136
137 G4double minPairEnergy;
138 G4double lowestKinEnergy;
139
140 G4double emin;
141 G4double emax;
142 G4double ymin = -5.0;
143 G4double dy = 0.005;
144
145 // Random numbers for sampling
146 G4double randNumbs[9];
147
148 G4int currentZ = 0;
149 G4int nYBinPerDecade = 4;
150 std::size_t nbiny = 1000;
151 std::size_t nbine = 0;
152
153 G4bool fTableToFile = false;
154
155 // static members
156 static const G4int NZDATPAIR = 5;
157 static const G4int NINTPAIR = 8;
158 static const G4int ZDATPAIR[NZDATPAIR];
159 static const G4double xgi[NINTPAIR];
160 static const G4double wgi[NINTPAIR];
161
162 private:
163
164 G4RiGeAngularGenerator* fAngularGenerator;
165 G4ParticleDefinition* theElectron;
166 G4ParticleDefinition* thePositron;
167 G4String dataName{""};
168 };
169
170 //....oooOO0OOooo........oooOO0OOooo........oo
171
172 inline void G4RiGeMuPairProductionModel::SetLo
173 {
174 lowestKinEnergy = e;
175 }
176
177 //....oooOO0OOooo........oooOO0OOooo........oo
178
179 inline
180 void G4RiGeMuPairProductionModel::SetParticle(
181 {
182 if(nullptr == particle) {
183 particle = p;
184 particleMass = particle->GetPDGMass();
185 }
186 }
187
188 //....oooOO0OOooo........oooOO0OOooo........oo
189
190 inline G4double
191 G4RiGeMuPairProductionModel::MaxSecondaryEnerg
192 G4double ZZ)
193 {
194 G4int Z = G4lrint(ZZ);
195 if(Z != currentZ) {
196 currentZ = Z;
197 z13 = nist->GetZ13(Z);
198 z23 = z13*z13;
199 lnZ = nist->GetLOGZ(Z);
200 }
201 return kineticEnergy + particleMass*(1.0 - 0
202 }
203
204 //....oooOO0OOooo........oooOO0OOooo........oo
205
206 #endif
207