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