| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25
26 // -------------------------------------------
27 //
28 // GEANT4 Class file
29 //
30 // File name: G4NistElementBuilder
31 //
32 // Author: Vladimir Ivanchenko
33 //
34 // Creation date: 23.12.2004
35 //
36 // Modifications:
37 // 02.05.2006 Subtract mass of atomic electron
38 // 11.05.2006 Do not subtract mass of atomic e
39 // 17.10.2006 Add natiral abundances flag to e
40 // use G4 units for isotope mass ve
41 // 10.05.2007 Add protection agains Z>101 (VI)
42 // 26.07.2007 Create one and only one Nist ele
43 // allow users to create there own
44 //
45 // -------------------------------------------
46 //
47 // Class Description:
48 //
49 // Element data from the NIST DB on Atomic Wei
50 // http://physics.nist.gov/PhysRefData/Composi
51
52 #include "G4NistElementBuilder.hh"
53
54 #include "G4AutoLock.hh"
55 #include "G4PhysicalConstants.hh"
56 #include "G4SystemOfUnits.hh"
57
58 #include <sstream>
59
60 namespace
61 {
62 G4Mutex nistElementMutex = G4MUTEX_INITIALIZER
63 }
64
65 //....oooOO0OOooo........oooOO0OOooo........oo
66
67 G4NistElementBuilder::G4NistElementBuilder(G4i
68 {
69 atomicMass[0] = 0.0;
70 nFirstIsotope[0] = 0;
71 nIsotopes[0] = 0;
72 idxIsotopes[0] = 0;
73 elmSymbol.push_back(" ");
74 Initialise();
75 for (G4int i=0; i<maxNumElements; ++i) {
76 elmIndex[i] = -1;
77 }
78 }
79
80 //....oooOO0OOooo........oooOO0OOooo........oo
81
82 G4int G4NistElementBuilder::GetZ(const G4Strin
83 {
84 G4int Z = maxNumElements;
85 // Loop checking, 07-Aug-2015, Vladimir Ivan
86 do {
87 --Z;
88 } while (Z > 0 && elmSymbol[Z] != name);
89 return Z;
90 }
91
92 //....oooOO0OOooo........oooOO0OOooo........oo
93
94 G4double G4NistElementBuilder::GetAtomicMassAm
95 {
96 G4int Z = maxNumElements;
97 // Loop checking, 07-Aug-2015, Vladimir Ivan
98 do {
99 --Z;
100 } while (Z > 0 && elmSymbol[Z] != name);
101 return GetAtomicMassAmu(Z);
102 }
103
104 //....oooOO0OOooo........oooOO0OOooo........oo
105
106 G4Element* G4NistElementBuilder::FindOrBuildEl
107 {
108 G4Element* anElement = FindElement(Z);
109 if (anElement == nullptr && Z > 0 && Z < max
110 anElement = BuildElement(Z);
111 }
112 return anElement;
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oo
116
117 G4Element* G4NistElementBuilder::FindOrBuildEl
118 {
119 G4Element* elm = nullptr;
120 const G4ElementTable* theElementTable = G4El
121 std::size_t nelm = theElementTable->size();
122 for (std::size_t i = 0; i < nelm; ++i) {
123 if (symb == ((*theElementTable)[i])->GetSy
124 elm = (*theElementTable)[i];
125 break;
126 }
127 }
128 if (nullptr == elm) {
129 for (G4int Z = 1; Z < maxNumElements; ++Z)
130 if (symb == elmSymbol[Z]) {
131 elm = BuildElement(Z);
132 break;
133 }
134 }
135 }
136 return elm;
137 }
138
139 //....oooOO0OOooo........oooOO0OOooo........oo
140
141 G4Element* G4NistElementBuilder::BuildElement(
142 {
143 G4Element* theElement = nullptr;
144 if (Z < 1 || Z >= maxNumElements) {
145 return theElement;
146 }
147 const G4ElementTable* theTable = G4Element::
148 // Nist element was already built
149 if (0 <= elmIndex[Z]) {
150 return (*theTable)[elmIndex[Z]];
151 }
152
153 G4AutoLock l(&nistElementMutex);
154 if (0 <= elmIndex[Z]) {
155 theElement = (*theTable)[elmIndex[Z]];
156 }
157 else {
158 G4double Aeff = atomicMass[Z];
159 if (verbose > 1) {
160 G4cout << "G4NistElementBuilder: Build E
161 << " Aeff= " << Aeff;
162 G4cout << " with natural isotope compos
163 }
164
165 // build Element with its Isotopes
166 //
167 G4int nc = nIsotopes[Z];
168 G4int n0 = nFirstIsotope[Z];
169 G4int idx = idxIsotopes[Z];
170 std::vector<G4Isotope*> iso;
171 for (G4int i = 0; i < nc; ++i) {
172 if (relAbundance[idx + i] > 0.0) {
173 G4int n1 = n0 + i;
174 std::ostringstream os;
175 os << n1;
176 const G4String& ss = elmSymbol[Z] + os
177 iso.push_back(new G4Isotope(ss, Z, n1,
178 CLHEP::g / (CLHEP::mole * CLHEP::a
179 /*
180 G4cout << " Z= " << Z << " N= " << n
181 << " miso(amu)= " << GetIsotopeMass
182 << " matom(amu)= " << GetAtomicMass(
183 */
184 }
185 }
186 auto ni = (G4int)iso.size();
187 theElement = new G4Element(elmSymbol[Z], e
188 for (G4int j = 0; j < ni; ++j) {
189 G4double w = relAbundance[idx + (iso[j])
190 theElement->AddIsotope(iso[j], w);
191 }
192 theElement->SetNaturalAbundanceFlag(true);
193 elmIndex[Z] = (G4int)theElement->GetIndex(
194 }
195 l.unlock();
196 return theElement;
197 }
198
199 //....oooOO0OOooo........oooOO0OOooo........oo
200
201 void G4NistElementBuilder::PrintElement(G4int
202 {
203 G4int imin = Z;
204 G4int imax = Z + 1;
205 if (Z == 0) {
206 imin = 1;
207 imax = maxNumElements;
208 }
209 if (imax > maxNumElements) {
210 imax = maxNumElements;
211 }
212
213 for (G4int i = imin; i < imax; ++i) {
214 G4int nc = nIsotopes[i];
215 G4cout << "Nist Element: <" << elmSymbol[i
216 << " " << nc << " isotopes:" << G4
217 G4int j;
218 G4int idx = idxIsotopes[i];
219 G4int n0 = nFirstIsotope[i];
220 G4cout << " N: ";
221 for (j = 0; j < nc; ++j) {
222 G4cout << n0 + j << " ";
223 }
224 G4cout << G4endl;
225 G4cout << " mass(amu): ";
226 for (j = 0; j < nc; ++j) {
227 G4cout << GetAtomicMass(i, n0 + j) << "
228 }
229 G4cout << G4endl;
230 G4cout << " abundance: ";
231 for (j = 0; j < nc; ++j) {
232 G4cout << relAbundance[idx + j] << " ";
233 }
234 G4cout << G4endl;
235 }
236 }
237
238 //....oooOO0OOooo........oooOO0OOooo........oo
239
240 void G4NistElementBuilder::AddElement(const G4
241 const G4double& A, const G4double& sA, const
242 {
243 if (verbose > 1) {
244 G4cout << "AddElement " << name << " Z= "
245 }
246 if (Z >= maxNumElements) {
247 G4cout << "G4NistElementBuilder::AddElemen
248 return;
249 }
250
251 if (index + nc >= maxAbundance) {
252 G4cout << "G4NistElementBuilder::AddElemen
253 << " is above array size " << maxAb
254 return;
255 }
256
257 elmSymbol.push_back(name);
258 atomicMass[Z] = 0.0;
259 idxIsotopes[Z] = index;
260 nIsotopes[Z] = nc;
261
262 nFirstIsotope[Z] = (&N)[0];
263 G4double ww = 0.0;
264 G4double www;
265
266 for (G4int i = 0; i < nc; ++i) {
267 www = 0.01 * (&W)[i];
268 // mass of the isotope in G4 units
269 massIsotopes[index] = (&A)[i] * amu_c2 - Z
270 sigMass[index] = (&sA)[i] * amu_c2;
271 relAbundance[index] = www;
272
273 // computation of mean atomic mass of the
274 atomicMass[Z] += www * (&A)[i];
275 ww += www;
276 ++index;
277 }
278
279 if (ww > 0.0) {
280 G4int idx = idxIsotopes[Z];
281 atomicMass[Z] /= ww;
282 for (G4int j = 0; j < nc; ++j) {
283 relAbundance[idx + j] /= ww;
284 }
285 }
286
287 if (verbose > 1) {
288 PrintElement(Z);
289 }
290 else if (1 == verbose) {
291 G4cout << "Nist Element " << elmSymbol[Z]
292 << " " << nIsotopes[Z] << " isotop
293 }
294 }
295
296 //....oooOO0OOooo........oooOO0OOooo........oo
297
298 void G4NistElementBuilder::Initialise()
299 {
300 // Parameterisation from D.Lunney,J.M.Pearso
301 // Rev.Mod.Phys. 75 (2003) 1021
302 bindingEnergy[0] = 0.0;
303 for (G4int i = 1; i < maxNumElements; ++i) {
304 auto Z = G4double(i);
305 bindingEnergy[i] = (14.4381 * std::pow(Z,
306 }
307
308 // NIST data
309 index = 0;
310
311 // Z = 1 -----------------------------------
312 G4int HN[6] = {1, 2, 3, 4, 5, 6};
313
314 G4double HA[6] = {1.00783, 2.0141, 3.01605,
315
316 // Garantee consistence with G4 masses
317 HA[0] = (proton_mass_c2 + electron_mass_c2 -
318 HA[1] = (1.875613 * GeV + electron_mass_c2 -
319 HA[2] = (2.80925 * GeV + electron_mass_c2 -
320
321 G4double HS[6] = {4, 4, 11, 12, 102, 28};
322
323 G4double HW[6] = {99.9885, 0.0115, 0, 0, 0,
324
325 AddElement("H", 1, 6, *HN, *HA, *HS, *HW);
326
327 // Z = 2 -----------------------------------
328 G4int HeN[8] = {3, 4, 5, 6, 7, 8, 9, 10};
329
330 G4double HeA[8] = {3.01603, 4.0026, 5.01222,
331
332 // Garantee consistence with G4 masses
333 HeA[0] = (2.80923 * GeV + 2.0 * electron_mas
334 HeA[1] = (3.727417 * GeV + 2.0 * electron_ma
335
336 G4double HeS[8] = {9, 10, 50, 11, 30, 8, 70,
337
338 G4double HeW[8] = {0.000137, 99.9999, 0, 0,
339
340 AddElement("He", 2, 8, *HeN, *HeA, *HeS, *He
341
342 // Z = 3 -----------------------------------
343 G4int LiN[9] = {4, 5, 6, 7, 8, 9, 10, 11, 12
344 G4double LiA[9] = {4.02718, 5.01254, 6.01512
345
346 G4double LiS[9] = {23, 50, 5, 5, 5, 21, 16,
347
348 G4double LiW[9] = {0, 0, 7.59, 92.41, 0, 0,
349
350 AddElement("Li", 3, 9, *LiN, *LiA, *LiS, *Li
351
352 // Z = 4 -----------------------------------
353 G4int BeN[10] = {5, 6, 7, 8, 9, 10, 11, 12,
354
355 G4double BeA[10] = {
356 5.04079, 6.01973, 7.01693, 8.00531, 9.0121
357
358 G4double BeS[10] = {429, 6, 5, 4, 4, 4, 7, 1
359
360 G4double BeW[10] = {0, 0, 0, 0, 100, 0, 0, 0
361
362 AddElement("Be", 4, 10, *BeN, *BeA, *BeS, *B
363
364 // Z = 5 -----------------------------------
365 G4int BN[13] = {7, 8, 9, 10, 11, 12, 13, 14,
366
367 G4double BA[13] = {7.02992, 8.02461, 9.01333
368 15.0311, 16.0398, 17.0469, 18.0562, 19.063
369
370 G4double BS[13] = {80, 12, 11, 4, 5, 15, 12,
371
372 G4double BW[13] = {0, 0, 0, 19.9, 80.1, 0, 0
373
374 AddElement("B", 5, 13, *BN, *BA, *BS, *BW);
375
376 // Z = 6 -----------------------------------
377 G4int CN[15] = {8, 9, 10, 11, 12, 13, 14, 15
378
379 G4double CA[15] = {8.03768, 9.03104, 10.0169
380 17.0226, 18.0268, 19.0353, 20.0403, 21.049
381
382 G4double CS[15] = {25, 23, 4, 10, 0, 10, 4,
383
384 G4double CW[15] = {0, 0, 0, 0, 98.93, 1.07,
385
386 AddElement("C", 6, 15, *CN, *CA, *CS, *CW);
387
388 // Z = 7 -----------------------------------
389 G4int NN[15] = {10, 11, 12, 13, 14, 15, 16,
390
391 G4double NA[15] = {10.0426, 11.0268, 12.0186
392 18.0141, 19.017, 20.0234, 21.0271, 22.0344
393
394 G4double NS[15] = {43, 19, 11, 29, 9, 9, 28,
395
396 G4double NW[15] = {0, 0, 0, 0, 99.632, 0.368
397
398 AddElement("N", 7, 15, *NN, *NA, *NS, *NW);
399
400 // Z = 8 -----------------------------------
401 G4int ON[15] = {12, 13, 14, 15, 16, 17, 18,
402
403 G4double OA[15] = {12.0344, 13.0248, 14.0086
404 20.0041, 21.0087, 22.01, 23.0157, 24.0204,
405
406 G4double Os[15] = {20, 10, 8, 5, 15, 22, 9,
407
408 G4double OW[15] = {0, 0, 0, 0, 99.757, 0.038
409
410 AddElement("O", 8, 15, *ON, *OA, *Os, *OW);
411
412 // Z = 9 -----------------------------------
413 G4int FN[16] = {14, 15, 16, 17, 18, 19, 20,
414
415 G4double FA[16] = {14.0361, 15.018, 16.0115,
416 23.0036, 24.0081, 25.0121, 26.0196, 27.026
417
418 G4double FS[16] = {43, 14, 9, 27, 6, 7, 9, 1
419
420 G4double FW[16] = {0, 0, 0, 0, 0, 100, 0, 0,
421
422 AddElement("F", 9, 16, *FN, *FA, *FS, *FW);
423
424 // Z = 10 ----------------------------------
425 G4int NeN[17] = {16, 17, 18, 19, 20, 21, 22,
426
427 G4double NeA[17] = {16.0258, 17.0177, 18.005
428 23.9936, 24.9978, 26.0005, 27.0076, 28.012
429
430 G4double NeS[17] = {22, 50, 16, 6, 20, 4, 23
431
432 G4double NeW[17] = {0, 0, 0, 0, 90.48, 0.27,
433
434 AddElement("Ne", 10, 17, *NeN, *NeA, *NeS, *
435
436 // Z = 11 ----------------------------------
437 G4int NaN[18] = {18, 19, 20, 21, 22, 23, 24,
438
439 G4double NaA[18] = {18.0272, 19.0139, 20.007
440 26.994, 27.9989, 29.0028, 30.0092, 31.0136
441
442 G4double NaS[18] = {43, 13, 7, 8, 5, 23, 23,
443
444 G4double NaW[18] = {0, 0, 0, 0, 0, 100, 0, 0
445
446 AddElement("Na", 11, 18, *NaN, *NaA, *NaS, *
447
448 // Z = 12 ----------------------------------
449 G4int MgN[18] = {20, 21, 22, 23, 24, 25, 26,
450
451 G4double MgA[18] = {20.0189, 21.0117, 21.999
452 27.9839, 28.9886, 29.9905, 30.9965, 31.999
453
454 G4double MgS[18] = {29, 18, 15, 13, 20, 20,
455
456 G4double MgW[18] = {0, 0, 0, 0, 78.99, 10, 1
457
458 AddElement("Mg", 12, 18, *MgN, *MgA, *MgS, *
459
460 // Z = 13 ----------------------------------
461 G4int AlN[19] = {21, 22, 23, 24, 25, 26, 27,
462
463 G4double AlA[19] = {21.028, 22.0195, 23.0073
464 28.9804, 29.983, 30.9839, 31.9881, 32.9909
465 39.0219};
466
467 G4double AlS[19] = {32, 10, 27, 4, 7, 21, 14
468
469 G4double AlW[19] = {0, 0, 0, 0, 0, 0, 100, 0
470
471 AddElement("Al", 13, 19, *AlN, *AlA, *AlS, *
472
473 // Z = 14 ----------------------------------
474 G4int SiN[21] = {
475 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32
476
477 G4double SiA[21] = {22.0345, 23.0255, 24.011
478 29.9738, 30.9754, 31.9741, 32.978, 33.9786
479 41.0127, 42.0161};
480
481 G4double SiS[21] = {
482 22, 21, 21, 11, 3, 17, 20, 3, 5, 7, 23, 17
483
484 G4double SiW[21] = {
485 0, 0, 0, 0, 0, 0, 92.2297, 4.6832, 3.0872,
486
487 AddElement("Si", 14, 21, *SiN, *SiA, *SiS, *
488
489 // Z = 15 ----------------------------------
490 G4int PN[23] = {
491 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34
492
493 G4double PA[23] = {24.0343, 25.0203, 26.0118
494 31.9739, 32.9717, 33.9736, 34.9733, 35.978
495 42.0001, 43.0033, 44.0099, 45.0151, 46.023
496
497 G4double PS[23] = {
498 54, 21, 21, 40, 4, 8, 4, 20, 20, 12, 5, 20
499
500 G4double PW[23] = {0, 0, 0, 0, 0, 0, 0, 100,
501
502 AddElement("P", 15, 23, *PN, *PA, *PS, *PW);
503
504 // Z = 16 ----------------------------------
505 G4int SN[24] = {
506 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36
507
508 G4double SA[24] = {26.0279, 27.0188, 28.0044
509 33.9679, 34.969, 35.9671, 36.9711, 37.9712
510 44.9948, 45.9996, 47.0076, 48.013, 49.022}
511
512 G4double SS[24] = {
513 32, 22, 17, 50, 3, 16, 12, 12, 11, 10, 25,
514
515 G4double SW[24] = {
516 0, 0, 0, 0, 0, 0, 94.93, 0.76, 4.29, 0, 0.
517
518 AddElement("S", 16, 24, *SN, *SA, *SS, *SW);
519
520 // Z = 17 ----------------------------------
521 G4int ClN[24] = {
522 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38
523
524 G4double ClA[24] = {28.0285, 29.0141, 30.004
525 35.9683, 36.9659, 37.968, 38.968, 39.9704,
526 46.9879, 47.9948, 48.9999, 50.0077, 51.013
527
528 G4double ClS[24] = {
529 54, 21, 21, 50, 7, 6, 13, 4, 8, 5, 12, 19,
530
531 G4double ClW[24] = {
532 0, 0, 0, 0, 0, 0, 0, 75.78, 0, 24.22, 0, 0
533
534 AddElement("Cl", 17, 24, *ClN, *ClA, *ClS, *
535
536 // Z = 18 ----------------------------------
537 G4int ArN[24] = {
538 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40
539
540 G4double ArA[24] = {30.0216, 31.0121, 31.997
541 37.9627, 38.9643, 39.9624, 40.9645, 41.963
542 47.9751, 48.9822, 49.9859, 50.9932, 51.998
543
544 G4double ArS[24] = {
545 32, 22, 50, 30, 3, 8, 27, 3, 5, 5, 3, 7, 4
546
547 G4double ArW[24] = {
548 0, 0, 0, 0, 0, 0, 0.3365, 0, 0.0632, 0, 99
549
550 AddElement("Ar", 18, 24, *ArN, *ArA, *ArS, *
551
552 // Z = 19 ----------------------------------
553 G4int KN[24] = {
554 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42
555
556 G4double KA[24] = {32.0219, 33.0073, 33.9984
557 40.9618, 41.9624, 42.9607, 43.9616, 44.960
558 50.9764, 51.9826, 52.9871, 53.994, 54.9994
559
560 G4double KS[24] = {
561 54, 21, 32, 21, 8, 29, 8, 3, 29, 28, 3, 10
562
563 G4double KW[24] = {
564 0, 0, 0, 0, 0, 0, 0, 93.2581, 0.0117, 6.73
565
566 AddElement("K", 19, 24, *KN, *KA, *KS, *KW);
567
568 // Z = 20 ----------------------------------
569 G4int CaN[24] = {
570 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44
571
572 G4double CaA[24] = {34.0141, 35.0048, 35.993
573 41.9586, 42.9588, 43.9555, 44.9562, 45.953
574 51.9651, 52.9701, 53.9747, 54.9806, 55.985
575
576 G4double CaS[24] = {
577 32, 70, 40, 24, 5, 19, 3, 4, 4, 5, 9, 10,
578
579 G4double CaW[24] = {0, 0, 0, 0, 0, 0, 96.941
580 0, 0, 0, 0, 0, 0};
581
582 AddElement("Ca", 20, 24, *CaN, *CaA, *CaS, *
583
584 // Z = 21 ----------------------------------
585 G4int ScN[24] = {
586 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46
587
588 G4double ScA[24] = {36.0149, 37.0031, 37.994
589 43.9594, 44.9559, 45.9552, 46.9524, 47.952
590 54.9674, 55.9727, 56.977, 57.9831, 58.988}
591
592 G4double ScS[24] = {
593 54, 32, 32, 26, 4, 3, 4, 20, 19, 12, 12, 2
594
595 G4double ScW[24] = {0, 0, 0, 0, 0, 0, 0, 0,
596
597 AddElement("Sc", 21, 24, *ScN, *ScA, *ScS, *
598
599 // Z = 22 ----------------------------------
600 G4int TiN[24] = {
601 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48
602
603 G4double TiA[24] = {38.0098, 39.0013, 39.990
604 45.9526, 46.9518, 47.9479, 48.9479, 49.944
605 55.958, 56.9629, 57.9661, 58.972, 59.9756,
606
607 G4double TiS[24] = {
608 27, 11, 17, 40, 6, 7, 8, 13, 12, 10, 10, 1
609
610 G4double TiW[24] = {
611 0, 0, 0, 0, 0, 0, 0, 0, 8.25, 7.44, 73.72,
612
613 AddElement("Ti", 22, 24, *TiN, *TiA, *TiS, *
614
615 // Z = 23 ----------------------------------
616 G4int VN[24] = {
617 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50
618
619 G4double VA[24] = {40.0111, 40.9997, 41.9912
620 47.9523, 48.9485, 49.9472, 50.944, 51.9448
621 57.9567, 58.9593, 59.9645, 60.9674, 61.973
622
623 G4double VS[24] = {
624 54, 27, 21, 25, 90, 18, 16, 12, 28, 14, 14
625
626 G4double VW[24] = {0, 0, 0, 0, 0, 0, 0, 0, 0
627
628 AddElement("V", 23, 24, *VN, *VA, *VS, *VW);
629
630 // Z = 24 ----------------------------------
631 G4int CrN[24] = {
632 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52
633
634 G4double CrA[24] = {42.0064, 42.9977, 43.985
635 50.9448, 51.9405, 52.9407, 53.9389, 54.940
636 60.9541, 61.9558, 62.9619, 63.9642, 64.970
637
638 G4double CrS[24] = {
639 32, 90, 14, 11, 22, 15, 8, 28, 14, 14, 15,
640
641 G4double CrW[24] = {
642 0, 0, 0, 0, 0, 0, 0, 0, 4.345, 0, 83.789,
643
644 AddElement("Cr", 24, 24, *CrN, *CrA, *CrS, *
645
646 // Z = 25 ----------------------------------
647 G4int MnN[24] = {
648 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54
649
650 G4double MnA[24] = {44.0069, 44.9945, 45.986
651 51.9456, 52.9413, 53.9404, 54.938, 55.9389
652 62.9498, 63.9537, 64.9561, 65.9608, 66.963
653
654 G4double MnS[24] = {
655 54, 32, 12, 17, 80, 26, 15, 14, 25, 15, 18
656
657 G4double MnW[24] = {0, 0, 0, 0, 0, 0, 0, 0,
658
659 AddElement("Mn", 25, 24, *MnN, *MnA, *MnS, *
660
661 // Z = 26 ----------------------------------
662 G4int FeN[25] = {45, 46, 47, 48, 49, 50, 51,
663 65, 66, 67, 68, 69};
664
665 G4double FeA[25] = {45.0146, 46.0008, 46.992
666 52.9453, 53.9396, 54.9383, 55.9349, 56.935
667 62.9401, 63.9409, 64.9449, 65.946, 66.95,
668
669 G4double FeS[25] = {43, 38, 28, 11, 17, 60,
670 30, 35, 50, 75, 86};
671
672 G4double FeW[25] = {
673 0, 0, 0, 0, 0, 0, 0, 0, 0, 5.845, 0, 91.75
674
675 AddElement("Fe", 26, 25, *FeN, *FeA, *FeS, *
676
677 // Z = 27 ----------------------------------
678 G4int CoN[25] = {48, 49, 50, 51, 52, 53, 54,
679 68, 69, 70, 71, 72};
680
681 G4double CoA[25] = {48.0018, 48.9897, 49.981
682 55.9398, 56.9363, 57.9358, 58.9332, 59.933
683 65.9398, 66.9406, 67.9444, 68.9452, 69.949
684
685 G4double CoS[25] = {43, 28, 18, 16, 70, 19,
686 30, 35, 40, 75, 86, 86};
687
688 G4double CoW[25] = {0, 0, 0, 0, 0, 0, 0, 0,
689
690 AddElement("Co", 27, 25, *CoN, *CoA, *CoS, *
691
692 // Z = 28 ----------------------------------
693 G4int NiN[29] = {50, 51, 52, 53, 54, 55, 56,
694 70, 71, 72, 73, 74, 75, 76, 77, 78};
695
696 G4double NiA[29] = {49.9959, 50.9877, 51.975
697 57.9353, 58.9344, 59.9308, 60.9311, 61.928
698 67.9318, 68.9352, 69.9361, 70.94, 71.9413,
699
700 G4double NiS[29] = {28, 28, 90, 17, 50, 12,
701 35, 40, 50, 64, 75, 86, 97, 107, 118};
702
703 G4double NiW[29] = {0, 0, 0, 0, 0, 0, 0, 0,
704 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
705
706 AddElement("Ni", 28, 29, *NiN, *NiA, *NiS, *
707
708 // Z = 29 ----------------------------------
709 G4int CuN[29] = {52, 53, 54, 55, 56, 57, 58,
710 72, 73, 74, 75, 76, 77, 78, 79, 80};
711
712 G4double CuA[29] = {51.9972, 52.9856, 53.976
713 59.9374, 60.9335, 61.9326, 62.9296, 63.929
714 69.9324, 70.9326, 71.9355, 72.9365, 73.940
715 79.9619};
716
717 G4double CuS[29] = {28, 28, 23, 32, 15, 17,
718 21, 32, 43, 54, 64, 75, 86, 97, 97};
719
720 G4double CuW[29] = {
721 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 69.17, 0,
722
723 AddElement("Cu", 29, 29, *CuN, *CuA, *CuS, *
724
725 // Z = 30 ----------------------------------
726 G4int ZnN[29] = {54, 55, 56, 57, 58, 59, 60,
727 74, 75, 76, 77, 78, 79, 80, 81, 82};
728
729 G4double ZnA[29] = {53.993, 54.984, 55.9724,
730 62.9332, 63.9291, 64.9292, 65.926, 66.9271
731 72.9298, 73.9295, 74.9329, 75.9334, 76.937
732
733 G4double ZnS[29] = {43, 27, 28, 15, 50, 40,
734 50, 80, 13, 14, 17, 29, 18, 43, 43};
735
736 G4double ZnW[29] = {0, 0, 0, 0, 0, 0, 0, 0,
737 0, 0, 0, 0, 0, 0, 0, 0};
738
739 AddElement("Zn", 30, 29, *ZnN, *ZnA, *ZnS, *
740
741 // Z = 31 ----------------------------------
742 G4int GaN[29] = {56, 57, 58, 59, 60, 61, 62,
743 76, 77, 78, 79, 80, 81, 82, 83, 84};
744
745 G4double GaA[29] = {55.9949, 56.9829, 57.974
746 63.9368, 64.9327, 65.9316, 66.9282, 67.928
747 74.9265, 75.9289, 76.9293, 77.9317, 78.932
748
749 G4double GaS[29] = {28, 28, 23, 18, 12, 21,
750 60, 90, 13, 13, 21, 32, 54, 64};
751
752 G4double GaW[29] = {0, 0, 0, 0, 0, 0, 0, 0,
753 0, 0, 0, 0, 0, 0};
754
755 AddElement("Ga", 31, 29, *GaN, *GaA, *GaS, *
756
757 // Z = 32 ----------------------------------
758 G4int GeN[29] = {58, 59, 60, 61, 62, 63, 64,
759 78, 79, 80, 81, 82, 83, 84, 85, 86};
760
761 G4double GeA[29] = {57.991, 58.9817, 59.9702
762 65.9339, 66.9327, 67.9281, 68.928, 69.9243
763 76.9235, 77.9229, 78.9254, 79.9254, 80.928
764
765 G4double GeS[29] = {34, 30, 25, 32, 15, 21,
766 4, 10, 25, 13, 26, 32, 43, 54, 64};
767
768 G4double GeW[29] = {0, 0, 0, 0, 0, 0, 0, 0,
769 0, 0, 0, 0, 0, 0, 0, 0, 0};
770
771 AddElement("Ge", 32, 29, *GeN, *GeA, *GeS, *
772
773 // Z = 33 ----------------------------------
774 G4int AsN[30] = {60, 61, 62, 63, 64, 65, 66,
775 80, 81, 82, 83, 84, 85, 86, 87, 88, 89};
776
777 G4double AsA[30] = {59.9931, 60.9806, 61.973
778 67.9368, 68.9323, 69.9309, 70.9271, 71.926
779 77.9218, 78.9209, 79.9226, 80.9221, 81.924
780 87.9446, 88.9492};
781
782 G4double AsS[30] = {64, 64, 32, 54, 38, 42,
783 23, 6, 21, 24, 32, 32, 43, 54, 64, 64};
784
785 G4double AsW[30] = {
786 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
787
788 AddElement("As", 33, 30, *AsN, *AsA, *AsS, *
789
790 // Z = 34 ----------------------------------
791 G4int SeN[28] = {65, 66, 67, 68, 69, 70, 71,
792 85, 86, 87, 88, 89, 90, 91, 92};
793
794 G4double SeA[28] = {64.9647, 65.9552, 66.950
795 72.9268, 73.9225, 74.9225, 75.9192, 76.919
796 82.9191, 83.9185, 84.9222, 85.9243, 86.928
797
798 G4double SeS[28] = {64, 32, 21, 32, 40, 22,
799 30, 17, 40, 50, 32, 43, 54, 64};
800
801 G4double SeW[28] = {0, 0, 0, 0, 0, 0, 0, 0,
802 0, 0, 0, 0, 0, 0, 0, 0, 0};
803
804 AddElement("Se", 34, 28, *SeN, *SeA, *SeS, *
805
806 // Z = 35 ----------------------------------
807 G4int BrN[28] = {67, 68, 69, 70, 71, 72, 73,
808 87, 88, 89, 90, 91, 92, 93, 94};
809
810 G4double BrA[28] = {66.9648, 67.9583, 68.950
811 74.9258, 75.9245, 76.9214, 77.9211, 78.918
812 84.9156, 85.9188, 86.9207, 87.9241, 88.926
813
814 G4double BrS[28] = {54, 58, 34, 39, 32, 28,
815 40, 60, 80, 80, 50, 32, 43};
816
817 G4double BrW[28] = {
818 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 50.69,
819
820 AddElement("Br", 35, 28, *BrN, *BrA, *BrS, *
821
822 // Z = 36 ----------------------------------
823 G4int KrN[29] = {69, 70, 71, 72, 73, 74, 75,
824 89, 90, 91, 92, 93, 94, 95, 96, 97};
825
826 G4double KrA[29] = {68.9653, 69.956, 70.9505
827 77.9204, 78.9201, 79.9164, 80.9166, 81.913
828 87.9144, 88.9176, 89.9195, 90.9234, 91.926
829
830 G4double KrS[29] = {54, 43, 32, 29, 15, 60,
831 20, 60, 13, 11, 32, 43, 54, 54};
832
833 G4double KrW[29] = {0, 0, 0, 0, 0, 0, 0, 0,
834 0, 0, 0, 0, 0, 0, 0, 0, 0};
835
836 AddElement("Kr", 36, 29, *KrN, *KrA, *KrS, *
837
838 // Z = 37 ----------------------------------
839 G4int RbN[32] = {71, 72, 73, 74, 75, 76, 77,
840 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 1
841
842 G4double RbA[32] = {70.9653, 71.9591, 72.950
843 78.924, 79.9225, 80.919, 81.9182, 82.9151,
844 89.9148, 90.9165, 91.9197, 92.922, 93.9264
845 99.9499, 100.953, 101.959};
846
847 G4double RbS[32] = {54, 54, 52, 77, 8, 8, 8,
848 21, 27, 30, 40, 16, 32, 18, 54};
849
850 G4double RbW[32] = {0, 0, 0, 0, 0, 0, 0, 0,
851 0, 0, 0, 0, 0, 0, 0, 0, 0};
852
853 AddElement("Rb", 37, 32, *RbN, *RbA, *RbS, *
854
855 // Z = 38 ----------------------------------
856 G4int SrN[32] = {73, 74, 75, 76, 77, 78, 79,
857 93, 94, 95, 96, 97, 98, 99, 100, 101, 102,
858
859 G4double SrA[32] = {72.966, 73.9563, 74.9499
860 80.9232, 81.9184, 82.9176, 83.9134, 84.912
861 90.9102, 91.911, 92.914, 93.9154, 94.9194,
862 101.943, 102.949, 103.952};
863
864 G4double SrS[32] = {64, 54, 32, 32, 16, 8, 9
865 26, 20, 27, 15, 14, 13, 12, 54, 75};
866
867 G4double SrW[32] = {0, 0, 0, 0, 0, 0, 0, 0,
868 0, 0, 0, 0, 0, 0, 0, 0, 0};
869
870 AddElement("Sr", 38, 32, *SrN, *SrA, *SrS, *
871
872 // Z = 39 ----------------------------------
873 G4int YN[30] = {77, 78, 79, 80, 81, 82, 83,
874 97, 98, 99, 100, 101, 102, 103, 104, 105,
875
876 G4double YA[30] = {76.9496, 77.9435, 78.9373
877 84.9164, 85.9149, 86.9109, 87.9095, 88.905
878 94.9128, 95.9159, 96.9181, 97.9222, 98.924
879 104.945, 105.95};
880
881 G4double YS[30] = {32, 43, 48, 43, 70, 11, 5
882 13, 26, 26, 80, 10, 90, 32, 43, 54, 75};
883
884 G4double YW[30] = {
885 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0
886
887 AddElement("Y", 39, 30, *YN, *YA, *YS, *YW);
888
889 // Z = 40 ----------------------------------
890 G4int ZrN[30] = {79, 80, 81, 82, 83, 84, 85,
891 99, 100, 101, 102, 103, 104, 105, 106, 107
892
893 G4double ZrA[30] = {78.9492, 79.9406, 80.936
894 86.9148, 87.9102, 88.9089, 89.9047, 90.905
895 97.9127, 98.9165, 99.9178, 100.921, 101.92
896 107.944};
897
898 G4double ZrS[30] = {43, 32, 32, 55, 10, 21,
899 21, 40, 30, 50, 12, 43, 43, 54, 64, 75};
900
901 G4double ZrW[30] = {0, 0, 0, 0, 0, 0, 0, 0,
902 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
903
904 AddElement("Zr", 40, 30, *ZrN, *ZrA, *ZrS, *
905
906 // Z = 41 ----------------------------------
907 G4int NbN[30] = {81, 82, 83, 84, 85, 86, 87,
908 101, 102, 103, 104, 105, 106, 107, 108, 10
909
910 G4double NbA[30] = {80.949, 81.9431, 82.9367
911 89.9113, 90.907, 91.9072, 92.9064, 93.9073
912 99.9142, 100.915, 101.918, 102.919, 103.92
913 109.943};
914
915 G4double NbS[30] = {43, 32, 34, 32, 24, 90,
916 20, 40, 70, 12, 11, 32, 43, 54, 54, 64};
917
918 G4double NbW[30] = {
919 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0
920
921 AddElement("Nb", 41, 30, *NbN, *NbA, *NbS, *
922
923 // Z = 42 ----------------------------------
924 G4int MoN[31] = {83, 84, 85, 86, 87, 88, 89,
925 102, 103, 104, 105, 106, 107, 108, 109, 11
926
927 G4double MoA[31] = {82.9487, 83.9401, 84.936
928 90.9118, 91.9068, 92.9068, 93.9051, 94.905
929 101.91, 102.913, 103.914, 104.917, 105.918
930 112.942};
931
932 G4double MoS[31] = {54, 43, 43, 47, 24, 22,
933 70, 80, 23, 17, 21, 32, 43, 54, 64, 64};
934
935 G4double MoW[31] = {0, 0, 0, 0, 0, 0, 0, 0,
936 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
937
938 AddElement("Mo", 42, 31, *MoN, *MoA, *MoS, *
939
940 // Z = 43 ----------------------------------
941 G4int TcN[31] = {85, 86, 87, 88, 89, 90, 91,
942 104, 105, 106, 107, 108, 109, 110, 111, 11
943
944 G4double TcA[31] = {84.9489, 85.9429, 86.936
945 92.9102, 93.9097, 94.9077, 95.9079, 96.906
946 102.909, 103.911, 104.912, 105.914, 106.91
947 112.931, 113.936, 114.938};
948
949 G4double TcS[31] = {54, 32, 32, 32, 23, 26,
950 15, 16, 14, 23, 43, 43, 54, 64, 64, 75};
951
952 G4double TcW[31] = {
953 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100
954
955 AddElement("Tc", 43, 31, *TcN, *TcA, *TcS, *
956
957 // Z = 44 ----------------------------------
958 G4int RuN[32] = {87, 88, 89, 90, 91, 92, 93,
959 106, 107, 108, 109, 110, 111, 112, 113, 11
960
961 G4double RuA[32] = {86.9492, 87.9404, 88.936
962 94.9104, 95.9076, 96.9076, 97.9053, 98.905
963 104.908, 105.907, 106.91, 107.91, 108.913,
964 115.93, 116.935, 117.937};
965
966 G4double RuS[32] = {64, 54, 54, 43, 54, 32,
967 13, 70, 25, 32, 58, 54, 39, 64, 75, 86, 97
968
969 G4double RuW[32] = {0, 0, 0, 0, 0, 0, 0, 0,
970 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
971
972 AddElement("Ru", 44, 32, *RuN, *RuA, *RuS, *
973
974 // Z = 45 ----------------------------------
975 G4int RhN[33] = {89, 90, 91, 92, 93, 94, 95,
976 107, 108, 109, 110, 111, 112, 113, 114, 11
977
978 G4double RhA[33] = {88.9494, 89.9429, 90.936
979 96.9113, 97.9107, 98.9081, 99.9081, 100.90
980 106.907, 107.909, 108.909, 109.911, 110.91
981 116.925, 117.929, 118.931, 119.936, 120.93
982
983 G4double RhS[33] = {54, 54, 43, 43, 43, 48,
984 24, 22, 54, 43, 32, 54, 54, 64, 75, 86, 86
985 G4double RhW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
986 0, 0, 0, 0, 0, 0, 0, 0};
987
988 AddElement("Rh", 45, 33, *RhN, *RhA, *RhS, *
989
990 // Z = 46 ----------------------------------
991 G4int PdN[33] = {91, 92, 93, 94, 95, 96, 97,
992 109, 110, 111, 112, 113, 114, 115, 116, 11
993
994 G4double PdA[33] = {90.9495, 91.9404, 92.935
995 98.9118, 99.9085, 100.908, 101.906, 102.90
996 108.906, 109.905, 110.908, 111.907, 112.91
997 119.924, 120.928, 121.93, 122.934};
998
999 G4double PdS[33] = {64, 54, 43, 43, 43, 16,
1000 19, 40, 26, 70, 60, 32, 23, 32, 43, 54, 5
1001
1002 G4double PdW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1003 11.72, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
1004
1005 AddElement("Pd", 46, 33, *PdN, *PdA, *PdS,
1006
1007 // Z = 47 ---------------------------------
1008 G4int AgN[34] = {94, 95, 96, 97, 98, 99, 10
1009 111, 112, 113, 114, 115, 116, 117, 118, 1
1010
1011 G4double AgA[34] = {93.9428, 94.9355, 95.93
1012 101.912, 102.909, 103.909, 104.907, 105.9
1013 111.907, 112.907, 113.909, 114.909, 115.9
1014 121.923, 122.925, 123.929, 124.931, 125.9
1015
1016 G4double AgS[34] = {54, 43, 43, 43, 16, 16,
1017 40, 50, 50, 70, 10, 80, 16, 22, 32, 43, 4
1018
1019 G4double AgW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1020 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1021
1022 AddElement("Ag", 47, 34, *AgN, *AgA, *AgS,
1023
1024 // Z = 48 ---------------------------------
1025 G4int CdN[35] = {96, 97, 98, 99, 100, 101,
1026 113, 114, 115, 116, 117, 118, 119, 120, 1
1027
1028 G4double CdA[35] = {95.9398, 96.9349, 97.92
1029 104.909, 105.906, 106.907, 107.904, 108.9
1030 114.905, 115.905, 116.907, 117.907, 118.9
1031 125.922, 126.926, 127.928, 128.932, 129.9
1032
1033 G4double CdS[35] = {54, 43, 22, 22, 10, 16,
1034 22, 90, 20, 90, 22, 40, 70, 70, 60, 80, 3
1035
1036 G4double CdW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1037 28.73, 0, 7.49, 0, 0, 0, 0, 0, 0, 0, 0, 0
1038
1039 AddElement("Cd", 48, 35, *CdN, *CdA, *CdS,
1040
1041 // Z = 49 ---------------------------------
1042 G4int InN[37] = {98, 99, 100, 101, 102, 103
1043 114, 115, 116, 117, 118, 119, 120, 121, 1
1044 133, 134};
1045
1046 G4double InA[37] = {97.9422, 98.9346, 99.93
1047 105.913, 106.91, 107.91, 108.907, 109.907
1048 116.905, 117.906, 118.906, 119.908, 120.9
1049 127.92, 128.922, 129.925, 130.927, 131.93
1050
1051 G4double InS[37] = {54, 54, 41, 32, 41, 27,
1052 40, 29, 50, 26, 50, 30, 40, 40, 50, 14, 5
1053
1054 G4double InW[37] = {0, 0, 0, 0, 0, 0, 0, 0,
1055 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1056
1057 AddElement("In", 49, 37, *InN, *InA, *InS,
1058
1059 // Z = 50 ---------------------------------
1060 G4int SnN[38] = {100, 101, 102, 103, 104, 1
1061 116, 117, 118, 119, 120, 121, 122, 123, 1
1062 135, 136, 137};
1063
1064 G4double SnA[38] = {99.939, 100.936, 101.93
1065 108.911, 109.908, 110.908, 111.905, 112.9
1066 118.903, 119.902, 120.904, 121.903, 122.9
1067 128.913, 129.914, 130.917, 131.918, 132.9
1068
1069 G4double SnS[38] = {46, 54, 43, 32, 16, 10,
1070 29, 29, 15, 16, 11, 27, 29, 13, 30, 80, 2
1071
1072 G4double SnW[38] = {0, 0, 0, 0, 0, 0, 0, 0,
1073 8.59, 32.58, 0, 4.63, 0, 5.79, 0, 0, 0, 0
1074
1075 AddElement("Sn", 50, 38, *SnN, *SnA, *SnS,
1076
1077 // Z = 51 ---------------------------------
1078 G4int SbN[37] = {103, 104, 105, 106, 107, 1
1079 119, 120, 121, 122, 123, 124, 125, 126, 1
1080 138, 139};
1081
1082 G4double SbA[37] = {102.94, 103.936, 104.93
1083 110.913, 111.912, 112.909, 113.909, 114.9
1084 120.904, 121.905, 122.904, 123.906, 124.9
1085 130.912, 131.914, 132.915, 133.921, 134.9
1086
1087 G4double SbS[37] = {54, 39, 17, 34, 32, 22,
1088 22, 22, 3, 30, 6, 27, 23, 27, 80, 25, 80,
1089
1090 G4double SbW[37] = {0, 0, 0, 0, 0, 0, 0, 0,
1091 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
1092
1093 AddElement("Sb", 51, 37, *SbN, *SbA, *SbS,
1094
1095 // Z = 52 ---------------------------------
1096 G4int TeN[37] = {106, 107, 108, 109, 110, 1
1097 122, 123, 124, 125, 126, 127, 128, 129, 1
1098 141, 142};
1099
1100 G4double TeA[37] = {105.938, 106.935, 107.9
1101 113.912, 114.912, 115.908, 116.909, 117.9
1102 123.903, 124.904, 125.903, 126.905, 127.9
1103 133.912, 134.916, 135.92, 136.925, 137.92
1104
1105 G4double TeS[37] = {43, 32, 16, 80, 60, 80,
1106 20, 4, 19, 3, 21, 22, 12, 80, 40, 10, 50,
1107
1108 G4double TeW[37] = {0, 0, 0, 0, 0, 0, 0, 0,
1109 18.84, 0, 31.74, 0, 34.08, 0, 0, 0, 0, 0,
1110
1111 AddElement("Te", 52, 37, *TeN, *TeA, *TeS,
1112
1113 // Z = 53 ---------------------------------
1114 G4int IN[37] = {108, 109, 110, 111, 112, 11
1115 124, 125, 126, 127, 128, 129, 130, 131, 1
1116 143, 144};
1117
1118 G4double IA[37] = {107.943, 108.938, 109.93
1119 116.914, 117.913, 118.91, 119.91, 120.907
1120 127.906, 128.905, 129.907, 130.906, 131.9
1121 138.926, 139.931, 140.935, 141.94, 142.94
1122
1123 G4double IS[37] = {39, 16, 33, 32, 23, 60,
1124 4, 4, 12, 11, 28, 16, 25, 50, 30, 90, 30,
1125
1126 G4double IW[37] = {0, 0, 0, 0, 0, 0, 0, 0,
1127 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1128
1129 AddElement("I", 53, 37, *IN, *IA, *IS, *IW)
1130
1131 // Z = 54 ---------------------------------
1132 G4int XeN[38] = {110, 111, 112, 113, 114, 1
1133 126, 127, 128, 129, 130, 131, 132, 133, 1
1134 145, 146, 147};
1135
1136 G4double XeA[38] = {109.944, 110.942, 111.9
1137 117.917, 118.916, 119.912, 120.911, 121.9
1138 127.904, 128.905, 129.904, 130.905, 131.9
1139 137.914, 138.919, 139.922, 140.927, 141.9
1140
1141 G4double XeS[38] = {43, 33, 16, 10, 22, 26,
1142 10, 10, 12, 4, 9, 11, 8, 8, 40, 23, 70, 1
1143
1144 G4double XeW[38] = {0, 0, 0, 0, 0, 0, 0, 0,
1145 21.18, 26.89, 0, 10.44, 0, 8.87, 0, 0, 0,
1146
1147 AddElement("Xe", 54, 38, *XeN, *XeA, *XeS,
1148
1149 // Z = 55 ---------------------------------
1150 G4int CsN[40] = {112, 113, 114, 115, 116, 1
1151 128, 129, 130, 131, 132, 133, 134, 135, 1
1152 147, 148, 149, 150, 151};
1153
1154 G4double CsA[40] = {111.95, 112.945, 113.94
1155 119.921, 120.917, 121.916, 122.913, 123.9
1156 129.907, 130.905, 131.906, 132.905, 133.9
1157 139.917, 140.92, 141.924, 142.927, 143.93
1158 150.962};
1159
1160 G4double CsS[40] = {33, 16, 33, 46, 38, 60,
1161 3, 3, 3, 4, 3, 10, 5, 9, 11, 11, 24, 30,
1162
1163 G4double CsW[40] = {0, 0, 0, 0, 0, 0, 0, 0,
1164 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1165
1166 AddElement("Cs", 55, 40, *CsN, *CsA, *CsS,
1167
1168 // Z = 56 ---------------------------------
1169 G4int BaN[40] = {114, 115, 116, 117, 118, 1
1170 130, 131, 132, 133, 134, 135, 136, 137, 1
1171 149, 150, 151, 152, 153};
1172
1173 G4double BaA[40] = {113.951, 114.948, 115.9
1174 121.92, 122.919, 123.915, 124.915, 125.91
1175 131.905, 132.906, 133.905, 134.906, 135.9
1176 141.916, 142.921, 143.923, 144.927, 145.9
1177 151.954, 152.96};
1178
1179 G4double BaS[40] = {48, 64, 54, 70, 54, 109
1180 3, 3, 3, 3, 3, 3, 9, 9, 7, 14, 15, 60, 80
1181
1182 G4double BaW[40] = {0, 0, 0, 0, 0, 0, 0, 0,
1183 6.592, 7.854, 11.232, 71.698, 0, 0, 0, 0,
1184
1185 AddElement("Ba", 56, 40, *BaN, *BaA, *BaS,
1186
1187 // Z = 57 ---------------------------------
1188 G4int LaN[39] = {117, 118, 119, 120, 121, 1
1189 133, 134, 135, 136, 137, 138, 139, 140, 1
1190 152, 153, 154, 155};
1191
1192 G4double LaA[39] = {116.95, 117.947, 118.94
1193 124.921, 125.919, 126.916, 127.915, 128.9
1194 135.908, 136.906, 137.907, 138.906, 139.9
1195 145.926, 146.928, 147.932, 148.934, 149.9
1196
1197 G4double LaS[39] = {96, 86, 75, 64, 54, 54,
1198 80, 50, 4, 3, 3, 5, 6, 16, 60, 70, 80, 80
1199
1200 G4double LaW[39] = {0, 0, 0, 0, 0, 0, 0, 0,
1201 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1202
1203 AddElement("La", 57, 39, *LaN, *LaA, *LaS,
1204
1205 // Z = 58 ---------------------------------
1206 G4int CeN[39] = {119, 120, 121, 122, 123, 1
1207 135, 136, 137, 138, 139, 140, 141, 142, 1
1208 154, 155, 156, 157};
1209
1210 G4double CeA[39] = {118.953, 119.947, 120.9
1211 126.923, 127.919, 128.918, 129.915, 130.9
1212 136.908, 137.906, 138.907, 139.905, 140.9
1213 146.923, 147.924, 148.928, 149.93, 150.93
1214 156.956};
1215
1216 G4double CeS[39] = {97, 86, 75, 64, 54, 54,
1217 11, 8, 3, 3, 4, 4, 4, 40, 70, 60, 13, 80,
1218
1219 G4double CeW[39] = {0, 0, 0, 0, 0, 0, 0, 0,
1220 0, 11.114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1221
1222 AddElement("Ce", 58, 39, *CeN, *CeA, *CeS,
1223
1224 // Z = 59 ---------------------------------
1225 G4int PrN[39] = {121, 122, 123, 124, 125, 1
1226 137, 138, 139, 140, 141, 142, 143, 144, 1
1227 156, 157, 158, 159};
1228
1229 G4double PrA[39] = {120.955, 121.952, 122.9
1230 128.925, 129.923, 130.92, 131.919, 132.91
1231 138.909, 139.909, 140.908, 141.91, 142.91
1232 148.924, 149.927, 150.928, 151.932, 152.9
1233 158.955};
1234
1235 G4double PrS[39] = {86, 86, 75, 64, 54, 54,
1236 3, 3, 3, 4, 8, 60, 40, 10, 11, 90, 40, 32
1237
1238 G4double PrW[39] = {0, 0, 0, 0, 0, 0, 0, 0,
1239 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
1240
1241 AddElement("Pr", 59, 39, *PrN, *PrA, *PrS,
1242
1243 // Z = 60 ---------------------------------
1244 G4int NdN[36] = {126, 127, 128, 129, 130, 1
1245 142, 143, 144, 145, 146, 147, 148, 149, 1
1246 161};
1247
1248 G4double NdA[36] = {125.943, 126.941, 127.9
1249 133.919, 134.918, 135.915, 136.915, 137.9
1250 144.913, 145.913, 146.916, 147.917, 148.9
1251 154.933, 155.935, 156.939, 157.942, 158.9
1252
1253 G4double NdS[36] = {75, 64, 64, 39, 54, 50,
1254 3, 3, 3, 4, 4, 30, 29, 12, 16, 43, 54, 64
1255
1256 G4double NdW[36] = {0, 0, 0, 0, 0, 0, 0, 0,
1257 0, 5.7, 0, 5.6, 0, 0, 0, 0, 0, 0, 0, 0, 0
1258
1259 AddElement("Nd", 60, 36, *NdN, *NdA, *NdS,
1260
1261 // Z = 61 ---------------------------------
1262 G4int PmN[36] = {128, 129, 130, 131, 132, 1
1263 144, 145, 146, 147, 148, 149, 150, 151, 1
1264 163};
1265
1266 G4double PmA[36] = {127.948, 128.943, 129.9
1267 136.921, 137.919, 138.917, 139.916, 140.9
1268 146.915, 147.917, 148.918, 149.921, 150.9
1269 156.933, 157.937, 158.939, 159.943, 160.9
1270
1271 G4double PmS[36] = {97, 86, 75, 64, 54, 54,
1272 5, 22, 6, 80, 12, 80, 30, 40, 32, 43, 54,
1273
1274 G4double PmW[36] = {0, 0, 0, 0, 0, 0, 0, 0,
1275 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1276
1277 AddElement("Pm", 61, 36, *PmN, *PmA, *PmS,
1278
1279 // Z = 62 ---------------------------------
1280 G4int SmN[36] = {130, 131, 132, 133, 134, 1
1281 146, 147, 148, 149, 150, 151, 152, 153, 1
1282 165};
1283
1284 G4double SmA[36] = {129.949, 130.946, 131.9
1285 137.924, 138.922, 139.919, 140.918, 141.9
1286 147.915, 148.917, 149.917, 150.92, 151.92
1287 158.933, 159.935, 160.939, 161.941, 162.9
1288
1289 G4double SmS[36] = {97, 97, 75, 64, 54, 54,
1290 3, 3, 3, 3, 10, 50, 80, 32, 43, 54, 64, 7
1291
1292 G4double SmW[36] = {0, 0, 0, 0, 0, 0, 0, 0,
1293 7.38, 0, 26.75, 0, 22.75, 0, 0, 0, 0, 0,
1294
1295 AddElement("Sm", 62, 36, *SmN, *SmA, *SmS,
1296
1297 // Z = 63 ---------------------------------
1298 G4int EuN[36] = {132, 133, 134, 135, 136, 1
1299 148, 149, 150, 151, 152, 153, 154, 155, 1
1300 167};
1301
1302 G4double EuA[36] = {131.954, 132.949, 133.9
1303 140.925, 141.923, 142.92, 143.919, 144.91
1304 151.922, 152.921, 153.923, 154.923, 155.9
1305 161.937, 162.939, 163.943, 164.946, 165.9
1306
1307 G4double EuS[36] = {97, 97, 75, 64, 54, 54,
1308 3, 3, 3, 6, 7, 80, 9, 22, 32, 43, 54, 64,
1309
1310 G4double EuW[36] = {0, 0, 0, 0, 0, 0, 0, 0,
1311 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1312
1313 AddElement("Eu", 63, 36, *EuN, *EuA, *EuS,
1314
1315 // Z = 64 ---------------------------------
1316 G4int GdN[34] = {136, 137, 138, 139, 140, 1
1317 152, 153, 154, 155, 156, 157, 158, 159, 1
1318
1319 G4double GdA[34] = {135.947, 136.945, 137.9
1320 143.923, 144.922, 145.918, 146.919, 147.9
1321 154.923, 155.922, 156.924, 157.924, 158.9
1322 164.939, 165.942, 166.946, 167.948, 168.9
1323
1324 G4double GdS[34] = {75, 64, 54, 54, 43, 32,
1325 3, 3, 3, 5, 32, 43, 54, 64, 64, 75, 86};
1326
1327 G4double GdW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1328 15.65, 24.84, 0, 21.86, 0, 0, 0, 0, 0, 0,
1329
1330 AddElement("Gd", 64, 34, *GdN, *GdA, *GdS,
1331
1332 // Z = 65 ---------------------------------
1333 G4int TbN[34] = {138, 139, 140, 141, 142, 1
1334 154, 155, 156, 157, 158, 159, 160, 161, 1
1335
1336 G4double TbA[34] = {137.953, 138.948, 139.9
1337 145.927, 146.924, 147.924, 148.923, 149.9
1338 155.925, 156.924, 157.925, 158.925, 159.9
1339 165.938, 166.94, 167.944, 168.946, 169.95
1340
1341 G4double TbS[34] = {86, 75, 96, 64, 82, 43,
1342 3, 3, 3, 40, 5, 11, 21, 32, 43, 54, 64, 7
1343
1344 G4double TbW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1345 0, 0, 0, 0, 0, 0, 0, 0, 0};
1346
1347 AddElement("Tb", 65, 34, *TbN, *TbA, *TbS,
1348
1349 // Z = 66 ---------------------------------
1350 G4int DyN[34] = {140, 141, 142, 143, 144, 1
1351 156, 157, 158, 159, 160, 161, 162, 163, 1
1352
1353 G4double DyA[34] = {139.954, 140.951, 141.9
1354 147.927, 148.927, 149.926, 150.926, 151.9
1355 157.924, 158.926, 159.925, 160.927, 161.9
1356 167.937, 168.94, 169.943, 170.946, 171.94
1357
1358 G4double DyS[34] = {97, 75, 85, 54, 43, 32,
1359 3, 3, 3, 3, 3, 60, 32, 32, 43, 54, 64, 75
1360
1361 G4double DyW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1362 25.51, 24.9, 28.18, 0, 0, 0, 0, 0, 0, 0,
1363
1364 AddElement("Dy", 66, 34, *DyN, *DyA, *DyS,
1365
1366 // Z = 67 ---------------------------------
1367 G4int HoN[34] = {142, 143, 144, 145, 146, 1
1368 158, 159, 160, 161, 162, 163, 164, 165, 1
1369
1370 G4double HoA[34] = {141.96, 142.955, 143.95
1371 150.932, 151.932, 152.93, 153.931, 154.92
1372 161.929, 162.929, 163.93, 164.93, 165.932
1373 172.947, 173.951, 174.954};
1374
1375 G4double HoS[34] = {107, 75, 64, 64, 54, 43
1376 5, 3, 3, 3, 3, 6, 30, 22, 50, 64, 43, 43,
1377
1378 G4double HoW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1379 0, 0, 0, 0, 0, 0, 0, 0, 0};
1380
1381 AddElement("Ho", 67, 34, *HoN, *HoA, *HoS,
1382
1383 // Z = 68 ---------------------------------
1384 G4int ErN[34] = {144, 145, 146, 147, 148, 1
1385 160, 161, 162, 163, 164, 165, 166, 167, 1
1386
1387 G4double ErA[34] = {143.961, 144.957, 145.9
1388 151.935, 152.935, 153.933, 154.933, 155.9
1389 162.93, 163.929, 164.931, 165.93, 166.932
1390 173.944, 174.948, 175.95, 176.954};
1391
1392 G4double ErS[34] = {86, 75, 64, 54, 43, 51,
1393 4, 3, 3, 3, 3, 3, 3, 5, 21, 32, 43, 43, 6
1394
1395 G4double ErW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1396 22.93, 26.78, 0, 14.93, 0, 0, 0, 0, 0, 0,
1397
1398 AddElement("Er", 68, 34, *ErN, *ErA, *ErS,
1399
1400 // Z = 69 ---------------------------------
1401 G4int TmN[34] = {146, 147, 148, 149, 150, 1
1402 162, 163, 164, 165, 166, 167, 168, 169, 1
1403
1404 G4double TmA[34] = {145.966, 146.961, 147.9
1405 153.941, 154.939, 155.939, 156.937, 157.9
1406 163.933, 164.932, 165.934, 166.933, 167.9
1407 173.942, 174.944, 175.947, 176.949, 177.9
1408
1409 G4double TmS[34] = {75, 64, 75, 64, 54, 15,
1410 12, 3, 4, 3, 3, 3, 7, 6, 50, 50, 11, 32,
1411
1412 G4double TmW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1413 0, 0, 0, 0, 0, 0, 0, 0, 0};
1414
1415 AddElement("Tm", 69, 34, *TmN, *TmA, *TmS,
1416
1417 // Z = 70 ---------------------------------
1418 G4int YbN[34] = {148, 149, 150, 151, 152, 1
1419 164, 165, 166, 167, 168, 169, 170, 171, 1
1420
1421 G4double YbA[34] = {147.967, 148.963, 149.9
1422 155.943, 156.943, 157.94, 158.94, 159.938
1423 166.935, 167.934, 168.935, 169.935, 170.9
1424 176.945, 177.947, 178.95, 179.952, 180.95
1425
1426 G4double YbS[34] = {86, 75, 64, 34, 38, 32,
1427 5, 5, 3, 3, 30, 30, 30, 30, 3, 3, 11, 32,
1428
1429 G4double YbW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1430 14.28, 21.83, 16.13, 31.83, 0, 12.76, 0,
1431
1432 AddElement("Yb", 70, 34, *YbN, *YbA, *YbS,
1433
1434 // Z = 71 ---------------------------------
1435 G4int LuN[35] = {150, 151, 152, 153, 154, 1
1436 166, 167, 168, 169, 170, 171, 172, 173, 1
1437
1438 G4double LuA[35] = {149.973, 150.967, 151.9
1439 157.949, 158.947, 159.946, 160.944, 161.9
1440 168.938, 169.938, 170.938, 171.939, 172.9
1441 178.947, 179.95, 180.952, 181.955, 182.95
1442
1443 G4double LuS[35] = {75, 65, 75, 64, 54, 14,
1444 20, 3, 4, 3, 30, 28, 28, 28, 3, 6, 80, 32
1445
1446 G4double LuW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1447 97.41, 2.59, 0, 0, 0, 0, 0, 0, 0, 0};
1448
1449 AddElement("Lu", 71, 35, *LuN, *LuA, *LuS,
1450
1451 // Z = 72 ---------------------------------
1452 G4int HfN[33] = {154, 155, 156, 157, 158, 1
1453 170, 171, 172, 173, 174, 175, 176, 177, 1
1454
1455 G4double HfA[33] = {153.964, 154.963, 155.9
1456 161.947, 162.947, 163.944, 164.945, 165.9
1457 172.941, 173.94, 174.942, 175.941, 176.94
1458 182.954, 183.955, 184.959, 185.961};
1459
1460 G4double HfS[33] = {75, 64, 38, 32, 11, 32,
1461 11, 3, 3, 29, 27, 27, 27, 27, 28, 7, 30,
1462
1463 G4double HfW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1464 18.6, 27.28, 13.62, 35.08, 0, 0, 0, 0, 0,
1465
1466 AddElement("Hf", 72, 33, *HfN, *HfA, *HfS,
1467
1468 // Z = 73 ---------------------------------
1469 G4int TaN[33] = {156, 157, 158, 159, 160, 1
1470 172, 173, 174, 175, 176, 177, 178, 179, 1
1471
1472 G4double TaA[33] = {155.972, 156.968, 157.9
1473 163.954, 164.951, 165.95, 166.948, 167.94
1474 173.944, 174.944, 175.945, 176.944, 177.9
1475 183.954, 184.956, 185.959, 186.96, 187.96
1476
1477 G4double TaS[33] = {64, 64, 54, 13, 33, 60,
1478 11, 11, 4, 11, 6, 3, 3, 3, 3, 28, 15, 60,
1479
1480 G4double TaW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1481 99.988, 0, 0, 0, 0, 0, 0, 0};
1482
1483 AddElement("Ta", 73, 33, *TaN, *TaA, *TaS,
1484
1485 // Z = 74 ---------------------------------
1486 G4int WN[33] = {158, 159, 160, 161, 162, 16
1487 174, 175, 176, 177, 178, 179, 180, 181, 1
1488
1489 G4double WA[33] = {157.974, 158.972, 159.96
1490 165.955, 166.955, 167.952, 168.952, 169.9
1491 175.946, 176.947, 177.946, 178.947, 179.9
1492 185.954, 186.957, 187.958, 188.962, 189.9
1493
1494 G4double WS[33] = {75, 64, 38, 33, 11, 33,
1495 11, 17, 5, 6, 3, 29, 29, 30, 3, 3, 4, 21,
1496
1497 G4double WW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1498 26.5, 14.31, 30.64, 0, 28.43, 0, 0, 0, 0}
1499
1500 AddElement("W", 74, 33, *WN, *WA, *WS, *WW)
1501
1502 // Z = 75 ---------------------------------
1503 G4int ReN[33] = {160, 161, 162, 163, 164, 1
1504 176, 177, 178, 179, 180, 181, 182, 183, 1
1505
1506 G4double ReA[33] = {159.981, 160.978, 161.9
1507 167.962, 168.959, 169.958, 170.956, 171.9
1508 177.951, 178.95, 179.951, 180.95, 181.951
1509 188.959, 189.962, 190.963, 191.966};
1510
1511 G4double ReS[33] = {64, 64, 55, 12, 33, 80,
1512 60, 40, 15, 11, 9, 6, 30, 3, 30, 30, 9, 2
1513
1514 G4double ReW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1515 37.4, 0, 62.6, 0, 0, 0, 0, 0};
1516
1517 AddElement("Re", 75, 33, *ReN, *ReA, *ReS,
1518
1519 // Z = 76 ---------------------------------
1520 G4int OsN[35] = {162, 163, 164, 165, 166, 1
1521 178, 179, 180, 181, 182, 183, 184, 185, 1
1522
1523 G4double OsA[35] = {161.984, 162.982, 163.9
1524 169.964, 170.963, 171.96, 172.96, 173.957
1525 180.953, 181.952, 182.953, 183.952, 184.9
1526 190.961, 191.961, 192.964, 193.965, 194.9
1527
1528 G4double OsS[35] = {75, 64, 38, 33, 11, 33,
1529 22, 27, 11, 3, 3, 3, 30, 30, 30, 3, 3, 4,
1530
1531 G4double OsW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1532 1.59, 1.96, 13.24, 16.15, 26.26, 0, 40.78
1533
1534 AddElement("Os", 76, 35, *OsN, *OsA, *OsS,
1535
1536 // Z = 77 ---------------------------------
1537 G4int IrN[35] = {165, 166, 167, 168, 169, 1
1538 181, 182, 183, 184, 185, 186, 187, 188, 1
1539
1540 G4double IrA[35] = {164.988, 165.986, 166.9
1541 172.968, 173.967, 174.964, 175.964, 176.9
1542 182.957, 183.957, 184.957, 185.958, 186.9
1543 192.963, 193.965, 194.966, 195.968, 196.9
1544
1545 G4double IrS[35] = {43, 55, 11, 35, 10, 16,
1546 29, 21, 22, 7, 8, 14, 21, 3, 3, 3, 3, 3,
1547
1548 G4double IrW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1549 37.3, 0, 62.7, 0, 0, 0, 0, 0, 0};
1550
1551 AddElement("Ir", 77, 35, *IrN, *IrA, *IrS,
1552
1553 // Z = 78 ---------------------------------
1554 G4int PtN[35] = {168, 169, 170, 171, 172, 1
1555 184, 185, 186, 187, 188, 189, 190, 191, 1
1556
1557 G4double PtA[35] = {167.988, 168.986, 169.9
1558 175.969, 176.968, 177.966, 178.965, 179.9
1559 185.959, 186.961, 187.959, 188.961, 189.9
1560 195.965, 196.967, 197.968, 198.971, 199.9
1561
1562 G4double PtS[35] = {38, 34, 11, 34, 40, 11,
1563 20, 6, 12, 7, 5, 4, 3, 3, 3, 3, 3, 4, 5,
1564
1565 G4double PtW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1566 0.782, 0, 32.967, 33.832, 25.242, 0, 7.16
1567
1568 AddElement("Pt", 78, 35, *PtN, *PtA, *PtS,
1569
1570 // Z = 79 ---------------------------------
1571 G4int AuN[35] = {171, 172, 173, 174, 175, 1
1572 187, 188, 189, 190, 191, 192, 193, 194, 1
1573
1574 G4double AuA[35] = {170.992, 171.99, 172.98
1575 179.972, 180.97, 181.97, 182.968, 183.967
1576 190.964, 191.965, 192.964, 193.965, 194.9
1577 200.972, 201.974, 202.975, 203.978, 204.9
1578
1579 G4double AuS[35] = {27, 36, 11, 16, 26, 43,
1580 17, 50, 17, 10, 12, 3, 4, 3, 3, 3, 60, 5,
1581
1582 G4double AuW[35] = {0, 0, 0, 0, 0, 0, 0, 0,
1583 100, 0, 0, 0, 0, 0, 0, 0, 0};
1584
1585 AddElement("Au", 79, 35, *AuN, *AuA, *AuS,
1586
1587 // Z = 80 ---------------------------------
1588 G4int HgN[34] = {175, 176, 177, 178, 179, 1
1589 191, 192, 193, 194, 195, 196, 197, 198, 1
1590
1591 G4double HgA[34] = {174.991, 175.987, 176.9
1592 182.975, 183.972, 184.972, 185.969, 186.9
1593 192.967, 193.965, 194.967, 195.966, 196.9
1594 202.973, 203.973, 204.976, 205.977, 206.9
1595
1596 G4double HgS[34] = {34, 40, 12, 16, 33, 21,
1597 25, 50, 4, 5, 3, 3, 3, 3, 3, 3, 3, 5, 22,
1598
1599 G4double HgW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1600 16.87, 23.1, 13.18, 29.86, 0, 6.87, 0, 0,
1601
1602 AddElement("Hg", 80, 34, *HgN, *HgA, *HgS,
1603
1604 // Z = 81 ---------------------------------
1605 G4int TlN[34] = {177, 178, 179, 180, 181, 1
1606 193, 194, 195, 196, 197, 198, 199, 200, 2
1607
1608 G4double TlA[34] = {176.997, 177.995, 178.9
1609 184.979, 185.979, 186.976, 187.976, 188.9
1610 194.97, 195.971, 196.97, 197.97, 198.97,
1611 205.976, 206.977, 207.982, 208.985, 209.9
1612
1613 G4double TlS[34] = {24, 23, 15, 48, 41, 43,
1614 15, 30, 90, 11, 7, 16, 16, 3, 3, 3, 3, 6,
1615
1616 G4double TlW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1617 29.524, 0, 70.476, 0, 0, 0, 0, 0};
1618
1619 AddElement("Tl", 81, 34, *TlN, *TlA, *TlS,
1620
1621 // Z = 82 ---------------------------------
1622 G4int PbN[34] = {181, 182, 183, 184, 185, 1
1623 197, 198, 199, 200, 201, 202, 203, 204, 2
1624
1625 G4double PbA[34] = {180.997, 181.993, 182.9
1626 188.981, 189.978, 190.978, 191.976, 192.9
1627 198.973, 199.972, 200.973, 201.972, 202.9
1628 208.981, 209.984, 210.989, 211.992, 212.9
1629
1630 G4double PbS[34] = {17, 18, 33, 21, 33, 50,
1631 14, 30, 11, 7, 3, 3, 3, 3, 3, 3, 3, 3, 29
1632
1633 G4double PbW[34] = {0, 0, 0, 0, 0, 0, 0, 0,
1634 24.1, 22.1, 52.4, 0, 0, 0, 0, 0, 0};
1635
1636 AddElement("Pb", 82, 34, *PbN, *PbA, *PbS,
1637
1638 // Z = 83 ---------------------------------
1639 G4int BiN[32] = {185, 186, 187, 188, 189, 1
1640 201, 202, 203, 204, 205, 206, 207, 208, 2
1641
1642 G4double BiA[32] = {184.998, 185.996, 186.9
1643 192.983, 193.983, 194.981, 195.981, 196.9
1644 202.977, 203.978, 204.977, 205.978, 206.9
1645 213.999, 215.002, 216.006};
1646
1647 G4double BiS[32] = {24, 48, 41, 32, 43, 39,
1648 28, 8, 9, 4, 4, 3, 3, 6, 3, 8, 12, 10, 11
1649
1650 G4double BiW[32] = {0, 0, 0, 0, 0, 0, 0, 0,
1651 0, 0, 0, 0, 0, 0, 0};
1652
1653 AddElement("Bi", 83, 32, *BiN, *BiA, *BiS,
1654
1655 // Z = 84 ---------------------------------
1656 G4int PoN[29] = {190, 191, 192, 193, 194, 1
1657 206, 207, 208, 209, 210, 211, 212, 213, 2
1658
1659 G4double PoA[29] = {189.995, 190.995, 191.9
1660 197.983, 198.984, 199.982, 200.982, 201.9
1661 208.982, 209.983, 210.987, 211.989, 212.9
1662
1663 G4double PoS[29] = {51, 32, 22, 30, 22, 24,
1664 3, 3, 3, 4, 3, 3, 29, 11, 27};
1665
1666 G4double PoW[29] = {
1667 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1668
1669 AddElement("Po", 84, 29, *PoN, *PoA, *PoS,
1670
1671 // Z = 85 ---------------------------------
1672 G4int AtN[31] = {193, 194, 195, 196, 197, 1
1673 209, 210, 211, 212, 213, 214, 215, 216, 2
1674
1675 G4double AtA[31] = {193.000, 193.999, 194.9
1676 200.988, 201.988, 202.987, 203.987, 204.9
1677 210.987, 211.991, 212.993, 213.996, 214.9
1678 221.018, 222.022, 223.025};
1679
1680 G4double AtS[31] = {43, 43, 43, 25, 38, 46,
1681 5, 8, 5, 8, 13, 90, 12, 32, 32, 43};
1682
1683 G4double AtW[31] = {
1684 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1685
1686 AddElement("At", 85, 31, *AtN, *AtA, *AtS,
1687
1688 // Z = 86 ---------------------------------
1689 G4int RnN[33] = {196, 197, 198, 199, 200, 2
1690 212, 213, 214, 215, 216, 217, 218, 219, 2
1691
1692 G4double RnA[33] = {196.002, 197.002, 197.9
1693 203.991, 204.992, 205.99, 206.991, 207.99
1694 214.999, 216.000, 217.004, 218.006, 219.0
1695 225.028, 226.031, 227.035, 228.038};
1696
1697 G4double RnS[33] = {22, 30, 22, 24, 19, 21,
1698 5, 4, 3, 29, 11, 27, 32, 32, 32, 43, 45,
1699
1700 G4double RnW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1701 100, 0, 0, 0, 0, 0, 0};
1702
1703 AddElement("Rn", 86, 33, *RnN, *RnA, *RnS,
1704
1705 // Z = 87 ---------------------------------
1706 G4int FrN[33] = {200, 201, 202, 203, 204, 2
1707 216, 217, 218, 219, 220, 221, 222, 223, 2
1708
1709 G4double FrA[33] = {200.006, 201.004, 202.0
1710 207.997, 208.996, 209.996, 210.996, 211.9
1711 218.008, 219.009, 220.012, 221.014, 222.0
1712 228.036, 229.038, 230.043, 231.045, 232.0
1713
1714 G4double FrS[33] = {25, 38, 46, 25, 22, 25,
1715 8, 23, 29, 50, 11, 10, 10, 22, 39, 48, 56
1716
1717 G4double FrW[33] = {0, 0, 0, 0, 0, 0, 0, 0,
1718 0, 0, 0, 0, 0, 0, 0, 0};
1719
1720 AddElement("Fr", 87, 33, *FrN, *FrA, *FrS,
1721
1722 // Z = 88 ---------------------------------
1723 G4int RaN[32] = {203, 204, 205, 206, 207, 2
1724 219, 220, 221, 222, 223, 224, 225, 226, 2
1725
1726 G4double RaA[32] = {203.009, 204.006, 205.0
1727 212, 213, 214, 215.003, 216.004, 217.006,
1728 224.02, 225.024, 226.025, 227.029, 228.03
1729 234.051};
1730
1731 G4double RaS[32] = {25, 19, 22, 16, 45, 15,
1732 29, 3, 27, 27, 27, 70, 40, 32, 39, 50, 58
1733
1734 G4double RaW[32] = {0, 0, 0, 0, 0, 0, 0, 0,
1735 0, 0, 0, 0, 0, 0, 0};
1736
1737 AddElement("Ra", 88, 32, *RaN, *RaA, *RaS,
1738
1739 // Z = 89 ---------------------------------
1740 G4int AcN[30] = {207, 208, 209, 210, 211, 2
1741 223, 224, 225, 226, 227, 228, 229, 230, 2
1742
1743 G4double AcA[30] = {207.012, 208.011, 209.0
1744 215.006, 216.009, 217.009, 218.012, 219.0
1745 225.023, 226.026, 227.028, 228.031, 229.0
1746 235.051, 236.055};
1747
1748 G4double AcS[30] = {25, 22, 26, 20, 14, 10,
1749 28, 50, 11, 11, 11, 32, 43, 45, 54};
1750
1751 G4double AcW[30] = {
1752 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1753
1754 AddElement("Ac", 89, 30, *AcN, *AcA, *AcS,
1755
1756 // Z = 90 ---------------------------------
1757 G4int ThN[29] = {210, 211, 212, 213, 214, 2
1758 226, 227, 228, 229, 230, 231, 232, 233, 2
1759
1760 G4double ThA[29] = {210.015, 211.015, 212.0
1761 218.013, 219.016, 220.016, 221.018, 222.0
1762 228.029, 229.032, 230.033, 231.036, 232.0
1763 238.056};
1764
1765 G4double ThS[29] = {17, 45, 15, 14, 10, 70,
1766 22, 22, 22, 22, 4, 50, 32, 39, 39};
1767
1768 G4double ThW[29] = {
1769 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1770
1771 AddElement("Th", 90, 29, *ThN, *ThA, *ThS,
1772
1773 // Z = 91 ---------------------------------
1774 G4int PaN[28] = {213, 214, 215, 216, 217, 2
1775 229, 230, 231, 232, 233, 234, 235, 236, 2
1776
1777 G4double PaA[28] = {213.021, 214.021, 215.0
1778 222.024, 223.024, 224.026, 225.026, 226.0
1779 232.039, 233.04, 234.043, 235.045, 236.04
1780
1781 G4double PaS[28] = {27, 21, 15, 11, 80, 80,
1782 24, 5, 50, 21, 11, 60, 32, 32};
1783
1784 G4double PaW[28] = {
1785 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1786
1787 AddElement("Pa", 91, 28, *PaN, *PaA, *PaS,
1788
1789 // Z = 92 ---------------------------------
1790 G4int UN[25] = {218, 219, 220, 221, 222, 22
1791 234, 235, 236, 237, 238, 239, 240, 241, 2
1792
1793 G4double UA[25] = {218.023, 219.025, 220.02
1794 226.029, 227.031, 228.031, 229.033, 230.0
1795 236.046, 237.049, 238.051, 239.054, 240.0
1796
1797 G4double US[25] = {
1798 10, 90, 22, 11, 11, 80, 27, 50, 20, 18, 1
1799
1800 G4double UW[25] = {
1801 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1802
1803 AddElement("U", 92, 25, *UN, *UA, *US, *UW)
1804
1805 // Z = 93 ---------------------------------
1806 G4int NpN[20] = {225, 226, 227, 228, 229, 2
1807 241, 242, 243, 244};
1808
1809 G4double NpA[20] = {225.034, 226.035, 227.0
1810 233.041, 234.043, 235.044, 236.047, 237.0
1811 243.064, 244.068};
1812
1813 G4double NpS[20] = {
1814 80, 10, 80, 21, 90, 60, 50, 11, 60, 9, 23
1815
1816 G4double NpW[20] = {0, 0, 0, 0, 0, 0, 0, 0,
1817
1818 AddElement("Np", 93, 20, *NpN, *NpA, *NpS,
1819
1820 // Z = 94 ---------------------------------
1821 G4int PuN[20] = {228, 229, 230, 231, 232, 2
1822 244, 245, 246, 247};
1823
1824 G4double PuA[20] = {228.039, 229.04, 230.04
1825 237.048, 238.05, 239.052, 240.054, 241.05
1826 247.074};
1827
1828 G4double PuS[20] = {30, 80, 26, 11, 20, 50,
1829
1830 G4double PuW[20] = {0, 0, 0, 0, 0, 0, 0, 0,
1831
1832 AddElement("Pu", 94, 20, *PuN, *PuA, *PuS,
1833
1834 // Z = 95 ---------------------------------
1835 G4int AmN[19] = {
1836 231, 232, 233, 234, 235, 236, 237, 238, 2
1837
1838 G4double AmA[19] = {231.046, 232.047, 233.0
1839 240.055, 241.057, 242.06, 243.061, 244.06
1840
1841 G4double AmS[19] = {32, 32, 23, 22, 22, 11,
1842
1843 G4double AmW[19] = {0, 0, 0, 0, 0, 0, 0, 0,
1844
1845 AddElement("Am", 95, 19, *AmN, *AmA, *AmS,
1846
1847 // Z = 96 ---------------------------------
1848 G4int CmN[20] = {233, 234, 235, 236, 237, 2
1849 249, 250, 251, 252};
1850
1851 G4double CmA[20] = {233.051, 234.05, 235.05
1852 241.058, 242.059, 243.061, 244.063, 245.0
1853 251.082, 252.085};
1854
1855 G4double CmS[20] = {43, 32, 24, 22, 22, 40,
1856
1857 G4double CmW[20] = {0, 0, 0, 0, 0, 0, 0, 0,
1858
1859 AddElement("Cm", 96, 20, *CmN, *CmA, *CmS,
1860
1861 // Z = 97 ---------------------------------
1862 G4int BkN[20] = {235, 236, 237, 238, 239, 2
1863 251, 252, 253, 254};
1864
1865 G4double BkA[20] = {235.057, 236.057, 237.0
1866 244.065, 245.066, 246.069, 247.07, 248.07
1867 254.091};
1868
1869 G4double BkS[20] = {43, 43, 32, 31, 31, 16,
1870
1871 G4double BkW[20] = {0, 0, 0, 0, 0, 0, 0, 0,
1872
1873 AddElement("Bk", 97, 20, *BkN, *BkA, *BkS,
1874
1875 // Z = 98 ---------------------------------
1876 G4int CfN[20] = {237, 238, 239, 240, 241, 2
1877 253, 254, 255, 256};
1878
1879 G4double CfA[20] = {237.062, 238.061, 239.0
1880 245.068, 246.069, 247.071, 248.072, 249.0
1881 255.091, 256.093};
1882
1883 G4double CfS[20] = {54, 43, 25, 22, 27, 40,
1884
1885 G4double CfW[20] = {0, 0, 0, 0, 0, 0, 0, 0,
1886
1887 AddElement("Cf", 98, 20, *CfN, *CfA, *CfS,
1888
1889 // Z = 99 ---------------------------------
1890 G4int EsN[18] = {
1891 240, 241, 242, 243, 244, 245, 246, 247, 2
1892
1893 G4double EsA[18] = {240.069, 241.069, 242.0
1894 249.076, 250.079, 251.08, 252.083, 253.08
1895
1896 G4double EsS[18] = {43, 32, 35, 31, 20, 22,
1897
1898 G4double EsW[18] = {0, 0, 0, 0, 0, 0, 0, 0,
1899
1900 AddElement("Es", 99, 18, *EsN, *EsA, *EsS,
1901
1902 // Z = 100 --------------------------------
1903 G4int FmN[18] = {
1904 242, 243, 244, 245, 246, 247, 248, 249, 2
1905
1906 G4double FmA[18] = {242.073, 243.075, 244.0
1907 250.08, 251.082, 252.082, 253.085, 254.08
1908
1909 G4double FmS[18] = {43, 25, 31, 30, 40, 16,
1910
1911 G4double FmW[18] = {0, 0, 0, 0, 0, 0, 0, 0,
1912
1913 AddElement("Fm", 100, 18, *FmN, *FmA, *FmS,
1914
1915 // Z = 101 --------------------------------
1916 G4int MdN[16] = {245, 246, 247, 248, 249, 2
1917
1918 G4double MdA[16] = {245.081, 246.082, 247.0
1919 253.087, 254.09, 255.091, 256.094, 257.09
1920
1921 G4double MdS[16] = {40, 41, 39, 26, 24, 32,
1922
1923 G4double MdW[16] = {0, 0, 0, 0, 0, 0, 0, 0,
1924
1925 AddElement("Md", 101, 16, *MdN, *MdA, *MdS,
1926
1927 // Z = 102 --------------------------------
1928 G4int NoN[14] = {249, 250, 251, 252, 253, 2
1929
1930 G4double NoA[14] = {249.088, 250.087, 251.0
1931 257.097, 258.098, 259.101, 260.103, 261.1
1932
1933 G4double NoS[14] = {37, 22, 19, 14, 26, 19,
1934
1935 G4double NoW[14] = {0, 0, 0, 0, 0, 0, 0, 0,
1936
1937 AddElement("No", 102, 14, *NoN, *NoA, *NoS,
1938
1939 // Z = 103 --------------------------------
1940 G4int LrN[13] = {251, 252, 253, 254, 255, 2
1941
1942 G4double LrA[13] = {251.094, 252.095, 253.0
1943 260.106, 261.107, 262.11, 263.111};
1944
1945 G4double LrS[13] = {32, 32, 24, 36, 22, 24,
1946
1947 G4double LrW[13] = {0, 0, 0, 0, 0, 0, 0, 0,
1948
1949 AddElement("Lr", 103, 13, *LrN, *LrA, *LrS,
1950
1951 // Z = 104 --------------------------------
1952 G4int RfN[12] = {253, 254, 255, 256, 257, 2
1953
1954 G4double RfA[12] = {253.101, 254.1, 255.101
1955 262.11, 263.113, 264.114};
1956
1957 G4double RfS[12] = {49, 31, 22, 29, 29, 22,
1958
1959 G4double RfW[12] = {0, 0, 0, 0, 0, 0, 0, 0,
1960
1961 AddElement("Rf", 104, 12, *RfN, *RfA, *RfS,
1962
1963 // Z = 105 --------------------------------
1964 G4int DbN[11] = {255, 256, 257, 258, 259, 2
1965
1966 G4double DbA[11] = {255.107, 256.108, 257.1
1967 263.115, 264.117, 265.119};
1968
1969 G4double DbS[11] = {45, 39, 25, 37, 31, 25,
1970
1971 G4double DbW[11] = {0, 0, 0, 0, 0, 0, 0, 10
1972
1973 AddElement("Db", 105, 11, *DbN, *DbA, *DbS,
1974
1975 // Z = 106 --------------------------------
1976 G4int SgN[9] = {258, 259, 260, 261, 262, 26
1977
1978 G4double SgA[9] = {
1979 258.113, 259.115, 260.114, 261.116, 262.1
1980
1981 G4double SgS[9] = {45, 23, 40, 30, 30, 13,
1982
1983 G4double SgW[9] = {0, 0, 0, 0, 0, 0, 0, 0,
1984
1985 AddElement("Sg", 106, 9, *SgN, *SgA, *SgS,
1986
1987 // Z = 107 --------------------------------
1988 G4int BhN[8] = {260, 261, 262, 263, 264, 26
1989
1990 G4double BhA[8] = {260.122, 261.122, 262.12
1991
1992 G4double BhS[8] = {66, 26, 40, 45, 30, 41,
1993
1994 G4double BhW[8] = {0, 0, 0, 0, 100, 0, 0, 0
1995
1996 AddElement("Bh", 107, 8, *BhN, *BhA, *BhS,
1997
1998 if (0 < verbose) {
1999 G4cout << "G4NistElementBuilder: " << max
2000 << " Isotopes" << G4endl;
2001 }
2002 }
2003