Context Navigation

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

Last change on this file since 1323 was 1315, checked in by garnier, 15 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
46	const G4int
47	G4AtomicShells::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
57	11 , 11 , 11 , 11 , 11 , 12 , 13 , 13 , 14 , 14 , // 31 - 40
58
59	14 , 14 , 14 , 14 , 14 , 15 , 15 , 15 , 16 , 16 , // 41 - 50
60
61	// ----------------------------------------------------------
62
63	16 , 16 , 16 , 17 , 18 , 18 , 19 , 19 , 19 , 19 , // 51 - 60
64
65	19 , 19 , 19 , 20 , 19 , 19 , 19 , 19 , 19 , 20 , // 61 - 70
66
67	21 , 21 , 21 , 21 , 21 , 21 , 21 , 21 , 22 , 22 , // 71 - 80
68
69	23 , 23 , 23 , 23 , 24 , 24 , 25 , 25 , 26 , 26 , // 81 - 90
70
71	27 , 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
80	const G4int
81	G4AtomicShells::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
97	const G4double
98	G4AtomicShells::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
538	const G4int
539	G4AtomicShells::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
746	G4int
747	G4AtomicShells::GetNumberOfShells(G4int Z)
748	{
749	assert (Z>0 && Z<101);
750	return fNumberOfShells[Z];
751	}
752
753	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
754
755	G4double
756	G4AtomicShells::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
764	G4int
765	G4AtomicShells::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
773	G4double 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.

Download in other formats: