ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4AtomicShells.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4AtomicShells.cc
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 
41 const G4int
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 
52 11 , 11 , 11 , 11 , 11 , 12 , 13 , 13 , 14 , 14 , // 31 - 40
53 
54 14 , 14 , 14 , 14 , 14 , 15 , 15 , 15 , 16 , 16 , // 41 - 50
55 
56 // ----------------------------------------------------------
57 
58 16 , 16 , 16 , 17 , 18 , 18 , 19 , 19 , 19 , 19 , // 51 - 60
59 
60 19 , 19 , 19 , 20 , 19 , 19 , 19 , 19 , 19 , 20 , // 61 - 70
61 
62 21 , 21 , 21 , 21 , 21 , 21 , 21 , 21 , 22 , 22 , // 71 - 80
63 
64 23 , 23 , 23 , 23 , 24 , 24 , 25 , 25 , 26 , 26 , // 81 - 90
65 
66 27 , 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 
75 const G4int
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 
92 const G4double
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 
535 const G4int
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 
743 G4int
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 
754 G4double
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 
768 G4int
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 
782 G4int
784 {
785 #ifdef G4VERBOSE
786  if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetNumberOfFreeElectrons"); }
787 #endif
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 
800 {
801 #ifdef G4VERBOSE
802  if(Z<0 || Z>100) { Z = PrintErrorZ(Z, "GetTotalBindingEnergy"); }
803 #endif
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 
815 G4int G4AtomicShells::PrintErrorZ(G4int Z, const G4String& ss)
816 {
817  G4String sss = "G4AtomicShells::"+ss+"()";
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 
826 G4int
827 G4AtomicShells::PrintErrorShell(G4int Z, G4int n, const G4String& ss)
828 {
829  G4String sss = "G4AtomicShells::"+ss+"()";
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....