source: trunk/source/materials/src/G4AtomicShells.cc @ 1315

Last change on this file since 1315 was 1315, checked in by garnier, 14 years ago

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