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