Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4AtomicShells.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//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
27
28// 24-04-98, implementation of the first version, V. Grichine
29// 16-11-98, GetBindingEnergy(Z,ShellNb), M.Maire
30// 19-07-04, add a protection in GetNumberOfShells(), mma
31// 11-02-05, GetNumberOfElectrons(Z,ShellNb), V.Ivanchenko
32// 30-04-10, added fIndexOfShells, V.Ivanchenko
33
34//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
35
36#include "G4AtomicShells.hh"
37#include "G4SystemOfUnits.hh"
38
39//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
40
41const G4int
42G4AtomicShells::fNumberOfShells[101] =
43{
44 0 , // nonexisting zero element
45
46 1 , 1 , 2 , 2 , 3 , 3 , 4 , 4 , 3 , 4 , // 1 - 10
47
48 5 , 5 , 6 , 6 , 6 , 6 , 6 , 7 , 8 , 8 , // 11 - 20
49
50 9 , 9 , 9 , 9 , 9 , 9 , 9 , 10 , 10 , 10 , // 21 - 30
51
5211 , 11 , 11 , 11 , 11 , 12 , 13 , 13 , 14 , 14 , // 31 - 40
53
5414 , 14 , 14 , 14 , 14 , 15 , 15 , 15 , 16 , 16 , // 41 - 50
55
56// ----------------------------------------------------------
57
5816 , 16 , 16 , 17 , 18 , 18 , 19 , 19 , 19 , 19 , // 51 - 60
59
6019 , 19 , 19 , 20 , 19 , 19 , 19 , 19 , 19 , 20 , // 61 - 70
61
6221 , 21 , 21 , 21 , 21 , 21 , 21 , 21 , 22 , 22 , // 71 - 80
63
6423 , 23 , 23 , 23 , 24 , 24 , 25 , 25 , 26 , 26 , // 81 - 90
65
6627 , 27 , 27 , 26 , 26 , 27 , 27 , 26 , 26 , 26 // 91 - 100
67
68};
69
70// The total shell number is:
71// 1 + G4AtomicShells::TotalNumberOfShells(100) = 1 + 1539 = 1540
72
73//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
74
75const G4int
76G4AtomicShells::fIndexOfShells[101] =
77 { 0,
78 1, 2, 3, 5, 7, 10, 13, 17, 21, 24,
79 28, 33, 38, 44, 50, 56, 62, 68, 75, 83,
80 91, 100, 109, 118, 127, 136, 145, 154, 164, 174,
81 184, 195, 206, 217, 228, 239, 251, 264, 277, 291,
82 305, 319, 333, 347, 361, 375, 390, 405, 420, 436,
83 452, 468, 484, 500, 517, 535, 553, 572, 591, 610,
84 629, 648, 667, 686, 706, 725, 744, 763, 782, 801,
85 821, 842, 863, 884, 905, 926, 947, 968, 989, 1011,
86 1033, 1056, 1079, 1102, 1125, 1149, 1173, 1198, 1223, 1249,
87 1275, 1302, 1329, 1356, 1382, 1408, 1435, 1462, 1488, 1514
88 };
89
90//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
91
92const G4double
93G4AtomicShells::fBindingEnergies[1540] =
94{
95 0.0 , // Nonexisting zero element
96
97 // H ---------------------------------------------------------
98 13.60 ,
99
100 // He ---------------------------------------------------------
101 24.59 ,
102
103 // Li --------------------------------------------------------
104 58.0 , 5.39 ,
105
106 // Be --------------------------------------------------------
107 115.0 , 9.32 ,
108
109 // B --------------------------------------------------------
110 192.0 , 12.93 , 8.3 ,
111
112 // C ---------------------------------------------------------
113 288.0 , 16.59 , 11.26 ,
114
115 // N ---------------------------------------------------------
116 403.0 , 37.3 , 20.33 , 14.53 ,
117
118 // O ---------------------------------------------------------
119 543.1 , 41.6 , 28.48 , 13.62 ,
120
121 // F ---------------------------------------------------------
122 696.7 , 37.85 , 17.42 ,
123
124 // Ne ---------------------------------------------------------
125 870.1 , 48.47 , 21.66 , 21.56 ,
126
127 // Na ---------------------------------------------------------
128 1075.0 , 66.0 , 34.0 , 34.0 , 5.14 ,
129
130 // Mg ---------------------------------------------------------
131 1308.0 , 92.0 , 54.0 , 54.0 , 7.65 ,
132
133 // Al ---------------------------------------------------------
134 1564.0 , 121. , 77.0 , 77.0 , 10.62 , 5.99 ,
135
136 // Si ---------------------------------------------------------
137 1844.0 , 154.0 , 104.0 , 104.0 , 13.46 , 8.15 ,
138
139 // P ---------------------------------------------------------
140 2148.0 , 191.0 , 135.0 , 134.0 , 16.15 , 10.49 ,
141
142 // S ---------------------------------------------------------
143 2476.0 , 232.0 , 170.0 , 168.0 , 20.20 , 10.36 ,
144
145 // Cl ---------------------------------------------------------
146 2829.0 , 277.0 , 208.0 , 206.0 , 24.54 , 12.97 ,
147
148 // Ar ---------------------------------------------------------
149 3206.3 , 326.5 , 250.6 , 248.5 , 29.24 , 15.94 , 15.76 ,
150
151 // K ---------------------------------------------------------
152 3610.0 , 381.0 , 299.0 , 296.0 , 37.0 , 19.0 , 18.7 , 4.34 ,
153
154 // Ca ---------------------------------------------------------
155 4041.0 , 441.0 , 353.0 , 349.0 , 46.0 , 28.0 , 28.0 , 6.11 ,
156
157 // Sc ---------------------------------------------------------
158 4494.0 , 503.0 , 408.0 , 403.0 , 55.0 , 33.0 , 33.0 , 8.0 , 6.54 ,
159
160 // Ti ---------------------------------------------------------
161 4966.0 , 567.0 , 465.0 , 459.0 , 64.0 , 39.0 , 38.0 , 8.0 , 6.82 ,
162
163 // V ---------------------------------------------------------
164 5465.0 , 633.0 , 525.0 , 518.0 , 72.0 , 44.0 , 43.0 , 8.0 , 6.74 ,
165
166 // Cr ---------------------------------------------------------
167 5989.0 , 702.0 , 589.0 , 580.0 , 80.0 , 49.0 , 48.0 , 8.25 , 6.77 ,
168
169 // Mn ---------------------------------------------------------
170 6539.0 , 755.0 , 656.0 , 645.0 , 89.0 , 55.0 , 53.0 , 9.0 , 7.43 ,
171
172 // Fe ---------------------------------------------------------
173 7112.0 , 851.0 , 726.0 , 713.0 , 98.0 , 61.0 , 59.0 , 9.0 , 7.87 ,
174
175 // Co ---------------------------------------------------------
176 7709.0 , 931.0 , 800.0 , 785.0 , 107.0 , 68.0 , 66.0 , 9.0 , 7.86 ,
177
178 // Ni ---------------------------------------------------------
179 8333.0 , 1015.0 , 877.0 , 860.0 , 117.0 , 75.0 , 73.0 , 10.0 , 10.0 ,
180 7.64 ,
181
182 // Cu ---------------------------------------------------------
183 8979.0 , 1103.0 , 958.0 , 938.0 , 127.0 , 82.0 , 80.0 , 11.0 , 10.4 ,
184 7.73 ,
185
186 // Zn ---------------------------------------------------------
187 9659.0 , 1198.0 , 1047.0 , 1024.0 , 141.0 , 94.0 , 91.0 , 12.0 , 11.2 ,
188 9.39 ,
189
190 // Ga ---------------------------------------------------------
191 10367.0 , 1302.0 , 1146.0 , 1119.0 , 162.0 , 111.0 , 107.0 , 21.0 ,
192 20.0 , 11.0 , 6.0 ,
193
194 // Ge ---------------------------------------------------------
195 11103.0 , 1413.0 , 1251.0 , 1220.0 , 184.0 , 130.0 , 125.0 , 33.0 ,
196 32.0 , 14.3 , 7.9 ,
197
198 // As ---------------------------------------------------------
199 11867.0 , 1531.0 , 1362.0 , 1327.0 , 208.0 , 151.0 , 145.0 , 46.0 ,
200 45.0 , 17.0 , 9.81 ,
201
202 // Se ---------------------------------------------------------
203 12658.0 , 1656.0 , 1479.0 , 1439.0 , 234.0 , 173.0 , 166.0 , 61.0 ,
204 60.0 , 20.15 , 9.75 ,
205
206 // Br ---------------------------------------------------------
207 13474.0 , 1787.0 , 1602.0 , 1556.0 , 262.0 , 197.0 , 189.0 , 77.0 ,
208 76.0 , 23.8 , 11.85 ,
209
210 // Kr ---------------------------------------------------------
211 14326.0 , 1924.6 , 1730.9 , 1678.4 , 292.8 , 222.2 , 214.4 , 95.0 ,
212 93.8 , 27.51 , 14.65 , 14.0 ,
213
214 // Rb ---------------------------------------------------------
215 15200.0 , 2068.0 , 1867.0 , 1807.0 , 325.0 , 251.0 , 242.0 , 116.0 ,
216 114.0 , 32.0 , 16.0 , 15.3 , 4.18 ,
217
218 // Sr ---------------------------------------------------------
219 16105.0 , 2219.0 , 2010.0 , 1943.0 , 361.0 , 283.0 , 273.0 , 139.0 ,
220 137.0 , 40.0 , 23.0 , 22.0 , 5.69 ,
221
222 // Y ---------------------------------------------------------
223 17038.0 , 2375.0 , 2158.0 , 2083.0 , 397.0 , 315.0 , 304.0 , 163.0 ,
224 161.0 , 48.0 , 30.0 , 29.0 , 6.48, 6.38 ,
225
226 // Zr ---------------------------------------------------------
227 17998.0 , 2536.0 , 2311.0 , 2227.0 , 434.0 , 348.0 , 335.0 , 187.0 ,
228 185.0 , 56.0 , 35.0 , 33.0 , 8.61 , 6.84 ,
229
230 // Nb ----------------------------------------------------------
231 18986.0 , 2702.0 , 2469.0 , 2375.0 , 472.0 , 382.0 , 367.0 , 212.0 ,
232 209.0 , 62.0 , 40.0 , 38.0 , 7.17 , 6.88 ,
233
234 // Mo -----------------------------------------------------------
235 20000.0 , 2872.0 , 2632.0 , 2527.0 , 511.0 , 416.0 , 399.0 , 237.0 ,
236 234.0 , 68.0 , 45.0 , 42.0 , 8.56 , 7.10 ,
237
238 // Tc -----------------------------------------------------------
239 21044.0 , 3048.0 , 2800.0 , 2683.0 , 551.0 , 451.0 , 432.0 , 263.0 ,
240 259.0 , 74.0 , 49.0 , 45.0 , 8.6 , 7.28 ,
241
242 // Ru -----------------------------------------------------------
243 22117.0 , 3230.0 , 2973.0 , 2844.0 , 592.0 , 488.0 , 466.0 , 290.0 ,
244 286.0 , 81.0 , 53.0 , 49.0 , 8.5 , 7.37 ,
245
246 // Rh -----------------------------------------------------------
247 23220.0 , 3418.0 , 3152.0 , 3010.0 , 634.0 , 526.0 , 501.0 , 318.0 ,
248 313.0 , 87.0 , 58.0 , 53.0 , 9.56 , 7.46 ,
249
250 // Pd -----------------------------------------------------------
251 24350.0 , 3611.0 , 3337.0 , 3180.0 , 677.0 , 565.0 , 537.0 , 347.0 ,
252 342.0 , 93.0 , 63.0 , 57.0 , 8.78 , 8.34 , 7.52 ,
253
254 // Ag -----------------------------------------------------------
255 25514.0 , 3812.0 , 3530.0 , 3357.0 , 724.0 , 608.0 , 577.0 , 379.0 ,
256 373.0 , 101.0 , 69.0 , 63.0 , 11.0 , 10.0 , 7.58 ,
257
258 // Cd ------------------------------------------------------------
259 26711.0 , 4022.0 , 3732.0 , 3542.0 , 775.0 , 655.0 , 621.0 , 415.0 ,
260 408.0 , 112.0 , 78.0 , 71.0 , 14.0 , 13.0 , 8.99 ,
261
262 // In ------------------------------------------------------------
263 27940.0 , 4242.0 , 3943.0 , 3735.0 , 830.0 , 707.0 , 669.0 , 455.0 ,
264 447.0 , 126.0 , 90.0 , 82.0 , 21.0 , 20.0 , 10.0 , 5.79 ,
265
266 // Sn ------------------------------------------------------------
267 29200.0 , 4469.0 , 4160.0 , 3933.0 , 888.0 , 761.0 , 719.0 , 497.0 ,
268 489.0 , 141.0 , 102.0 , 93.0 , 29.0 , 28.0 , 12.0 , 7.34 ,
269
270 // Sb ------------------------------------------------------------
271 30419.0 , 4698.0 , 4385.0 , 4137.0 , 949.0 , 817.0 , 771.0 , 542.0 ,
272 533.0 , 157.0 , 114.0 , 104.0 , 38.0 , 37.0 , 15.0 , 8.64 ,
273
274 // Te ------------------------------------------------------------
275 31814.0 , 4939.0 , 4612.0 , 4347.0 , 1012.0 , 876.0 , 825.0 , 589.0 ,
276 578.0 , 174.0 , 127.0 , 117.0 , 48.0 , 46.0 , 17.84 , 9.01 ,
277
278 // I ------------------------------------------------------------
279 33169.0 , 5188.0 , 4852.0 , 4557.0 , 1078.0 , 937.0 , 881.0 , 638.0 ,
280 626.0 , 193.0 , 141.0 , 131.0 , 58.0 , 56.0 , 20.61 , 10.45 ,
281
282 // Xe ------------------------------------------------------------
283 34570.0 , 5460.0 , 5110.0 , 4790.0 , 1148.7 , 1002.1 , 940.6 , 689.0 ,
284 676.4 , 213.2 , 157.0 , 145.5 , 69.5 , 67.5 , 23.39 , 13.43 ,
285 12.13 ,
286
287 // Cs ------------------------------------------------------------
288 35985.0 , 5714.0 , 5359.0 , 5012.0 , 1220.0 , 1068.0 , 1000.0 , 742.0 ,
289 728.0 , 233.0 , 174.0 , 164.0 , 81.0 , 79.0 , 25.0 , 14.0 ,
290 12.3 , 3.89 ,
291
292 // Ba -------------------------------------------------------------
293 37441.0 , 5989.0 , 5624.0 , 5247.0 , 1293.0 , 1138.0 , 1063.0 , 797.0 ,
294 782.0 , 254.0 , 193.0 , 181.0 , 94.0 , 92.0 , 31.0 , 18.0 ,
295 16.0 , 5.21 ,
296
297 // La -------------------------------------------------------------
298 38925.0 , 6266.0 , 5891.0 , 5483.0 , 1365.0 , 1207.0 , 1124.0 , 851.0 ,
299 834.0 , 273.0 , 210.0 , 196.0 , 105.0 , 103.0 , 36.0 , 22.0 ,
300 19.0 , 5.75 , 5.58 ,
301
302 // Ce -------------------------------------------------------------
303 40443.0 , 6548.0 , 6164.0 , 5723.0 , 1437.0 , 1275.0 , 1184.0 , 903.0 ,
304 885.0 , 291.0 , 225.0 , 209.0 , 114.0 , 111.0 , 39.0 , 25.0 ,
305 22.0 , 6.0 , 5.65 ,
306
307 // Pr -------------------------------------------------------------
308 41991.0 , 6835.0 , 6440.0 , 5964.0 , 1509.0 , 1342.0 , 1244.0 , 954.0 ,
309 934.0 , 307.0 , 238.0 , 220.0 , 121.0 , 117.0 , 41.0 ,
310 27.0 , 24.0 , 6.0 , 5.42 ,
311
312 // Nd -------------------------------------------------------------
313 43569.0 , 7126.0 , 6722.0 , 6208.0 , 1580.0 , 1408.0 , 1303.0 , 1005.0 ,
314 983.0 , 321.0 , 250.0 , 230.0 , 126.0 , 122.0 , 42.0 ,
315 28.0 , 25.0 , 6.0 , 5.49 ,
316
317 // Pm -------------------------------------------------------------
318 45184.0 , 7428.0 , 7013.0 , 6459.0 , 1653.0 , 1476.0 , 1362.0 , 1057.0 ,
319 1032.0 , 325.0 , 261.0 , 240.0 , 131.0 , 127.0 , 43.0 ,
320 28.0 , 25.0 , 6.0 , 5.55 ,
321
322 // Sm -------------------------------------------------------------
323 46834.0 , 7737.0 , 7312.0 , 6716.0 , 1728.0 , 1546.0 , 1422.0 , 1110.0 ,
324 1083.0 , 349.0 , 273.0 , 251.0 , 137.0 , 132.0 , 44.0 ,
325 29.0 , 25.0 , 6.0 , 5.63 ,
326
327 // Eu -------------------------------------------------------------
328 48519.0 , 8052.0 , 7617.0 , 6977.0 , 1805.0 , 1618.0 , 1484.0 , 1164.0 ,
329 1135.0 , 364.0 , 286.0 , 262.0 , 143.0 , 137.0 , 45.0 ,
330 30.0 , 26.0 , 6.0 , 5.68 ,
331
332 // Gd -------------------------------------------------------------
333 50239.0 , 8376.0 , 7930.0 , 7243.0 , 1884.0 , 1692.0 , 1547.0 , 1220.0 ,
334 1189.0 , 380.0 , 300.0 , 273.0 , 150.0 , 143.0 , 46.0 ,
335 31.0 , 27.0 , 6.16 , 6.0 , 6.0 ,
336
337 // Tb -------------------------------------------------------------
338 51996.0 , 8708.0 , 8252.0 , 7514.0 , 1965.0 , 1768.0 , 1612.0 , 1277.0 ,
339 1243.0 , 398.0 , 315.0 , 285.0 , 157.0 , 150.0 , 48.0 ,
340 32.0 , 28.0 , 6.0 , 5.85 ,
341
342 // Dy -------------------------------------------------------------
343 53789.0 , 9046.0 , 8581.0 , 7790.0 , 2048.0 , 1846.0 , 1678.0 , 1335.0 ,
344 1298.0 , 416.0 , 331.0 , 297.0 , 164.0 , 157.0 , 50.0 ,
345 33.0 , 28.0 , 6.0 , 5.93 ,
346
347 // Ho -------------------------------------------------------------
348 55618.0 , 9394.0 , 8918.0 , 8071.0 , 2133.0 , 1926.0 , 1746.0 , 1395.0 ,
349 1354.0 , 434.0 , 348.0 , 310.0 , 172.0 , 164.0 , 52.0 ,
350 34.0 , 29.0 , 6.02 , 6.0 ,
351
352 // Er -------------------------------------------------------------
353 57486.0 , 9751.0 , 9264.0 , 8358.0 , 2220.0 , 2008.0 , 1815.0 , 1456.0 ,
354 1412.0 , 452.0 , 365.0 , 323.0 , 181.0 , 172.0 , 54.0 ,
355 35.0 , 30.0 , 6.10 , 6.0 ,
356
357 // Tu -------------------------------------------------------------
358 59390.0 , 10116.0 , 9617.0 , 8648.0 , 2309.0 , 2092.0 , 1885.0 , 1518.0 ,
359 1471.0 , 471.0 , 382.0 , 336.0 , 190.0 , 181.0 , 56.0 ,
360 36.0 , 30.0 , 7.0 , 6.18 ,
361
362 // Yb -------------------------------------------------------------
363 61332.0 , 10486.0 , 9978.0 , 8944.0 , 2401.0 , 2178.0 , 1956.0 , 1580.0 ,
364 1531.0 , 490.0 , 399.0 , 349.0 , 200.0 , 190.0 ,
365 58.0 , 37.0 , 31.0 , 8.0 , 7.0 , 6.25 ,
366
367 // Lu -------------------------------------------------------------
368 63314.0 , 10870.0 , 10349.0 , 9244.0 , 2499.0 , 2270.0 , 2032.0 , 1647.0 ,
369 1596.0 , 514.0 , 420.0 , 366.0 , 213.0 , 202.0 ,
370 62.0 , 39.0 , 32.0 , 13.0 , 12.0 , 7.0 , 6.6 ,
371
372 // Hf -------------------------------------------------------------
373 65351.0 , 11271.0 , 10739.0 , 9561.0 , 2604.0 , 2369.0 , 2113.0 , 1720.0 ,
374 1665.0 , 542.0 , 444.0 , 386.0 , 229.0 , 217.0 ,
375 68.0 , 43.0 , 35.0 , 21.0 , 20.0 , 7.5 , 7.0 ,
376
377 // Ta -------------------------------------------------------------
378 67416.0 , 11682.0 , 11136.0 , 9881.0 , 2712.0 , 2472.0 , 2197.0 , 1796.0 ,
379 1737.0 , 570.0 , 469.0 , 407.0 , 245.0 , 232.0 ,
380 74.0 , 47.0 , 38.0 , 30.0 , 28.0 , 8.3 , 7.9 ,
381
382 // W -------------------------------------------------------------
383 69525.0 , 12100.0 , 11544.0 , 10207.0 , 2823.0 , 2577.0 , 2283.0 , 1874.0 ,
384 1811.0 , 599.0 , 495.0 , 428.0 , 261.0 , 248.0 ,
385 80.0 , 51.0 , 41.0 , 38.0 , 36.0 , 9.0 , 8.0 ,
386
387 // Re -------------------------------------------------------------
388 71676.0 , 12527.0 , 11959.0 , 10535.0 , 2937.0 , 2686.0 , 2371.0 , 1953.0 ,
389 1887.0 , 629.0 , 522.0 , 450.0 , 278.0 , 264.0 ,
390 86.0 , 56.0 , 47.0 , 45.0 , 45.0 , 9.6 , 7.9 ,
391
392 // Os -------------------------------------------------------------
393 73871.0 , 12968.0 , 12385.0 , 10871.0 , 3054.0 , 2797.0 , 2461.0 , 2035.0 ,
394 1964.0 , 660.0 , 551.0 , 473.0 , 295.0 , 280.0 ,
395 92.0 , 61.0 , 56.0 , 54.0 , 49.0 , 9.6 , 8.5 ,
396
397 // Ir -------------------------------------------------------------
398 76111.0 , 13419.0 , 12824.0 , 11215.0 , 3175.0 , 2912.0 , 2554.0 , 2119.0 ,
399 2044.0 , 693.0 , 581.0 , 497.0 , 314.0 , 298.0 ,
400 99.0 , 67.0 , 66.0 , 64.0 , 53.0 , 9.6 , 9.1 ,
401
402 // Pt -------------------------------------------------------------
403 78395.0 , 13880.0 , 13273.0 , 11564.0 , 3300.0 , 3030.0 , 2649.0 , 2206.0 ,
404 2126.0 , 727.0 , 612.0 , 522.0 , 335.0 , 318.0 ,
405 106.0 , 78.0 , 75.0 , 71.0 , 57.0 , 9.6 , 9.0 ,
406
407 // Au -------------------------------------------------------------
408 80725.0 , 14353.0 , 13734.0 , 11919.0 , 3430.0 , 3153.0 , 2748.0 , 2295.0 ,
409 2210.0 , 764.0 , 645.0 , 548.0 , 357.0 , 339.0 ,
410 114.0 , 91.0 , 87.0 , 76.0 , 61.0 , 12.5 , 11.1 , 9.23 ,
411
412 // Hg -------------------------------------------------------------
413 83102.0 , 14839.0 , 14209.0 , 12284.0 , 3567.0 , 3283.0 , 2852.0 , 2390.0 ,
414 2300.0 , 806.0 , 683.0 , 579.0 , 382.0 , 363.0 ,
415 125.0 , 107.0 , 103.0 , 85.0 , 68.0 , 14.0 , 12.0 , 10.4 ,
416
417 // Tl -------------------------------------------------------------
418 85530.0 , 15347.0 , 14698.0 , 12658.0 , 3710.0 , 3420.0 , 2961.0 , 2490.0 ,
419 2394.0 , 852.0 , 726.0 , 615.0 , 411.0 , 391.0 ,
420 139.0 , 127.0 , 123.0 , 98.0 , 79.0 , 21.0 , 19.0 , 8.0 ,
421 6.11 ,
422
423 // Pb -------------------------------------------------------------
424 88005.0 , 15861.0 , 15200.0 , 13055.0 , 3857.0 , 3560.0 , 3072.0 , 2592.0 ,
425 2490.0 , 899.0 , 769.0 , 651.0 , 441.0 , 419.0 ,
426 153.0 , 148.0 , 144.0 , 111.0 , 90.0 , 27.0 , 25.0 , 10.0 ,
427 7.42 ,
428
429 // Bi -------------------------------------------------------------
430 90526.0 , 16388.0 , 15711.0 , 13419.0 , 4007.0 , 3704.0 , 3185.0 , 2696.0 ,
431 2588.0 , 946.0 , 813.0 , 687.0 , 472.0 , 448.0 , 170.0 ,
432 167.0 , 165.0 , 125.0 , 101.0 , 34.0 , 32.0 , 12.0 , 7.29 ,
433
434 // Po -------------------------------------------------------------
435 93105.0 , 16939.0 , 16244.0 , 13814.0 , 4161.0 , 3852.0 , 3301.0 , 2802.0 ,
436 2687.0 , 994.0 , 858.0 , 724.0 , 503.0 , 478.0 , 193.0 , 187.0 ,
437 181.0 , 139.0 , 112.0 , 41.0 , 38.0 , 15.0 , 8.43 ,
438
439 // At -------------------------------------------------------------
440 95730.0 , 17493.0 , 16785.0 , 14214.0 , 4320.0 , 4005.0 , 3420.0 , 2910.0 ,
441 2788.0 , 1044.0 , 904.0 , 761.0 , 535.0 , 508.0 , 217.0 , 211.0 ,
442 196.0 , 153.0 , 123.0 , 48.0 , 44.0 , 19.0 , 11.0 , 9.3 ,
443
444 // Rn -------------------------------------------------------------
445 98404.0 , 18049.0 , 17337.0 , 14619.0 , 4483.0 , 4162.0 , 3452.0 , 3109.0 ,
446 2890.0 , 1096.0 , 951.0 , 798.0 , 567.0 , 538.0 , 242.0 , 235.0 ,
447 212.0 , 167.0 , 134.0 , 55.0 , 51.0 , 24.0 , 14.0 , 10.7 ,
448
449 // Fr -------------------------------------------------------------
450 101137.0 , 18639.0 , 17907.0 , 15031.0 , 4652.0 , 4324.0 , 3666.0 , 3134.0 ,
451 2998.0 , 1153.0 , 1003.0 , 839.0 , 603.0 , 572.0 , 268.0 , 260.0 ,
452 231.0 , 183.0 , 147.0 , 65.0 , 61.0 , 33.0 , 19.0 , 14.0 ,
453 4.0 ,
454
455 // Ra -------------------------------------------------------------
456 103922.0 , 19237.0 , 18484.0 , 15444.0 , 4822.0 , 4491.0 , 3793.0 , 3254.0 ,
457 3111.0 , 1214.0 , 1060.0 , 884.0 , 642.0 , 609.0 , 296.0 , 287.0 ,
458 253.0 , 201.0 , 161.0 , 77.0 , 73.0 , 40.0 , 25.0 , 19.0 ,
459 5.28 ,
460
461 // Ac -------------------------------------------------------------
462 106755.0 , 19840.0 , 19083.0 , 15871.0 , 5002.0 , 4656.0 , 3921.0 , 3374.0 ,
463 3223.0 , 1274.0 , 1116.0 , 928.0 , 680.0 , 645.0 , 322.0 , 313.0 ,
464 274.0 , 218.0 , 174.0 , 88.0 , 83.0 , 45.0 , 29.0 , 22.0 ,
465 6.3 , 5.7 ,
466
467 // Th -------------------------------------------------------------
468 109651.0 , 20472.0 , 19693.0 , 16300.0 , 5182.0 , 4830.0 , 4049.0 , 3494.0 ,
469 3335.0 , 1333.0 , 1171.0 , 970.0 , 717.0 , 679.0 , 347.0 , 338.0 ,
470 293.0 , 233.0 , 185.0 , 97.0 , 91.0 , 50.0 , 33.0 , 25.0 ,
471 6.0 , 6.0 ,
472
473 // Pa -------------------------------------------------------------
474 112601.0 , 21105.0 , 20314.0 , 16733.0 , 5367.0 , 5001.0 , 4178.0 , 3613.0 ,
475 3446.0 , 1390.0 , 1225.0 , 1011.0 , 752.0 , 712.0 , 372.0 , 362.0 ,
476 312.0 , 248.0 , 195.0 , 104.0 , 97.0 , 50.0 , 32.0 ,
477 24.0 , 6.0 , 6.0 , 6.0 ,
478
479 // U -------------------------------------------------------------
480 115606.0 , 21757.0 , 20948.0 , 17166.0 , 5548.0 , 5182.0 , 4308.0 , 3733.0 ,
481 3557.0 , 1446.0 , 1278.0 , 1050.0 , 785.0 , 743.0 , 396.0 , 386.0 ,
482 329.0 , 261.0 , 203.0 , 110.0 , 101.0 , 52.0 , 34.0 ,
483 24.0 , 6.1 , 6.0 , 6.0 ,
484
485 // Np -------------------------------------------------------------
486 118678.0 , 22426.0 , 21600.0 , 17610.0 , 5723.0 , 5366.0 , 4440.0 , 3854.0 ,
487 3669.0 , 1504.0 , 1331.0 , 1089.0 , 819.0 , 774.0 , 421.0 , 410.0 ,
488 346.0 , 274.0 , 211.0 , 116.0 , 106.0 , 54.0 , 35.0 ,
489 25.0 , 6.0 , 6.0 , 6.0 ,
490
491 // Pu -------------------------------------------------------------
492 121818.0 , 23097.0 , 22266.0 , 18056.0 , 5933.0 , 5541.0 , 4557.0 , 3977.0 ,
493 3783.0 , 1563.0 , 1384.0 , 1128.0 , 853.0 , 805.0 , 446.0 , 434.0 ,
494 356.0 , 287.0 , 219.0 , 122.0 , 111.0 , 53.0 , 34.0 ,
495 23.0 , 6.0 , 6.0 ,
496
497 // Am -------------------------------------------------------------
498 125027.0 , 23773.0 , 22944.0 , 18504.0 , 6121.0 , 5710.0 , 4667.0 , 4102.0 ,
499 3898.0 , 1623.0 , 1439.0 , 1167.0 , 887.0 , 836.0 , 467.0 , 452.0 ,
500 355.0 , 301.0 , 220.0 , 123.0 , 112.0 , 54.0 , 44.0 ,
501 36.0 , 6.0 , 6.0 ,
502
503 // Cm -------------------------------------------------------------
504 128220.0 , 24460.0 , 23779.0 , 18930.0 , 6288.0 , 5895.0 , 4797.0 , 4236.0 ,
505 4014.0 , 1664.0 , 1493.0 , 1194.0 , 919.0 , 864.0 , 494.0 , 479.0 ,
506 384.0 , 314.0 , 239.0 , 126.0 , 119.0 , 60.0 , 39.0 ,
507 27.0 , 11.0 , 5.0 , 6.0 ,
508
509 // Bk -------------------------------------------------------------
510 131590.0 , 25275.0 , 24385.0 , 19452.0 , 6556.0 , 6147.0 , 4977.0 , 4366.0 ,
511 4133.0 , 1729.0 , 1554.0 , 1236.0 , 955.0 , 898.0 , 520.0 , 504.0 ,
512 401.0 , 329.0 , 248.0 , 142.0 , 124.0 , 63.0 , 41.0 ,
513 27.0 , 12.0 , 6.0 , 4.0 ,
514
515 // Cf -------------------------------------------------------------
516 135960.0 , 26110.0 , 25250.0 , 19930.0 , 6754.0 , 6359.0 , 5109.0 , 4492.0 ,
517 4247.0 , 1789.0 , 1610.0 , 1273.0 , 987.0 , 925.0 , 546.0 , 529.0 ,
518 412.0 , 338.0 , 251.0 , 142.0 , 129.0 , 61.0 , 39.0 ,
519 25.0 , 9.0 , 6.0 ,
520
521 // Es -------------------------------------------------------------
522 139490.0 , 26900.0 , 26020.0 , 20410.0 , 6977.0 , 6754.0 , 5252.0 , 4630.0 ,
523 4369.0 , 1857.0 , 1674.0 , 1316.0 , 1024.0 , 959.0 , 573.0 , 554.0 ,
524 429.0 , 353.0 , 260.0 , 148.0 , 135.0 , 63.0 , 40.0 ,
525 25.0 , 9.0 , 6.0 ,
526
527 // Fm -------------------------------------------------------------
528 143090.0 , 27700.0 , 26810.0 , 20900.0 , 7205.0 , 6793.0 , 5397.0 , 4766.0 ,
529 4498.0 , 1933.0 , 1746.0 , 1366.0 , 1068.0 , 1000.0 , 606.0 , 587.0 ,
530 453.0 , 375.0 , 275.0 , 160.0 , 145.0 , 69.0 , 45.0 ,
531 29.0 , 15.0 , 7.0 ,
532
533};
534
535const G4int
536G4AtomicShells::fNumberOfElectrons[1540] =
537{
538 0, // Nonexisting zero element
539 // H ---------------------------------------------------------
540 1,
541 // He ---------------------------------------------------------
542 2,
543 // Li --------------------------------------------------------
544 2, 1,
545 // Be --------------------------------------------------------
546 2, 2,
547 // B --------------------------------------------------------
548 2, 2, 1,
549 // C ---------------------------------------------------------
550 2, 2, 2,
551 // N ---------------------------------------------------------
552 2, 2, 2, 1,
553 // O ---------------------------------------------------------
554 2, 2, 2, 2,
555 // F ---------------------------------------------------------
556 2, 2, 5,
557 // Ne ---------------------------------------------------------
558 2, 2, 2, 4,
559 // Na ---------------------------------------------------------
560 2, 2, 2, 4, 1,
561 // Mg ---------------------------------------------------------
562 2, 2, 2, 4, 2,
563 // Al ---------------------------------------------------------
564 2, 2, 2, 4, 2, 1,
565 // Si ---------------------------------------------------------
566 2, 2, 2, 4, 2, 2,
567 // P ---------------------------------------------------------
568 2, 2, 2, 4, 2, 3,
569 // S ---------------------------------------------------------
570 2, 2, 2, 4, 2, 4,
571 // Cl ---------------------------------------------------------
572 2, 2, 2, 4, 2, 5,
573 // Ar ---------------------------------------------------------
574 2, 2, 2, 4, 2, 2, 4,
575 // K ---------------------------------------------------------
576 2, 2, 2, 4, 2, 2, 4, 1,
577 // Ca ---------------------------------------------------------
578 2, 2, 2, 4, 2, 2, 4, 2,
579 // Sc ---------------------------------------------------------
580 2, 2, 2, 4, 2, 2, 4, 1, 2,
581 // Ti ---------------------------------------------------------
582 2, 2, 2, 4, 2, 2, 4, 2, 2,
583 // V ---------------------------------------------------------
584 2, 2, 2, 4, 2, 2, 4, 3, 2,
585 // Cr ---------------------------------------------------------
586 2, 2, 2, 4, 2, 2, 4, 4, 2,
587 // Mn ---------------------------------------------------------
588 2, 2, 2, 4, 2, 2, 4, 5, 2,
589 // Fe ---------------------------------------------------------
590 2, 2, 2, 4, 2, 2, 4, 6, 2,
591 // Co ---------------------------------------------------------
592 2, 2, 2, 4, 2, 2, 4, 7, 2,
593 // Ni ---------------------------------------------------------
594 2, 2, 2, 4, 2, 2, 4, 4, 4, 2,
595 // Cu ---------------------------------------------------------
596 2, 2, 2, 4, 2, 2, 4, 4, 5, 2,
597 // Zn ---------------------------------------------------------
598 2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
599 // Ga ---------------------------------------------------------
600 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 1,
601 // Ge ---------------------------------------------------------
602 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2,
603 // As ---------------------------------------------------------
604 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 3,
605 // Se ---------------------------------------------------------
606 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 4,
607 // Br ---------------------------------------------------------
608 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 5,
609 // Kr ---------------------------------------------------------
610 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4,
611 // Rb ---------------------------------------------------------
612 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 1,
613 // Sr ---------------------------------------------------------
614 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2,
615 // Y ---------------------------------------------------------
616 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 1,
617 // Zr ---------------------------------------------------------
618 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,
619 // Nb ----------------------------------------------------------
620 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 3, 2,
621 // Mo -----------------------------------------------------------
622 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 2,
623 // Tc -----------------------------------------------------------
624 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 5, 2,
625 // Ru -----------------------------------------------------------
626 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,
627 // Rh -----------------------------------------------------------
628 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,
629 // Pd -----------------------------------------------------------
630 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 4, 2,
631 // Ag -----------------------------------------------------------
632 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 5, 2,
633 // Cd ------------------------------------------------------------
634 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2,
635 // In ------------------------------------------------------------
636 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 1,
637 // Sn ------------------------------------------------------------
638 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2,
639 // Sb ------------------------------------------------------------
640 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 3,
641 // Te ------------------------------------------------------------
642 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 4,
643 // I ------------------------------------------------------------
644 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 5,
645 // Xe ------------------------------------------------------------
646 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4,
647 // Cs ------------------------------------------------------------
648 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1,
649 // Ba -------------------------------------------------------------
650 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2,
651 // La -------------------------------------------------------------
652 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 1, 2,
653 // Ce -------------------------------------------------------------
654 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2, 2,
655 // Pr -------------------------------------------------------------
656 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 3, 2,
657 // Nd -------------------------------------------------------------
658 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 4, 2,
659 // Pm -------------------------------------------------------------
660 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 5, 2,
661 // Sm -------------------------------------------------------------
662 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,
663 // Eu -------------------------------------------------------------
664 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,
665 // Gd -------------------------------------------------------------
666 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2, 7, 1,
667 // Tb -------------------------------------------------------------
668 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 9, 2,
669 // Dy -------------------------------------------------------------
670 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 10,2,
671 // Ho -------------------------------------------------------------
672 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2,11,
673 // Er -------------------------------------------------------------
674 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 2,12,
675 // Tu -------------------------------------------------------------
676 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4,13, 2,
677 // Yb -------------------------------------------------------------
678 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 8, 2,
679 // Lu -------------------------------------------------------------
680 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 8, 2, 1,
681 // Hf -------------------------------------------------------------
682 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 8, 2, 2,
683 // Ta -------------------------------------------------------------
684 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 8, 3, 2,
685 // W -------------------------------------------------------------
686 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 4, 6, 8, 4, 2,
687 // Re -------------------------------------------------------------
688 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 6, 8, 4, 5, 2,
689 // Os -------------------------------------------------------------
690 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 2, 6, 8, 4, 6, 2,
691 // Ir -------------------------------------------------------------
692 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 2, 8, 4, 7, 2,
693 // Pt -------------------------------------------------------------
694 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 8, 2, 4, 9, 1,
695 // Au -------------------------------------------------------------
696 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 8, 2, 4, 4, 6, 1,
697 // Hg -------------------------------------------------------------
698 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 8, 2, 4, 4, 6, 2,
699 // Tl -------------------------------------------------------------
700 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 8, 2, 4, 4, 6, 2, 1,
701 // Pb -------------------------------------------------------------
702 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 2, 6, 8, 2, 4, 4, 6, 2, 2,
703 // Bi -------------------------------------------------------------
704 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 2, 8, 2, 4, 4, 6, 2, 3,
705 // Po -------------------------------------------------------------
706 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 4,
707 // At -------------------------------------------------------------
708 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 3,
709 // Rn -------------------------------------------------------------
710 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,
711 // Fr -------------------------------------------------------------
712 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,
713 // Ra -------------------------------------------------------------
714 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,
715 // Ac -------------------------------------------------------------
716 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, 1,
717 // Th -------------------------------------------------------------
718 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,
719 // Pa -------------------------------------------------------------
720 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, 1, 2,
721 // U -------------------------------------------------------------
722 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, 3, 2,
723 // Np -------------------------------------------------------------
724 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 4, 1, 2,
725 // Pu -------------------------------------------------------------
726 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 6, 2,
727 // Am -------------------------------------------------------------
728 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2,
729 // Cm -------------------------------------------------------------
730 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 7, 2, 1,
731 // Bk -------------------------------------------------------------
732 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4, 8, 2, 1,
733 // Cf -------------------------------------------------------------
734 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,10, 2,
735 // Es -------------------------------------------------------------
736 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,11, 2,
737 // Fm -------------------------------------------------------------
738 2, 2, 2, 4, 2, 2, 4, 4, 6, 2, 2, 4, 4, 6, 6, 8, 2, 2, 4, 4, 6, 2, 2, 4,12, 2
739};
740
741//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
742
743G4int
744G4AtomicShells::GetNumberOfShells(G4int Z)
745{
746#ifdef G4VERBOSE
747 if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetNumberOfShells"); }
748#endif
749 return fNumberOfShells[Z];
750}
751
752//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
753
754G4double
755G4AtomicShells::GetBindingEnergy(G4int Z, G4int ShellNb)
756{
757#ifdef G4VERBOSE
758 if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetBindingEnergy"); }
759 if(ShellNb < 0 || ShellNb>=fNumberOfShells[Z]) {
760 ShellNb = PrintErrorShell(Z, ShellNb, "GetBindingEnergy");
761 }
762#endif
763 return fBindingEnergies[fIndexOfShells[Z] + ShellNb]*CLHEP::eV;
764}
765
766//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
767
768G4int
769G4AtomicShells::GetNumberOfElectrons(G4int Z, G4int ShellNb)
770{
771#ifdef G4VERBOSE
772 if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetNumberOfElectrons"); }
773 if(ShellNb < 0 || ShellNb>=fNumberOfShells[Z]) {
774 ShellNb = PrintErrorShell(Z, ShellNb, "GetNumberOfElectrons");
775 }
776#endif
777 return fNumberOfElectrons[fIndexOfShells[Z] + ShellNb];
778}
779
780//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
781
782G4int
783G4AtomicShells::GetNumberOfFreeElectrons(G4int Z, G4double th)
784{
785#ifdef G4VERBOSE
786 if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetNumberOfFreeElectrons"); }
787#endif
788 G4int idx = fIndexOfShells[Z];
789 G4int idxmax = idx + fNumberOfShells[Z];
790 G4int n = 0;
791 for (G4int i=idx; i<idxmax; ++i) {
792 if(fBindingEnergies[i] <= th) { n += fNumberOfElectrons[i]; }
793 }
794 return n;
795}
796
797//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
798
799G4double G4AtomicShells::GetTotalBindingEnergy(G4int Z)
800{
801#ifdef G4VERBOSE
802 if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetTotalBindingEnergy"); }
803#endif
804 G4int idx = fIndexOfShells[Z];
805 G4int idxmax = idx + fNumberOfShells[Z];
806 G4double energy = 0.0;
807 for (G4int i=idx; i<idxmax; ++i) {energy += fBindingEnergies[i];}
808 return energy*CLHEP::eV;
809}
810
811//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
812
813#ifdef G4VERBOSE
814
815G4int G4AtomicShells::PrintErrorZ(G4int Z, const G4String& ss)
816{
817 G4String sss = "G4AtomicShells::"+ss+"()";
818 G4ExceptionDescription ed;
819 ed << "Atomic number out of range Z= " << Z;
820 G4Exception(sss,"mat060",FatalException,ed,"");
821 return 1;
822}
823
824//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
825
826G4int
827G4AtomicShells::PrintErrorShell(G4int Z, G4int n, const G4String& ss)
828{
829 G4String sss = "G4AtomicShells::"+ss+"()";
830 G4ExceptionDescription ed;
831 ed << "Shell number out of range Nshell= " << n << " Z= " << Z;
832 G4Exception(sss,"mat061",FatalException,ed,"");
833 return 0;
834}
835#endif
836
837//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4AtomicShells::GetTotalBindingEnergy
static G4double GetTotalBindingEnergy(G4int Z)
Definition: G4AtomicShells.cc:799
G4AtomicShells::GetNumberOfElectrons
static G4int GetNumberOfElectrons(G4int Z, G4int SubshellNb)
Definition: G4AtomicShells.cc:769
G4AtomicShells::GetBindingEnergy
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
Definition: G4AtomicShells.cc:755
G4AtomicShells::GetNumberOfShells
static G4int GetNumberOfShells(G4int Z)
Definition: G4AtomicShells.cc:744
G4AtomicShells::GetNumberOfFreeElectrons
static G4int GetNumberOfFreeElectrons(G4int Z, G4double th)
Definition: G4AtomicShells.cc:783