Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4AtomicShells_XDB_EADL.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27////////////////////////////////////////////////////////////////////////////////
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#include "G4SystemOfUnits.hh"
43
44
45const G4int
46G4AtomicShells_XDB_EADL::fNumberOfShells[121] =
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
76const G4int
77G4AtomicShells_XDB_EADL::fIndexOfShells[121] =
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
107const G4double
108G4AtomicShells_XDB_EADL::fBindingEnergies[2171] =
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
715const G4int
716G4AtomicShells_XDB_EADL::fNumberOfElectrons[2171] =
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
1232G4int
1233G4AtomicShells_XDB_EADL::GetNumberOfShells(G4int Z)
1234{
1235#ifdef G4VERBOSE
1236 if (Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetNumberOfShells"); }
1237#endif
1238 return fNumberOfShells[Z];
1239}
1240
1241
1242G4double
1243G4AtomicShells_XDB_EADL::GetBindingEnergy(G4int Z, G4int ShellNb)
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
1255G4int
1256G4AtomicShells_XDB_EADL::GetNumberOfElectrons(G4int Z, G4int ShellNb)
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
1269G4int
1270G4AtomicShells_XDB_EADL::GetNumberOfFreeElectrons(G4int Z, G4double th)
1271{
1272#ifdef G4VERBOSE
1273 if(Z<0 || Z>120) { Z = PrintErrorZ(Z, "GetNumberOfFreeElectrons"); }
1274#endif
1275 G4int idx = fIndexOfShells[Z];
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
1284G4double G4AtomicShells_XDB_EADL::GetTotalBindingEnergy(G4int Z)
1285{
1286#ifdef G4VERBOSE
1287 if ( Z < 0 || Z > 120) { Z = PrintErrorZ(Z, "GetTotalBindingEnergy"); }
1288#endif
1289 G4int idx = fIndexOfShells[Z];
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
1299G4int G4AtomicShells_XDB_EADL::PrintErrorZ(G4int Z, const G4String& ss)
1300{
1301 G4String sss = "G4AtomicShells::"+ss+"()";
1302 G4ExceptionDescription ed;
1303 ed << "Atomic number out of range Z= " << Z;
1304 G4Exception(sss,"mat060",FatalException,ed,"");
1305 return 1;
1306}
1307
1308
1309G4int
1310G4AtomicShells_XDB_EADL::PrintErrorShell(G4int Z, G4int n, const G4String& ss)
1311{
1312 G4String sss = "G4AtomicShells::"+ss+"()";
1313 G4ExceptionDescription ed;
1314 ed << "Shell number out of range Nshell= " << n << " Z= " << Z;
1315 G4Exception(sss,"mat061",FatalException,ed,"");
1316 return 0;
1317}
1318#endif
1319
G4AtomicShells_XDB_EADL.hh
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4AtomicShells_XDB_EADL::GetNumberOfShells
static G4int GetNumberOfShells(G4int Z)
Definition: G4AtomicShells_XDB_EADL.cc:1233
G4AtomicShells_XDB_EADL::GetNumberOfElectrons
static G4int GetNumberOfElectrons(G4int Z, G4int SubshellNb)
Definition: G4AtomicShells_XDB_EADL.cc:1256
G4AtomicShells_XDB_EADL::GetBindingEnergy
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
Definition: G4AtomicShells_XDB_EADL.cc:1243
G4AtomicShells_XDB_EADL::GetTotalBindingEnergy
static G4double GetTotalBindingEnergy(G4int Z)
Definition: G4AtomicShells_XDB_EADL.cc:1284
G4AtomicShells_XDB_EADL::GetNumberOfFreeElectrons
static G4int GetNumberOfFreeElectrons(G4int Z, G4double th)
Definition: G4AtomicShells_XDB_EADL.cc:1270