ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4AtomicShells_XDB_EADL.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4AtomicShells_XDB_EADL.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 //
28 // //
29 // File: G4AtomicShells_XDB_EADL.cc //
30 // Author: D.H. Wright (SLAC) //
31 // Date: 9 August 2018 //
32 // //
33 // Description: Class containing number of shells, electron configurations //
34 // and binding energies for atoms from Z = 1 to Z = 120. //
35 // Most entries are taken from the X-ray Data Book, with //
36 // unmeasured values supplied by Kibedi. See documentation //
37 // file in G4EMLOW7.3 or later, directory fluor/ //
38 // //
40 
42 #include "G4SystemOfUnits.hh"
43 
44 
45 const G4int
47 {
48  1 , // nonexistent zero element
49 
50  1, 1, 2, 2, 3, 3, 4, 4, 4, 4, // 1 - 10
51 
52  5, 5, 6, 6, 7, 7, 7, 7, 8, 8, // 11 - 20
53 
54  9, 9, 9, 10, 10, 10, 10, 10, 10, 10, // 21 - 30
55 
56  11, 11, 12, 12, 12, 12, 13, 13, 14, 14, // 31 - 40
57 
58  14, 15, 15, 15, 15, 14, 15, 15, 16, 16, // 41 - 50
59 
60  17, 17, 17, 17, 18, 18, 19, 21, 21, 20, // 51 - 60
61 
62  19, 20, 20, 21, 21, 20, 20, 20, 20, 20, // 61 - 70
63 
64  21, 21, 21, 21, 22, 22, 22, 22, 22, 22, // 71 - 80
65 
66  23, 23, 24, 24, 24, 24, 25, 25, 26, 26, // 81 - 90
67 
68  27, 27, 27, 26, 27, 28, 28, 27, 27, 27, // 91 - 100
69 
70  27, 27, 28, 28, 28, 28, 29, 29, 29, 29, // 101 - 110
71 
72  29, 29, 30, 30, 31, 31, 31, 31, 32, 32
73 };
74 
75 
76 const G4int
78 {
79  0, // nonexistent zero element
80 
81  1, 2, 3, 5, 7, 10, 13, 17, 21, 25, // 1 - 10
82 
83  29, 34, 39, 45, 51, 58, 65, 72, 79, 87, // 11 - 20
84 
85  95, 104, 113, 122, 132, 142, 152, 162, 172, 182, // 21 - 30
86 
87  192, 203, 214, 226, 238, 250, 262, 275, 288, 302, // 31 - 40
88 
89  316, 330, 345, 360, 375, 390, 404, 419, 434, 450, // 41 - 50
90 
91  466, 483, 500, 517, 534, 552, 570, 589, 610, 631, // 51 - 60
92 
93  651, 670, 690, 710, 731, 752, 772, 792, 812, 832, // 61 - 70
94 
95  852, 873, 894, 915, 936, 958, 980, 1002, 1024, 1046, // 71 - 80
96 
97  1068, 1091, 1114, 1138, 1162, 1186, 1210, 1235, 1260, 1286, // 81 - 90
98 
99  1312, 1339, 1366, 1393, 1419, 1446, 1474, 1502, 1529, 1556, // 91 - 100
100 
101  1583, 1610, 1637, 1665, 1693, 1721, 1749, 1778, 1807, 1836, // 101 - 110
102 
103  1865, 1894, 1923, 1953, 1983, 2014, 2045, 2076, 2107, 2139
104 };
105 
106 
107 const G4double
109 {
110  // 0 ---------------------------------------------------------------------
111  0.0000, // nonexistent zero element
112 
113  // H ---------------------------------------------------------------------
114  0.0136,
115 
116  // He ---------------------------------------------------------------------
117  0.0246,
118 
119  // Li ---------------------------------------------------------------------
120  0.0547, 0.0054,
121 
122  // Be ---------------------------------------------------------------------
123  0.1115, 0.0093,
124 
125  // B ---------------------------------------------------------------------
126  0.1880, 0.0126, 0.0047,
127 
128  // C ---------------------------------------------------------------------
129  0.2842, 0.0181, 0.0090,
130 
131  // N ---------------------------------------------------------------------
132  0.4099, 0.0373, 0.0089, 0.0097,
133 
134  // O ---------------------------------------------------------------------
135  0.5431, 0.0416, 0.0068, 0.0074,
136 
137  // F ---------------------------------------------------------------------
138  0.6967, 0.0340, 0.0084, 0.0087,
139 
140  // Ne ---------------------------------------------------------------------
141  0.8702, 0.0485, 0.0217, 0.0216,
142 
143  // Na ---------------------------------------------------------------------
144  1.0708, 0.0635, 0.03065, 0.03081, 0.0007,
145 
146  // Mg ---------------------------------------------------------------------
147  1.3030, 0.0887, 0.04978, 0.0495, 0.0021,
148 
149  // Al ---------------------------------------------------------------------
150  1.5596, 0.1178, 0.07295, 0.07255, 0.0007, 0.0060,
151 
152  // Si ---------------------------------------------------------------------
153  1.8390, 0.1497, 0.09982, 0.09942, 0.0076, 0.0030,
154 
155  // P ---------------------------------------------------------------------
156  2.1455, 0.1890, 0.1360, 0.1350, 0.0162, 0.0096, 0.0101,
157 
158  // S ---------------------------------------------------------------------
159  2.4720, 0.2309, 0.1636, 0.1625, 0.0158, 0.0078, 0.0082,
160 
161  // Cl ---------------------------------------------------------------------
162  2.8224, 0.2700, 0.2020, 0.2000, 0.0175, 0.0067, 0.0067,
163 
164  // Ar ---------------------------------------------------------------------
165  3.2059, 0.3263, 0.2506, 0.2484, 0.0293, 0.0159, 0.0157,
166 
167  // K ---------------------------------------------------------------------
168  3.6084, 0.3786, 0.2973, 0.2946, 0.0348, 0.0183, 0.0183,
169  0.0001,
170 
171  // Ca ---------------------------------------------------------------------
172  4.0385, 0.4384, 0.3497, 0.3462, 0.0443, 0.0254, 0.0254,
173  0.0018,
174 
175  // Sc ---------------------------------------------------------------------
176  4.4920, 0.4980, 0.4036, 0.3987, 0.0511, 0.0283, 0.0283,
177  0.0066, 0.0017,
178 
179  // Ti ---------------------------------------------------------------------
180  4.9660, 0.5609, 0.4602, 0.4538, 0.0587, 0.0326, 0.0326,
181  0.0037, 0.0016,
182 
183  // V ---------------------------------------------------------------------
184  5.4650, 0.6267, 0.5198, 0.5121, 0.0663, 0.0372, 0.0372,
185  0.0022, 0.0017,
186 
187  // Cr ---------------------------------------------------------------------
188  5.9890, 0.6960, 0.5838, 0.5741, 0.0741, 0.0422, 0.0422,
189  0.0029, 0.0022, 0.0010,
190 
191  // Mn ---------------------------------------------------------------------
192  6.5390, 0.7691, 0.6499, 0.6387, 0.0823, 0.0472, 0.0472,
193  0.0035, 0.0027, 0.0019,
194 
195  // Fe ---------------------------------------------------------------------
196  7.1120, 0.8446, 0.7199, 0.7068, 0.0913, 0.0527, 0.0527,
197  0.0039, 0.0031, 0.0021,
198 
199  // Co ---------------------------------------------------------------------
200  7.7090, 0.9251, 0.7932, 0.7781, 0.1010, 0.0589, 0.0599,
201  0.0033, 0.0027, 0.0019,
202 
203  // Ni ---------------------------------------------------------------------
204  8.3330, 1.0086, 0.8700, 0.8527, 0.1108, 0.0680, 0.0662,
205  0.0039, 0.0033, 0.0022,
206 
207  // Cu ---------------------------------------------------------------------
208  8.9790, 1.0967, 0.9523, 0.9327, 0.1225, 0.0773, 0.0751,
209  0.0018, 0.0015, 0.0012,
210 
211  // Zn ---------------------------------------------------------------------
212  9.6590, 1.1962, 1.0449, 1.0218, 0.1398, 0.0914, 0.0886,
213  0.0102, 0.0101, 0.0013,
214 
215  // Ga ---------------------------------------------------------------------
216  10.3670, 1.2990, 1.1432, 1.1164, 0.1595, 0.1035, 0.1000,
217  0.0187, 0.0187, 0.0056, 0.0008,
218 
219  // Ge ---------------------------------------------------------------------
220  11.1030, 1.4146, 1.2481, 1.2170, 0.1801, 0.1249, 0.1208,
221  0.0298, 0.0292, 0.0090, 0.0023,
222 
223  // As ---------------------------------------------------------------------
224  11.8670, 1.5270, 1.3591, 1.3236, 0.2047, 0.1462, 0.1412,
225  0.0417, 0.0417, 0.0125, 0.0025, 0.0025,
226 
227  // Se ---------------------------------------------------------------------
228  12.6580, 1.6520, 1.4743, 1.4339, 0.2296, 0.1665, 0.1607,
229  0.0555, 0.0546, 0.0162, 0.0056, 0.0056,
230 
231  // Br ---------------------------------------------------------------------
232  13.4740, 1.7820, 1.5960, 1.5500, 0.2570, 0.1890, 0.1820,
233  0.0700, 0.0690, 0.0273, 0.0052, 0.0046,
234 
235  // Kr ---------------------------------------------------------------------
236  14.3260, 1.9210, 1.7309, 1.6784, 0.2928, 0.2222, 0.2144,
237  0.0950, 0.0938, 0.0275, 0.0141, 0.0141,
238 
239  // Rb ---------------------------------------------------------------------
240  15.2000, 2.0650, 1.8640, 1.8040, 0.3267, 0.2487, 0.2391,
241  0.1130, 0.1120, 0.0305, 0.0163, 0.0153, 0.0040,
242 
243  // Sr ---------------------------------------------------------------------
244  16.1050, 2.2160, 2.0070, 1.9400, 0.3587, 0.2803, 0.2700,
245  0.1360, 0.1342, 0.0389, 0.0213, 0.0201, 0.0050,
246 
247  // Y ---------------------------------------------------------------------
248  17.0380, 2.3730, 2.1560, 2.0800, 0.3920, 0.3106, 0.2988,
249  0.1577, 0.1558, 0.0438, 0.0244, 0.0231, 0.0024, 0.0060,
250 
251  // Zr ---------------------------------------------------------------------
252  17.9980, 2.5320, 2.3070, 2.2230, 0.4303, 0.3435, 0.3298,
253  0.1811, 0.1788, 0.0506, 0.0285, 0.0271, 0.0030, 0.0070,
254 
255  // Nb ---------------------------------------------------------------------
256  18.9860, 2.6980, 2.4650, 2.3710, 0.4666, 0.3761, 0.3606,
257  0.2050, 0.2023, 0.0564, 0.0326, 0.0308, 0.0032, 0.0070,
258 
259  // Mo ---------------------------------------------------------------------
260  20.0000, 2.8660, 2.6250, 2.5200, 0.5063, 0.4116, 0.3940,
261  0.2311, 0.2279, 0.0632, 0.0376, 0.0355, 0.0019, 0.0012,
262  0.0070,
263 
264  // Tc ---------------------------------------------------------------------
265  21.0440, 3.0430, 2.7930, 2.6770, 0.5440, 0.4476, 0.4177,
266  0.2576, 0.2539, 0.0695, 0.0423, 0.0399, 0.0020, 0.0012,
267  0.0070,
268 
269  // Ru ---------------------------------------------------------------------
270  22.1170, 3.2240, 2.9670, 2.8380, 0.5861, 0.4835, 0.4614,
271  0.2842, 0.2800, 0.0750, 0.0463, 0.0432, 0.0024, 0.0018,
272  0.0070,
273 
274  // Rh ---------------------------------------------------------------------
275  23.2200, 3.4120, 3.1460, 3.0040, 0.6281, 0.5213, 0.4965,
276  0.3119, 0.3072, 0.0814, 0.0505, 0.0473, 0.0028, 0.0022,
277  0.0070,
278 
279  // Pd ---------------------------------------------------------------------
280  24.3500, 3.6040, 3.3300, 3.1730, 0.6716, 0.5599, 0.5323,
281  0.3405, 0.3352, 0.0871, 0.0557, 0.0509, 0.0017, 0.0013,
282 
283  // Ag ---------------------------------------------------------------------
284  25.5140, 3.8060, 3.5240, 3.3510, 0.7190, 0.6038, 0.5730,
285  0.3740, 0.3683, 0.0970, 0.0637, 0.0583, 0.0036, 0.0031,
286  0.0080,
287 
288  // Cd ---------------------------------------------------------------------
289  26.7110, 4.0180, 3.7270, 3.5380, 0.7720, 0.6526, 0.6184,
290  0.4119, 0.4052, 0.1098, 0.0639, 0.0639, 0.0117, 0.0107,
291  0.0090,
292 
293  // In ---------------------------------------------------------------------
294  27.9400, 4.2380, 3.9380, 3.7300, 0.8272, 0.7032, 0.6653,
295  0.4514, 0.4439, 0.1229, 0.0735, 0.0735, 0.0177, 0.0169,
296  0.0001, 0.0008,
297 
298  // Sn ---------------------------------------------------------------------
299  29.2000, 4.4650, 4.1560, 3.9290, 0.8847, 0.7565, 0.7146,
300  0.4932, 0.4849, 0.1371, 0.0836, 0.0836, 0.0249, 0.0239,
301  0.0009, 0.0011,
302 
303  // Sb ---------------------------------------------------------------------
304  30.4910, 4.6980, 4.3800, 4.1320, 0.9460, 0.8127, 0.7664,
305  0.5375, 0.5282, 0.1532, 0.0956, 0.0956, 0.0333, 0.0321,
306  0.0067, 0.0022, 0.0020,
307 
308  // Te ---------------------------------------------------------------------
309  31.8140, 4.9390, 4.6120, 4.3410, 1.0060, 0.8708, 0.8200,
310  0.5834, 0.5730, 0.1694, 0.1033, 0.1033, 0.0419, 0.0404,
311  0.0116, 0.0026, 0.0020,
312 
313  // I ---------------------------------------------------------------------
314  33.1690, 5.1880, 4.8520, 4.5570, 1.0720, 0.9310, 0.8750,
315  0.6308, 0.6193, 0.1860, 0.1230, 0.1230, 0.0506, 0.0489,
316  0.0136, 0.0038, 0.0029,
317 
318  // Xe ---------------------------------------------------------------------
319  34.5610, 5.4530, 5.1070, 4.7860, 1.1487, 1.0021, 0.9406,
320  0.6890, 0.6764, 0.2132, 0.1467, 0.1455, 0.0695, 0.0675,
321  0.0233, 0.0134, 0.0121,
322 
323  // Cs ---------------------------------------------------------------------
324  35.9850, 5.7140, 5.3590, 5.0120, 1.2110, 1.0710, 1.0030,
325  0.7405, 0.7266, 0.2323, 0.1724, 0.1613, 0.0798, 0.0775,
326  0.0227, 0.0142, 0.0121, 0.004,
327 
328  // Ba ---------------------------------------------------------------------
329  37.4410, 5.9890, 5.6240, 5.2470, 1.2930, 1.1370, 1.0630,
330  0.7957, 0.7805, 0.2535, 0.1920, 0.1786, 0.0926, 0.0899,
331  0.0303, 0.0170, 0.0148, 0.0050,
332 
333  // La ---------------------------------------------------------------------
334  38.9250, 6.2660, 5.8910, 5.4830, 1.3620, 1.2090, 1.1280,
335  0.8530, 0.8360, 0.2747, 0.2058, 0.1960, 0.1053, 0.1025,
336  0.0343, 0.0193, 0.0168, 0.0060, 0.0050,
337 
338  // Ce ---------------------------------------------------------------------
339  40.4430, 6.5490, 6.1640, 5.7230, 1.4360, 1.2740, 1.1870,
340  0.9024, 0.8838, 0.2910, 0.2232, 0.2065, 0.1090, 0.1076,
341  0.0001, 0.0001, 0.0378, 0.0198, 0.0170, 0.0060, 0.0050,
342 
343  // Pr ---------------------------------------------------------------------
344  41.9910, 6.8350, 6.4400, 5.9640, 1.5110, 1.3370, 1.2429,
345  0.9483, 0.9288, 0.3045, 0.2363, 0.2176, 0.1151, 0.1151,
346  0.0020, 0.0020, 0.0374, 0.0223, 0.0223, 0.0060, 0.0050,
347 
348  // Nd ---------------------------------------------------------------------
349  43.5690, 7.1260, 6.7220, 6.2080, 1.5750, 1.4030, 1.2970,
350  1.0033, 0.9804, 0.3192, 0.2433, 0.2246, 0.1205, 0.1205,
351  0.0015, 0.0015, 0.0375, 0.0211, 0.0211, 0.0050,
352 
353  // Pm ---------------------------------------------------------------------
354  45.1840, 7.4280, 7.0130, 6.4590, 1.6486, 1.4710, 1.3570,
355  1.0520, 1.0270, 0.3314, 0.2420, 0.2420, 0.1200, 0.1200,
356  0.0035, 0.0360, 0.0245, 0.0201, 0.0050,
357 
358  // Sm ---------------------------------------------------------------------
359  46.8340, 7.7370, 7.3120, 6.7160, 1.7230, 1.5410, 1.4200,
360  1.1109, 1.0834, 0.3472, 0.2656, 0.2474, 0.1290, 0.1290,
361  0.0052, 0.0052, 0.0374, 0.0213, 0.0213, 0.0050,
362 
363  // Eu ---------------------------------------------------------------------
364  48.5190, 8.0520, 7.6170, 6.9770, 1.8000, 1.6140, 1.4810,
365  1.1586, 1.1275, 0.3600, 0.2840, 0.2570, 0.1330, 0.1277,
366  0.0000, 0.0000, 0.0320, 0.0220, 0.0220, 0.0060,
367 
368  // Gd ---------------------------------------------------------------------
369  50.2390, 8.3760, 7.9300, 7.2430, 1.8810, 1.6880, 1.5440,
370  1.2219, 1.1896, 0.3786, 0.2860, 0.2710, 0.1495, 0.1426,
371  0.0086, 0.0086, 0.0360, 0.0280, 0.0210, 0.0060, 0.0060,
372 
373  // Tb ---------------------------------------------------------------------
374  51.9960, 8.7080, 8.2520, 7.5140, 1.9680, 1.7680, 1.6110,
375  1.2769, 1.2411, 0.3960, 0.3224, 0.2841, 0.1505, 0.1505,
376  0.0077, 0.0024, 0.0456, 0.0287, 0.0226, 0.0060, 0.0060,
377 
378  // Dy ---------------------------------------------------------------------
379  53.7890, 9.0460, 8.5810, 7.7900, 2.0470, 1.8420, 1.6760,
380  1.3330, 1.2926, 0.4142, 0.3335, 0.2932, 0.1536, 0.1536,
381  0.0080, 0.0043, 0.0499, 0.0263, 0.0263, 0.0060,
382 
383  // Ho ---------------------------------------------------------------------
384  55.6180, 9.3940, 8.9180, 8.0710, 2.1280, 1.9230, 1.7410,
385  1.3920, 1.3510, 0.4324, 0.3435, 0.3082, 0.1600, 0.1600,
386  0.0086, 0.0052, 0.0493, 0.0308, 0.0241, 0.0060,
387 
388  // Er ---------------------------------------------------------------------
389  57.4860, 9.7510, 9.2640, 8.3580, 2.2070, 2.0060, 1.8120,
390  1.4530, 1.4090, 0.4498, 0.3662, 0.3202, 0.1676, 0.1676,
391  0.0053, 0.0047, 0.0506, 0.0314, 0.0247, 0.0060,
392 
393  // Tm ---------------------------------------------------------------------
394  59.3900, 10.1160, 9.6170, 8.6480, 2.3070, 2.0900, 1.8850,
395  1.5150, 1.4680, 0.4709, 0.3859, 0.3326, 0.1755, 0.1755,
396  0.0062, 0.0046, 0.0547, 0.0318, 0.0250, 0.0060,
397 
398  // Yb ---------------------------------------------------------------------
399  61.3320, 10.4860, 9.9780, 8.9440, 2.3980, 2.1730, 1.9500,
400  1.5760, 1.5280, 0.4805, 0.3887, 0.3397, 0.1912, 0.1824,
401  0.0025, 0.0013, 0.0520, 0.0303, 0.0241, 0.0060,
402 
403  // Lu ---------------------------------------------------------------------
404  63.3140, 10.8700, 10.3490, 9.2440, 2.4910, 2.2640, 2.0240,
405  1.6390, 1.5890, 0.5068, 0.4124, 0.3592, 0.2061, 0.1963,
406  0.0089, 0.0075, 0.0573, 0.0336, 0.0267, 0.0046, 0.0070,
407 
408  // Hf ---------------------------------------------------------------------
409  65.3510, 11.2710, 10.7390, 9.5610, 2.6010, 2.3650, 2.1080,
410  1.7160, 1.6620, 0.5380, 0.4382, 0.3807, 0.2200, 0.2115,
411  0.0159, 0.0142, 0.0642, 0.0380, 0.0299, 0.0066, 0.0070,
412 
413  // Ta ---------------------------------------------------------------------
414  67.4160, 11.6820, 11.1360, 9.8810, 2.7080, 2.4690, 2.1940,
415  1.7930, 1.7350, 0.5634, 0.4634, 0.4009, 0.2379, 0.2264,
416  0.0235, 0.0216, 0.0697, 0.0422, 0.0327, 0.0057, 0.0080,
417 
418  // W ---------------------------------------------------------------------
419  69.5250, 12.1000, 11.5440, 10.2070, 2.8200, 2.5750, 2.2810,
420  1.8720, 1.8090, 0.5941, 0.4904, 0.4236, 0.2559, 0.2435,
421  0.0336, 0.0314, 0.0756, 0.0453, 0.0368, 0.0061, 0.0080,
422 
423  // Re ---------------------------------------------------------------------
424  71.6760, 12.5270, 11.9590, 10.5350, 2.9320, 2.6820, 2.3670,
425  1.9490, 1.8830, 0.6254, 0.5187, 0.4468, 0.2739, 0.2605,
426  0.0429, 0.0405, 0.0830, 0.0456, 0.0346, 0.0038, 0.0025,
427  0.0080,
428 
429  // Os ---------------------------------------------------------------------
430  73.8710, 12.9680, 12.3850, 10.8710, 3.0490, 2.7920, 2.4570,
431  2.0310, 1.9600, 0.6582, 0.5491, 0.4707, 0.2931, 0.2785,
432  0.0534, 0.0507, 0.0840, 0.0580, 0.0445, 0.0009, 0.0004,
433  0.0080,
434 
435  // Ir ---------------------------------------------------------------------
436  76.1110, 13.4190, 12.8240, 11.2150, 3.1740, 2.9090, 2.5510,
437  2.1160, 2.0400, 0.6911, 0.5778, 0.4958, 0.3119, 0.2963,
438  0.0638, 0.0608, 0.0952, 0.0630, 0.0480, 0.0042, 0.0032,
439  0.0090,
440 
441  // Pt ---------------------------------------------------------------------
442  78.3950, 13.8800, 13.2730, 11.5640, 3.2960, 3.0270, 2.6450,
443  2.2020, 2.1220, 0.7254, 0.6091, 0.5194, 0.3316, 0.3146,
444  0.0745, 0.0712, 0.1017, 0.0653, 0.0517, 0.0028, 0.0014,
445  0.0090,
446 
447  // Au ---------------------------------------------------------------------
448  80.7250, 14.3530, 13.7340, 11.9190, 3.4250, 3.1480, 2.7430,
449  2.2910, 2.2060, 0.7621, 0.6427, 0.5463, 0.3532, 0.3351,
450  0.0876, 0.0840, 0.1072, 0.0742, 0.0572, 0.0033, 0.0018,
451  0.0090,
452 
453  // Hg ---------------------------------------------------------------------
454  83.1020, 14.8390, 14.2090, 12.2840, 3.5620, 3.2790, 2.8470,
455  2.3850, 2.2950, 0.8022, 0.6802, 0.5766, 0.3782, 0.3588,
456  0.1040, 0.0999, 0.1270, 0.0831, 0.0645, 0.0096, 0.0078,
457  0.0090,
458 
459  // Tl --------------------------------------------------------------------
460  85.5300, 15.3470, 14.6980, 12.6580, 3.7040, 3.4160, 2.9570,
461  2.4850, 2.3890, 0.8462, 0.7205, 0.6095, 0.4057, 0.3850,
462  0.1222, 0.1178, 0.1360, 0.0946, 0.0735, 0.0147, 0.0125,
463  0.0060, 0.0008,
464 
465  // Pb --------------------------------------------------------------------
466  88.0050, 15.8610, 15.200, 13.0350, 3.8510, 3.5540, 3.0660,
467  2.5860, 2.4840, 0.8918, 0.7619, 0.6435, 0.4343, 0.4122,
468  0.1417, 0.1369, 0.1470, 0.1064, 0.0833, 0.0207, 0.0181,
469  0.0031, 0.0007,
470 
471  // Bi --------------------------------------------------------------------
472  90.5240, 16.3880, 15.7110, 13.4190, 3.9990, 3.6960, 3.1770,
473  2.6880, 2.5800, 0.9390, 0.8052, 0.6788, 0.4640, 0.4401,
474  0.1623, 0.1570, 0.1593, 0.1190, 0.0926, 0.0269, 0.0238,
475  0.0075, 0.0012, 0.0002,
476 
477  // Po --------------------------------------------------------------------
478  93.1050, 16.9390, 16.2440, 13.8140, 4.1490, 3.8540, 3.3020,
479  2.7980, 2.6830, 0.9950, 0.8510, 0.7050, 0.5000, 0.4730,
480  0.1840, 0.1840, 0.1770, 0.1320, 0.1040, 0.0310, 0.0310,
481  0.0110, 0.0032, 0.0014,
482 
483  // At --------------------------------------------------------------------
484  95.7300, 17.4930, 16.7850, 14.2140, 4.3170, 4.0080, 3.4260,
485  2.9090, 2.7870, 1.0420, 0.8860, 0.7400, 0.5330, 0.5070,
486  0.2100, 0.2100, 0.1950, 0.1480, 0.1150, 0.0400, 0.0400,
487  0.0150, 0.0057, 0.0028,
488 
489  // Rn --------------------------------------------------------------------
490  98.4040, 18.0490, 17.3370, 14.6190, 4.4820, 4.1590, 3.5380,
491  3.0220, 2.8920, 1.0970, 0.9290, 0.7680, 0.5670, 0.5410,
492  0.2380, 0.2380, 0.2140, 0.1640, 0.1270, 0.0480, 0.0480,
493  0.0260, 0.0076, 0.0041,
494 
495  // Fr --------------------------------------------------------------------
496  101.1370, 18.6390, 17.9070, 15.0310, 4.6520, 4.3270, 3.6630,
497  3.1360, 3.0000, 1.1530, 0.9800, 0.8100, 0.6030, 0.5770,
498  0.2680, 0.2680, 0.2340, 0.1820, 0.1400, 0.0580, 0.0580,
499  0.0340, 0.0150, 0.0150, 0.0040,
500 
501  // Ra --------------------------------------------------------------------
502  103.9220, 19.2370, 18.4840, 15.4440, 4.8220, 4.4900, 3.7920,
503  3.2480, 3.1050, 1.2080, 1.0580, 0.8790, 0.6360, 0.6030,
504  0.2990, 0.2990, 0.2540, 0.2000, 0.1530, 0.0680, 0.0680,
505  0.0440, 0.0190, 0.0190, 0.0050,
506 
507  // Ac --------------------------------------------------------------------
508  106.7550, 19.8400, 19.0830, 15.8710, 5.0020, 4.6560, 3.9090,
509  3.3700, 3.2190, 1.2690, 1.0800, 0.8900, 0.6750, 0.6390,
510  0.3190, 0.3190, 0.2720, 0.2150, 0.1670, 0.0800, 0.0800,
511  0.0398, 0.0241, 0.0170, 0.0060, 0.0060,
512 
513  // Th --------------------------------------------------------------------
514  109.6510, 20.4720, 19.6930, 16.3000, 5.1820, 4.8300, 4.0460,
515  3.4910, 3.3320, 1.3300, 1.1680, 0.9664, 0.7121, 0.6752,
516  0.3424, 0.3331, 0.2900, 0.2290, 0.1820, 0.0925, 0.0854,
517  0.0414, 0.0245, 0.0166, 0.0060, 0.0060,
518 
519  // Pa --------------------------------------------------------------------
520  112.6010, 21.1050, 20.3140, 16.7330, 5.3670, 5.0010, 4.1740,
521  3.6110, 3.4420, 1.3870, 1.2240, 1.0070, 0.7430, 0.7080,
522  0.3710, 0.3600, 0.3100, 0.2320, 0.2320, 0.0940, 0.0940,
523  0.0060, 0.0467, 0.0281, 0.0189, 0.0060, 0.0060,
524 
525  // U --------------------------------------------------------------------
526  115.6060, 21.7570, 20.9480, 17.1660, 5.5480, 5.1820, 4.3030,
527  3.7280, 3.5520, 1.4390, 1.2710, 1.0430, 0.7783, 0.7362,
528  0.3882, 0.3774, 0.3210, 0.2570, 0.1920, 0.1028, 0.0942,
529  0.0060, 0.0439, 0.0268, 0.0168, 0.0060, 0.0060,
530 
531  // Np --------------------------------------------------------------------
532  118.6690, 22.4281, 21.6013, 17.6106, 5.7396, 5.3667, 4.4334,
533  3.8498, 3.6652, 1.5001, 1.3277, 1.0860, 0.8161, 0.7708,
534  0.4143, 0.4034, 0.3398, 0.2834, 0.2062, 0.1093, 0.1012,
535  0.0060, 0.0500, 0.0293, 0.0175, 0.0060, 0.0060,
536 
537  // Pu --------------------------------------------------------------------
538  121.7910, 23.0978, 22.2668, 18.0570, 5.9329, 5.5412, 4.5566,
539  3.9726, 3.7781, 1.5593, 1.3774, 1.1209, 0.8489, 0.8015,
540  0.4374, 0.4252, 0.3519, 0.2825, 0.2153, 0.1160, 0.1052,
541  0.0060, 0.0486, 0.0306, 0.0184, 0.0060,
542 
543  // Am --------------------------------------------------------------------
544  124.9810, 23.8051, 22.9494, 18.5062, 6.1326, 5.7470, 4.7060,
545  4.0921, 3.8869, 1.6192, 1.4351, 1.1680, 0.8804, 0.8300,
546  0.4633, 0.4490, 0.3730, 0.3030, 0.2164, 0.1180, 0.1079,
547  0.0060, 0.0060, 0.0504, 0.0311, 0.0181, 0.0060,
548 
549  // Cm --------------------------------------------------------------------
550  128.2000, 24.5200, 23.6600, 18.9700, 6.3350, 5.9400, 4.8500,
551  4.2300, 4.0150, 1.6850, 1.4950, 1.2050, 0.9250, 0.8700,
552  0.4950, 0.4850, 0.3970, 0.2350, 0.1350, 0.1200, 0.0920,
553  0.0060, 0.0060, 0.0670, 0.0500, 0.0060, 0.0060, 0.0970,
554 
555  // Bk --------------------------------------------------------------------
556  131.6000, 25.2600, 24.3800, 19.4400, 6.5450, 6.1450, 4.9900,
557  4.3600, 4.1350, 1.7450, 1.5550, 1.2450, 0.9600, 0.9000,
558  0.5250, 0.5100, 0.4150, 0.3300, 0.2430, 0.1400, 0.1250,
559  0.0060, 0.0060, 0.0960, 0.0700, 0.0500, 0.0044, 0.0060,
560 
561  // Cf --------------------------------------------------------------------
562  135.0000, 26.0100, 25.1200, 19.9100, 6.7550, 6.3500, 5.1350,
563  4.4950, 4.2600, 1.8100, 1.6150, 1.2850, 0.9950, 0.9300,
564  0.5550, 0.5400, 0.4300, 0.3450, 0.2500, 0.1450, 0.1300,
565  0.0060, 0.0060, 0.1000, 0.0700, 0.0500, 0.006,
566 
567  // Es --------------------------------------------------------------------
568  138.4000, 26.7800, 25.8800, 20.3900, 6.9750, 6.5600, 5.2800,
569  4.6300, 4.3850, 1.8750, 1.6750, 1.3250, 1.0300, 0.9650,
570  0.5850, 0.5700, 0.4450, 0.3600, 0.2600, 0.1550, 0.1350,
571  0.0060, 0.0060, 0.1000, 0.0700, 0.0500, 0.0060,
572 
573  // Fm --------------------------------------------------------------------
574  142.0000, 27.5700, 26.6600, 20.8700, 7.2000, 6.7800, 5.4250,
575  4.7650, 4.5100, 1.9400, 1.7400, 1.3700, 1.0650, 1.0000,
576  0.6150, 0.6000, 0.4600, 0.3750, 0.2700, 0.1650, 0.1400,
577  0.0060, 0.0060, 0.1000, 0.0800, 0.0500, 0.0060,
578 
579  // Md --------------------------------------------------------------------
580  145.6000, 28.3800, 27.4600, 21.3600, 7.4300, 7.0000, 5.5750,
581  4.9050, 4.6330, 2.0050, 1.8050, 1.4150, 1.1050, 1.0350,
582  0.6450, 0.6300, 0.4800, 0.3900, 0.2800, 0.1750, 0.1450,
583  0.0070, 0.0070, 0.1100, 0.0800, 0.0500, 0.0060,
584 
585  // No --------------------------------------------------------------------
586  149.3000, 29.2220, 28.2800, 21.8500, 7.6700, 7.2300, 5.7300,
587  5.0500, 4.7650, 2.0750, 1.8750, 1.4600, 1.1450, 1.0700,
588  0.6800, 0.6500, 0.5000, 0.4100, 0.2900, 0.1850, 0.1550,
589  0.0100, 0.0100, 0.1100, 0.0800, 0.0500, 0.0060,
590 
591  // Lr --------------------------------------------------------------------
592  153.0000, 30.0800, 29.1300, 22.3500, 7.9200, 7.4700, 5.8900,
593  5.2000, 4.9000, 2.1500, 1.9500, 1.5050, 1.1900, 1.1100,
594  0.7200, 0.6900, 0.5200, 0.4300, 0.3000, 0.2000, 0.1600,
595  0.0150, 0.0150, 0.1200, 0.0900, 0.0600, 0.0070, 0.0070,
596 
597  // Rf --------------------------------------------------------------------
598  156.9000, 30.9600, 30.0000, 22.8600, 8.1700, 7.7200, 6.0550,
599  5.3500, 5.0400, 2.2250, 2.0250, 1.5550, 1.2350, 1.1500,
600  0.7600, 0.7200, 0.5400, 0.4550, 0.3150, 0.2150, 0.1700,
601  0.0200, 0.0200, 0.1200, 0.0900, 0.0600, 0.0080, 0.0080,
602 
603  // Db --------------------------------------------------------------------
604  160.8000, 31.8700, 30.9000, 23.3800, 8.4300, 7.9800, 6.6250,
605  5.5050, 5.1850, 2.3050, 2.1050, 1.6050, 1.2800, 1.1900,
606  0.8000, 0.7500, 0.5600, 0.4800, 0.3300, 0.2300, 0.1800,
607  0.0250, 0.0250, 0.1300, 0.1000, 0.0600, 0.0090, 0.0080,
608 
609  // Sg --------------------------------------------------------------------
610  164.9000, 32.8000, 31.8200, 23.9100, 8.7000, 8.2400, 6.4000,
611  5.6650, 5.3300, 2.3850, 2.1900, 1.6600, 1.3250, 1.2350,
612  0.8400, 0.7800, 0.5850, 0.5050, 0.3450, 0.2450, 0.1900,
613  0.0300, 0.0300, 0.1400, 0.1100, 0.0600, 0.0090, 0.0080,
614 
615  // Bh --------------------------------------------------------------------
616  169.0000, 33.7600, 32.7700, 24.4400, 8.9800, 8.5100, 6.5800,
617  5.8300, 5.4750, 2.4700, 2.2750, 1.7150, 1.3700, 1.2800,
618  0.8800, 0.8200, 0.6100, 0.5300, 0.3600, 0.2600, 0.2000,
619  0.0400, 0.0400, 0.1400, 0.1100, 0.0700, 0.0100, 0.0100,
620  0.0090,
621 
622  // Hs --------------------------------------------------------------------
623  173.2000, 34.7500, 33.7400, 24.9200, 9.2700, 8.7900, 6.7600,
624  5.9950, 5.6200, 2.5550, 2.3650, 1.7700, 1.4200, 1.3250,
625  0.9200, 0.8600, 0.6350, 0.5550, 0.3750, 0.2750, 0.2100,
626  0.0500, 0.0500, 0.1500, 0.1200, 0.0700, 0.0100, 0.0100,
627  0.0090,
628 
629  // Mt --------------------------------------------------------------------
630  177.6000, 35.7700, 34.7400, 25.5200, 9.5700, 9.0800, 6.9400,
631  6.1600, 5.7700, 2.6450, 2.4600, 1.8250, 1.4700, 1.3700,
632  0.9600, 0.9000, 0.6600, 0.5800, 0.3900, 0.2900, 0.2200,
633  0.0600, 0.0600, 0.1600, 0.1200, 0.0700, 0.0100, 0.0100,
634  0.0100,
635 
636  // Ds --------------------------------------------------------------------
637  182.0000, 36.8100, 35.7700, 26.0700, 9.8700, 9.3700, 7.1200,
638  6.3300, 5.9200, 2.7350, 2.5600, 1.8800, 1.5200, 1.4150,
639  1.0000, 0.9400, 0.6900, 0.6100, 0.4050, 0.3050, 0.2300,
640  0.0700, 0.0700, 0.1600, 0.1300, 0.0800, 0.0100, 0.0100,
641  0.0100,
642 
643  // Rg --------------------------------------------------------------------
644  186.1360, 37.8810, 36.7990, 26.6480, 10.1870, 9.6870, 7.2240,
645  6.4470, 6.0250, 2.9260, 2.6730, 1.9780, 1.5970, 1.4830,
646  0.9604, 0.9267, 0.7283, 0.6450, 0.2810, 0.2510, 0.2445,
647  0.0671, 0.0614, 0.1267, 0.0901, 0.0498, 0.0094, 0.0066,
648  0.0094,
649 
650  // Cp --------------------------------------------------------------------
651  190.7170, 38.9840, 37.8960, 27.2130, 10.5080, 10.0020, 7.4060,
652  6.6180, 6.1780, 3.0330, 2.7760, 2.0390, 1.6520, 1.5320,
653  1.0016, 0.9662, 0.7624, 0.6770, 0.2960, 0.2650, 0.2577,
654  0.0753, 0.0690, 0.1349, 0.0969, 0.0528, 0.0104, 0.0073,
655  0.0098,
656 
657  // Uut --------------------------------------------------------------------
658  195.4100, 40.1210, 39.0340, 27.7880, 10.8420, 10.3310, 7.5930,
659  6.7930, 6.3350, 3.1460, 2.8860, 2.1070, 1.7120, 1.5870,
660  1.0464, 1.0092, 0.8006, 0.7140, 0.3150, 0.2820, 0.2740,
661  0.0866, 0.0798, 0.1464, 0.1069, 0.0586, 0.0138, 0.0101,
662  0.0116, 0.0049,
663 
664  // Fl --------------------------------------------------------------------
665  200.2140, 41.2940, 40.2070, 28.3700, 11.1840, 10.6720, 7.7830,
666  6.9720, 6.4950, 3.2640, 2.9990, 2.1740, 1.7740, 1.6420,
667  1.0923, 1.0533, 0.8409, 0.7530, 0.3350, 0.3000, 0.2910,
668  0.0984, 0.0911, 0.1583, 0.1174, 0.0646, 0.0173, 0.0130,
669  0.0135, 0.0057,
670 
671  // Uup --------------------------------------------------------------------
672  205.1360, 42.5010, 41.4200, 28.9590, 11.5410, 11.0220, 7.9790,
673  7.1550, 6.6580, 3.3860, 3.1180, 2.2450, 1.8380, 1.7010,
674  1.1410, 1.1000, 0.8831, 0.7940, 0.3570, 0.3200, 0.3100,
675  0.1124, 0.1045, 0.1726, 0.1301, 0.0724, 0.0227, 0.0177,
676  0.0169, 0.0078, 0.0038,
677 
678  // Lv --------------------------------------------------------------------
679  210.1770, 43.7500, 42.6740, 29.5560, 11.9070, 11.3870, 8.1760,
680  7.3410, 6.8240, 3.5130, 3.2410, 2.3160, 1.9030, 1.7590,
681  1.1907, 1.1477, 0.9275, 0.8380, 0.3790, 0.3400, 0.3296,
682  0.1269, 0.1184, 0.1874, 0.1434, 0.0805, 0.0282, 0.0226,
683  0.0203, 0.0101, 0.0044,
684 
685  // Uus --------------------------------------------------------------------
686  215.3370, 44.0310, 43.9750, 30.1580, 12.2830, 11.7620, 8.3760,
687  7.5300, 6.9910, 3.6430, 3.4690, 2.3900, 1.9700, 1.8190,
688  1.2413, 1.1962, 0.9733, 0.8820, 0.4020, 0.3610, 0.3496,
689  0.1417, 0.1327, 0.2026, 0.1571, 0.0886, 0.0337, 0.0275,
690  0.0237, 0.0124, 0.0052,
691 
692  // Uuo --------------------------------------------------------------------
693  220.6300, 46.3540, 45.3540, 30.7670, 12.6700, 12.1510, 8.5790,
694  7.7200, 7.1610, 3.7780, 3.5000, 2.4640, 2.0370, 1.8800,
695  1.2928, 1.2455, 1.0205, 0.9280, 0.4250, 0.3820, 0.3699,
696  0.1570, 0.1472, 0.2183, 0.1714, 0.0968, 0.0394, 0.0324,
697  0.0272, 0.0149, 0.0060,
698 
699  // Uue --------------------------------------------------------------------
700  226.0630, 47.7240, 46.7140, 31.3860, 13.0740, 12.5560, 8.7890,
701  7.9180, 7.3360, 3.9200, 3.6390, 2.5430, 2.1090, 1.9440,
702  1.3476, 1.2982, 1.0720, 0.9800, 0.4520, 0.4060, 0.3931,
703  0.1751, 0.1646, 0.2371, 0.1887, 0.1075, 0.0476, 0.0399,
704  0.0331, 0.0197, 0.0084, 0.0036,
705 
706  // Ubn --------------------------------------------------------------------
707  231.6340, 49.1390, 48.1580, 32.0130, 13.4890, 12.9750, 9.0000,
708  8.1170, 7.5130, 4.0660, 3.7840, 2.6220, 2.1820, 2.0100,
709  1.4040, 1.3523, 1.1256, 1.0330, 0.4790, 0.4310, 0.4173,
710  0.1941, 0.1830, 0.2570, 0.2071, 0.1190, 0.0565, 0.0479,
711  0.0396, 0.0251, 0.0113, 0.0042
712 };
713 
714 
715 const G4int
717 {
718  0, // nonexistent zero element
719 
720  // H ---------------------------------------------------------
721  1,
722 
723  // He ---------------------------------------------------------
724  2,
725 
726  // Li --------------------------------------------------------
727  2, 1,
728 
729  // Be --------------------------------------------------------
730  2, 2,
731 
732  // B --------------------------------------------------------
733  2, 2, 1,
734 
735  // C ---------------------------------------------------------
736  2, 2, 2,
737 
738  // N ---------------------------------------------------------
739  2, 2, 2, 1,
740 
741  // O ---------------------------------------------------------
742  2, 2, 2, 2,
743 
744  // F ---------------------------------------------------------
745  2, 2, 2, 3,
746 
747  // Ne ---------------------------------------------------------
748  2, 2, 2, 4,
749 
750  // Na ----------------------------------------------
751  2, 2, 2, 4, 1,
752 
753  // Mg ----------------------------------------------
754  2, 2, 2, 4, 2,
755 
756  // Al ----------------------------------------------
757  2, 2, 2, 4, 2, 1,
758 
759  // Si ----------------------------------------------
760  2, 2, 2, 4, 2, 2,
761 
762  // P ----------------------------------------------
763  2, 2, 2, 4, 2, 2, 1,
764 
765  // S ----------------------------------------------
766  2, 2, 2, 4, 2, 2, 2,
767 
768  // Cl ----------------------------------------------
769  2, 2, 2, 4, 2, 2, 3,
770 
771  // Ar ----------------------------------------------
772  2, 2, 2, 4, 2, 2, 4,
773 
774  // K ----------------------------------------------
775  2, 2, 2, 4, 2, 2, 4, 1,
776 
777  // Ca ----------------------------------------------
778  2, 2, 2, 4, 2, 2, 4, 2,
779 
780  // Sc ----------------------------------------------
781  2, 2, 2, 4, 2, 2, 4, 1, 2,
782 
783  // Ti ----------------------------------------------
784  2, 2, 2, 4, 2, 2, 4, 2, 2,
785 
786  // V ----------------------------------------------
787  2, 2, 2, 4, 2, 2, 4, 3, 2,
788 
789  // Cr ----------------------------------------------
790  2, 2, 2, 4, 2, 2, 4, 4, 1, 1,
791 
792  // Mn ----------------------------------------------
793  2, 2, 2, 4, 2, 2, 4, 4, 1, 2,
794 
795  // Fe ----------------------------------------------
796  2, 2, 2, 4, 2, 2, 4, 4, 2, 2,
797 
798  // Co ----------------------------------------------
799  2, 2, 2, 4, 2, 2, 4, 4, 3, 2,
800 
801  // Ni ----------------------------------------------
802  2, 2, 2, 4, 2, 2, 4, 4, 4, 2,
803 
804  // Cu ----------------------------------------------
805  2, 2, 2, 4, 2, 2, 4, 4, 6, 1,
806 
807  // Zn ----------------------------------------------
808  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
809 
810  // Ga ----------------------------------------------
811  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
812  1,
813 
814  // Ge ----------------------------------------------
815  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
816  2,
817 
818  // As ----------------------------------------------
819  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
820  2, 1,
821 
822  // Se ----------------------------------------------
823  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
824  2, 2,
825 
826  // Br ----------------------------------------------
827  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
828  2, 3,
829 
830  // Kr ----------------------------------------------
831  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
832  2, 4,
833 
834  // Rb ----------------------------------------------
835  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
836  2, 4, 1,
837 
838  // Sr ----------------------------------------------
839  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
840  2, 4, 2,
841 
842  // Y ----------------------------------------------
843  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
844  2, 4, 1, 2,
845 
846  // Zr ----------------------------------------------
847  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
848  2, 4, 2, 2,
849 
850  // Nb ----------------------------------------------
851  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
852  2, 4, 4, 1,
853 
854  // Mo ----------------------------------------------
855  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
856  2, 4, 4, 1, 1,
857 
858  // Tc ----------------------------------------------
859  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
860  2, 4, 4, 1, 2,
861 
862  // Ru ----------------------------------------------
863  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
864  2, 4, 4, 3, 1,
865 
866  // Rh ----------------------------------------------
867  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
868  2, 4, 4, 4, 1,
869 
870  // Pd ----------------------------------------------
871  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
872  2, 4, 4, 6,
873 
874  // Ag ----------------------------------------------
875  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
876  2, 4, 4, 6, 1,
877 
878  // Cd ----------------------------------------------
879  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
880  2, 4, 4, 6, 2,
881 
882  // In ----------------------------------------------
883  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
884  2, 4, 4, 6, 2, 1,
885 
886  // Sn ----------------------------------------------
887  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
888  2, 4, 4, 6, 2, 2,
889 
890  // Sb ----------------------------------------------
891  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
892  2, 4, 4, 6, 2, 2, 1,
893 
894  // Te ----------------------------------------------
895  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
896  2, 4, 4, 6, 2, 2, 2,
897 
898  // I ----------------------------------------------
899  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
900  2, 4, 4, 6, 2, 2, 3,
901 
902  // Xe ----------------------------------------------
903  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
904  2, 4, 4, 6, 2, 2, 4,
905 
906  // Cs ----------------------------------------------
907  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
908  2, 4, 4, 6, 2, 2, 4, 1,
909 
910  // Ba ----------------------------------------------
911  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
912  2, 4, 4, 6, 2, 2, 4, 2,
913 
914  // La ----------------------------------------------
915  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
916  2, 4, 4, 6, 2, 2, 4, 1, 2,
917 
918  // Ce ----------------------------------------------
919  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
920  2, 4, 4, 6, 1, 1, 2, 2, 4, 1,
921  1,
922 
923  // Pr ----------------------------------------------
924  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
925  2, 4, 4, 6, 2, 1, 2, 2, 4, 1,
926  1,
927 
928  // Nd ----------------------------------------------
929  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
930  2, 4, 4, 6, 4, 1, 2, 2, 4, 1,
931 
932  // Pm ----------------------------------------------
933  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
934  2, 4, 4, 6, 5, 2, 2, 4, 2,
935 
936  // Sm ----------------------------------------------
937  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
938  2, 4, 4, 6, 6, 1, 2, 2, 4, 1,
939 
940  // Eu ----------------------------------------------
941  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
942  2, 4, 4, 6, 6, 1, 2, 2, 4, 2,
943 
944  // Gd ----------------------------------------------
945  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
946  2, 4, 4, 6, 6, 1, 2, 2, 4, 1,
947  2,
948 
949  // Tb ----------------------------------------------
950  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
951  2, 4, 4, 6, 6, 2, 2, 2, 4, 1,
952  2,
953 
954  // Dy ----------------------------------------------
955  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
956  2, 4, 4, 6, 6, 4, 2, 2, 4, 2,
957 
958  // Ho ----------------------------------------------
959  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
960  2, 4, 4, 6, 6, 5, 2, 2, 4, 2,
961 
962  // Er ----------------------------------------------
963  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
964  2, 4, 4, 6, 6, 6, 2, 2, 4, 2,
965 
966  // Tm ----------------------------------------------
967  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
968  2, 4, 4, 6, 6, 7, 2, 2, 4, 2,
969 
970  // Yb ----------------------------------------------
971  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
972  2, 4, 4, 6, 6, 8, 2, 2, 4, 2,
973 
974  // Lu ----------------------------------------------
975  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
976  2, 4, 4, 6, 6, 8, 2, 2, 4, 1,
977  2,
978 
979  // Hf ----------------------------------------------
980  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
981  2, 4, 4, 6, 6, 8, 2, 2, 4, 2,
982  2,
983 
984  // Ta ----------------------------------------------
985  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
986  2, 4, 4, 6, 6, 8, 2, 2, 4, 3,
987  2,
988 
989  // W ----------------------------------------------
990  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
991  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
992  2,
993 
994  // Re ----------------------------------------------
995  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
996  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
997  1, 2,
998 
999  // Os ----------------------------------------------
1000  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1001  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1002  2, 2,
1003 
1004  // Ir ----------------------------------------------
1005  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1006  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1007  3, 2,
1008 
1009  // Pt ----------------------------------------------
1010  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1011  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1012  5, 1,
1013 
1014  // Au ----------------------------------------------
1015  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1016  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1017  6, 1,
1018 
1019  // Hg ----------------------------------------------
1020  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1021  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1022  6, 2,
1023 
1024  // Tl ----------------------------------------------
1025  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1026  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1027  6, 2, 1,
1028 
1029  // Pb ----------------------------------------------
1030  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1031  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1032  6, 2, 2,
1033 
1034  // Bi ----------------------------------------------
1035  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1036  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1037  6, 2, 2, 1,
1038 
1039  // Po ----------------------------------------------
1040  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1041  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1042  6, 2, 2, 2,
1043 
1044  // At ----------------------------------------------
1045  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1046  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1047  6, 2, 2, 3,
1048 
1049  // Rn ----------------------------------------------
1050  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1051  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1052  6, 2, 2, 4,
1053 
1054  // Fr ----------------------------------------------
1055  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1056  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1057  6, 2, 2, 4, 1,
1058 
1059  // Ra ----------------------------------------------
1060  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1061  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1062  6, 2, 2, 4, 2,
1063 
1064  // Ac ----------------------------------------------
1065  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1066  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1067  6, 2, 2, 4, 1, 2,
1068 
1069  // Th ----------------------------------------------
1070  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1071  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1072  6, 2, 2, 4, 2, 2,
1073 
1074  // Pa ----------------------------------------------
1075  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1076  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1077  6, 2, 2, 2, 4, 1, 2,
1078 
1079  // U ----------------------------------------------
1080  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1081  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1082  6, 3, 2, 2, 4, 1, 2,
1083 
1084  // Np ----------------------------------------------
1085  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1086  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1087  6, 4, 2, 2, 4, 1, 2,
1088 
1089  // Pu ----------------------------------------------
1090  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1091  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1092  6, 6, 2, 2, 4, 2,
1093 
1094  // Am ----------------------------------------------
1095  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1096  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1097  6, 6, 1, 2, 2, 4, 2,
1098 
1099  // Cm ----------------------------------------------
1100  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1101  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1102  6, 6, 1, 2, 2, 4, 1, 2,
1103 
1104  // Bk ----------------------------------------------
1105  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1106  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1107  6, 6, 2, 2, 2, 4, 1, 2,
1108 
1109  // Cf ----------------------------------------------
1110  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1111  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1112  6, 6, 4, 2, 2, 4, 2,
1113 
1114  // Es ----------------------------------------------
1115  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1116  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1117  6, 6, 5, 2, 2, 4, 2,
1118 
1119  // Fm ----------------------------------------------
1120  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1121  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1122  6, 6, 6, 2, 2, 4, 2,
1123 
1124  // Md ----------------------------------------------
1125  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1126  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1127  6, 6, 7, 2, 2, 4, 2,
1128 
1129  // No ----------------------------------------------
1130  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1131  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1132  6, 6, 8, 2, 2, 4, 2,
1133 
1134  // Lr ----------------------------------------------
1135  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1136  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1137  6, 6, 8, 2, 2, 4, 1, 2,
1138 
1139  // Rf ----------------------------------------------
1140  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1141  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1142  6, 6, 8, 2, 2, 4, 2, 2,
1143 
1144  // Db ----------------------------------------------
1145  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1146  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1147  6, 6, 8, 2, 2, 4, 3, 2,
1148 
1149  // Sg ----------------------------------------------
1150  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1151  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1152  6, 6, 8, 2, 2, 4, 4, 2,
1153 
1154  // Bh ----------------------------------------------
1155  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1156  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1157  6, 6, 8, 2, 2, 4, 4, 1, 2,
1158 
1159  // Hs ----------------------------------------------
1160  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1161  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1162  6, 6, 8, 2, 2, 4, 4, 2, 2,
1163 
1164  // Mt ----------------------------------------------
1165  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1166  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1167  6, 6, 8, 2, 2, 4, 4, 3, 2,
1168 
1169  // Ds ----------------------------------------------
1170  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1171  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1172  6, 6, 8, 2, 2, 4, 4, 5, 1,
1173 
1174  // Rg ----------------------------------------------
1175  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1176  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1177  6, 6, 8, 2, 2, 4, 4, 6, 1,
1178 
1179  // Cp ----------------------------------------------
1180  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1181  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1182  6, 6, 8, 2, 2, 4, 4, 6, 2,
1183 
1184  // Uut ---------------------------------------------
1185  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1186  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1187  6, 6, 8, 2, 2, 4, 4, 6, 2, 1,
1188 
1189  // Fl ----------------------------------------------
1190  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1191  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1192  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1193 
1194  // Uup ---------------------------------------------
1195  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1196  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1197  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1198  1,
1199 
1200  // Lv ----------------------------------------------
1201  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1202  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1203  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1204  2,
1205 
1206  // Uus ---------------------------------------------
1207  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1208  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1209  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1210  3,
1211 
1212  // Uuo ---------------------------------------------
1213  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1214  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1215  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1216  4,
1217 
1218  // Uue ---------------------------------------------
1219  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1220  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1221  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1222  4, 1,
1223 
1224  // Ubn ---------------------------------------------
1225  2, 2, 2, 4, 2, 2, 4, 4, 6, 2,
1226  2, 4, 4, 6, 6, 8, 2, 2, 4, 4,
1227  6, 6, 8, 2, 2, 4, 4, 6, 2, 2,
1228  4, 2
1229 };
1230 
1231 
1232 G4int
1234 {
1235 #ifdef G4VERBOSE
1236  if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetNumberOfShells"); }
1237 #endif
1238  return fNumberOfShells[Z];
1239 }
1240 
1241 
1242 G4double
1244 {
1245 #ifdef G4VERBOSE
1246  if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetBindingEnergy"); }
1247  if(ShellNb < 0 || ShellNb >= fNumberOfShells[Z]) {
1248  ShellNb = PrintErrorShell(Z, ShellNb, "GetBindingEnergy");
1249  }
1250 #endif
1251  return fBindingEnergies[fIndexOfShells[Z] + ShellNb]*CLHEP::keV;
1252 }
1253 
1254 
1255 G4int
1257 {
1258 #ifdef G4VERBOSE
1259  if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetNumberOfElectrons"); }
1260  if(ShellNb < 0 || ShellNb >= fNumberOfShells[Z]) {
1261  ShellNb = PrintErrorShell(Z, ShellNb, "GetNumberOfElectrons");
1262  }
1263 #endif
1264  return fNumberOfElectrons[fIndexOfShells[Z] + ShellNb];
1265 }
1266 
1267 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
1268 
1269 G4int
1271 {
1272 #ifdef G4VERBOSE
1273  if(Z<0 || Z>120) { Z = PrintErrorZ(Z, "GetNumberOfFreeElectrons"); }
1274 #endif
1276  G4int idxmax = idx + fNumberOfShells[Z];
1277  G4int n = 0;
1278  for (G4int i=idx; i<idxmax; ++i) {
1279  if(fBindingEnergies[i] <= th) { n += fNumberOfElectrons[i]; }
1280  }
1281  return n;
1282 }
1283 
1285 {
1286 #ifdef G4VERBOSE
1287  if ( Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetTotalBindingEnergy"); }
1288 #endif
1290  G4int idxmax = idx + fNumberOfShells[Z];
1291  G4double energy = 0.0;
1292  for (G4int i=idx; i<idxmax; ++i) {energy += fBindingEnergies[i];}
1293  return energy*CLHEP::keV;
1294 }
1295 
1296 
1297 #ifdef G4VERBOSE
1298 
1299 G4int G4AtomicShells_XDB_EADL::PrintErrorZ(G4int Z, const G4String& ss)
1300 {
1301  G4String sss = "G4AtomicShells::"+ss+"()";
1303  ed << "Atomic number out of range Z= " << Z;
1304  G4Exception(sss,"mat060",FatalException,ed,"");
1305  return 1;
1306 }
1307 
1308 
1309 G4int
1310 G4AtomicShells_XDB_EADL::PrintErrorShell(G4int Z, G4int n, const G4String& ss)
1311 {
1312  G4String sss = "G4AtomicShells::"+ss+"()";
1314  ed << "Shell number out of range Nshell= " << n << " Z= " << Z;
1315  G4Exception(sss,"mat061",FatalException,ed,"");
1316  return 0;
1317 }
1318 #endif
1319