Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4hQAOModel.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: G4hQAOModel
33//
34// Author: V.Ivanchenko ([email protected])
35//
36// Creation date: 27 April 2004
37//
38// Modifications:
39//
40// Class Description:
41//
42// Electronic stopping power for negative heavy partiles
43//
44// Class Description: End
45//
46// -------------------------------------------------------------------
47//
48//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
49
50#include "G4hQAOModel.hh"
51#include "G4PhysicalConstants.hh"
52#include "G4SystemOfUnits.hh"
53#include "G4UnitsTable.hh"
54#include "G4Element.hh"
55#include "G4Material.hh"
56
57//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
58
59G4hQAOModel::G4hQAOModel():G4VhElectronicStoppingPower(),
60 numberOfMaterials(6),
61 sizeL0(67),
62 sizeL1(22),
63 sizeL2(14),
64 currentMaterial(0),
65 currentElement(0)
66{
67 theZieglerFactor = eV*cm2*1.0e-15;
68 thePlasmonFactor = 28.816 * 28.816 * 1e-6;
69}
70
71//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
72
73G4hQAOModel::~G4hQAOModel()
74{}
75
76//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
77
78G4double G4hQAOModel::StoppingPower(const G4Material* material,
79 G4double kineticEnergy)
80{
81 G4double eloss = 0.0 ;
82 G4int numberOfElements = material->GetNumberOfElements() ;
83 const G4ElementVector* theElementVector = material->GetElementVector();
84 const G4double* theAtomicNumDensityVector = material->GetAtomicNumDensityVector() ;
85 currentMaterial = material;
86
87 // loop for the elements in the material
88 for (G4int i=0; i<numberOfElements; i++)
89 {
90 currentElement = (*theElementVector)[i] ;
91 G4double z = currentElement->GetZ() ;
92 eloss += ElectronicStoppingPower(z,kineticEnergy)*theAtomicNumDensityVector[i];
93 }
94
95 return eloss;
96}
97
98//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
99
100G4double G4hQAOModel::ElectronicStoppingPower(G4double z,
101 G4double kineticEnergy) const
102{
103 G4int Z = (G4int)z;
104 G4int nbOfShell = GetNumberOfShell(Z);
105 if(nbOfShell < 1) nbOfShell = 1;
106
107 G4double dedx=0.0;
108
109 G4double v = c_light * std::sqrt( 2.0 * kineticEnergy / proton_mass_c2 );
110 G4double coeff = twopi*proton_mass_c2*z / (electron_mass_c2*theZieglerFactor) ;
111 G4double fBetheVelocity = fine_structure_const * c_light / v;
112 coeff *= fine_structure_const * fine_structure_const * hbarc_squared /
113 kineticEnergy ;
114
115 G4double l0Term = 0, l1Term = 0, l2Term = 0;
116
117 for (G4int nos = 0 ; nos < nbOfShell ; nos++){
118
119 G4double l0 = 0, l1 = 0, l2 = 0;
120 G4double NormalizedEnergy = ( 2.0 * electron_mass_c2 * v * v ) /
121 ( c_squared * GetShellEnergy(Z, nos) );
122 G4double shStrength = GetShellStrength(Z, nos);
123
124 l0 = GetL0(NormalizedEnergy);
125 l0Term += shStrength * l0;
126
127 l1 = GetL1(NormalizedEnergy);
128 l1Term += shStrength * l1;
129
130 l2 = GetL2(NormalizedEnergy);
131 l2Term += shStrength * l2;
132
133/*
134 //if(Z == 6){
135 G4cout << nos << ". "
136 << " E(MeV)= " << kineticEnergy/MeV
137 << " normE= " << NormalizedEnergy
138 << " sh en= " << GetShellEnergy(Z,nos)
139 << " str= " << shStrength
140 << " v0/v= " << fBetheVelocity
141 << " l0= " << l0Term
142 << " l1= " << l1Term
143 << " l2= " << l2Term
144 << G4endl;
145 // }
146*/
147 }
148
149 dedx = coeff * (l0Term - fBetheVelocity*l1Term + fBetheVelocity*fBetheVelocity*l2Term);
150
151 //G4cout << " E(MeV)= " << kineticEnergy/MeV
152 // << " dedx(Mev/mm)= " << dedx*mm/MeV << G4endl;
153
154 if(dedx < 0.0) dedx = 0.0;
155 return dedx;
156}
157
158
159G4int G4hQAOModel::GetNumberOfShell(G4int Z) const
160{
161 // Set return value from table
162 G4int nShell = 0;
163
164 // Set return value if in material available from Aahrus
165 for(G4int i=0; i<numberOfMaterials; i++) {
166
167 if(materialAvailable[i] == Z){
168 nShell = nbofShellForMaterial[i];
169 break;
170 }
171 else nShell = fNumberOfShells[Z];
172 }
173
174 return nShell;
175}
176
177
178
179G4double G4hQAOModel::GetShellEnergy(G4int Z, G4int nbOfTheShell) const
180{
181 G4double shellEnergy = alShellEnergy[0];
182 if (Z == 13) shellEnergy = alShellEnergy[nbOfTheShell];
183 else if(Z == 14) shellEnergy = siShellEnergy[nbOfTheShell];
184 else if(Z == 29) shellEnergy = cuShellEnergy[nbOfTheShell];
185 else if(Z == 73) shellEnergy = taShellEnergy[nbOfTheShell];
186 else if(Z == 79) shellEnergy = auShellEnergy[nbOfTheShell];
187 else if(Z == 78) shellEnergy = ptShellEnergy[nbOfTheShell];
188 else shellEnergy = GetOscillatorEnergy(Z, nbOfTheShell);
189
190 return shellEnergy;
191}
192
193
194G4double G4hQAOModel::GetOscillatorEnergy(G4int Z, G4int nbOfTheShell) const
195{
196 G4double squaredPlasmonEnergy = thePlasmonFactor;
197 G4double zeff = (G4double)Z;
198 if(currentMaterial)
199 squaredPlasmonEnergy *= (currentMaterial->GetDensity()*zeff*cm3)
200 / (g * currentElement->GetN());
201
202 G4double occn = GetOccupationNumber(Z,nbOfTheShell);
203 G4double plasmonTerm = 0.66667 * occn * squaredPlasmonEnergy/(zeff*zeff);
204
205 G4double ionTerm = std::exp(0.5) * currentElement->GetAtomicShell(nbOfTheShell);
206
207 ionTerm = ionTerm*ionTerm ;
208
209 G4double oscShellEnergy = std::sqrt( ionTerm + plasmonTerm );
210
211/* if(material->GetName()=="Graphite"){
212 G4cout << "\t" << Z
213 << "\t" << nbOfTheShell
214 << "\t" << squaredPlasmonEnergy
215 << "\t" << plasmonTerm
216 << "\t" << ionTerm
217 << "\t" << GetOccupationNumber(Z,nbOfTheShell)
218 << "\t" << element->GetAtomicShell(nbOfTheShell)
219 << "\t" << oscShellEnergy << G4endl;}
220*/
221 return oscShellEnergy;
222}
223
224
225G4double G4hQAOModel::GetShellStrength(G4int Z, G4int nbOfTheShell) const
226{
227 G4double shellStrength = alShellStrength[0];
228
229 if(Z == 13) shellStrength = alShellStrength[nbOfTheShell];
230 else if(Z == 14) shellStrength = siShellStrength[nbOfTheShell];
231 else if(Z == 29) shellStrength = cuShellStrength[nbOfTheShell];
232 else if(Z == 73) shellStrength = taShellStrength[nbOfTheShell];
233 else if(Z == 79) shellStrength = auShellStrength[nbOfTheShell];
234 else if(Z == 78) shellStrength = ptShellStrength[nbOfTheShell];
235 else shellStrength = GetOccupationNumber(Z,nbOfTheShell) / (G4double)Z;
236
237 return shellStrength;
238}
239
240G4double G4hQAOModel::GetOccupationNumber(G4int Z, G4int ShellNb) const
241{
242
243 G4int indice = ShellNb ;
244 for (G4int z = 1 ; z < Z ; z++) {indice += fNumberOfShells[z];}
245 G4double x = (G4double)(nbOfElectronPerSubShell[indice+1]);
246 return x;
247}
248
249
250G4double G4hQAOModel::GetL0(G4double normEnergy) const
251{
252 G4int n;
253
254 for(n = 0; n < sizeL0; n++) {
255 if( normEnergy < L0[n][0] ) break;
256 }
257 if(0 == n) n = 1 ;
258 if(n >= sizeL0) n = sizeL0 - 1 ;
259
260 G4double l0 = L0[n][1];
261 G4double l0p = L0[n-1][1];
262 G4double bethe = l0p + (l0 - l0p) * ( normEnergy - L0[n-1][0]) /
263 (L0[n][0] - L0[n-1][0]);
264 return bethe ;
265
266}
267
268G4double G4hQAOModel::GetL1(G4double normEnergy) const
269{
270 G4int n;
271
272 for(n = 0; n < sizeL1; n++) {
273 if( normEnergy < L1[n][0] ) break;
274 }
275 if(0 == n) n = 1 ;
276 if(n >= sizeL1) n = sizeL1 - 1 ;
277
278 G4double l1 = L1[n][1];
279 G4double l1p = L1[n-1][1];
280 G4double barkas= l1p + (l1 - l1p) * ( normEnergy - L1[n-1][0]) /
281 (L1[n][0] - L1[n-1][0]);
282
283 return barkas;
284
285}
286
287
288G4double G4hQAOModel::GetL2(G4double normEnergy) const
289{
290 G4int n;
291 for(n = 0; n < sizeL2; n++) {
292 if( normEnergy < L2[n][0] ) break;
293 }
294 if(0 == n) n = 1 ;
295 if(n >= sizeL2) n = sizeL2 - 1 ;
296
297 G4double l2 = L2[n][1];
298 G4double l2p = L2[n-1][1];
299 G4double bloch = l2p + (l2 - l2p) * ( normEnergy - L2[n-1][0]) /
300 (L2[n][0] - L2[n-1][0]);
301
302 return bloch;
303}
304
305
306const G4int G4hQAOModel::materialAvailable[6] = {13,14,29,73,79,78};
307 /*
308 "Aluminum",
309 "Silicon",
310 "Copper",
311 "Tantalum",
312 "Gold",
313 "Platinum"};
314 */
315const G4int G4hQAOModel::nbofShellForMaterial[6] = {3,3,4,6,6,6};
316
317const G4double G4hQAOModel::alShellEnergy[3] ={ 2795e-6, 202e-6, 16.9e-6};
318const G4double G4hQAOModel::alShellStrength[3]={ 0.1349, 0.6387, 0.2264};
319const G4double G4hQAOModel::siShellEnergy[3] ={ 3179e-6, 249e-6, 20.3e-6 };
320const G4double G4hQAOModel::siShellStrength[3]={ 0.1222, 0.5972, 0.2806};
321const G4double G4hQAOModel::cuShellEnergy[4] ={ 16931e-6, 1930e-6, 199e-6, 39.6e-6};
322const G4double G4hQAOModel::cuShellStrength[4]={ 0.0505, 0.2561, 0.4913, 0.2021};
323const G4double G4hQAOModel::taShellEnergy[6] ={ 88926e-6, 18012e-6, 3210e-6, 575e-6, 108.7e-6, 30.8e-6};
324const G4double G4hQAOModel::taShellStrength[6]={ 0.0126, 0.0896, 0.2599, 0.3413, 0.2057, 0.0908};
325const G4double G4hQAOModel::auShellEnergy[6]={ 96235e-6, 25918e-6, 4116e-6, 599e-6, 87.3e-6, 36.9e-6};
326const G4double G4hQAOModel::auShellStrength[6]={ 0.0139, 0.0803, 0.2473, 0.423, 0.1124, 0.1231};
327const G4double G4hQAOModel::ptShellEnergy[6]={ 95017e-6, 25590e-6, 4063e-6, 576e-6, 81.9e-6, 31.4e-6};
328const G4double G4hQAOModel::ptShellStrength[6]={ 0.0129, 0.0745, 0.2295, 0.4627, 0.1324, 0.0879};
329
330
331const G4double G4hQAOModel::L0[67][2] =
332{
333 {0.00, 0.000001},
334 {0.10, 0.000001},
335 {0.12, 0.00001},
336 {0.14, 0.00005},
337 {0.16, 0.00014},
338 {0.18, 0.00030},
339 {0.20, 0.00057},
340 {0.25, 0.00189},
341 {0.30, 0.00429},
342 {0.35, 0.00784},
343 {0.40, 0.01248},
344 {0.45, 0.01811},
345 {0.50, 0.02462},
346 {0.60, 0.03980},
347 {0.70, 0.05731},
348 {0.80, 0.07662},
349 {0.90, 0.09733},
350 {1.00, 0.11916},
351 {1.20, 0.16532},
352 {1.40, 0.21376},
353 {1.60, 0.26362},
354 {1.80, 0.31428},
355 {2.00, 0.36532},
356 {2.50, 0.49272},
357 {3.00, 0.61765},
358 {3.50, 0.73863},
359 {4.00, 0.85496},
360 {4.50, 0.96634},
361 {5.00, 1.07272},
362 {6.00, 1.27086},
363 {7.00, 1.45075},
364 {8.00, 1.61412},
365 {9.00, 1.76277},
366 {10.00, 1.89836},
367 {12.00, 2.13625},
368 {14.00, 2.33787},
369 {16.00, 2.51093},
370 {18.00, 2.66134},
371 {20.00, 2.79358},
372 {25.00, 3.06539},
373 {30.00, 3.27902},
374 {35.00, 3.45430},
375 {40.00, 3.60281},
376 {45.00, 3.73167},
377 {50.00, 3.84555},
378 {60.00, 4.04011},
379 {70.00, 4.20264},
380 {80.00, 4.34229},
381 {90.00, 4.46474},
382 {100.00, 4.57378},
383 {120.00, 4.76155},
384 {140.00, 4.91953},
385 {160.00, 5.05590},
386 {180.00, 5.17588},
387 {200.00, 5.28299},
388 {250.00, 5.50925},
389 {300.00, 5.69364},
390 {350.00, 5.84926},
391 {400.00, 5.98388},
392 {450.00, 6.10252},
393 {500.00, 6.20856},
394 {600.00, 6.39189},
395 {700.00, 6.54677},
396 {800.00, 6.68084},
397 {900.00, 6.79905},
398 {1000.00, 6.90474}
399};
400
401
402const G4double G4hQAOModel::L1[22][2] =
403{
404 {0.00, -0.000001},
405 {0.10, -0.00001},
406 {0.20, -0.00049},
407 {0.30, -0.00084},
408 {0.40, 0.00085},
409 {0.50, 0.00519},
410 {0.60, 0.01198},
411 {0.70, 0.02074},
412 {0.80, 0.03133},
413 {0.90, 0.04369},
414 {1.00, 0.06035},
415 {2.00, 0.24023},
416 {3.00, 0.44284},
417 {4.00, 0.62012},
418 {5.00, 0.77031},
419 {6.00, 0.90390},
420 {7.00, 1.02705},
421 {8.00, 1.10867},
422 {9.00, 1.17546},
423 {10.00, 1.21599},
424 {15.00, 1.24349},
425 {20.00, 1.16752}
426};
427
428
429const G4double G4hQAOModel::L2[14][2] =
430{
431 {0.00, 0.000001},
432 {0.10, 0.00001},
433 {0.20, 0.00000},
434 {0.40, -0.00120},
435 {0.60, -0.00036},
436 {0.80, 0.00372},
437 {1.00, 0.01298},
438 {2.00, 0.08296},
439 {4.00, 0.21953},
440 {6.00, 0.23903},
441 {8.00, 0.20893},
442 {10.00, 0.10879},
443 {20.00, -0.88409},
444 {40.00, -1.13902}
445};
446
447
448const G4int G4hQAOModel::nbOfElectronPerSubShell[1540] =
449{
450 0, // consistency with G4AtomicShells
451 1,//------ H
452 2,//------ He
453 2, 1,//------ Li
454 2, 2,//------ Be
455 2, 2, 1,//------ B
456 2, 2, 2,//------ C
457 2, 2, 2, 1,//------ N
458 2, 2, 2, 2,//------ O
459 2, 2, 5,//------ F
460 2, 2, 2, 4,//------ Ne
461 2, 2, 2, 4, 1,//------ Na
462 2, 2, 2, 4, 2,//------ Mg
463 2, 2, 2, 4, 2, 1,//------ Al
464 2, 2, 2, 4, 2, 2,//------ Si
465 2, 2, 2, 4, 2, 3,//------ P
466 2, 2, 2, 4, 2, 4,//------
467 2, 2, 2, 4, 2, 5,//------
468 2, 2, 2, 4, 2, 2, 4,//------
469 2, 2, 2, 4, 2, 2, 4, 1,//------
470 2, 2, 2, 4, 2, 2, 4, 2,//------
471 2, 2, 2, 4, 2, 2, 4, 1, 2,//------
472 2, 2, 2, 4, 2, 2, 4, 2, 2,//------
473 2, 2, 2, 4, 2, 2, 4, 3, 2,//------
474 2, 2, 2, 4, 2, 2, 4, 4, 2,//------
475 2, 2, 2, 4, 2, 2, 4, 5, 2,//------
476 2, 2, 2, 4, 2, 2, 4, 6, 2,//------
477 2, 2, 2, 4, 2, 2, 4, 7, 2,//------
478 2, 2, 2, 4, 2, 2, 4, 4, 4, 2,//------
479 2, 2, 2, 4, 2, 2, 4, 4, 5, 2,//------
480 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,//------
481 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 1,//------
482 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2,//------
483 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 3,//------
484 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 4,//------
485 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 5,//------
486 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4,//------
487 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1,//------
488 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2,//------
489 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,//------
490 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,//------
491 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 3, 2,//------
492 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 2,//------
493 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 5, 2,//------
494 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,//------
495 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,//------
496 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 4, 2,//------
497 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 5, 2,//------
498 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,//------
499 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 1,//------
500 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2,//------
501 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 3,//------
502 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 4,//------
503 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 5,//------
504 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4,//------
505 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1,//------
506 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2,//------
507 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,//------
508 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 2, 4, 2,//------
509 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 3, 2, 2, 4, 2,//------
510 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 4, 2, 2, 4, 2,//------
511 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 5, 2, 2, 4, 2,//------
512 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 2, 2, 4, 2,//------
513 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 7, 2, 2, 4, 2,//------
514 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 7, 2, 2, 4, 1, 2,//------
515 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 9, 2, 2, 4, 2,//------
516 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 10, 2, 2, 4, 2,//------
517 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 11, 2, 2, 4, 2,//------
518 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 12, 2, 2, 4, 2,//------
519 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 13, 2, 2, 4, 2,//------
520 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 2,//------
521 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 1, 2,//------
522 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 2, 2,//------
523 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 3, 2,//------
524 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 2,//------
525 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 5, 2,//------
526 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 6, 2,//------
527 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 7, 2,//------
528 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 9, 1,//------
529 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 1,//------
530 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2,//------
531 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 1,//------
532 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2,//------
533 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 3,//------
534 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 4,//------
535 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 3,//------
536 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,//------
537 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 1,//------
538 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 2,//------
539 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,//------
540 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,//------
541 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 2, 4, 1, 2,//------
542 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 3, 2, 2, 4, 1, 2,//------
543 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 4, 2, 2, 4, 1, 2,//------
544 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 6, 2, 2, 4, 2,//------
545 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 7, 2, 2, 4, 2,//------
546 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 7, 2, 2, 4, 1, 2,//------
547 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 8, 2, 2, 4, 1, 2,//------
548 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 10, 2, 2, 4, 2,//------
549 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 11, 2, 2, 4, 2,//------
550 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 12, 2, 2, 4, 2 //-----
551};
552
553const G4int G4hQAOModel::fNumberOfShells[101] =
554{
555 0 , // nonexisting zero element
556
557 1 , 1 , 2 , 2 , 3 , 3 , 4 , 4 , 3 , 4 , // 1 - 10
558
559 5 , 5 , 6 , 6 , 6 , 6 , 6 , 7 , 8 , 8 , // 11 - 20
560
561 9 , 9 , 9 , 9 , 9 , 9 , 9 , 10 , 10 , 10 , // 21 - 30
562
56311 , 11 , 11 , 11 , 11 , 12 , 13 , 13 , 14 , 14 , // 31 - 40
564
56514 , 14 , 14 , 14 , 14 , 15 , 15 , 15 , 16 , 16 , // 41 - 50
566
567// ----------------------------------------------------------
568
56916 , 16 , 16 , 17 , 18 , 18 , 19 , 19 , 19 , 19 , // 51 - 60
570
57119 , 19 , 19 , 20 , 19 , 19 , 19 , 19 , 19 , 20 , // 61 - 70
572
57321 , 21 , 21 , 21 , 21 , 21 , 21 , 21 , 22 , 22 , // 71 - 80
574
57523 , 23 , 23 , 23 , 24 , 24 , 25 , 25 , 26 , 26 , // 81 - 90
576
57727 , 27 , 27 , 26 , 26 , 27 , 27 , 26 , 26 , 26 // 91 - 100
578
579};
580
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:44
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4Element::GetN
G4double GetN() const
Definition: G4Element.hh:134
G4Element::GetAtomicShell
G4double GetAtomicShell(G4int index) const
Definition: G4Element.cc:367
G4Material::GetDensity
G4double GetDensity() const
Definition: G4Material.hh:179
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:189
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:185
G4Material::GetAtomicNumDensityVector
const G4double * GetAtomicNumDensityVector() const
Definition: G4Material.hh:215
G4hQAOModel::ElectronicStoppingPower
G4double ElectronicStoppingPower(G4double z, G4double kineticEnergy) const
Definition: G4hQAOModel.cc:100
G4hQAOModel::StoppingPower
G4double StoppingPower(const G4Material *material, G4double kineticEnergy)
Definition: G4hQAOModel.cc:78