Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4EmCalculator.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// $Id$
27//
28// -------------------------------------------------------------------
29//
30// GEANT4 Class file
31//
32//
33// File name: G4EmCalculator
34//
35// Author: Vladimir Ivanchenko
36//
37// Creation date: 28.06.2004
38//
39// Modifications:
40// 12.09.2004 Add verbosity (V.Ivanchenko)
41// 17.11.2004 Change signature of methods, add new methods (V.Ivanchenko)
42// 08.04.2005 Major optimisation of internal interfaces (V.Ivantchenko)
43// 08.05.2005 Use updated interfaces (V.Ivantchenko)
44// 23.10.2005 Fix computations for ions (V.Ivantchenko)
45// 11.01.2006 Add GetCSDARange (V.Ivantchenko)
46// 26.01.2006 Rename GetRange -> GetRangeFromRestricteDEDX (V.Ivanchenko)
47// 14.03.2006 correction in GetCrossSectionPerVolume (mma)
48// suppress GetCrossSectionPerAtom
49// elm->GetA() in ComputeCrossSectionPerAtom
50// 22.03.2006 Add ComputeElectronicDEDX and ComputeTotalDEDX (V.Ivanchenko)
51// 13.05.2006 Add Corrections for ion stopping (V.Ivanchenko)
52// 29.09.2006 Uncomment computation of smoothing factor (V.Ivanchenko)
53// 27.10.2006 Change test energy to access lowEnergy model from
54// 10 keV to 1 keV (V. Ivanchenko)
55// 15.03.2007 Add ComputeEnergyCutFromRangeCut methods (V.Ivanchenko)
56// 21.04.2008 Updated computations for ions (V.Ivanchenko)
57//
58// Class Description:
59//
60// -------------------------------------------------------------------
61//
62//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
63//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
64
65#include "G4EmCalculator.hh"
66#include "G4SystemOfUnits.hh"
67#include "G4LossTableManager.hh"
68#include "G4VEmProcess.hh"
69#include "G4VEnergyLossProcess.hh"
70#include "G4VMultipleScattering.hh"
71#include "G4Material.hh"
72#include "G4MaterialCutsCouple.hh"
73#include "G4ParticleDefinition.hh"
74#include "G4ParticleTable.hh"
75#include "G4PhysicsTable.hh"
76#include "G4ProductionCutsTable.hh"
77#include "G4ProcessManager.hh"
78#include "G4ionEffectiveCharge.hh"
79#include "G4RegionStore.hh"
80#include "G4Element.hh"
81#include "G4EmCorrections.hh"
82#include "G4GenericIon.hh"
83#include "G4ProcessVector.hh"
84#include "G4Gamma.hh"
85
86//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
87
88G4EmCalculator::G4EmCalculator()
89{
90 manager = G4LossTableManager::Instance();
91 corr = manager->EmCorrections();
92 nLocalMaterials = 0;
93 verbose = 0;
94 currentCoupleIndex = 0;
95 currentCouple = 0;
96 currentMaterial = 0;
97 currentParticle = 0;
98 lambdaParticle = 0;
99 baseParticle = 0;
100 currentLambda = 0;
101 currentModel = 0;
102 currentProcess = 0;
103 loweModel = 0;
104 chargeSquare = 1.0;
105 massRatio = 1.0;
106 mass = 0.0;
107 currentCut = 0.0;
108 currentParticleName= "";
109 currentMaterialName= "";
110 currentName = "";
111 lambdaName = "";
112 theGenericIon = G4GenericIon::GenericIon();
113 ionEffCharge = new G4ionEffectiveCharge();
114 isIon = false;
115 isApplicable = false;
116}
117
118//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
119
120G4EmCalculator::~G4EmCalculator()
121{
122 delete ionEffCharge;
123 for (G4int i=0; i<nLocalMaterials; ++i) {
124 delete localCouples[i];
125 }
126}
127
128//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
129
130G4double G4EmCalculator::GetDEDX(G4double kinEnergy, const G4ParticleDefinition* p,
131 const G4Material* mat, const G4Region* region)
132{
133 G4double res = 0.0;
134 const G4MaterialCutsCouple* couple = FindCouple(mat, region);
135 if(couple && UpdateParticle(p, kinEnergy) ) {
136 res = manager->GetDEDX(p, kinEnergy, couple);
137
138 if(isIon) {
139 if(FindEmModel(p, currentProcessName, kinEnergy)) {
140 G4double length = CLHEP::nm;
141 G4double eloss = res*length;
142 //G4cout << "### GetDEDX: E= " << kinEnergy << " dedx0= " << res
143 // << " de= " << eloss << G4endl;;
144 G4double niel = 0.0;
145 dynParticle.SetKineticEnergy(kinEnergy);
146 currentModel->GetChargeSquareRatio(p, mat, kinEnergy);
147 currentModel->CorrectionsAlongStep(couple,&dynParticle,eloss,niel,length);
148 res = eloss/length;
149 //G4cout << " de1= " << eloss << " res1= " << res
150 // << " " << p->GetParticleName() <<G4endl;;
151 }
152 }
153
154 if(verbose>0) {
155 G4cout << "G4EmCalculator::GetDEDX: E(MeV)= " << kinEnergy/MeV
156 << " DEDX(MeV/mm)= " << res*mm/MeV
157 << " DEDX(MeV*cm^2/g)= " << res*gram/(MeV*cm2*mat->GetDensity())
158 << " " << p->GetParticleName()
159 << " in " << mat->GetName()
160 << " isIon= " << isIon
161 << G4endl;
162 }
163 }
164 return res;
165}
166
167//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
168
169G4double G4EmCalculator::GetDEDX(G4double kinEnergy, const G4String& particle,
170 const G4String& material, const G4String& reg)
171{
172 return GetDEDX(kinEnergy,FindParticle(particle),
173 FindMaterial(material),FindRegion(reg));
174}
175
176//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
177
178G4double G4EmCalculator::GetRangeFromRestricteDEDX(G4double kinEnergy,
179 const G4ParticleDefinition* p,
180 const G4Material* mat,
181 const G4Region* region)
182{
183 G4double res = 0.0;
184 const G4MaterialCutsCouple* couple = FindCouple(mat,region);
185 if(couple && UpdateParticle(p, kinEnergy)) {
186 res = manager->GetRangeFromRestricteDEDX(p, kinEnergy, couple);
187 if(verbose>0) {
188 G4cout << "G4EmCalculator::GetRange: E(MeV)= " << kinEnergy/MeV
189 << " range(mm)= " << res/mm
190 << " " << p->GetParticleName()
191 << " in " << mat->GetName()
192 << G4endl;
193 }
194 }
195 return res;
196}
197
198//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
199
200G4double G4EmCalculator::GetCSDARange(G4double kinEnergy,
201 const G4ParticleDefinition* p,
202 const G4Material* mat,
203 const G4Region* region)
204{
205 G4double res = 0.0;
206 const G4MaterialCutsCouple* couple = FindCouple(mat,region);
207 if(couple && UpdateParticle(p, kinEnergy)) {
208 res = manager->GetCSDARange(p, kinEnergy, couple);
209 if(verbose>0) {
210 G4cout << "G4EmCalculator::GetRange: E(MeV)= " << kinEnergy/MeV
211 << " range(mm)= " << res/mm
212 << " " << p->GetParticleName()
213 << " in " << mat->GetName()
214 << G4endl;
215 }
216 }
217 return res;
218}
219
220//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
221
222G4double G4EmCalculator::GetRange(G4double kinEnergy,
223 const G4ParticleDefinition* p,
224 const G4Material* mat,
225 const G4Region* region)
226{
227 G4double res = 0.0;
228 const G4MaterialCutsCouple* couple = FindCouple(mat,region);
229 if(couple && UpdateParticle(p, kinEnergy)) {
230 res = manager->GetRange(p, kinEnergy, couple);
231 if(verbose>0) {
232 G4cout << "G4EmCalculator::GetRange: E(MeV)= " << kinEnergy/MeV
233 << " range(mm)= " << res/mm
234 << " " << p->GetParticleName()
235 << " in " << mat->GetName()
236 << G4endl;
237 }
238 }
239 return res;
240}
241
242//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
243
244G4double G4EmCalculator::GetRangeFromRestricteDEDX(G4double kinEnergy,
245 const G4String& particle,
246 const G4String& material,
247 const G4String& reg)
248{
249 return GetRangeFromRestricteDEDX(kinEnergy,FindParticle(particle),
250 FindMaterial(material),FindRegion(reg));
251}
252
253//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
254
255G4double G4EmCalculator::GetCSDARange(G4double kinEnergy,
256 const G4String& particle,
257 const G4String& material,
258 const G4String& reg)
259{
260 return GetCSDARange(kinEnergy,FindParticle(particle),
261 FindMaterial(material),FindRegion(reg));
262}
263
264//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
265
266G4double G4EmCalculator::GetRange(G4double kinEnergy,
267 const G4String& particle,
268 const G4String& material,
269 const G4String& reg)
270{
271 return GetRange(kinEnergy,FindParticle(particle),
272 FindMaterial(material),FindRegion(reg));
273}
274
275//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
276
277G4double G4EmCalculator::GetKinEnergy(G4double range,
278 const G4ParticleDefinition* p,
279 const G4Material* mat,
280 const G4Region* region)
281{
282 G4double res = 0.0;
283 const G4MaterialCutsCouple* couple = FindCouple(mat,region);
284 if(couple && UpdateParticle(p, 1.0*GeV)) {
285 res = manager->GetEnergy(p, range, couple);
286 if(verbose>0) {
287 G4cout << "G4EmCalculator::GetKinEnergy: Range(mm)= " << range/mm
288 << " KinE(MeV)= " << res/MeV
289 << " " << p->GetParticleName()
290 << " in " << mat->GetName()
291 << G4endl;
292 }
293 }
294 return res;
295}
296
297//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
298
299G4double G4EmCalculator::GetKinEnergy(G4double range, const G4String& particle,
300 const G4String& material, const G4String& reg)
301{
302 return GetKinEnergy(range,FindParticle(particle),
303 FindMaterial(material),FindRegion(reg));
304}
305
306//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
307
308G4double G4EmCalculator::GetCrossSectionPerVolume(G4double kinEnergy,
309 const G4ParticleDefinition* p,
310 const G4String& processName,
311 const G4Material* mat,
312 const G4Region* region)
313{
314 G4double res = 0.0;
315 const G4MaterialCutsCouple* couple = FindCouple(mat,region);
316
317 if(couple && UpdateParticle(p, kinEnergy)) {
318 G4int idx = couple->GetIndex();
319 FindLambdaTable(p, processName, kinEnergy);
320
321 if(currentLambda) {
322 G4double e = kinEnergy*massRatio;
323 res = (((*currentLambda)[idx])->Value(e))*chargeSquare;
324 } else {
325 res = ComputeCrossSectionPerVolume(kinEnergy, p, processName, mat,
326 kinEnergy);
327 }
328 if(verbose>0) {
329 G4cout << "G4EmCalculator::GetXSPerVolume: E(MeV)= " << kinEnergy/MeV
330 << " cross(cm-1)= " << res*cm
331 << " " << p->GetParticleName()
332 << " in " << mat->GetName();
333 if(verbose>1)
334 G4cout << " idx= " << idx << " Escaled((MeV)= "
335 << kinEnergy*massRatio
336 << " q2= " << chargeSquare;
337 G4cout << G4endl;
338 }
339 }
340 return res;
341}
342
343//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
344
345G4double G4EmCalculator::GetCrossSectionPerVolume(G4double kinEnergy,
346 const G4String& particle,
347 const G4String& processName,
348 const G4String& material,
349 const G4String& reg)
350{
351 return GetCrossSectionPerVolume(kinEnergy,FindParticle(particle),processName,
352 FindMaterial(material),FindRegion(reg));
353}
354
355//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
356
357G4double G4EmCalculator::GetShellIonisationCrossSectionPerAtom(
358 const G4String& particle,
359 G4int Z,
360 G4AtomicShellEnumerator shell,
361 G4double kinEnergy)
362{
363 G4double res = 0.0;
364 const G4ParticleDefinition* p = FindParticle(particle);
365 G4VAtomDeexcitation* ad = manager->AtomDeexcitation();
366 if(p && ad) {
367 res = ad->GetShellIonisationCrossSectionPerAtom(p, Z, shell, kinEnergy);
368 }
369 return res;
370}
371
372//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
373
374G4double G4EmCalculator::GetMeanFreePath(G4double kinEnergy,
375 const G4ParticleDefinition* p,
376 const G4String& processName,
377 const G4Material* mat,
378 const G4Region* region)
379{
380 G4double res = DBL_MAX;
381 G4double x = GetCrossSectionPerVolume(kinEnergy,p, processName, mat,region);
382 if(x > 0.0) { res = 1.0/x; }
383 if(verbose>1) {
384 G4cout << "G4EmCalculator::GetMeanFreePath: E(MeV)= " << kinEnergy/MeV
385 << " MFP(mm)= " << res/mm
386 << " " << p->GetParticleName()
387 << " in " << mat->GetName()
388 << G4endl;
389 }
390 return res;
391}
392
393//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
394
395G4double G4EmCalculator::GetMeanFreePath(G4double kinEnergy,
396 const G4String& particle,
397 const G4String& processName,
398 const G4String& material,
399 const G4String& reg)
400{
401 return GetMeanFreePath(kinEnergy,FindParticle(particle),processName,
402 FindMaterial(material),FindRegion(reg));
403}
404
405//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
406
407void G4EmCalculator::PrintDEDXTable(const G4ParticleDefinition* p)
408{
409 const G4VEnergyLossProcess* elp = FindEnergyLossProcess(p);
410 G4cout << "##### DEDX Table for " << p->GetParticleName() << G4endl;
411 if(elp) G4cout << *(elp->DEDXTable()) << G4endl;
412}
413
414//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
415
416void G4EmCalculator::PrintRangeTable(const G4ParticleDefinition* p)
417{
418 const G4VEnergyLossProcess* elp = FindEnergyLossProcess(p);
419 G4cout << "##### Range Table for " << p->GetParticleName() << G4endl;
420 if(elp) G4cout << *(elp->RangeTableForLoss()) << G4endl;
421}
422
423//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
424
425void G4EmCalculator::PrintInverseRangeTable(const G4ParticleDefinition* p)
426{
427 const G4VEnergyLossProcess* elp = FindEnergyLossProcess(p);
428 G4cout << "### G4EmCalculator: Inverse Range Table for "
429 << p->GetParticleName() << G4endl;
430 if(elp) G4cout << *(elp->InverseRangeTable()) << G4endl;
431}
432
433//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
434
435G4double G4EmCalculator::ComputeDEDX(G4double kinEnergy,
436 const G4ParticleDefinition* p,
437 const G4String& processName,
438 const G4Material* mat,
439 G4double cut)
440{
441 currentMaterial = mat;
442 currentMaterialName = mat->GetName();
443 G4double res = 0.0;
444 if(verbose > 1) {
445 G4cout << "### G4EmCalculator::ComputeDEDX: " << p->GetParticleName()
446 << " in " << currentMaterialName
447 << " e(MeV)= " << kinEnergy/MeV << " cut(MeV)= " << cut/MeV
448 << G4endl;
449 }
450 if(UpdateParticle(p, kinEnergy)) {
451 if(FindEmModel(p, processName, kinEnergy)) {
452 G4double escaled = kinEnergy*massRatio;
453 if(baseParticle) {
454 res = currentModel->ComputeDEDXPerVolume(
455 mat, baseParticle, escaled, cut) * chargeSquare;
456 if(verbose > 1) {
457 G4cout << baseParticle->GetParticleName()
458 << " Escaled(MeV)= " << escaled;
459 }
460 } else {
461 res = currentModel->ComputeDEDXPerVolume(mat, p, kinEnergy, cut);
462 if(verbose > 1) { G4cout << " no basePart E(MeV)= " << kinEnergy << " "; }
463 }
464 if(verbose > 1) {
465 G4cout << currentModel->GetName() << ": DEDX(MeV/mm)= " << res*mm/MeV
466 << " DEDX(MeV*cm^2/g)= "
467 << res*gram/(MeV*cm2*mat->GetDensity())
468 << G4endl;
469 }
470
471 // emulate smoothing procedure
472 G4double eth = currentModel->LowEnergyLimit();
473 // G4cout << "massRatio= " << massRatio << " eth= " << eth << G4endl;
474 if(loweModel) {
475 G4double res0 = 0.0;
476 G4double res1 = 0.0;
477 if(baseParticle) {
478 res1 = currentModel->ComputeDEDXPerVolume(mat, baseParticle, eth, cut)
479 * chargeSquare;
480 res0 = loweModel->ComputeDEDXPerVolume(mat, baseParticle, eth, cut)
481 * chargeSquare;
482 } else {
483 res1 = currentModel->ComputeDEDXPerVolume(mat, p, eth, cut);
484 res0 = loweModel->ComputeDEDXPerVolume(mat, p, eth, cut);
485 }
486 if(verbose > 1) {
487 G4cout << "At boundary energy(MeV)= " << eth/MeV
488 << " DEDX(MeV/mm)= " << res1*mm/MeV
489 << G4endl;
490 }
491
492 //G4cout << "eth= " << eth << " escaled= " << escaled
493 // << " res0= " << res0 << " res1= "
494 // << res1 << " q2= " << chargeSquare << G4endl;
495
496 if(res1 > 0.0 && escaled > 0.0) {
497 res *= (1.0 + (res0/res1 - 1.0)*eth/escaled);
498 }
499 }
500
501 // low energy correction for ions
502 if(isIon) {
503 G4double length = CLHEP::nm;
504 const G4Region* r = 0;
505 const G4MaterialCutsCouple* couple = FindCouple(mat, r);
506 G4double eloss = res*length;
507 G4double niel = 0.0;
508 dynParticle.SetKineticEnergy(kinEnergy);
509 currentModel->GetChargeSquareRatio(p, mat, kinEnergy);
510 currentModel->CorrectionsAlongStep(couple,&dynParticle,eloss,niel,length);
511 res = eloss/length;
512
513 if(verbose > 1) {
514 G4cout << "After Corrections: DEDX(MeV/mm)= " << res*mm/MeV
515 << " DEDX(MeV*cm^2/g)= " << res*gram/(MeV*cm2*mat->GetDensity())
516 << G4endl;
517 }
518 }
519 }
520
521 if(verbose > 0) {
522 G4cout << "Sum: E(MeV)= " << kinEnergy/MeV
523 << " DEDX(MeV/mm)= " << res*mm/MeV
524 << " DEDX(MeV*cm^2/g)= " << res*gram/(MeV*cm2*mat->GetDensity())
525 << " cut(MeV)= " << cut/MeV
526 << " " << p->GetParticleName()
527 << " in " << currentMaterialName
528 << " Zi^2= " << chargeSquare
529 << " isIon=" << isIon
530 << G4endl;
531 }
532 }
533 return res;
534}
535
536//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
537
538G4double G4EmCalculator::ComputeElectronicDEDX(G4double kinEnergy,
539 const G4ParticleDefinition* part,
540 const G4Material* mat,
541 G4double cut)
542{
543 currentMaterial = mat;
544 currentMaterialName = mat->GetName();
545 G4double dedx = 0.0;
546 if(UpdateParticle(part, kinEnergy)) {
547
548 G4LossTableManager* lManager = G4LossTableManager::Instance();
549 const std::vector<G4VEnergyLossProcess*> vel =
550 lManager->GetEnergyLossProcessVector();
551 G4int n = vel.size();
552
553 //G4cout << "ComputeElectronicDEDX for " << part->GetParticleName()
554 // << " n= " << n << G4endl;
555
556 for(G4int i=0; i<n; ++i) {
557 if(vel[i]) {
558 G4VProcess* p = reinterpret_cast<G4VProcess*>(vel[i]);
559 if(ActiveForParticle(part, p)) {
560 //G4cout << "idx= " << i << " " << (vel[i])->GetProcessName()
561 // << " " << (vel[i])->Particle()->GetParticleName() << G4endl;
562 dedx += ComputeDEDX(kinEnergy,part,(vel[i])->GetProcessName(),mat,cut);
563 }
564 }
565 }
566 }
567 return dedx;
568}
569
570//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
571
572G4double G4EmCalculator::ComputeElectronicDEDX(G4double kinEnergy, const G4String& part,
573 const G4String& mat, G4double cut)
574{
575 return ComputeElectronicDEDX(kinEnergy,FindParticle(part),FindMaterial(mat),cut);
576}
577
578//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
579
580G4double G4EmCalculator::ComputeTotalDEDX(G4double kinEnergy,
581 const G4ParticleDefinition* part,
582 const G4Material* mat,
583 G4double cut)
584{
585 G4double dedx = ComputeElectronicDEDX(kinEnergy,part,mat,cut);
586 if(mass > 700.*MeV) { dedx += ComputeNuclearDEDX(kinEnergy,part,mat); }
587 return dedx;
588}
589
590//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
591
592G4double G4EmCalculator::ComputeTotalDEDX(G4double kinEnergy,
593 const G4String& part,
594 const G4String& mat,
595 G4double cut)
596{
597 return ComputeTotalDEDX(kinEnergy,FindParticle(part),FindMaterial(mat),cut);
598}
599
600//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
601
602G4double G4EmCalculator::ComputeDEDX(G4double kinEnergy,
603 const G4String& particle,
604 const G4String& processName,
605 const G4String& material,
606 G4double cut)
607{
608 return ComputeDEDX(kinEnergy,FindParticle(particle),processName,
609 FindMaterial(material),cut);
610}
611
612//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
613
614G4double G4EmCalculator::ComputeNuclearDEDX(G4double kinEnergy,
615 const G4ParticleDefinition* p,
616 const G4Material* mat)
617{
618
619 G4double res = corr->NuclearDEDX(p, mat, kinEnergy, false);
620
621 if(verbose > 1) {
622 G4cout << p->GetParticleName() << " E(MeV)= " << kinEnergy/MeV
623 << " NuclearDEDX(MeV/mm)= " << res*mm/MeV
624 << " NuclearDEDX(MeV*cm^2/g)= "
625 << res*gram/(MeV*cm2*mat->GetDensity())
626 << G4endl;
627 }
628 return res;
629}
630
631//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
632
633G4double G4EmCalculator::ComputeNuclearDEDX(G4double kinEnergy,
634 const G4String& particle,
635 const G4String& material)
636{
637 return ComputeNuclearDEDX(kinEnergy,FindParticle(particle),
638 FindMaterial(material));
639}
640
641//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
642
643G4double G4EmCalculator::ComputeCrossSectionPerVolume(
644 G4double kinEnergy,
645 const G4ParticleDefinition* p,
646 const G4String& processName,
647 const G4Material* mat,
648 G4double cut)
649{
650 currentMaterial = mat;
651 currentMaterialName = mat->GetName();
652 G4double res = 0.0;
653 if(UpdateParticle(p, kinEnergy)) {
654 if(FindEmModel(p, processName, kinEnergy)) {
655 G4double e = kinEnergy;
656 if(baseParticle) {
657 e *= kinEnergy*massRatio;
658 res = currentModel->CrossSectionPerVolume(
659 mat, baseParticle, e, cut, e) * chargeSquare;
660 } else {
661 res = currentModel->CrossSectionPerVolume(mat, p, e, cut, e);
662 }
663 if(verbose>0) {
664 G4cout << "G4EmCalculator::ComputeXSPerVolume: E(MeV)= " << kinEnergy/MeV
665 << " cross(cm-1)= " << res*cm
666 << " cut(keV)= " << cut/keV
667 << " " << p->GetParticleName()
668 << " in " << mat->GetName()
669 << G4endl;
670 }
671 }
672 }
673 return res;
674}
675
676//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
677
678G4double G4EmCalculator::ComputeCrossSectionPerVolume(
679 G4double kinEnergy,
680 const G4String& particle,
681 const G4String& processName,
682 const G4String& material,
683 G4double cut)
684{
685 return ComputeCrossSectionPerVolume(kinEnergy,FindParticle(particle),
686 processName,
687 FindMaterial(material),cut);
688}
689
690//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
691
692G4double G4EmCalculator::ComputeCrossSectionPerAtom(
693 G4double kinEnergy,
694 const G4ParticleDefinition* p,
695 const G4String& processName,
696 G4double Z, G4double A,
697 G4double cut)
698{
699 G4double res = 0.0;
700 if(UpdateParticle(p, kinEnergy)) {
701 if(FindEmModel(p, processName, kinEnergy)) {
702 G4double e = kinEnergy;
703 if(baseParticle) {
704 e *= kinEnergy*massRatio;
705 res = currentModel->ComputeCrossSectionPerAtom(
706 baseParticle, e, Z, A, cut) * chargeSquare;
707 } else {
708 res = currentModel->ComputeCrossSectionPerAtom(p, e, Z, A, cut);
709 }
710 if(verbose>0) {
711 G4cout << "E(MeV)= " << kinEnergy/MeV
712 << " cross(barn)= " << res/barn
713 << " " << p->GetParticleName()
714 << " Z= " << Z << " A= " << A/(g/mole) << " g/mole"
715 << G4endl;
716 }
717 }
718 }
719 return res;
720}
721
722//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
723
724G4double G4EmCalculator::ComputeCrossSectionPerAtom(G4double kinEnergy,
725 const G4String& particle,
726 const G4String& processName,
727 const G4Element* elm,
728 G4double cut)
729{
730 return ComputeCrossSectionPerAtom(kinEnergy,FindParticle(particle),
731 processName,
732 elm->GetZ(),elm->GetN(),cut);
733}
734
735//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
736
737G4double
738G4EmCalculator::ComputeGammaAttenuationLength(G4double kinEnergy,
739 const G4Material* mat)
740{
741 G4double res = 0.0;
742 const G4ParticleDefinition* gamma = G4Gamma::Gamma();
743 res += ComputeCrossSectionPerVolume(kinEnergy, gamma, "conv", mat, 0.0);
744 res += ComputeCrossSectionPerVolume(kinEnergy, gamma, "compt", mat, 0.0);
745 res += ComputeCrossSectionPerVolume(kinEnergy, gamma, "phot", mat, 0.0);
746 res += ComputeCrossSectionPerVolume(kinEnergy, gamma, "Rayl", mat, 0.0);
747 if(res > 0.0) { res = 1.0/res; }
748 return res;
749}
750
751//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
752
753G4double G4EmCalculator::ComputeShellIonisationCrossSectionPerAtom(
754 const G4String& particle,
755 G4int Z,
756 G4AtomicShellEnumerator shell,
757 G4double kinEnergy,
758 const G4Material* mat)
759{
760 G4double res = 0.0;
761 const G4ParticleDefinition* p = FindParticle(particle);
762 G4VAtomDeexcitation* ad = manager->AtomDeexcitation();
763 if(p && ad) {
764 res = ad->ComputeShellIonisationCrossSectionPerAtom(p, Z, shell,
765 kinEnergy, mat);
766 }
767 return res;
768}
769
770//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
771
772G4double G4EmCalculator::ComputeMeanFreePath(G4double kinEnergy,
773 const G4ParticleDefinition* p,
774 const G4String& processName,
775 const G4Material* mat,
776 G4double cut)
777{
778 G4double mfp = DBL_MAX;
779 G4double x = ComputeCrossSectionPerVolume(kinEnergy, p, processName, mat, cut);
780 if(x > 0.0) { mfp = 1.0/x; }
781 if(verbose>1) {
782 G4cout << "E(MeV)= " << kinEnergy/MeV
783 << " MFP(mm)= " << mfp/mm
784 << " " << p->GetParticleName()
785 << " in " << mat->GetName()
786 << G4endl;
787 }
788 return mfp;
789}
790
791//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
792
793G4double G4EmCalculator::ComputeMeanFreePath(G4double kinEnergy,
794 const G4String& particle,
795 const G4String& processName,
796 const G4String& material,
797 G4double cut)
798{
799 return ComputeMeanFreePath(kinEnergy,FindParticle(particle),processName,
800 FindMaterial(material),cut);
801}
802
803//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
804
805G4double G4EmCalculator::ComputeEnergyCutFromRangeCut(
806 G4double range,
807 const G4ParticleDefinition* part,
808 const G4Material* mat)
809{
810 return G4ProductionCutsTable::GetProductionCutsTable()->
811 ConvertRangeToEnergy(part, mat, range);
812}
813
814//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
815
816G4double G4EmCalculator::ComputeEnergyCutFromRangeCut(
817 G4double range,
818 const G4String& particle,
819 const G4String& material)
820{
821 return ComputeEnergyCutFromRangeCut(range,FindParticle(particle),
822 FindMaterial(material));
823}
824
825
826//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
827
828G4bool G4EmCalculator::UpdateParticle(const G4ParticleDefinition* p,
829 G4double kinEnergy)
830{
831 if(p != currentParticle) {
832
833 // new particle
834 currentParticle = p;
835 dynParticle.SetDefinition(const_cast<G4ParticleDefinition*>(p));
836 dynParticle.SetKineticEnergy(kinEnergy);
837 baseParticle = 0;
838 currentParticleName = p->GetParticleName();
839 massRatio = 1.0;
840 mass = p->GetPDGMass();
841 chargeSquare = 1.0;
842 currentProcess = FindEnergyLossProcess(p);
843 currentProcessName = "";
844 isIon = false;
845
846 // ionisation process exist
847 if(currentProcess) {
848 currentProcessName = currentProcess->GetProcessName();
849 baseParticle = currentProcess->BaseParticle();
850
851 // base particle is used
852 if(baseParticle) {
853 massRatio = baseParticle->GetPDGMass()/p->GetPDGMass();
854 G4double q = p->GetPDGCharge()/baseParticle->GetPDGCharge();
855 chargeSquare = q*q;
856 }
857
858 if(p->GetParticleType() == "nucleus"
859 && currentParticleName != "deuteron"
860 && currentParticleName != "triton"
861 && currentParticleName != "alpha+"
862 && currentParticleName != "helium"
863 && currentParticleName != "hydrogen"
864 ) {
865 isIon = true;
866 massRatio = theGenericIon->GetPDGMass()/p->GetPDGMass();
867 baseParticle = theGenericIon;
868 // G4cout << p->GetParticleName()
869 // << " in " << currentMaterial->GetName()
870 // << " e= " << kinEnergy << G4endl;
871 }
872 }
873 }
874
875 // Effective charge for ions
876 if(isIon) {
877 chargeSquare =
878 corr->EffectiveChargeSquareRatio(p, currentMaterial, kinEnergy)
879 * corr->EffectiveChargeCorrection(p,currentMaterial,kinEnergy);
880 if(currentProcess) {
881 currentProcess->SetDynamicMassCharge(massRatio,chargeSquare);
882 //G4cout << "NewP: massR= " << massRatio << " q2= " << chargeSquare << G4endl;
883 }
884 }
885 return true;
886}
887
888//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
889
890const G4ParticleDefinition* G4EmCalculator::FindParticle(const G4String& name)
891{
892 const G4ParticleDefinition* p = 0;
893 if(name != currentParticleName) {
894 p = G4ParticleTable::GetParticleTable()->FindParticle(name);
895 if(!p) {
896 G4cout << "### WARNING: G4EmCalculator::FindParticle fails to find "
897 << name << G4endl;
898 }
899 } else {
900 p = currentParticle;
901 }
902 return p;
903}
904
905//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
906
907const G4ParticleDefinition* G4EmCalculator::FindIon(G4int Z, G4int A)
908{
909 const G4ParticleDefinition* p =
910 G4ParticleTable::GetParticleTable()->FindIon(Z,A,0,Z);
911 return p;
912}
913
914//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
915
916const G4Material* G4EmCalculator::FindMaterial(const G4String& name)
917{
918 if(name != currentMaterialName) {
919 currentMaterial = G4Material::GetMaterial(name);
920 currentMaterialName = name;
921 if(!currentMaterial)
922 G4cout << "### WARNING: G4EmCalculator::FindMaterial fails to find "
923 << name << G4endl;
924 }
925 return currentMaterial;
926}
927
928//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
929
930const G4Region* G4EmCalculator::FindRegion(const G4String& reg)
931{
932 const G4Region* r = 0;
933 if(reg != "" && reg != "world") {
934 r = G4RegionStore::GetInstance()->GetRegion(reg);
935 } else {
936 r = G4RegionStore::GetInstance()->GetRegion("DefaultRegionForTheWorld");
937 }
938 return r;
939}
940
941//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
942
943const G4MaterialCutsCouple* G4EmCalculator::FindCouple(
944 const G4Material* material,
945 const G4Region* region)
946{
947 if(!material) return 0;
948 currentMaterial = material;
949 currentMaterialName = material->GetName();
950 // Access to materials
951 const G4ProductionCutsTable* theCoupleTable=
952 G4ProductionCutsTable::GetProductionCutsTable();
953 const G4Region* r = region;
954 if(!r) {
955 r = G4RegionStore::GetInstance()->GetRegion("DefaultRegionForTheWorld");
956 }
957 return theCoupleTable->GetMaterialCutsCouple(material,r->GetProductionCuts());
958
959}
960
961//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
962
963G4bool G4EmCalculator::UpdateCouple(const G4Material* material, G4double cut)
964{
965 if(!material) return false;
966 currentMaterial = material;
967 currentMaterialName = material->GetName();
968 for (G4int i=0; i<nLocalMaterials; ++i) {
969 if(material == localMaterials[i] && cut == localCuts[i]) {
970 currentCouple = localCouples[i];
971 currentCoupleIndex = currentCouple->GetIndex();
972 currentCut = cut;
973 return true;
974 }
975 }
976 const G4MaterialCutsCouple* cc = new G4MaterialCutsCouple(material);
977 localMaterials.push_back(material);
978 localCouples.push_back(cc);
979 localCuts.push_back(cut);
980 nLocalMaterials++;
981 currentCouple = cc;
982 currentCoupleIndex = currentCouple->GetIndex();
983 currentCut = cut;
984 return true;
985}
986
987//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
988
989void G4EmCalculator::FindLambdaTable(const G4ParticleDefinition* p,
990 const G4String& processName,
991 G4double kinEnergy)
992{
993 // Search for the process
994 if (!currentLambda || p != lambdaParticle || processName != lambdaName) {
995 lambdaName = processName;
996 currentLambda = 0;
997 lambdaParticle = p;
998
999 const G4ParticleDefinition* part = p;
1000 if(isIon) { part = theGenericIon; }
1001
1002 // Search for energy loss process
1003 currentName = processName;
1004 currentModel = 0;
1005 loweModel = 0;
1006
1007 G4VEnergyLossProcess* elproc = FindEnLossProcess(part, processName);
1008 if(elproc) {
1009 currentLambda = elproc->LambdaTable();
1010 if(currentLambda) {
1011 isApplicable = true;
1012 if(verbose>1) {
1013 G4cout << "G4VEnergyLossProcess is found out: " << currentName
1014 << G4endl;
1015 }
1016 }
1017 return;
1018 }
1019
1020 // Search for discrete process
1021 G4VEmProcess* proc = FindDiscreteProcess(part, processName);
1022 if(proc) {
1023 currentLambda = proc->LambdaTable();
1024 if(currentLambda) {
1025 isApplicable = true;
1026 if(verbose>1) {
1027 G4cout << "G4VEmProcess is found out: " << currentName << G4endl;
1028 }
1029 }
1030 return;
1031 }
1032
1033 // Search for msc process
1034 G4VMultipleScattering* msc = FindMscProcess(part, processName);
1035 if(msc) {
1036 currentModel = msc->SelectModel(kinEnergy,0);
1037 /*
1038 if(currentModel) {
1039 currentLambda = currentModel->GetCrossSectionTable();
1040 if(currentLambda) {
1041 isApplicable = true;
1042 if(verbose>1) {
1043 G4cout << "G4VMultipleScattering is found out: " << currentName
1044 << G4endl;
1045 }
1046 }
1047 }
1048 */
1049 }
1050 }
1051}
1052
1053//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1054
1055G4bool G4EmCalculator::FindEmModel(const G4ParticleDefinition* p,
1056 const G4String& processName,
1057 G4double kinEnergy)
1058{
1059 isApplicable = false;
1060 if(!p) {
1061 G4cout << "G4EmCalculator::FindEmModel WARNING: no particle defined"
1062 << G4endl;
1063 return isApplicable;
1064 }
1065 G4String partname = p->GetParticleName();
1066 const G4ParticleDefinition* part = p;
1067 G4double scaledEnergy = kinEnergy*massRatio;
1068 if(isIon) { part = theGenericIon; }
1069
1070 if(verbose > 1) {
1071 G4cout << "## G4EmCalculator::FindEmModel for " << partname
1072 << " (type= " << p->GetParticleType()
1073 << ") and " << processName << " at E(MeV)= " << scaledEnergy
1074 << G4endl;
1075 if(p != part) { G4cout << " GenericIon is the base particle" << G4endl; }
1076 }
1077
1078 // Search for energy loss process
1079 currentName = processName;
1080 currentModel = 0;
1081 loweModel = 0;
1082 size_t idx = 0;
1083
1084 G4VEnergyLossProcess* elproc = FindEnLossProcess(part, processName);
1085 if(elproc) {
1086 currentModel = elproc->SelectModelForMaterial(scaledEnergy, idx);
1087 G4double eth = currentModel->LowEnergyLimit();
1088 if(eth > 0.0) {
1089 loweModel = elproc->SelectModelForMaterial(eth - CLHEP::eV, idx);
1090 if(loweModel == currentModel) { loweModel = 0; }
1091 }
1092 }
1093
1094 // Search for discrete process
1095 if(!currentModel) {
1096 G4VEmProcess* proc = FindDiscreteProcess(part, processName);
1097 if(proc) {
1098 currentModel = proc->SelectModelForMaterial(kinEnergy, idx);
1099 G4double eth = currentModel->LowEnergyLimit();
1100 if(eth > 0.0) {
1101 loweModel = proc->SelectModelForMaterial(eth - CLHEP::eV, idx);
1102 if(loweModel == currentModel) { loweModel = 0; }
1103 }
1104 }
1105 }
1106
1107 // Search for msc process
1108 if(!currentModel) {
1109 G4VMultipleScattering* proc = FindMscProcess(part, processName);
1110 if(proc) {
1111 currentModel = proc->SelectModel(kinEnergy, idx);
1112 loweModel = 0;
1113 }
1114 }
1115 if(currentModel) {
1116 if(loweModel == currentModel) { loweModel = 0; }
1117 isApplicable = true;
1118 if(verbose > 1) {
1119 G4cout << " Model <" << currentModel->GetName()
1120 << "> Emin(MeV)= " << currentModel->LowEnergyLimit()/MeV
1121 << " for " << part->GetParticleName();
1122 if(elproc) {
1123 G4cout << " and " << elproc->GetProcessName() << " " << elproc
1124 << G4endl;
1125 }
1126 if(loweModel) {
1127 G4cout << " LowEnergy model <" << loweModel->GetName() << ">";
1128 }
1129 G4cout << G4endl;
1130 }
1131 }
1132 return isApplicable;
1133}
1134
1135//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1136
1137G4VEnergyLossProcess* G4EmCalculator::FindEnergyLossProcess(
1138 const G4ParticleDefinition* p)
1139{
1140 G4VEnergyLossProcess* elp = 0;
1141 G4String partname = p->GetParticleName();
1142 const G4ParticleDefinition* part = p;
1143
1144 if(p->GetParticleType() == "nucleus"
1145 && currentParticleName != "deuteron"
1146 && currentParticleName != "triton"
1147 && currentParticleName != "alpha+"
1148 && currentParticleName != "helium"
1149 && currentParticleName != "hydrogen"
1150 ) { part = theGenericIon; }
1151
1152 elp = G4LossTableManager::Instance()->GetEnergyLossProcess(part);
1153 return elp;
1154}
1155
1156//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1157
1158G4VEnergyLossProcess*
1159G4EmCalculator::FindEnLossProcess(const G4ParticleDefinition* part,
1160 const G4String& processName)
1161{
1162 G4VEnergyLossProcess* proc = 0;
1163 const std::vector<G4VEnergyLossProcess*> v =
1164 G4LossTableManager::Instance()->GetEnergyLossProcessVector();
1165 G4int n = v.size();
1166 for(G4int i=0; i<n; ++i) {
1167 if((v[i])->GetProcessName() == processName) {
1168 G4VProcess* p = reinterpret_cast<G4VProcess*>(v[i]);
1169 if(ActiveForParticle(part, p)) {
1170 proc = v[i];
1171 break;
1172 }
1173 }
1174 }
1175 return proc;
1176}
1177
1178//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1179
1180G4VEmProcess*
1181G4EmCalculator::FindDiscreteProcess(const G4ParticleDefinition* part,
1182 const G4String& processName)
1183{
1184 G4VEmProcess* proc = 0;
1185 const std::vector<G4VEmProcess*> v =
1186 G4LossTableManager::Instance()->GetEmProcessVector();
1187 G4int n = v.size();
1188 for(G4int i=0; i<n; ++i) {
1189 if((v[i])->GetProcessName() == processName) {
1190 G4VProcess* p = reinterpret_cast<G4VProcess*>(v[i]);
1191 if(ActiveForParticle(part, p)) {
1192 proc = v[i];
1193 break;
1194 }
1195 }
1196 }
1197 return proc;
1198}
1199
1200//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1201
1202G4VMultipleScattering*
1203G4EmCalculator::FindMscProcess(const G4ParticleDefinition* part,
1204 const G4String& processName)
1205{
1206 G4VMultipleScattering* proc = 0;
1207 const std::vector<G4VMultipleScattering*> v =
1208 G4LossTableManager::Instance()->GetMultipleScatteringVector();
1209 G4int n = v.size();
1210 for(G4int i=0; i<n; ++i) {
1211 if((v[i])->GetProcessName() == processName) {
1212 G4VProcess* p = reinterpret_cast<G4VProcess*>(v[i]);
1213 if(ActiveForParticle(part, p)) {
1214 proc = v[i];
1215 break;
1216 }
1217 }
1218 }
1219 return proc;
1220}
1221
1222//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1223
1224G4bool G4EmCalculator::ActiveForParticle(const G4ParticleDefinition* part,
1225 G4VProcess* proc)
1226{
1227 G4ProcessManager* pm = part->GetProcessManager();
1228 G4ProcessVector* pv = pm->GetProcessList();
1229 G4int n = pv->size();
1230 G4bool res = false;
1231 for(G4int i=0; i<n; ++i) {
1232 if((*pv)[i] == proc) {
1233 if(pm->GetProcessActivation(i)) { res = true; }
1234 break;
1235 }
1236 }
1237 return res;
1238}
1239
1240//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1241
1242void G4EmCalculator::SetVerbose(G4int verb)
1243{
1244 verbose = verb;
1245}
1246
1247//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
1248
G4AtomicShellEnumerator
G4AtomicShellEnumerator
Definition: G4AtomicShellEnumerator.hh:52
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4bool
bool G4bool
Definition: G4Types.hh:67
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4DynamicParticle::SetDefinition
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
Definition: G4DynamicParticle.cc:271
G4DynamicParticle::SetKineticEnergy
void SetKineticEnergy(G4double aEnergy)
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4Element::GetN
G4double GetN() const
Definition: G4Element.hh:134
G4EmCalculator::FindParticle
const G4ParticleDefinition * FindParticle(const G4String &)
Definition: G4EmCalculator.cc:890
G4EmCalculator::ComputeMeanFreePath
G4double ComputeMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
Definition: G4EmCalculator.cc:772
G4EmCalculator::GetShellIonisationCrossSectionPerAtom
G4double GetShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy)
Definition: G4EmCalculator.cc:357
G4EmCalculator::GetRangeFromRestricteDEDX
G4double GetRangeFromRestricteDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:178
G4EmCalculator::ComputeTotalDEDX
G4double ComputeTotalDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:580
G4EmCalculator::ComputeNuclearDEDX
G4double ComputeNuclearDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *)
Definition: G4EmCalculator.cc:614
G4EmCalculator::FindIon
const G4ParticleDefinition * FindIon(G4int Z, G4int A)
Definition: G4EmCalculator.cc:907
G4EmCalculator::ComputeGammaAttenuationLength
G4double ComputeGammaAttenuationLength(G4double kinEnergy, const G4Material *)
Definition: G4EmCalculator.cc:738
G4EmCalculator::ComputeDEDX
G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:435
G4EmCalculator::SetVerbose
void SetVerbose(G4int val)
Definition: G4EmCalculator.cc:1242
G4EmCalculator::ComputeEnergyCutFromRangeCut
G4double ComputeEnergyCutFromRangeCut(G4double range, const G4ParticleDefinition *, const G4Material *)
Definition: G4EmCalculator.cc:805
G4EmCalculator::FindCouple
const G4MaterialCutsCouple * FindCouple(const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:943
G4EmCalculator::ComputeCrossSectionPerVolume
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
Definition: G4EmCalculator.cc:643
G4EmCalculator::ComputeCrossSectionPerAtom
G4double ComputeCrossSectionPerAtom(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4double Z, G4double A, G4double cut=0.0)
Definition: G4EmCalculator.cc:692
G4EmCalculator::GetCrossSectionPerVolume
G4double GetCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:308
G4EmCalculator::PrintDEDXTable
void PrintDEDXTable(const G4ParticleDefinition *)
Definition: G4EmCalculator.cc:407
G4EmCalculator::FindRegion
const G4Region * FindRegion(const G4String &)
Definition: G4EmCalculator.cc:930
G4EmCalculator::GetCSDARange
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:200
G4EmCalculator::GetKinEnergy
G4double GetKinEnergy(G4double range, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:277
G4EmCalculator::ComputeShellIonisationCrossSectionPerAtom
G4double ComputeShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy, const G4Material *mat=0)
Definition: G4EmCalculator.cc:753
G4EmCalculator::GetDEDX
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:130
G4EmCalculator::PrintRangeTable
void PrintRangeTable(const G4ParticleDefinition *)
Definition: G4EmCalculator.cc:416
G4EmCalculator::PrintInverseRangeTable
void PrintInverseRangeTable(const G4ParticleDefinition *)
Definition: G4EmCalculator.cc:425
G4EmCalculator::GetMeanFreePath
G4double GetMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:374
G4EmCalculator::GetRange
G4double GetRange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
Definition: G4EmCalculator.cc:222
G4EmCalculator::ComputeElectronicDEDX
G4double ComputeElectronicDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *mat, G4double cut=DBL_MAX)
Definition: G4EmCalculator.cc:538
G4EmCalculator::FindMaterial
const G4Material * FindMaterial(const G4String &)
Definition: G4EmCalculator.cc:916
G4EmCorrections::EffectiveChargeSquareRatio
G4double EffectiveChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4EmCorrections.hh:335
G4EmCorrections::NuclearDEDX
G4double NuclearDEDX(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy, G4bool fluct=true)
Definition: G4EmCorrections.cc:669
G4EmCorrections::EffectiveChargeCorrection
G4double EffectiveChargeCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4EmCorrections.cc:742
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
G4GenericIon::GenericIon
static G4GenericIon * GenericIon()
Definition: G4GenericIon.cc:92
G4LossTableManager::Instance
static G4LossTableManager * Instance()
Definition: G4LossTableManager.cc:107
G4LossTableManager::GetEmProcessVector
const std::vector< G4VEmProcess * > & GetEmProcessVector()
Definition: G4LossTableManager.cc:992
G4LossTableManager::GetCSDARange
G4double GetCSDARange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
Definition: G4LossTableManager.cc:1171
G4LossTableManager::GetEnergyLossProcess
G4VEnergyLossProcess * GetEnergyLossProcess(const G4ParticleDefinition *)
Definition: G4LossTableManager.cc:1128
G4LossTableManager::GetEnergy
G4double GetEnergy(const G4ParticleDefinition *aParticle, G4double range, const G4MaterialCutsCouple *couple)
Definition: G4LossTableManager.cc:1212
G4LossTableManager::GetMultipleScatteringVector
const std::vector< G4VMultipleScattering * > & GetMultipleScatteringVector()
Definition: G4LossTableManager.cc:1000
G4LossTableManager::GetEnergyLossProcessVector
const std::vector< G4VEnergyLossProcess * > & GetEnergyLossProcessVector()
Definition: G4LossTableManager.cc:985
G4LossTableManager::AtomDeexcitation
G4VAtomDeexcitation * AtomDeexcitation()
Definition: G4LossTableManager.cc:1106
G4LossTableManager::GetRangeFromRestricteDEDX
G4double GetRangeFromRestricteDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
Definition: G4LossTableManager.cc:1184
G4LossTableManager::EmCorrections
G4EmCorrections * EmCorrections()
Definition: G4LossTableManager.cc:1078
G4LossTableManager::GetRange
G4double GetRange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
Definition: G4LossTableManager.cc:1198
G4LossTableManager::GetDEDX
G4double GetDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
Definition: G4LossTableManager.cc:1145
G4Material::GetDensity
G4double GetDensity() const
Definition: G4Material.hh:179
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:177
G4Material::GetMaterial
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
Definition: G4Material.cc:576
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
G4ParticleDefinition::GetParticleType
const G4String & GetParticleType() const
Definition: G4ParticleDefinition.hh:136
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:119
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:115
G4ParticleTable::FindIon
G4ParticleDefinition * FindIon(G4int atomicNumber, G4int atomicMass, G4double excitationEnergy)
Definition: G4ParticleTable.cc:399
G4ParticleTable::FindParticle
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
Definition: G4ParticleTable.cc:472
G4ParticleTable::GetParticleTable
static G4ParticleTable * GetParticleTable()
Definition: G4ParticleTable.cc:63
G4ProcessManager::GetProcessActivation
G4bool GetProcessActivation(G4VProcess *aProcess) const
Definition: G4ProcessManager.cc:1124
G4ProcessManager::GetProcessList
G4ProcessVector * GetProcessList() const
G4ProcessVector::size
G4int size() const
G4ProductionCutsTable::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
Definition: G4ProductionCutsTable.hh:262
G4ProductionCutsTable::GetProductionCutsTable
static G4ProductionCutsTable * GetProductionCutsTable()
Definition: G4ProductionCutsTable.cc:63
G4RegionStore::GetInstance
static G4RegionStore * GetInstance()
Definition: G4RegionStore.cc:162
G4RegionStore::GetRegion
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
Definition: G4RegionStore.cc:216
G4Region::GetProductionCuts
G4ProductionCuts * GetProductionCuts() const
G4VAtomDeexcitation::ComputeShellIonisationCrossSectionPerAtom
virtual G4double ComputeShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
G4VAtomDeexcitation::GetShellIonisationCrossSectionPerAtom
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
G4VEmModel::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:240
G4VEmModel::CorrectionsAlongStep
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
Definition: G4VEmModel.cc:279
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:529
G4VEmModel::CrossSectionPerVolume
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:186
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:262
G4VEmModel::GetName
const G4String & GetName() const
Definition: G4VEmModel.hh:655
G4VEmModel::ComputeDEDXPerVolume
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.cc:177
G4VEmProcess::SelectModelForMaterial
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idxRegion) const
Definition: G4VEmProcess.hh:454
G4VEmProcess::LambdaTable
const G4PhysicsTable * LambdaTable() const
Definition: G4VEmProcess.hh:597
G4VEnergyLossProcess::BaseParticle
const G4ParticleDefinition * BaseParticle() const
Definition: G4VEnergyLossProcess.hh:887
G4VEnergyLossProcess::RangeTableForLoss
G4PhysicsTable * RangeTableForLoss() const
Definition: G4VEnergyLossProcess.hh:1108
G4VEnergyLossProcess::LambdaTable
G4PhysicsTable * LambdaTable()
Definition: G4VEnergyLossProcess.hh:1122
G4VEnergyLossProcess::InverseRangeTable
G4PhysicsTable * InverseRangeTable() const
Definition: G4VEnergyLossProcess.hh:1115
G4VEnergyLossProcess::SelectModelForMaterial
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idx) const
Definition: G4VEnergyLossProcess.hh:606
G4VEnergyLossProcess::SetDynamicMassCharge
void SetDynamicMassCharge(G4double massratio, G4double charge2ratio)
Definition: G4VEnergyLossProcess.hh:630
G4VEnergyLossProcess::DEDXTable
G4PhysicsTable * DEDXTable() const
Definition: G4VEnergyLossProcess.hh:1062
G4VMultipleScattering::SelectModel
G4VEmModel * SelectModel(G4double kinEnergy, size_t idx)
Definition: G4VMultipleScattering.hh:307
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:379
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83