Geant4 11.4.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4RiGeAngularGenerator.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// -------------------------------------------------------------------
28//
29// GEANT4 Class file
30//
31//
32// File name: G4RiGeAngularGenerator
33//
34// Authors: Gerardo Depaola & Ricardo Pacheco
35//
36// Creation date: 27 October 2024
37//
38// -------------------------------------------------------------------
39//
40
41#include "G4RiGeAngularGenerator.hh"
42#include "Randomize.hh"
43#include "G4Log.hh"
44#include "G4Exp.hh"
45#include <CLHEP/Units/PhysicalConstants.h>
46
47//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
48
52
53//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
54
55G4ThreeVector&
56G4RiGeAngularGenerator::SampleDirection(const G4DynamicParticle* dp,
57 G4double gEnergy, G4int,
58 const G4Material*)
59{
60 // Sample gamma angle (Z - axis along the parent particle).
61 G4double cost = SampleCosTheta(dp->GetKineticEnergy(), gEnergy,
62 dp->GetDefinition()->GetPDGMass());
63 G4double sint = std::sqrt((1.0 - cost)*(1.0 + cost));
64 G4double phi = CLHEP::twopi*G4UniformRand();
65
66 fLocalDirection.set(sint*std::cos(phi), sint*std::sin(phi), cost);
67 fLocalDirection.rotateUz(dp->GetMomentumDirection());
68
69 return fLocalDirection;
70}
71
72//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
73
74G4double G4RiGeAngularGenerator::SampleCosTheta(G4double primKinEnergy,
75 G4double gEnergy,
76 G4double mass)
77{
78 G4double gam = 1.0 + primKinEnergy/mass;
79 G4double rmax = gam*CLHEP::halfpi*std::min(1.0, gam*mass/gEnergy - 1.0);
80 G4double rmax2= rmax*rmax;
81 G4double x = G4UniformRand()*rmax2/(1.0 + rmax2);
82
83 return std::cos(std::sqrt(x/(1.0 - x))/gam);
84}
85
86//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
87
88void G4RiGeAngularGenerator::Sample5DPairDirections(const G4DynamicParticle* dp,
89 G4ThreeVector& dirElectron,
90 G4ThreeVector& dirPositron,
91 const G4double gEnergy, const G4double q2,
92 const G4double gMomentum,
93 const G4double muFinalMomentum,
94 const G4double muFinalEnergy,
95 const G4double* randNumbs,
96 const G4double* W)
97{
98 G4double muKinEnergy = dp->GetKineticEnergy();
99 G4ThreeVector muMomentumVector = dp->GetMomentum();
100 G4double muMomentum = muMomentumVector.mag();
101
102 // Electron mass
103 G4double eMass = CLHEP::electron_mass_c2;
104 G4double eMass2 = eMass*eMass;
105
106 // Muon mass
107 G4double muMass = dp->GetDefinition()->GetPDGMass();
108
109 G4double muEnergy = muKinEnergy + muMass;
110 G4LorentzVector muFinalFourMomentum(muMomentumVector, muEnergy);
111
112 G4double mint3 = 0.;
113 G4double maxt3 = CLHEP::pi;
114 G4double Cmin = std::cos(maxt3);
115 G4double Cmax = std::cos(mint3);
116
117 if (randNumbs[7] < W[0]) {
118 G4double A1 = -(q2 - 2.*muEnergy*gEnergy);
119 G4double B1 = -(2.*gMomentum*muMomentum);
120 G4double tginterval = G4Log((A1 + B1)/(A1 - B1))/B1;
121
122 G4double costg = (-A1 + (A1 - B1)*G4Exp(B1*tginterval*randNumbs[1]))/B1;
123 G4double sintg = std::sqrt((1.0 - costg)*(1.0 + costg));
124 G4double phig = CLHEP::twopi*randNumbs[2];
125 G4double sinpg = std::sin(phig);
126 G4double cospg = std::cos(phig);
127
128 G4ThreeVector dirGamma;
129 dirGamma.set(sintg*cospg, sintg*sinpg, costg);
130 G4LorentzVector gFourMomentum(dirGamma*gMomentum, gEnergy);
131
132 G4double cost5 = -1. + 2.*randNumbs[6];
133 G4double phi5 = CLHEP::twopi*randNumbs[8];
134
135 G4LorentzVector eFourMomentumMQ = eDP2(q2, eMass2, eMass2, cost5, phi5);
136 G4LorentzVector pFourMomentumMQ = pDP2(eMass2, eFourMomentumMQ);
137
138 G4LorentzVector eFourMomentum = eFourMomentumMQ.boost(gFourMomentum.boostVector());
139 G4LorentzVector pFourMomentum = pFourMomentumMQ.boost(gFourMomentum.boostVector());
140
141 dirElectron = eFourMomentum.vect().unit();
142 dirPositron = pFourMomentum.vect().unit();
143
144 G4double phi = CLHEP::twopi*randNumbs[3];
145 PhiRotation(dirElectron, phi);
146 PhiRotation(dirPositron, phi);
147
148 } else if (randNumbs[7] >= W[0] && randNumbs[7] < W[1]) {
149 G4double A3 = q2 + 2.*gEnergy*muFinalEnergy;
150 G4double B3 = -2.*gMomentum*muFinalMomentum;
151
152 G4double tQ3interval = G4Log((A3 + B3)/(A3 - B3))/B3;
153 G4double tQMG = (-A3 + (A3 - B3)*G4Exp(B3*tQ3interval*randNumbs[0]))/B3;
154 G4double phiQP = CLHEP::twopi*randNumbs[2];
155
156 G4double sintQ3 = std::sqrt(1. - tQMG*tQMG);
157 G4double cospQP = std::cos(phiQP);
158 G4double sinpQP = std::sin(phiQP);
159
160 G4double Ap = muMomentum*muMomentum + muFinalMomentum*muFinalMomentum + gMomentum*gMomentum;
161 G4double A = Ap + 2.*muFinalMomentum*gMomentum*tQMG;
162 G4double B = -2.*muMomentum*gMomentum*sintQ3*cospQP;
163 G4double C = -2.*muMomentum*gMomentum*tQMG - 2.*muMomentum*muFinalMomentum;
164
165 G4double absB = std::abs(B);
166 G4double t3interval = (1./(A + C + absB*mint3) - 1./(A + C + absB*maxt3))/absB;
167 G4double t3 = (-(A + C) + 1./(1./(A + C + absB*mint3) - absB*t3interval*randNumbs[0]))/absB;
168 G4double sint3 = std::sin(t3);
169 G4double cost3 = std::cos(t3);
170
171 G4double cost = -sint3*sintQ3*cospQP + cost3*tQMG;
172 G4double sint = std::sqrt((1. + cost)*(1. - cost));
173 G4double cosp = (sintQ3*cospQP*cost3 + sint3*tQMG)/sint;
174 G4double sinp = sintQ3*sinpQP/sint;
175
176 G4ThreeVector dirGamma;
177 dirGamma.set(sint*cosp, sint*sinp, cost);
178 G4LorentzVector gFourMomentum(dirGamma*gMomentum, gEnergy);
179
180 G4double cost5 = -1. + 2.*randNumbs[6];
181 G4double phi5 = CLHEP::twopi*randNumbs[8];
182
183 G4LorentzVector eFourMomentumMQ = eDP2(q2, eMass2, eMass2, cost5, phi5);
184 G4LorentzVector pFourMomentumMQ = pDP2(eMass2, eFourMomentumMQ);
185
186 G4LorentzVector eFourMomentum = eFourMomentumMQ.boost(gFourMomentum.boostVector());
187 G4LorentzVector pFourMomentum = pFourMomentumMQ.boost(gFourMomentum.boostVector());
188
189 dirElectron = eFourMomentum.vect().unit();
190 dirPositron = pFourMomentum.vect().unit();
191
192 G4double phi = CLHEP::twopi*randNumbs[3];
193 PhiRotation(dirElectron, phi);
194 PhiRotation(dirPositron, phi);
195
196 } else if (randNumbs[7] >= W[1] && randNumbs[7] < W[2]) {
197 G4double phi3 = CLHEP::twopi*randNumbs[0];
198 G4double phi5 = CLHEP::twopi*randNumbs[1];
199 G4double phi6 = CLHEP::twopi*randNumbs[2];
200 G4double minmuFinalEnergy = muMass;
201 G4double muEnergyInterval = muEnergy - 2.*eMass - minmuFinalEnergy;
202 G4double muFEnergy = minmuFinalEnergy + muEnergyInterval*randNumbs[3];
203
204 G4double mineEnergy = eMass;
205 G4double maxeEnergy = muEnergy - muFEnergy - eMass;
206 G4double eEnergyInterval = maxeEnergy - mineEnergy;
207 G4double eEnergy = mineEnergy + eEnergyInterval*randNumbs[4];
208
209 G4double cosp3 = 1.;
210 G4double sinp3 = 0.;
211 G4double cosp5 = std::cos(phi5);
212 G4double sinp5 = std::sin(phi5);
213 G4double cosp6 = std::cos(phi6);
214 G4double sinp6 = std::sin(phi6);
215
216 G4double muFMomentum = std::sqrt(muFEnergy*muFEnergy - muMass*muMass);
217 G4double eMomentum = std::sqrt(eEnergy*eEnergy - eMass*eMass);
218 G4double pEnergy = muEnergy - muFEnergy - eEnergy;
219 G4double pMomentum = std::sqrt(pEnergy*pEnergy - eMass*eMass);
220
221 G4double A3 = -2.*muMass*muMass + 2.*muEnergy*muFEnergy;
222 G4double B3 = -2.*std::sqrt(muEnergy*muEnergy - muMass*muMass)*muFMomentum;
223 G4double cost3interval = G4Log((A3 + B3*Cmax)/(A3 + B3*Cmin))/B3;
224
225 G4double expanCost3r6 = G4Exp(B3*cost3interval*randNumbs[5]);
226 G4double cost3 = A3*(expanCost3r6 - 1.)/B3 + Cmin*expanCost3r6;
227 G4double sint3 = std::sqrt((1. - cost3)*(1. + cost3));
228
229 G4ThreeVector muFinalMomentumVector(muFMomentum*sint3, 0., muFMomentum*cost3);
230
231 G4LorentzVector muFourMomentum(muMomentumVector, muEnergy);
232 muFinalFourMomentum.set(muFinalMomentumVector, muFEnergy);
233 G4LorentzVector auxVec1 = muFourMomentum - muFinalFourMomentum;
234 G4double A5 = auxVec1.mag2() - 2.*eEnergy*(muEnergy - muFEnergy) +
235 2.*muMomentumVector[2]*eMomentum - 2.*muFMomentum*eMomentum*cost3;
236 G4double B5 = -2.*muFMomentum*eMomentum*(sint3*cosp3*cosp5 + sint3*sinp3*sinp5);
237 G4double absA5 = std::abs(A5);
238 G4double absB5 = std::abs(B5);
239 G4double mint5 = 0.;
240 G4double maxt5 = CLHEP::pi;
241 G4double t5interval = G4Log((absA5 + absB5*maxt5)/(absA5 + absB5*mint5))/absB5;
242 G4double argexp = absB5*t5interval*randNumbs[6] + G4Log(absA5 + absB5*mint5);
243 G4double t5 = -absA5/absB5 + G4Exp(argexp)/absB5;
244 G4double sint5 = std::sin(t5);
245 G4double cost5 = std::cos(t5);
246
247 dirElectron.set(sint5*cosp5, sint5*sinp5, cost5);
248 G4ThreeVector eMomentumVector = eMomentum*dirElectron;
249
250 G4ThreeVector auxVec2 = muMomentumVector - muFinalMomentumVector - eMomentumVector;
251 G4double p1mp3mp52 = auxVec2.dot(auxVec2);
252 G4double Bp = muFinalMomentum*(sint3*cosp3*cosp6 + sint3*sinp3*sinp6) +
253 eMomentum*(sint5*cosp5*cosp6 + sint5*sinp5*sinp6);
254 G4double Cp = -muMomentum + muFMomentum*cost3 + eMomentum*cost5;
255 G4double A6 = p1mp3mp52 + pMomentum*pMomentum;
256 G4double B6 = 2.*pMomentum*Bp;
257 G4double C6 = 2.*pMomentum*Cp;
258 G4double mint6 = 0.;
259 G4double maxt6 = CLHEP::pi;
260 G4double absA6C6 = std::abs(A6 + C6);
261 G4double absB6 = std::abs(B6);
262 G4double t6interval = (1./(absA6C6 + absB6*mint6) - 1./(absA6C6 + absB6*maxt6))/absB6;
263 G4double t6 = (-absA6C6 + 1./(1./(absA6C6 + absB6*mint6) - absB6*t6interval*randNumbs[8]))/absB6;
264 G4double sint6 = std::sin(t6);
265 G4double cost6 = std::cos(t6);
266
267 dirPositron.set(sint6*cosp6, sint6*sinp6, cost6);
268
269 PhiRotation(dirElectron, phi3);
270 PhiRotation(dirPositron, phi3);
271
272 } else if (randNumbs[7] >= W[2]) {
273 G4double phi3 = CLHEP::twopi*randNumbs[0];
274 G4double phi6 = CLHEP::twopi*randNumbs[1];
275 G4double phi5 = CLHEP::twopi*randNumbs[2];
276 G4double minmuFinalEnergy = muMass;
277 G4double muFinalEnergyinterval = muEnergy - 2.*eMass - minmuFinalEnergy;
278 G4double muFEnergy = minmuFinalEnergy + muFinalEnergyinterval*randNumbs[3];
279
280 G4double minpEnergy = eMass;
281 G4double maxpEnergy = muEnergy - muFEnergy - eMass;
282 G4double pEnergyinterval = maxpEnergy - minpEnergy;
283 G4double pEnergy = minpEnergy + pEnergyinterval*randNumbs[4];
284
285 G4double cosp3 = 1.;
286 G4double sinp3 = 0.;
287 G4double cosp5 = std::cos(phi5);
288 G4double sinp5 = std::sin(phi5);
289 G4double cosp6 = std::cos(phi6);
290 G4double sinp6 = std::sin(phi6);
291
292 G4double muFMomentum = std::sqrt(muFEnergy*muFEnergy - muMass*muMass);
293 G4double pMomentum = std::sqrt(pEnergy*pEnergy - eMass*eMass);
294 G4double eEnergy = muEnergy - muFEnergy - pEnergy;
295 G4double eMomentum = std::sqrt(eEnergy*eEnergy - eMass*eMass);
296
297 G4double A3 = -2.*muMass*muMass + 2.*muEnergy*muFEnergy;
298 G4double B3 = -2.*std::sqrt(muEnergy*muEnergy - muMass*muMass)*muFMomentum;
299 G4double cost3interval = G4Log((A3 + B3*Cmax)/(A3 + B3*Cmin))/B3;
300
301 G4double expanCost3r6 = G4Exp(B3*cost3interval*randNumbs[5]);
302 G4double cost3 = A3*(expanCost3r6 - 1.)/B3 + Cmin*expanCost3r6;
303 G4double sint3 = std::sqrt((1. - cost3)*(1. + cost3));
304
305 G4ThreeVector muFinalMomentumVector;
306 muFinalMomentumVector.set(muFMomentum*sint3*cosp3, muFMomentum*sint3*sinp3,
307 muFMomentum*cost3);
308
309 G4LorentzVector muFourMomentum(muMomentumVector, muEnergy);
310 muFinalFourMomentum.set(muFinalMomentumVector, muFEnergy);
311 G4LorentzVector auxVec1 = muFourMomentum - muFinalFourMomentum;
312 G4double A6 = auxVec1.mag2() - 2.*pEnergy*(muEnergy - muFEnergy) +
313 2.*muMomentumVector[2]*pMomentum - 2.*muFMomentum*pMomentum*cost3;
314 G4double B6 = -2.*muFMomentum*pMomentum*(sint3*cosp3*cosp6 + sint3*sinp3*sinp6);
315 G4double absA6 = std::abs(A6);
316 G4double absB6 = std::abs(B6);
317 G4double mint6 = 0.;
318 G4double maxt6 = CLHEP::pi;
319 G4double t6interval = G4Log((absA6 + absB6*maxt6)/(absA6 + absB6*mint6))/absB6;
320 G4double argexp = absB6*t6interval*randNumbs[6] + G4Log(absA6 + absB6*mint6);
321 G4double t6 = -absA6/absB6 + G4Exp(argexp)/absB6;
322 G4double sint6 = std::sin(t6);
323 G4double cost6 = std::cos(t6);
324
325 dirPositron.set(sint6*cosp6, sint6*sinp6, cost6);
326 G4ThreeVector pMomentumVector = pMomentum*dirPositron;
327
328 G4ThreeVector auxVec2 = muMomentumVector - muFinalMomentumVector - pMomentumVector;
329 G4double p1mp3mp62 = auxVec2.dot(auxVec2);
330 G4double Bp = muFMomentum*(sint3*cosp3*cosp5 + sint3*sinp3*sinp5) +
331 pMomentum*(sint6*cosp6*cosp5 + sint6*sinp6*sinp5);
332 G4double Cp = -muMomentum + muFMomentum*cost3 + pMomentum*cost6;
333 G4double A5 = p1mp3mp62 + eMomentum*eMomentum;
334 G4double B5 = 2.*eMomentum*Bp;
335 G4double C5 = 2.*eMomentum*Cp;
336 G4double mint5 = 0.;
337 G4double maxt5 = CLHEP::pi;
338 G4double absA5C5 = std::abs(A5 + C5);
339 G4double absB5 = std::abs(B5);
340 G4double t5interval = (1./(absA5C5 + absB5*mint5) - 1./(absA5C5 + absB5*maxt5))/absB5;
341 G4double t5 = (-absA5C5 + 1./(1./(absA5C5 + absB5*mint5) - absB5*t5interval*randNumbs[8]))/absB5;
342 G4double sint5 = std::sin(t5);
343 G4double cost5 = std::cos(t5);
344
345 dirElectron.set(sint5*cosp5, sint5*sinp5, cost5);
346
347 PhiRotation(dirElectron, phi3);
348 PhiRotation(dirPositron, phi3);
349 }
350}
351
352//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
353//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
354
355void G4RiGeAngularGenerator::PhiRotation(G4ThreeVector& dir, G4double phi)
356{
357 G4double sinp = std::sin(phi);
358 G4double cosp = std::cos(phi);
359
360 G4double newX = dir.x()*cosp + dir.y()*sinp;
361 G4double newY = -dir.x()*sinp + dir.y()*cosp;
362
363 dir.setX(newX);
364 dir.setY(newY);
365}
366
367//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
368
369G4LorentzVector G4RiGeAngularGenerator::eDP2(G4double x1, G4double x2,
370 G4double x3, G4double x4,
371 G4double x5)
372{
373 G4double sint = std::sqrt((1.0 - x4)*(1.0 + x4));
374 G4double cosp = std::cos(x5);
375 G4double sinp = std::sin(x5);
376
377 G4double QJM2 = (x1 + x3 - x2)*(x1 + x3 - x2)/(4.*x1) - x3;
378
379 if (QJM2 < 0.) {
380 QJM2 = 1.e-13;
381 }
382
383 G4double QJM = std::sqrt(QJM2);
384
385 G4LorentzVector x6(QJM*sint*cosp, QJM*sint*sinp, QJM*x4, std::sqrt(x2 + QJM2));
386
387 return x6;
388}
389
390//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
391
392G4LorentzVector G4RiGeAngularGenerator::pDP2(G4double x3, const G4LorentzVector& x6)
393{
394 G4LorentzVector x7(-x6.vect(), std::sqrt(x3 + x6.vect().dot(x6.vect())));
395 return x7;
396}
397
398//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
399
400void G4RiGeAngularGenerator::PrintGeneratorInformation() const
401{
402 G4cout << "\n" << G4endl;
403 G4cout << "Angular Generator by RiGe algorithm" << G4endl;
404}
405
406//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
C
G4double C(G4double temp)
Definition G4DNAElectronHoleRecombination.cc:88
B
G4double B(G4double temperature)
Definition G4DNAElectronHoleRecombination.cc:71
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition G4Exp.hh:132
G4Log
G4double G4Log(G4double x)
Definition G4Log.hh:169
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition G4LorentzVector.hh:33
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition G4ThreeVector.hh:36
G4double
double G4double
Definition G4Types.hh:83
G4int
int G4int
Definition G4Types.hh:85
A
const G4double A[17]
Definition G4WaterStopping.cc:53
G4endl
#define G4endl
Definition G4ios.hh:67
G4cout
G4GLOB_DLL std::ostream G4cout
C5
#define C5
G4UniformRand
#define G4UniformRand()
Definition Randomize.hh:52
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::x
double x() const
CLHEP::Hep3Vector::setY
void setY(double)
CLHEP::Hep3Vector::y
double y() const
CLHEP::Hep3Vector::dot
double dot(const Hep3Vector &) const
CLHEP::Hep3Vector::mag
double mag() const
CLHEP::Hep3Vector::set
void set(double x, double y, double z)
CLHEP::Hep3Vector::setX
void setX(double)
CLHEP::HepLorentzVector::boostVector
Hep3Vector boostVector() const
Definition LorentzVector.cc:173
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition LorentzVector.cc:55
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
CLHEP::HepLorentzVector::set
void set(double x, double y, double z, double t)
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4DynamicParticle::GetMomentum
G4ThreeVector GetMomentum() const
G4RiGeAngularGenerator::PhiRotation
void PhiRotation(G4ThreeVector &dir, G4double phi)
Definition G4RiGeAngularGenerator.cc:355
G4RiGeAngularGenerator::eDP2
G4LorentzVector eDP2(G4double x1, G4double x2, G4double x3, G4double x4, G4double x5)
Definition G4RiGeAngularGenerator.cc:369
G4RiGeAngularGenerator::pDP2
G4LorentzVector pDP2(G4double x3, const G4LorentzVector &x6)
Definition G4RiGeAngularGenerator.cc:392
G4RiGeAngularGenerator::SampleDirection
G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double gEnergy, G4int Z, const G4Material *mat=nullptr) override
Definition G4RiGeAngularGenerator.cc:56
G4RiGeAngularGenerator::PrintGeneratorInformation
void PrintGeneratorInformation() const override
Definition G4RiGeAngularGenerator.cc:400
G4RiGeAngularGenerator::Sample5DPairDirections
void Sample5DPairDirections(const G4DynamicParticle *dp, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, const G4double gEnergy, const G4double q2, const G4double gMomentum, G4double muFinalMomentum, G4double muFinalEnergy, const G4double *randNumbs, const G4double *W)
Definition G4RiGeAngularGenerator.cc:88
G4VEmAngularDistribution::G4VEmAngularDistribution
G4VEmAngularDistribution(const G4String &name)
Definition G4VEmAngularDistribution.cc:55
W
#define W
Definition crc32.c:85