Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4eBremParametrizedModel.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: G4eBremParametrizedModel
33//
34// Author: Andreas Schaelicke
35//
36// Creation date: 06.04.2011
37//
38// Modifications:
39//
40// Main References:
41// - based on G4eBremsstrahlungModel and G4eBremsstrahlungRelModel
42// -------------------------------------------------------------------
43//
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
45//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46
47#include "G4eBremParametrizedModel.hh"
48#include "G4PhysicalConstants.hh"
49#include "G4SystemOfUnits.hh"
50#include "G4Electron.hh"
51#include "G4Positron.hh"
52#include "G4Gamma.hh"
53#include "Randomize.hh"
54#include "G4Material.hh"
55#include "G4Element.hh"
56#include "G4ElementVector.hh"
57#include "G4ProductionCutsTable.hh"
58#include "G4ParticleChangeForLoss.hh"
59#include "G4LossTableManager.hh"
60#include "G4ModifiedTsai.hh"
61#include "G4Exp.hh"
62#include "G4Log.hh"
63
64//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65
66const G4double G4eBremParametrizedModel::xgi[]={ 0.0199, 0.1017, 0.2372, 0.4083,
67 0.5917, 0.7628, 0.8983, 0.9801 };
68const G4double G4eBremParametrizedModel::wgi[]={ 0.0506, 0.1112, 0.1569, 0.1813,
69 0.1813, 0.1569, 0.1112, 0.0506 };
70
71static const G4double tlow = 1.*CLHEP::MeV;
72
73//
74// GEANT4 internal units.
75//
76static const G4double
77 ah10 = 4.67733E+00, ah11 =-6.19012E-01, ah12 = 2.02225E-02,
78 ah20 =-7.34101E+00, ah21 = 1.00462E+00, ah22 =-3.20985E-02,
79 ah30 = 2.93119E+00, ah31 =-4.03761E-01, ah32 = 1.25153E-02;
80
81static const G4double
82 bh10 = 4.23071E+00, bh11 =-6.10995E-01, bh12 = 1.95531E-02,
83 bh20 =-7.12527E+00, bh21 = 9.69160E-01, bh22 =-2.74255E-02,
84 bh30 = 2.69925E+00, bh31 =-3.63283E-01, bh32 = 9.55316E-03;
85
86static const G4double
87 al00 =-2.05398E+00, al01 = 2.38815E-02, al02 = 5.25483E-04,
88 al10 =-7.69748E-02, al11 =-6.91499E-02, al12 = 2.22453E-03,
89 al20 = 4.06463E-02, al21 =-1.01281E-02, al22 = 3.40919E-04;
90
91static const G4double
92 bl00 = 1.04133E+00, bl01 =-9.43291E-03, bl02 =-4.54758E-04,
93 bl10 = 1.19253E-01, bl11 = 4.07467E-02, bl12 =-1.30718E-03,
94 bl20 =-1.59391E-02, bl21 = 7.27752E-03, bl22 =-1.94405E-04;
95
96using namespace std;
97
98G4eBremParametrizedModel::G4eBremParametrizedModel(const G4ParticleDefinition* p,
99 const G4String& nam)
100 : G4VEmModel(nam),
101 particle(0),
102 isElectron(true),
103 fMigdalConstant(classic_electr_radius*electron_Compton_length*electron_Compton_length*4.0*pi),
104 bremFactor(fine_structure_const*classic_electr_radius*classic_electr_radius*16./3.),
105 isInitialised(false)
106{
107 theGamma = G4Gamma::Gamma();
108
109 minThreshold = 0.1*keV;
110 lowKinEnergy = 10.*MeV;
111 SetLowEnergyLimit(lowKinEnergy);
112
113 nist = G4NistManager::Instance();
114
115 SetAngularDistribution(new G4ModifiedTsai());
116
117 particleMass = kinEnergy = totalEnergy = currentZ = z13 = z23 = lnZ = Fel = Finel
118 = densityFactor = densityCorr = fMax = fCoulomb = 0.;
119
120 InitialiseConstants();
121 if(p) { SetParticle(p); }
122}
123
124//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125
126void G4eBremParametrizedModel::InitialiseConstants()
127{
128 facFel = G4Log(184.15);
129 facFinel = G4Log(1194.);
130}
131
132//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
133
134G4eBremParametrizedModel::~G4eBremParametrizedModel()
135{
136}
137
138//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
139
140void G4eBremParametrizedModel::SetParticle(const G4ParticleDefinition* p)
141{
142 particle = p;
143 particleMass = p->GetPDGMass();
144 if(p == G4Electron::Electron()) { isElectron = true; }
145 else { isElectron = false;}
146}
147
148//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
149
150G4double G4eBremParametrizedModel::MinEnergyCut(const G4ParticleDefinition*,
151 const G4MaterialCutsCouple*)
152{
153 return minThreshold;
154}
155
156//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
157
158void G4eBremParametrizedModel::SetupForMaterial(const G4ParticleDefinition*,
159 const G4Material* mat,
160 G4double kineticEnergy)
161{
162 densityFactor = mat->GetElectronDensity()*fMigdalConstant;
163
164 // calculate threshold for density effect
165 kinEnergy = kineticEnergy;
166 totalEnergy = kineticEnergy + particleMass;
167 densityCorr = densityFactor*totalEnergy*totalEnergy;
168}
169
170
171//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
172
173void G4eBremParametrizedModel::Initialise(const G4ParticleDefinition* p,
174 const G4DataVector& cuts)
175{
176 if(p) { SetParticle(p); }
177
178 lowKinEnergy = LowEnergyLimit();
179
180 currentZ = 0.;
181
182 if(IsMaster()) { InitialiseElementSelectors(p, cuts); }
183
184 if(isInitialised) { return; }
185 fParticleChange = GetParticleChangeForLoss();
186 isInitialised = true;
187}
188
189//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
190
191void G4eBremParametrizedModel::InitialiseLocal(const G4ParticleDefinition*,
192 G4VEmModel* masterModel)
193{
194 SetElementSelectors(masterModel->GetElementSelectors());
195}
196
197//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
198
199G4double G4eBremParametrizedModel::ComputeDEDXPerVolume(
200 const G4Material* material,
201 const G4ParticleDefinition* p,
202 G4double kineticEnergy,
203 G4double cutEnergy)
204{
205 if(!particle) { SetParticle(p); }
206 if(kineticEnergy < lowKinEnergy) { return 0.0; }
207 G4double cut = std::min(cutEnergy, kineticEnergy);
208 if(cut == 0.0) { return 0.0; }
209
210 SetupForMaterial(particle, material,kineticEnergy);
211
212 const G4ElementVector* theElementVector = material->GetElementVector();
213 const G4double* theAtomicNumDensityVector = material->GetAtomicNumDensityVector();
214
215 G4double dedx = 0.0;
216
217 // loop for elements in the material
218 for (size_t i=0; i<material->GetNumberOfElements(); i++) {
219
220 G4VEmModel::SetCurrentElement((*theElementVector)[i]);
221 SetCurrentElement((*theElementVector)[i]->GetZ());
222
223 dedx += theAtomicNumDensityVector[i]*currentZ*currentZ*ComputeBremLoss(cut);
224 }
225 dedx *= bremFactor;
226
227 return dedx;
228}
229
230//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
231
232G4double G4eBremParametrizedModel::ComputeBremLoss(G4double cut)
233{
234 G4double loss = 0.0;
235
236 // number of intervals and integration step
237 G4double vcut = cut/totalEnergy;
238 G4int n = (G4int)(20*vcut) + 3;
239 G4double delta = vcut/G4double(n);
240
241 G4double e0 = 0.0;
242 G4double xs;
243
244 // integration
245 for(G4int l=0; l<n; l++) {
246
247 for(G4int i=0; i<8; i++) {
248
249 G4double eg = (e0 + xgi[i]*delta)*totalEnergy;
250
251 xs = ComputeDXSectionPerAtom(eg);
252
253 loss += wgi[i]*xs/(1.0 + densityCorr/(eg*eg));
254 }
255 e0 += delta;
256 }
257
258 loss *= delta*totalEnergy;
259
260 return loss;
261}
262
263//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
264
265G4double G4eBremParametrizedModel::ComputeCrossSectionPerAtom(
266 const G4ParticleDefinition* p,
267 G4double kineticEnergy,
268 G4double Z, G4double,
269 G4double cutEnergy,
270 G4double maxEnergy)
271{
272 if(!particle) { SetParticle(p); }
273 if(kineticEnergy < lowKinEnergy) { return 0.0; }
274 G4double cut = std::min(cutEnergy, kineticEnergy);
275 G4double tmax = std::min(maxEnergy, kineticEnergy);
276
277 if(cut >= tmax) { return 0.0; }
278
279 SetCurrentElement(Z);
280
281 G4double cross = ComputeXSectionPerAtom(cut);
282
283 // allow partial integration
284 if(tmax < kinEnergy) { cross -= ComputeXSectionPerAtom(tmax); }
285
286 cross *= Z*Z*bremFactor;
287
288 return cross;
289}
290
291//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
292
293G4double G4eBremParametrizedModel::ComputeXSectionPerAtom(G4double cut)
294{
295 G4double cross = 0.0;
296
297 // number of intervals and integration step
298 G4double vcut = G4Log(cut/totalEnergy);
299 G4double vmax = G4Log(kinEnergy/totalEnergy);
300 G4int n = (G4int)(0.45*(vmax - vcut)) + 4;
301 // n=1; // integration test
302 G4double delta = (vmax - vcut)/G4double(n);
303
304 G4double e0 = vcut;
305 G4double xs;
306
307 // integration
308 for(G4int l=0; l<n; l++) {
309
310 for(G4int i=0; i<8; i++) {
311
312 G4double eg = G4Exp(e0 + xgi[i]*delta)*totalEnergy;
313
314 xs = ComputeDXSectionPerAtom(eg);
315
316 cross += wgi[i]*xs/(1.0 + densityCorr/(eg*eg));
317 }
318 e0 += delta;
319 }
320
321 cross *= delta;
322
323 return cross;
324}
325
326//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
327
328// compute the value of the screening function 3*PHI1 - PHI2
329
330G4double G4eBremParametrizedModel::ScreenFunction1(G4double ScreenVariable)
331{
332 G4double screenVal;
333
334 if (ScreenVariable > 1.)
335 screenVal = 42.24 - 8.368*G4Log(ScreenVariable+0.952);
336 else
337 screenVal = 42.392 - ScreenVariable* (7.796 - 1.961*ScreenVariable);
338
339 return screenVal;
340}
341
342//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
343
344// compute the value of the screening function 1.5*PHI1 - 0.5*PHI2
345
346G4double G4eBremParametrizedModel::ScreenFunction2(G4double ScreenVariable)
347{
348 G4double screenVal;
349
350 if (ScreenVariable > 1.)
351 screenVal = 42.24 - 8.368*G4Log(ScreenVariable+0.952);
352 else
353 screenVal = 41.734 - ScreenVariable* (6.484 - 1.250*ScreenVariable);
354
355 return screenVal;
356}
357
358//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
359
360// Parametrized cross section
361G4double G4eBremParametrizedModel::ComputeParametrizedDXSectionPerAtom(
362 G4double kineticEnergy,
363 G4double gammaEnergy, G4double Z)
364{
365 SetCurrentElement(Z);
366 G4double FZ = lnZ* (4.- 0.55*lnZ);
367 G4double Z3 = z13;
368 G4double ZZ = z13*nist->GetZ13(G4lrint(Z)+1);
369
370 totalEnergy = kineticEnergy + electron_mass_c2;
371
372 // G4double x, epsil, greject, migdal, grejmax, q;
373 G4double epsil, greject;
374 G4double U = G4Log(kineticEnergy/electron_mass_c2);
375 G4double U2 = U*U;
376
377 // precalculated parameters
378 G4double ah, bh;
379
380 if (kineticEnergy > tlow) {
381
382 G4double ah1 = ah10 + ZZ* (ah11 + ZZ* ah12);
383 G4double ah2 = ah20 + ZZ* (ah21 + ZZ* ah22);
384 G4double ah3 = ah30 + ZZ* (ah31 + ZZ* ah32);
385
386 G4double bh1 = bh10 + ZZ* (bh11 + ZZ* bh12);
387 G4double bh2 = bh20 + ZZ* (bh21 + ZZ* bh22);
388 G4double bh3 = bh30 + ZZ* (bh31 + ZZ* bh32);
389
390 ah = 1. + (ah1*U2 + ah2*U + ah3) / (U2*U);
391 bh = 0.75 + (bh1*U2 + bh2*U + bh3) / (U2*U);
392
393 // limit of the screening variable
394 G4double screenfac =
395 136.*electron_mass_c2/(Z3*totalEnergy);
396
397 epsil = gammaEnergy/totalEnergy; // epsil = x*kineticEnergy/totalEnergy;
398 G4double screenvar = screenfac*epsil/(1.0-epsil);
399 G4double F1 = max(ScreenFunction1(screenvar) - FZ ,0.);
400 G4double F2 = max(ScreenFunction2(screenvar) - FZ ,0.);
401
402
403 greject = (F1 - epsil* (ah*F1 - bh*epsil*F2))/8.; // 1./(42.392 - FZ);
404
405 std::cout << " yy = "<<epsil<<std::endl;
406 std::cout << " F1/(...) "<<F1/(42.392 - FZ)<<std::endl;
407 std::cout << " F2/(...) "<<F2/(42.392 - FZ)<<std::endl;
408 std::cout << " (42.392 - FZ) " << (42.392 - FZ) <<std::endl;
409
410 } else {
411
412 G4double al0 = al00 + ZZ* (al01 + ZZ* al02);
413 G4double al1 = al10 + ZZ* (al11 + ZZ* al12);
414 G4double al2 = al20 + ZZ* (al21 + ZZ* al22);
415
416 G4double bl0 = bl00 + ZZ* (bl01 + ZZ* bl02);
417 G4double bl1 = bl10 + ZZ* (bl11 + ZZ* bl12);
418 G4double bl2 = bl20 + ZZ* (bl21 + ZZ* bl22);
419
420 ah = al0 + al1*U + al2*U2;
421 bh = bl0 + bl1*U + bl2*U2;
422
423 G4double x=gammaEnergy/kineticEnergy;
424 greject=(1. + x* (ah + bh*x));
425
426 /*
427 // Compute the maximum of the rejection function
428 grejmax = max(1. + xmin* (ah + bh*xmin), 1.+ah+bh);
429 G4double xm = -ah/(2.*bh);
430 if ( xmin < xm && xm < xmax) grejmax = max(grejmax, 1.+ xm* (ah + bh*xm));
431 */
432 }
433
434 return greject;
435}
436
437//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
438
439G4double G4eBremParametrizedModel::ComputeDXSectionPerAtom(G4double gammaEnergy)
440{
441
442 if(gammaEnergy < 0.0) { return 0.0; }
443
444 G4double y = gammaEnergy/totalEnergy;
445
446 G4double main=0.;
447 //secondTerm=0.;
448
449 // ** form factors complete screening case **
450 // only valid for high energies (and if LPM suppression does not play a role)
451 main = (3./4.*y*y - y + 1.) * ( (Fel-fCoulomb) + Finel/currentZ );
452 // secondTerm = (1.-y)/12.*(1.+1./currentZ);
453
454 std::cout<<" F1(0) "<<ScreenFunction1(0.) <<std::endl;
455 std::cout<<" F1(0) "<<ScreenFunction2(0.) <<std::endl;
456 std::cout<<"Ekin = "<<kinEnergy<<std::endl;
457 std::cout<<"Z = "<<currentZ<<std::endl;
458 std::cout<<"main = "<<main<<std::endl;
459 std::cout<<" y = "<<y<<std::endl;
460 std::cout<<" Fel-fCoulomb "<< (Fel-fCoulomb) <<std::endl;
461
462 G4double main2 = ComputeParametrizedDXSectionPerAtom(kinEnergy,gammaEnergy,currentZ);
463 std::cout<<"main2 = "<<main2<<std::endl;
464 std::cout<<"main2tot = "<<main2 * ( (Fel-fCoulomb) + Finel/currentZ )/(Fel-fCoulomb);
465
466 G4double cross = main2; //main+secondTerm;
467 return cross;
468}
469
470//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
471
472void G4eBremParametrizedModel::SampleSecondaries(
473 std::vector<G4DynamicParticle*>* vdp,
474 const G4MaterialCutsCouple* couple,
475 const G4DynamicParticle* dp,
476 G4double cutEnergy,
477 G4double maxEnergy)
478{
479 G4double kineticEnergy = dp->GetKineticEnergy();
480 if(kineticEnergy < lowKinEnergy) { return; }
481 G4double cut = std::min(cutEnergy, kineticEnergy);
482 G4double emax = std::min(maxEnergy, kineticEnergy);
483 if(cut >= emax) { return; }
484
485 SetupForMaterial(particle, couple->GetMaterial(),kineticEnergy);
486
487 const G4Element* elm = SelectTargetAtom(couple,particle,kineticEnergy,
488 dp->GetLogKineticEnergy(),cut,emax);
489 SetCurrentElement(elm->GetZ());
490
491 kinEnergy = kineticEnergy;
492 totalEnergy = kineticEnergy + particleMass;
493 densityCorr = densityFactor*totalEnergy*totalEnergy;
494
495 G4double xmin = G4Log(cut*cut + densityCorr);
496 G4double xmax = G4Log(emax*emax + densityCorr);
497 G4double gammaEnergy, f, x;
498
499 CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
500
501 do {
502 x = G4Exp(xmin + rndmEngine->flat()*(xmax - xmin)) - densityCorr;
503 if(x < 0.0) x = 0.0;
504 gammaEnergy = sqrt(x);
505 f = ComputeDXSectionPerAtom(gammaEnergy);
506
507 if ( f > fMax ) {
508 G4cout << "### G4eBremParametrizedModel Warning: Majoranta exceeded! "
509 << f << " > " << fMax
510 << " Egamma(MeV)= " << gammaEnergy
511 << " E(mEV)= " << kineticEnergy
512 << G4endl;
513 }
514
515 // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
516 } while (f < fMax*rndmEngine->flat());
517
518 //
519 // angles of the emitted gamma. ( Z - axis along the parent particle)
520 // use general interface
521 //
522 G4ThreeVector gammaDirection =
523 GetAngularDistribution()->SampleDirection(dp, totalEnergy-gammaEnergy,
524 G4lrint(currentZ),
525 couple->GetMaterial());
526
527 // create G4DynamicParticle object for the Gamma
528 G4DynamicParticle* gamma = new G4DynamicParticle(theGamma,gammaDirection,
529 gammaEnergy);
530 vdp->push_back(gamma);
531
532 G4double totMomentum = sqrt(kineticEnergy*(totalEnergy + electron_mass_c2));
533 G4ThreeVector direction = (totMomentum*dp->GetMomentumDirection()
534 - gammaEnergy*gammaDirection).unit();
535
536 // energy of primary
537 G4double finalE = kineticEnergy - gammaEnergy;
538
539 // stop tracking and create new secondary instead of primary
540 if(gammaEnergy > SecondaryThreshold()) {
541 fParticleChange->ProposeTrackStatus(fStopAndKill);
542 fParticleChange->SetProposedKineticEnergy(0.0);
543 G4DynamicParticle* el =
544 new G4DynamicParticle(const_cast<G4ParticleDefinition*>(particle),
545 direction, finalE);
546 vdp->push_back(el);
547
548 // continue tracking
549 } else {
550 fParticleChange->SetProposedMomentumDirection(direction);
551 fParticleChange->SetProposedKineticEnergy(finalE);
552 }
553}
554
555//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
556
557
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:43
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:179
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:226
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::HepRandomEngine::flat
virtual double flat()=0
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetLogKineticEnergy
G4double GetLogKineticEnergy() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:93
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:130
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:85
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:166
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:188
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:184
G4Material::GetAtomicNumDensityVector
const G4double * GetAtomicNumDensityVector() const
Definition: G4Material.hh:214
G4Material::GetElectronDensity
G4double GetElectronDensity() const
Definition: G4Material.hh:215
G4NistManager::GetZ13
G4double GetZ13(G4double Z) const
Definition: G4NistManager.hh:589
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:70
G4ParticleChangeForLoss::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForLoss.hh:135
G4ParticleChangeForLoss::SetProposedMomentumDirection
void SetProposedMomentumDirection(const G4ThreeVector &dir)
Definition: G4ParticleChangeForLoss.hh:185
G4VEmAngularDistribution::SampleDirection
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
G4VEmModel::SetElementSelectors
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
Definition: G4VEmModel.hh:842
G4VEmModel::GetAngularDistribution
G4VEmAngularDistribution * GetAngularDistribution()
Definition: G4VEmModel.hh:611
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:652
G4VEmModel::GetElementSelectors
std::vector< G4EmElementSelector * > * GetElementSelectors()
Definition: G4VEmModel.hh:834
G4VEmModel::IsMaster
G4bool IsMaster() const
Definition: G4VEmModel.hh:736
G4VEmModel::SetCurrentElement
void SetCurrentElement(const G4Element *)
Definition: G4VEmModel.hh:487
G4VEmModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:764
G4VEmModel::SetAngularDistribution
void SetAngularDistribution(G4VEmAngularDistribution *)
Definition: G4VEmModel.hh:618
G4VEmModel::SelectTargetAtom
const G4Element * SelectTargetAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:587
G4VEmModel::InitialiseElementSelectors
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
Definition: G4VEmModel.cc:148
G4VEmModel::SecondaryThreshold
G4double SecondaryThreshold() const
Definition: G4VEmModel.hh:680
G4VEmModel::GetParticleChangeForLoss
G4ParticleChangeForLoss * GetParticleChangeForLoss()
Definition: G4VEmModel.cc:118
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4eBremParametrizedModel::InitialiseLocal
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel) override
Definition: G4eBremParametrizedModel.cc:191
G4eBremParametrizedModel::wgi
static const G4double wgi[8]
Definition: G4eBremParametrizedModel.hh:132
G4eBremParametrizedModel::ComputeDEDXPerVolume
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
Definition: G4eBremParametrizedModel.cc:199
G4eBremParametrizedModel::MinEnergyCut
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *) override
Definition: G4eBremParametrizedModel.cc:150
G4eBremParametrizedModel::fParticleChange
G4ParticleChangeForLoss * fParticleChange
Definition: G4eBremParametrizedModel.hh:130
G4eBremParametrizedModel::particle
const G4ParticleDefinition * particle
Definition: G4eBremParametrizedModel.hh:128
G4eBremParametrizedModel::SampleSecondaries
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double cutEnergy, G4double maxEnergy) override
Definition: G4eBremParametrizedModel.cc:472
G4eBremParametrizedModel::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double tkin, G4double Z, G4double, G4double cutEnergy, G4double maxEnergy=DBL_MAX) override
Definition: G4eBremParametrizedModel.cc:265
G4eBremParametrizedModel::G4eBremParametrizedModel
G4eBremParametrizedModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="eBremParam")
Definition: G4eBremParametrizedModel.cc:98
G4eBremParametrizedModel::xgi
static const G4double xgi[8]
Definition: G4eBremParametrizedModel.hh:132
G4eBremParametrizedModel::SetupForMaterial
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double) override
Definition: G4eBremParametrizedModel.cc:158
G4eBremParametrizedModel::Initialise
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
Definition: G4eBremParametrizedModel.cc:173
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition: G4INCLGlobals.hh:112
G4lrint
int G4lrint(double ad)
Definition: templates.hh:134