Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4DNAPTBElasticModel.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// Authors: S. Meylan and C. Villagrasa (IRSN, France)
27// Models come from
28// M. Bug et al, Rad. Phys and Chem. 130, 459-479 (2017)
29//
30
31#include "G4DNAPTBElasticModel.hh"
32#include "G4DNAChampionElasticModel.hh"
33#include "G4PhysicalConstants.hh"
34#include "G4SystemOfUnits.hh"
35#include "G4DNAMolecularMaterial.hh"
36#include "G4Proton.hh"
37
38G4DNAPTBElasticModel::G4DNAPTBElasticModel(const G4String& applyToMaterial, const G4ParticleDefinition*,
39 const G4String& nam)
40 : G4VDNAModel(nam, applyToMaterial)
41{
42 fKillBelowEnergy = 10*eV; // will be override by the limits defined for each material
43
44 verboseLevel= 0;
45 // Verbosity scale:
46 // 0 = nothing
47 // 1 = warning for energy non-conservation
48 // 2 = details of energy budget
49 // 3 = calculation of cross sections, file openings, sampling of atoms
50 // 4 = entering in methods
51
52 if( verboseLevel>0 )
53 {
54 G4cout << "PTB Elastic model is constructed " << G4endl;
55 }
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
59
60G4DNAPTBElasticModel::~G4DNAPTBElasticModel()
61{
62
63}
64
65//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66
67void G4DNAPTBElasticModel::Initialise(const G4ParticleDefinition* particle,
68 const G4DataVector& /*cuts*/, G4ParticleChangeForGamma*)
69{
70 if (verboseLevel > 3)
71 G4cout << "Calling G4DNAPTBElasticModel::Initialise()" << G4endl;
72
73 G4double scaleFactor = 1e-16*cm*cm;
74
75 G4ParticleDefinition* electronDef = G4Electron::ElectronDefinition();
76
77 //*******************************************************
78 // Cross section data
79 //*******************************************************
80
81 if(particle == electronDef)
82 {
83 G4String particleName = particle->GetParticleName();
84
85 // MPietrzak, adding paths for N2
86 AddCrossSectionData("N2",
87 particleName,
88 "dna/sigma_elastic_e-_PTB_N2",
89 "dna/sigmadiff_cumulated_elastic_e-_PTB_N2",
90 scaleFactor);
91 SetLowELimit("N2", particleName, 10*eV);
92 SetHighELimit("N2", particleName, 1.02*MeV);
93 // MPietrzak
94
95 AddCrossSectionData("THF",
96 particleName,
97 "dna/sigma_elastic_e-_PTB_THF",
98 "dna/sigmadiff_cumulated_elastic_e-_PTB_THF",
99 scaleFactor);
100 SetLowELimit("THF", particleName, 10*eV);
101 SetHighELimit("THF", particleName, 1*keV);
102
103 AddCrossSectionData("PY",
104 particleName,
105 "dna/sigma_elastic_e-_PTB_PY",
106 "dna/sigmadiff_cumulated_elastic_e-_PTB_PY",
107 scaleFactor);
108 SetLowELimit("PY", particleName, 10*eV);
109 SetHighELimit("PY", particleName, 1*keV);
110
111 AddCrossSectionData("PU",
112 particleName,
113 "dna/sigma_elastic_e-_PTB_PU",
114 "dna/sigmadiff_cumulated_elastic_e-_PTB_PU",
115 scaleFactor);
116 SetLowELimit("PU", particleName, 10*eV);
117 SetHighELimit("PU", particleName, 1*keV);
118
119 AddCrossSectionData("TMP",
120 particleName,
121 "dna/sigma_elastic_e-_PTB_TMP",
122 "dna/sigmadiff_cumulated_elastic_e-_PTB_TMP",
123 scaleFactor);
124 SetLowELimit("TMP", particleName, 10*eV);
125 SetHighELimit("TMP", particleName, 1*keV);
126
127 AddCrossSectionData("G4_WATER",
128 particleName,
129 "dna/sigma_elastic_e_champion",
130 "dna/sigmadiff_cumulated_elastic_e_champion",
131 scaleFactor);
132 SetLowELimit("G4_WATER", particleName, 10*eV);
133 SetHighELimit("G4_WATER", particleName, 1*keV);
134
135 // DNA materials
136 //
137 AddCrossSectionData("backbone_THF",
138 particleName,
139 "dna/sigma_elastic_e-_PTB_THF",
140 "dna/sigmadiff_cumulated_elastic_e-_PTB_THF",
141 scaleFactor*33./30);
142 SetLowELimit("backbone_THF", particleName, 10*eV);
143 SetHighELimit("backbone_THF", particleName, 1*keV);
144
145 AddCrossSectionData("cytosine_PY",
146 particleName,
147 "dna/sigma_elastic_e-_PTB_PY",
148 "dna/sigmadiff_cumulated_elastic_e-_PTB_PY",
149 scaleFactor*42./30);
150 SetLowELimit("cytosine_PY", particleName, 10*eV);
151 SetHighELimit("cytosine_PY", particleName, 1*keV);
152
153 AddCrossSectionData("thymine_PY",
154 particleName,
155 "dna/sigma_elastic_e-_PTB_PY",
156 "dna/sigmadiff_cumulated_elastic_e-_PTB_PY",
157 scaleFactor*48./30);
158 SetLowELimit("thymine_PY", particleName, 10*eV);
159 SetHighELimit("thymine_PY", particleName, 1*keV);
160
161 AddCrossSectionData("adenine_PU",
162 particleName,
163 "dna/sigma_elastic_e-_PTB_PU",
164 "dna/sigmadiff_cumulated_elastic_e-_PTB_PU",
165 scaleFactor*50./44);
166 SetLowELimit("adenine_PU", particleName, 10*eV);
167 SetHighELimit("adenine_PU", particleName, 1*keV);
168
169 AddCrossSectionData("guanine_PU",
170 particleName,
171 "dna/sigma_elastic_e-_PTB_PU",
172 "dna/sigmadiff_cumulated_elastic_e-_PTB_PU",
173 scaleFactor*56./44);
174 SetLowELimit("guanine_PU", particleName, 10*eV);
175 SetHighELimit("guanine_PU", particleName, 1*keV);
176
177 AddCrossSectionData("backbone_TMP",
178 particleName,
179 "dna/sigma_elastic_e-_PTB_TMP",
180 "dna/sigmadiff_cumulated_elastic_e-_PTB_TMP",
181 scaleFactor*33./50);
182 SetLowELimit("backbone_TMP", particleName, 10*eV);
183 SetHighELimit("backbone_TMP", particleName, 1*keV);
184 }
185
186 //*******************************************************
187 // Load the data
188 //*******************************************************
189
190 LoadCrossSectionData(particle->GetParticleName() );
191
192 //*******************************************************
193 // Verbose output
194 //*******************************************************
195
196 if (verboseLevel > 2)
197 G4cout << "Loaded cross section files for PTB Elastic model" << G4endl;
198
199 if( verboseLevel>0 )
200 {
201 G4cout << "PTB Elastic model is initialized " << G4endl;
202 }
203}
204
205//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
206
207void G4DNAPTBElasticModel::ReadDiffCSFile(const G4String& materialName,
208 const G4String& particleName,
209 const G4String& file,
210 const G4double)
211{
212 // Method to read and save the information contained within the differential cross section files.
213 // This method is not yet standard.
214
215 // get the path of the G4LEDATA data folder
216 const char* path = G4FindDataDir("G4LEDATA");
217 // if it is not found then quit and print error message
218 if(!path)
219 {
220 G4Exception("G4DNAPTBElasticModel::ReadAllDiffCSFiles","em0006",
221 FatalException,"G4LEDATA environment variable not set.");
222 return;
223 }
224
225 // build the fullFileName path of the data file
226 std::ostringstream fullFileName;
227 fullFileName << path <<"/"<< file<<".dat";
228
229 // open the data file
230 std::ifstream diffCrossSection (fullFileName.str().c_str());
231 // error if file is not there
232 std::stringstream endPath;
233 if (!diffCrossSection)
234 {
235 endPath << "Missing data file: "<<file;
236 G4Exception("G4DNAPTBElasticModel::Initialise","em0003",
237 FatalException, endPath.str().c_str());
238 }
239
240 tValuesVec[materialName][particleName].push_back(0.);
241
242 G4String line;
243
244 // read the file line by line until we reach the end of file point
245 while(std::getline(diffCrossSection, line))
246 {
247 // check if the line is comment or empty
248 //
249 std::istringstream testIss(line);
250 G4String test;
251 testIss >> test;
252 // check first caracter to determine if following information is data or comments
253 if(test=="#")
254 {
255 // skip the line by beginning a new while loop.
256 continue;
257 }
258 // check if line is empty
259 else if(line.empty())
260 {
261 // skip the line by beginning a new while loop.
262 continue;
263 }
264 //
265 // end of the check
266
267 // transform the line into a iss
268 std::istringstream iss(line);
269
270 // Variables to be filled by the input file
271 double tDummy;
272 double eDummy;
273
274 // fill the variables with the content of the line
275 iss>>tDummy>>eDummy;
276
277 // SI : mandatory Vecm initialization
278
279 // Fill two vectors contained in maps of types:
280 // [materialName][particleName]=vector
281 // [materialName][particleName][T]=vector
282 // to list all the incident energies (tValues) and all the output energies (eValues) within the file
283 //
284 // Check if we already have the current T value in the vector.
285 // If not then add it
286 if (tDummy != tValuesVec[materialName][particleName].back())
287 {
288 // Add the current T value
289 tValuesVec[materialName][particleName].push_back(tDummy);
290
291 // Make it correspond to a default zero E value
292 eValuesVect[materialName][particleName][tDummy].push_back(0.);
293 }
294
295 // Put the differential cross section value of the input file within the diffCrossSectionData map
296 iss>>diffCrossSectionData[materialName][particleName][tDummy][eDummy];
297
298 // If the current E value (eDummy) is different from the one already registered in the eVector then add it to the vector
299 if (eDummy != eValuesVect[materialName][particleName][tDummy].back()) eValuesVect[materialName][particleName][tDummy].push_back(eDummy);
300 }
301}
302
303//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
304
305G4double G4DNAPTBElasticModel::CrossSectionPerVolume(const G4Material* /*material*/,
306 const G4String& materialName,
307 const G4ParticleDefinition* p,
308 G4double ekin,
309 G4double /*emin*/,
310 G4double /*emax*/)
311{
312 if (verboseLevel > 3)
313 G4cout << "Calling CrossSectionPerVolume() of G4DNAPTBElasticModel" << G4endl;
314
315 // Get the name of the current particle
316 const G4String& particleName = p->GetParticleName();
317
318 // set killBelowEnergy value for current material
319 fKillBelowEnergy = GetLowELimit(materialName, particleName);
320
321 // initialise the return value (cross section) to zero
322 G4double sigma(0);
323
324 // check if we are below the high energy limit
325 if (ekin < GetHighELimit(materialName, particleName) )
326 {
327 // This is used to kill the particle if its kinetic energy is below fKillBelowEnergy.
328 // If the energy is lower then we return a maximum cross section and thus the SampleSecondaries method will be called for sure.
329 // SampleSecondaries will remove the particle from the simulation.
330 //
331 //SI : XS must not be zero otherwise sampling of secondaries method ignored
332 if (ekin < fKillBelowEnergy) return DBL_MAX;
333
334 // Get the tables with the cross section data
335 TableMapData* tableData = GetTableData();
336
337 // Retrieve the cross section value
338 sigma = (*tableData)[materialName][particleName]->FindValue(ekin);
339 }
340
341 if (verboseLevel > 2)
342 {
343 G4cout << "__________________________________" << G4endl;
344 G4cout << "°°° G4DNAPTBElasticModel - XS INFO START" << G4endl;
345 G4cout << "°°° Kinetic energy(eV)=" << ekin/eV << " particle : " << particleName << G4endl;
346 G4cout << "°°° Cross section per molecule (cm^2)=" << sigma/cm/cm << G4endl;
347 G4cout << "°°° G4DNAPTBElasticModel - XS INFO END" << G4endl;
348 }
349
350 // Return the cross section
351 return sigma;
352}
353
354//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
355
356void G4DNAPTBElasticModel::SampleSecondaries(std::vector<G4DynamicParticle*>* /*fvect*/,
357 const G4MaterialCutsCouple* /*couple*/,
358 const G4String& materialName,
359 const G4DynamicParticle* aDynamicElectron,
360 G4ParticleChangeForGamma* particleChangeForGamma,
361 G4double /*tmin*/,
362 G4double /*tmax*/)
363{
364 if (verboseLevel > 3)
365 G4cout << "Calling SampleSecondaries() of G4DNAPTBElasticModel" << G4endl;
366
367 G4double electronEnergy0 = aDynamicElectron->GetKineticEnergy();
368
369 const G4String& particleName = aDynamicElectron->GetParticleDefinition()->GetParticleName();
370
371 // set killBelowEnergy value for material
372 fKillBelowEnergy = GetLowELimit(materialName, particleName);
373
374 // If the particle (electron here) energy is below the kill limit then we remove it from the simulation
375 if (electronEnergy0 < fKillBelowEnergy)
376 {
377 particleChangeForGamma->SetProposedKineticEnergy(0.);
378 particleChangeForGamma->ProposeTrackStatus(fStopAndKill);
379 particleChangeForGamma->ProposeLocalEnergyDeposit(electronEnergy0);
380 }
381 // If we are above the kill limite and below the high limit then we proceed
382 else if (electronEnergy0>= fKillBelowEnergy && electronEnergy0 < GetHighELimit(materialName, particleName) )
383 {
384 // Random sampling of the cosTheta
385 G4double cosTheta = RandomizeCosTheta(electronEnergy0, materialName);
386
387 // Random sampling of phi
388 G4double phi = 2. * pi * G4UniformRand();
389
390 G4ThreeVector zVers = aDynamicElectron->GetMomentumDirection();
391 G4ThreeVector xVers = zVers.orthogonal();
392 G4ThreeVector yVers = zVers.cross(xVers);
393
394 G4double xDir = std::sqrt(1. - cosTheta*cosTheta);
395 G4double yDir = xDir;
396 xDir *= std::cos(phi);
397 yDir *= std::sin(phi);
398
399 // Particle direction after ModelInterface
400 G4ThreeVector zPrikeVers((xDir*xVers + yDir*yVers + cosTheta*zVers));
401
402 // Give the new direction
403 particleChangeForGamma->ProposeMomentumDirection(zPrikeVers.unit()) ;
404
405 // Update the energy which does not change here
406 particleChangeForGamma->SetProposedKineticEnergy(electronEnergy0);
407 }
408}
409
410//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
411
412G4double G4DNAPTBElasticModel::Theta
413(G4ParticleDefinition * particleDefinition, G4double k, G4double integrDiff, const G4String& materialName)
414{
415 G4double theta = 0.;
416 G4double valueT1 = 0;
417 G4double valueT2 = 0;
418 G4double valueE21 = 0;
419 G4double valueE22 = 0;
420 G4double valueE12 = 0;
421 G4double valueE11 = 0;
422 G4double xs11 = 0;
423 G4double xs12 = 0;
424 G4double xs21 = 0;
425 G4double xs22 = 0;
426 G4String particleName = particleDefinition->GetParticleName();
427
428 if (particleDefinition == G4Electron::ElectronDefinition())
429 {
430 std::vector<double>::iterator t2 = std::upper_bound(tValuesVec[materialName][particleName].begin(),tValuesVec[materialName][particleName].end(), k);
431 std::vector<double>::iterator t1 = t2-1;
432
433 std::vector<double>::iterator e12 = std::upper_bound(eValuesVect[materialName][particleName][(*t1)].begin(),eValuesVect[materialName][particleName][(*t1)].end(), integrDiff);
434 std::vector<double>::iterator e11 = e12-1;
435
436 std::vector<double>::iterator e22 = std::upper_bound(eValuesVect[materialName][particleName][(*t2)].begin(),eValuesVect[materialName][particleName][(*t2)].end(), integrDiff);
437 std::vector<double>::iterator e21 = e22-1;
438
439 valueT1 =*t1;
440 valueT2 =*t2;
441 valueE21 =*e21;
442 valueE22 =*e22;
443 valueE12 =*e12;
444 valueE11 =*e11;
445
446 xs11 = diffCrossSectionData[materialName][particleName][valueT1][valueE11];
447 xs12 = diffCrossSectionData[materialName][particleName][valueT1][valueE12];
448 xs21 = diffCrossSectionData[materialName][particleName][valueT2][valueE21];
449 xs22 = diffCrossSectionData[materialName][particleName][valueT2][valueE22];
450 }
451
452 if (xs11==0 && xs12==0 && xs21==0 && xs22==0) return (0.);
453
454 theta = QuadInterpolator ( valueE11, valueE12,
455 valueE21, valueE22,
456 xs11, xs12,
457 xs21, xs22,
458 valueT1, valueT2,
459 k, integrDiff );
460
461 return theta;
462}
463
464//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
465
466G4double G4DNAPTBElasticModel::LinLogInterpolate(G4double e1,
467 G4double e2,
468 G4double e,
469 G4double xs1,
470 G4double xs2)
471{
472 G4double d1 = std::log(xs1);
473 G4double d2 = std::log(xs2);
474 G4double value = std::exp(d1 + (d2 - d1)*(e - e1)/ (e2 - e1));
475 return value;
476}
477
478//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
479
480G4double G4DNAPTBElasticModel::LinLinInterpolate(G4double e1,
481 G4double e2,
482 G4double e,
483 G4double xs1,
484 G4double xs2)
485{
486 G4double d1 = xs1;
487 G4double d2 = xs2;
488 G4double value = (d1 + (d2 - d1)*(e - e1)/ (e2 - e1));
489 return value;
490}
491
492//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
493
494G4double G4DNAPTBElasticModel::LogLogInterpolate(G4double e1,
495 G4double e2,
496 G4double e,
497 G4double xs1,
498 G4double xs2)
499{
500 G4double a = (std::log10(xs2)-std::log10(xs1)) / (std::log10(e2)-std::log10(e1));
501 G4double b = std::log10(xs2) - a*std::log10(e2);
502 G4double sigma = a*std::log10(e) + b;
503 G4double value = (std::pow(10.,sigma));
504 return value;
505}
506
507//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
508
509G4double G4DNAPTBElasticModel::QuadInterpolator(G4double e11, G4double e12,
510 G4double e21, G4double e22,
511 G4double xs11, G4double xs12,
512 G4double xs21, G4double xs22,
513 G4double t1, G4double t2,
514 G4double t, G4double e)
515{
516 // Log-Log
517 /*
518 G4double interpolatedvalue1 = LogLogInterpolate(e11, e12, e, xs11, xs12);
519 G4double interpolatedvalue2 = LogLogInterpolate(e21, e22, e, xs21, xs22);
520 G4double value = LogLogInterpolate(t1, t2, t, interpolatedvalue1, interpolatedvalue2);
521
522
523 // Lin-Log
524 G4double interpolatedvalue1 = LinLogInterpolate(e11, e12, e, xs11, xs12);
525 G4double interpolatedvalue2 = LinLogInterpolate(e21, e22, e, xs21, xs22);
526 G4double value = LinLogInterpolate(t1, t2, t, interpolatedvalue1, interpolatedvalue2);
527*/
528
529 // Lin-Lin
530 G4double interpolatedvalue1 = LinLinInterpolate(e11, e12, e, xs11, xs12);
531 G4double interpolatedvalue2 = LinLinInterpolate(e21, e22, e, xs21, xs22);
532 G4double value = LinLinInterpolate(t1, t2, t, interpolatedvalue1, interpolatedvalue2);
533
534 return value;
535}
536
537//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
538
539G4double G4DNAPTBElasticModel::RandomizeCosTheta(G4double k, const G4String& materialName)
540{
541 G4double integrdiff=0;
542 G4double uniformRand=G4UniformRand();
543 integrdiff = uniformRand;
544
545 G4double theta=0.;
546 G4double cosTheta=0.;
547 theta = Theta(G4Electron::ElectronDefinition(),k/eV,integrdiff, materialName);
548
549 cosTheta= std::cos(theta*pi/180);
550
551 return cosTheta;
552}
553
554
555
556
G4FindDataDir
const char * G4FindDataDir(const char *)
Definition: G4FindDataDir.cc:93
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:59
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
G4double
double G4double
Definition: G4Types.hh:83
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::orthogonal
Hep3Vector orthogonal() const
CLHEP::Hep3Vector::cross
Hep3Vector cross(const Hep3Vector &) const
G4DNAPTBElasticModel::CrossSectionPerVolume
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4String &materialName, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
CrossSectionPerVolume This method is mandatory for any model class. It finds and return the cross sec...
Definition: G4DNAPTBElasticModel.cc:305
G4DNAPTBElasticModel::G4DNAPTBElasticModel
G4DNAPTBElasticModel(const G4String &applyToMaterial="all", const G4ParticleDefinition *p=0, const G4String &nam="DNAPTBElasticModel")
G4DNAPTBElasticModel Constructor.
Definition: G4DNAPTBElasticModel.cc:38
G4DNAPTBElasticModel::SampleSecondaries
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4String &materialName, const G4DynamicParticle *, G4ParticleChangeForGamma *particleChangeForGamma, G4double tmin, G4double tmax)
SampleSecondaries Method called after CrossSectionPerVolume if the process is the one which is select...
Definition: G4DNAPTBElasticModel.cc:356
G4DNAPTBElasticModel::~G4DNAPTBElasticModel
virtual ~G4DNAPTBElasticModel()
~G4DNAPTBElasticModel Destructor
Definition: G4DNAPTBElasticModel.cc:60
G4DNAPTBElasticModel::Initialise
virtual void Initialise(const G4ParticleDefinition *particle, const G4DataVector &, G4ParticleChangeForGamma *fpChangeForGamme=nullptr)
Initialise Mandatory method for every model class. The material/particle for which the model can be u...
Definition: G4DNAPTBElasticModel.cc:67
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Electron::ElectronDefinition
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:88
G4ParticleChangeForGamma::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForGamma.hh:106
G4ParticleChangeForGamma::ProposeMomentumDirection
void ProposeMomentumDirection(const G4ThreeVector &Pfinal)
Definition: G4ParticleChangeForGamma.hh:120
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:109
G4VDNAModel
The G4VDNAModel class.
Definition: G4VDNAModel.hh:50
G4VDNAModel::GetTableData
TableMapData * GetTableData()
GetTableData.
Definition: G4VDNAModel.hh:193
G4VDNAModel::SetHighELimit
void SetHighELimit(const G4String &material, const G4String &particle, G4double lim)
SetHighEnergyLimit.
Definition: G4VDNAModel.hh:169
G4VDNAModel::TableMapData
std::map< G4String, std::map< G4String, G4DNACrossSectionDataSet *, std::less< G4String > > > TableMapData
Definition: G4VDNAModel.hh:183
G4VDNAModel::AddCrossSectionData
void AddCrossSectionData(G4String materialName, G4String particleName, G4String fileCS, G4String fileDiffCS, G4double scaleFactor)
AddCrossSectionData Method used during the initialization of the model class to add a new material....
Definition: G4VDNAModel.cc:58
G4VDNAModel::SetLowELimit
void SetLowELimit(const G4String &material, const G4String &particle, G4double lim)
SetLowEnergyLimit.
Definition: G4VDNAModel.hh:177
G4VDNAModel::GetHighELimit
G4double GetHighELimit(const G4String &material, const G4String &particle)
GetHighEnergyLimit.
Definition: G4VDNAModel.hh:153
G4VDNAModel::LoadCrossSectionData
void LoadCrossSectionData(const G4String &particleName)
LoadCrossSectionData Method to loop on all the registered materials in the model and load the corresp...
Definition: G4VDNAModel.cc:75
G4VDNAModel::GetLowELimit
G4double GetLowELimit(const G4String &material, const G4String &particle)
GetLowEnergyLimit.
Definition: G4VDNAModel.hh:161
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
DBL_MAX
#define DBL_MAX
Definition: templates.hh:62