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