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

Last change on this file since 850 was 850, checked in by garnier, 16 years ago

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