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