Geant4-11
G4NistElementBuilder.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// -------------------------------------------------------------------
28//
29// GEANT4 Class file
30//
31// File name: G4NistElementBuilder
32//
33// Author: Vladimir Ivanchenko
34//
35// Creation date: 23.12.2004
36//
37// Modifications:
38// 02.05.2006 Subtract mass of atomic electrons from NIST mass (VI)
39// 11.05.2006 Do not subtract mass of atomic electrons from NIST mass (VI)
40// 17.10.2006 Add natiral abundances flag to element and
41// use G4 units for isotope mass vector (VI)
42// 10.05.2007 Add protection agains Z>101 (VI)
43// 26.07.2007 Create one and only one Nist element with given Z and
44// allow users to create there own elements with the same Z (VI)
45//
46// -------------------------------------------------------------------
47//
48// Class Description:
49//
50// Element data from the NIST DB on Atomic Weights and Isotope Compositions
51// http://physics.nist.gov/PhysRefData/Compositions/index.html
52//
53//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55
56#include <sstream>
57
58#include "G4NistElementBuilder.hh"
59#include "G4PhysicalConstants.hh"
60#include "G4SystemOfUnits.hh"
61
62#ifdef G4MULTITHREADED
63G4Mutex G4NistElementBuilder::nistElementMutex = G4MUTEX_INITIALIZER;
64#endif
65
66//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
67
68G4NistElementBuilder::G4NistElementBuilder(G4int vb):
69 verbose(vb)
70{
71 nFirstIsotope[0] = 0;
72 nIsotopes[0] = 0;
73 relAbundance[0] = 0.0;
74 Initialise();
75 for(G4int i=0; i<maxNumElements; ++i) {elmIndex[i] = -1;}
76}
77
78//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
79
80G4NistElementBuilder::~G4NistElementBuilder()
81{}
82
83//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84
85G4int G4NistElementBuilder::GetZ(const G4String& name) const
86{
87 G4int Z = maxNumElements;
88 // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
89 do {--Z;} while( Z>0 && elmSymbol[Z] != name);
90 return Z;
91}
92
93//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
94
95G4double G4NistElementBuilder::GetAtomicMassAmu(const G4String& name) const
96{
97 G4int Z = maxNumElements;
98 // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
99 do {--Z;} while( Z>0 && elmSymbol[Z] != name);
100 return GetAtomicMassAmu(Z);
101}
102
103//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
104
105G4Element* G4NistElementBuilder::FindOrBuildElement(G4int Z, G4bool)
106{
107 G4Element* anElement = FindElement(Z);
108 if(anElement == nullptr && Z > 0 && Z < maxNumElements) {
109#ifdef G4MULTITHREADED
110 G4MUTEXLOCK(&nistElementMutex);
111#endif
112 if(elmIndex[Z] < 0) { anElement = BuildElement(Z); }
113 if(anElement) { elmIndex[Z] = anElement->GetIndex(); }
114#ifdef G4MULTITHREADED
115 G4MUTEXUNLOCK(&nistElementMutex);
116#endif
117 }
118 return anElement;
119}
120
121
122//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
123
124G4Element*
125G4NistElementBuilder::FindOrBuildElement(const G4String& symb, G4bool)
126{
127 G4Element* elm = nullptr;
128 const G4ElementTable* theElementTable = G4Element::GetElementTable();
129 size_t nelm = theElementTable->size();
130 for(size_t i=0; i<nelm; ++i) {
131 if(symb == ((*theElementTable)[i])->GetSymbol()) {
132 elm = (*theElementTable)[i];
133 break;
134 }
135 }
136 if(nullptr == elm) {
137 for(G4int Z = 1; Z<maxNumElements; ++Z) {
138 if(symb == elmSymbol[Z]) {
139#ifdef G4MULTITHREADED
140 G4MUTEXLOCK(&nistElementMutex);
141#endif
142 if(elmIndex[Z] < 0) { elm = BuildElement(Z); }
143 if(elm) { elmIndex[Z] = elm->GetIndex(); }
144#ifdef G4MULTITHREADED
145 G4MUTEXUNLOCK(&nistElementMutex);
146#endif
147 break;
148 }
149 }
150 }
151 return elm;
152}
153
154//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
155
156G4Element* G4NistElementBuilder::BuildElement(G4int Z)
157{
158 G4Element* theElement = nullptr;
159 if(Z<1 || Z>=maxNumElements) { return theElement; }
160 G4double Aeff = atomicMass[Z];
161 if (verbose > 1) {
162 G4cout << "G4NistElementBuilder: Build Element <" << elmSymbol[Z]
163 << "> Z= " << Z
164 << " Aeff= " << Aeff;
165 G4cout << " with natural isotope composition" << G4endl;
166 }
167
168 //build Element with its Isotopes
169 //
170 G4int nc = nIsotopes[Z];
171 G4int n0 = nFirstIsotope[Z];
172 G4int idx = idxIsotopes[Z];
173 std::vector<G4Isotope*> iso;
174 G4Isotope* ist;
175 for (G4int i=0; i<nc; ++i) {
176 if (relAbundance[idx + i] > 0.0) {
177 std::ostringstream os;
178 os << elmSymbol[Z] << n0 + i;
179 ist = new G4Isotope(os.str(), Z, n0 + i,
180 GetAtomicMass(Z, n0 + i)*g/(mole*amu_c2));
181 /*
182 G4cout << " Z= " << Z << " N= " << n0 + i
183 << " miso(amu)= " << GetIsotopeMass(Z, n0 + i)/amu_c2
184 << " matom(amu)= " << GetAtomicMass(Z, n0 + i)/amu_c2 << G4endl;
185 */
186 iso.push_back(ist);
187 }
188 }
189 G4int ni = iso.size();
190 G4double w;
191 theElement = new G4Element(elmSymbol[Z],elmSymbol[Z],ni);
192 for(G4int j=0; j<ni; ++j) {
193 w = relAbundance[idx + (iso[j])->GetN() - n0];
194 ist = iso[j];
195 theElement->AddIsotope(ist, w);
196 }
197 theElement->SetNaturalAbundanceFlag(true);
198 return theElement;
199}
200
201//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
202
203void G4NistElementBuilder::PrintElement(G4int Z) const
204{
205 G4int imin = Z;
206 G4int imax = Z+1;
207 if (Z == 0) {
208 imin = 1;
209 imax = maxNumElements;
210 }
211 if(imax > maxNumElements) { imax = maxNumElements; }
212
213 for(G4int i=imin; i<imax; ++i) {
214 G4int nc = nIsotopes[i];
215 G4cout << "Nist Element: <" << elmSymbol[i]
216 << "> Z= " << i
217 << " Aeff(amu)= " << atomicMass[i] << " "
218 << nc << " isotopes:"
219 << G4endl;
220 G4int j;
221 G4int idx = idxIsotopes[i];
222 G4int n0 = nFirstIsotope[i];
223 G4cout << " N: ";
224 for(j=0; j<nc; ++j) {G4cout << n0 + j << " ";}
225 G4cout << G4endl;
226 G4cout << " mass(amu): ";
227 for(j=0; j<nc; ++j) {G4cout << GetAtomicMass(i, n0 + j) << " ";}
228 G4cout << G4endl;
229 G4cout << " abundance: ";
230 for(j=0; j<nc; ++j) {G4cout << relAbundance[idx + j] << " ";}
231 G4cout << G4endl;
232 }
233}
234
235//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
236
237void G4NistElementBuilder::AddElement(const G4String& name, G4int Z, G4int nc,
238 const G4int& N, const G4double& A,
239 const G4double& sA, const G4double& W)
240{
241 if (verbose > 1) {
242 G4cout << "AddElement " << name << " Z= " << Z << " nc= " << nc << G4endl;
243 }
244 if (Z >= maxNumElements) {
245 G4cout << "G4NistElementBuilder::AddElement: Warning! Z= " << Z
246 << " is too big" << G4endl;
247 return;
248 }
249
250 if (index+nc >= maxAbundance) {
251 G4cout << "G4NistElementBuilder::AddElement: ERROR index= " << index
252 << " + nc= " << nc
253 << " is above array size " << maxAbundance << G4endl;
254 return;
255 }
256
257 elmSymbol[Z] = 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*electron_mass_c2 + bindingEnergy[Z];
270 sigMass[index] = (&sA)[i]*amu_c2;
271 relAbundance[index] = www;
272
273 // computation of mean atomic mass of the element in atomic units
274 atomicMass[Z] += www*(&A)[i];
275 ww += www;
276 ++index;
277 }
278
279 if(ww != 1.0) {
280 G4int idx = idxIsotopes[Z];
281 atomicMass[Z] /= ww;
282 for(G4int j=0; j<nc; ++j) {relAbundance[idx + j] /= ww;}
283 }
284
285 if (verbose > 1) { PrintElement(Z); }
286 else if(1 == verbose) {
287 G4cout << "Nist Element " << elmSymbol[Z]
288 << " Z= " << Z
289 << " Aeff(amu)= " << atomicMass[Z] << " "
290 << nIsotopes[Z] << " isotopes:"
291 << G4endl;
292 }
293}
294
295//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
296
297void G4NistElementBuilder::Initialise()
298{
299 // Parameterisation from D.Lunney,J.M.Pearson,C.Thibault,
300 // Rev.Mod.Phys. 75 (2003) 1021
301 bindingEnergy[0] = 0.0;
302 for(G4int i=1; i<maxNumElements; ++i) {
303 G4double Z = G4double(i);
304 bindingEnergy[i] = (14.4381*std::pow(Z,2.39) + 1.55468e-6*std::pow(Z,5.35))*eV;
305 }
306
307 // NIST data
308 index = 0;
309
310 // Z = 1 ---------------------------------------------------------------------
311 G4int HN[6] =
312 {1, 2, 3, 4, 5, 6};
313
314 G4double HA[6] =
315 {1.00783, 2.0141, 3.01605, 4.02783, 5.03954, 6.04494};
316
317 // Garantee consistence with G4 masses
318 HA[0] = (proton_mass_c2 + electron_mass_c2 - bindingEnergy[1])/amu_c2;
319 HA[1] = (1.875613*GeV + electron_mass_c2 - bindingEnergy[1])/amu_c2;
320 HA[2] = (2.80925*GeV + electron_mass_c2 - bindingEnergy[1])/amu_c2;
321
322 G4double HS[6] =
323 {4, 4, 11, 12, 102, 28};
324
325 G4double HW[6] =
326 {99.9885, 0.0115, 0, 0, 0, 0};
327
328 AddElement("H", 1, 6, *HN , *HA , *HS , *HW);
329
330 // Z = 2 ---------------------------------------------------------------------
331 G4int HeN[8] =
332 {3, 4, 5, 6, 7, 8, 9, 10};
333
334 G4double HeA[8] =
335 {3.01603, 4.0026, 5.01222, 6.01889, 7.02803, 8.03392, 9.04382, 10.0524};
336
337 // Garantee consistence with G4 masses
338 HeA[0] = (2.80923*GeV + 2.0*electron_mass_c2 - bindingEnergy[2])/amu_c2;
339 HeA[1] = (3.727417*GeV + 2.0*electron_mass_c2 - bindingEnergy[2])/amu_c2;
340
341 G4double HeS[8] =
342 {9, 10, 50, 11, 30, 8, 70, 80};
343
344 G4double HeW[8] =
345 {0.000137, 99.9999, 0, 0, 0, 0, 0, 0};
346
347 AddElement("He", 2, 8, *HeN , *HeA , *HeS , *HeW);
348
349 // Z = 3 ---------------------------------------------------------------------
350 G4int LiN[9] =
351 {4, 5, 6, 7, 8, 9, 10, 11, 12};
352 G4double LiA[9] =
353 {4.02718, 5.01254, 6.01512, 7.016, 8.02249, 9.02679, 10.0355, 11.0438,
354 12.0538};
355
356 G4double LiS[9] =
357 {23, 50, 5, 5, 5, 21, 16, 29, 107};
358
359 G4double LiW[9] =
360 {0, 0, 7.59, 92.41, 0, 0, 0, 0, 0};
361
362 AddElement("Li", 3, 9, *LiN , *LiA , *LiS , *LiW);
363
364 // Z = 4 ---------------------------------------------------------------------
365 G4int BeN[10] =
366 {5, 6, 7, 8, 9, 10, 11, 12, 13, 14};
367
368 G4double BeA[10] =
369 {5.04079, 6.01973, 7.01693, 8.00531, 9.01218, 10.0135, 11.0217, 12.0269,
370 13.0361, 14.0428};
371
372 G4double BeS[10] =
373 {429, 6, 5, 4, 4, 4, 7, 16, 54, 12};
374
375 G4double BeW[10] =
376 {0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
377
378 AddElement("Be", 4, 10, *BeN , *BeA , *BeS , *BeW);
379
380 // Z = 5 ---------------------------------------------------------------------
381 G4int BN[13] =
382 {7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19};
383
384 G4double BA[13] =
385 {7.02992, 8.02461, 9.01333, 10.0129, 11.0093, 12.0144, 13.0178, 14.0254,
386 15.0311, 16.0398, 17.0469, 18.0562, 19.0637};
387
388 G4double BS[13] =
389 {80, 12, 11, 4, 5, 15, 12, 23, 24, 60, 15, 86, 43};
390
391 G4double BW[13] =
392 {0, 0, 0, 19.9, 80.1, 0, 0, 0, 0, 0, 0, 0, 0};
393
394 AddElement("B", 5, 13, *BN , *BA , *BS , *BW);
395
396 // Z = 6 ---------------------------------------------------------------------
397 G4int CN[15] =
398 {8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22};
399
400 G4double CA[15] =
401 {8.03768, 9.03104, 10.0169, 11.0114, 12. , 13.0034, 14.0032, 15.0106,
402 16.0147, 17.0226, 18.0268, 19.0353, 20.0403, 21.0493, 22.0565};
403
404 G4double CS[15] =
405 {25, 23, 4, 10, 0, 10, 4, 9, 4, 19, 30, 12, 22, 54, 97};
406
407 G4double CW[15] =
408 {0, 0, 0, 0, 98.93, 1.07, 0, 0, 0, 0, 0, 0, 0, 0, 0};
409
410 AddElement("C", 6, 15, *CN , *CA , *CS , *CW);
411
412 // Z = 7 ---------------------------------------------------------------------
413 G4int NN[15] =
414 {10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};
415
416 G4double NA[15] =
417 {10.0426, 11.0268, 12.0186, 13.0057, 14.0031, 15.0001, 16.0061, 17.0084,
418 18.0141, 19.017 , 20.0234, 21.0271, 22.0344, 23.0405, 24.0505};
419
420 G4double NS[15] =
421 {43, 19, 11, 29, 9, 9, 28, 16, 21, 18, 60, 10, 21, 76, 54};
422
423 G4double NW[15] =
424 {0, 0, 0, 0, 99.632, 0.368, 0, 0, 0, 0, 0, 0, 0, 0, 0};
425
426 AddElement("N", 7, 15, *NN , *NA , *NS , *NW);
427
428 // Z = 8 ---------------------------------------------------------------------
429 G4int ON[15] =
430 {12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26};
431
432 G4double OA[15] =
433 {12.0344, 13.0248, 14.0086, 15.0031, 15.9949, 16.9991, 17.9992, 19.0036,
434 20.0041, 21.0087, 22.01 , 23.0157, 24.0204, 25.0291, 26.0377};
435
436 G4double Os[15] =
437 {20, 10, 8, 5, 15, 22, 9, 3, 13, 13, 60, 11, 33, 40, 46};
438
439 G4double OW[15] =
440 {0, 0, 0, 0, 99.757, 0.038, 0.205, 0, 0, 0, 0, 0, 0, 0, 0};
441
442 AddElement("O", 8, 15, *ON , *OA , *Os , *OW);
443
444 // Z = 9 ---------------------------------------------------------------------
445 G4int FN[16] =
446 {14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29};
447
448 G4double FA[16] =
449 {14.0361, 15.018 , 16.0115, 17.0021, 18.0009, 18.9984, 20 , 20.9999,
450 22.003 , 23.0036, 24.0081, 25.0121, 26.0196, 27.0269, 28.0357, 29.0433};
451
452 G4double FS[16] =
453 {43, 14, 9, 27, 6, 7, 9, 19, 13, 90, 70, 80, 13, 45, 55, 62};
454
455 G4double FW[16] =
456 {0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
457
458 AddElement("F", 9, 16, *FN , *FA , *FS , *FW);
459
460 // Z = 10 --------------------------------------------------------------------
461 G4int NeN[17] =
462 {16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32};
463
464 G4double NeA[17] =
465 {16.0258, 17.0177, 18.0057, 19.0019, 19.9924, 20.9938, 21.9914, 22.9945,
466 23.9936, 24.9978, 26.0005, 27.0076, 28.0121, 29.0193, 30.0239, 31.0331,
467 32.0399};
468
469 G4double NeS[17] =
470 {22, 50, 16, 6, 20, 4, 23, 26, 11, 50, 60, 10, 12, 32, 88, 97, 94};
471
472 G4double NeW[17] =
473 {0, 0, 0, 0, 90.48, 0.27, 9.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
474
475 AddElement("Ne", 10, 17, *NeN , *NeA , *NeS , *NeW);
476
477 // Z = 11 --------------------------------------------------------------------
478 G4int NaN[18] =
479 {18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35};
480
481 G4double NaA[18] =
482 {18.0272, 19.0139, 20.0073, 20.9977, 21.9944, 22.9898, 23.991 , 24.99 ,
483 25.9926, 26.994 , 27.9989, 29.0028, 30.0092, 31.0136, 32.0196, 33.0274,
484 34.0349, 35.0442};
485
486 G4double NaS[18] =
487 {43, 13, 7, 8, 5, 23, 23, 13, 15, 40, 80, 10, 10, 18, 52, 160, 113, 166};
488
489 G4double NaW[18] =
490 {0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
491
492 AddElement("Na", 11, 18, *NaN , *NaA , *NaS , *NaW);
493
494 // Z = 12 --------------------------------------------------------------------
495 G4int MgN[18] =
496 {20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37};
497
498 G4double MgA[18] =
499 {20.0189, 21.0117, 21.9996, 22.9941, 23.985 , 24.9858, 25.9826, 26.9843,
500 27.9839, 28.9886, 29.9905, 30.9965, 31.9992, 33.0056, 34.0091, 35.0175,
501 36.0224, 37.0312};
502
503 G4double MgS[18] =
504 {29, 18, 15, 13, 20, 20, 21, 21, 22, 30, 70, 80, 10, 16, 28, 47, 97, 97};
505
506 G4double MgW[18] =
507 {0, 0, 0, 0, 78.99, 10, 11.01, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
508
509 AddElement("Mg", 12, 18, *MgN , *MgA , *MgS , *MgW);
510
511 // Z = 13 --------------------------------------------------------------------
512 G4int AlN[19] =
513 {21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39};
514
515 G4double AlA[19] =
516 {21.028 , 22.0195, 23.0073, 23.9999, 24.9904, 25.9869, 26.9815, 27.9819,
517 28.9804, 29.983 , 30.9839, 31.9881, 32.9909, 33.9969, 34.9999, 36.0063,
518 37.0103, 38.0169, 39.0219};
519
520 G4double AlS[19] =
521 {32, 10, 27, 4, 7, 21, 14, 15, 13, 15, 22, 90, 70, 10, 15, 29, 58, 60, 64};
522
523 G4double AlW[19] =
524 {0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
525
526 AddElement("Al", 13, 19, *AlN , *AlA , *AlS , *AlW);
527
528 // Z = 14 --------------------------------------------------------------------
529 G4int SiN[21] =
530 {22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40,
531 41, 42};
532
533 G4double SiA[21] =
534 {22.0345, 23.0255, 24.0115, 25.0041, 25.9923, 26.9867, 27.9769, 28.9765,
535 29.9738, 30.9754, 31.9741, 32.978 , 33.9786, 34.9846, 35.9867, 36.993 ,
536 37.996 , 39.0023, 40.0058, 41.0127, 42.0161};
537
538 G4double SiS[21] =
539 {22, 21, 21, 11, 3, 17, 20, 3, 5, 7, 23, 17, 15, 40, 11, 13, 29, 43, 54, 64,
540 75};
541
542 G4double SiW[21] =
543 {0, 0, 0, 0, 0, 0, 92.2297, 4.6832, 3.0872, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
544 0};
545
546 AddElement("Si", 14, 21, *SiN , *SiA , *SiS , *SiW);
547
548 // Z = 15 --------------------------------------------------------------------
549 G4int PN[23] =
550 {24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42,
551 43, 44, 45, 46};
552
553 G4double PA[23] =
554 {24.0343, 25.0203, 26.0118, 26.9992, 27.9923, 28.9818, 29.9783, 30.9738,
555 31.9739, 32.9717, 33.9736, 34.9733, 35.9783, 36.9796, 37.9845, 38.9864,
556 39.9911, 40.9948, 42.0001, 43.0033, 44.0099, 45.0151, 46.0238};
557
558 G4double PS[23] =
559 {54, 21, 21, 40, 4, 8, 4, 20, 20, 12, 5, 20, 14, 40, 15, 16, 21, 50, 54, 54,
560 75, 86, 97};
561
562 G4double PW[23] =
563 {0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
564
565 AddElement("P", 15, 23, *PN , *PA , *PS , *PW);
566
567 // Z = 16 --------------------------------------------------------------------
568 G4int SN[24] =
569 {26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44,
570 45, 46, 47, 48, 49};
571
572 G4double SA[24] =
573 {26.0279, 27.0188, 28.0044, 28.9966, 29.9849, 30.9796, 31.9721, 32.9715,
574 33.9679, 34.969 , 35.9671, 36.9711, 37.9712, 38.9751, 39.9755, 40.98 ,
575 41.9815, 42.9866, 43.9883, 44.9948, 45.9996, 47.0076, 48.013 , 49.022 };
576
577 G4double SS[24] =
578 {32, 22, 17, 50, 3, 16, 12, 12, 11, 10, 25, 27, 8, 50, 25, 23, 35, 90, 54,
579 64, 75, 86, 97, 107};
580
581 G4double SW[24] =
582 {0, 0, 0, 0, 0, 0, 94.93, 0.76, 4.29, 0, 0.02, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
583 0, 0, 0};
584
585 AddElement("S", 16, 24, *SN , *SA , *SS , *SW);
586
587 // Z = 17 --------------------------------------------------------------------
588 G4int ClN[24] =
589 {28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46,
590 47, 48, 49, 50, 51};
591
592 G4double ClA[24] =
593 {28.0285, 29.0141, 30.0048, 30.9924, 31.9857, 32.9775, 33.9738, 34.9689,
594 35.9683, 36.9659, 37.968 , 38.968 , 39.9704, 40.9706, 41.9732, 42.9742,
595 43.9785, 44.9797, 45.9841, 46.9879, 47.9948, 48.9999, 50.0077, 51.0135};
596
597 G4double ClS[24] =
598 {54, 21, 21, 50, 7, 6, 13, 4, 8, 5, 12, 19, 30, 70, 12, 17, 24, 70, 54,
599 64, 75, 86, 97, 107};
600
601 G4double ClW[24] =
602 {0, 0, 0, 0, 0, 0, 0, 75.78, 0, 24.22, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
603 0};
604
605 AddElement("Cl", 17, 24, *ClN , *ClA , *ClS , *ClW);
606
607 // Z = 18 --------------------------------------------------------------------
608 G4int ArN[24] =
609 {30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48,
610 49, 50, 51, 52, 53};
611
612 G4double ArA[24] =
613 {30.0216, 31.0121, 31.9977, 32.9899, 33.9803, 34.9753, 35.9675, 36.9668,
614 37.9627, 38.9643, 39.9624, 40.9645, 41.9631, 42.9657, 43.9654, 44.9681,
615 45.9681, 46.9722, 47.9751, 48.9822, 49.9859, 50.9932, 51.9982, 53.0062};
616
617 G4double ArS[24] =
618 {32, 22, 50, 30, 3, 8, 27, 3, 5, 5, 3, 7, 40, 80, 22, 60, 40, 11, 32, 54, 75,
619 75, 97, 107};
620
621 G4double ArW[24] =
622 {0, 0, 0, 0, 0, 0, 0.3365, 0, 0.0632, 0, 99.6003, 0, 0, 0, 0, 0, 0, 0, 0, 0,
623 0, 0, 0, 0};
624
625 AddElement("Ar", 18, 24, *ArN , *ArA , *ArS , *ArW);
626
627 // Z = 19 --------------------------------------------------------------------
628 G4int KN[24] =
629 {32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50,
630 51, 52, 53, 54, 55};
631
632 G4double KA[24] =
633 {32.0219, 33.0073, 33.9984, 34.988 , 35.9813, 36.9734, 37.9691, 38.9637,
634 39.964 , 40.9618, 41.9624, 42.9607, 43.9616, 44.9607, 45.962 , 46.9617,
635 47.9655, 48.9674, 49.9728, 50.9764, 51.9826, 52.9871, 53.994 , 54.9994};
636
637 G4double KS[24] =
638 {54, 21, 32, 21, 8, 29, 8, 3, 29, 28, 3, 10, 40, 11, 17, 9, 26, 80, 30, 54,
639 75, 75, 97, 107};
640
641 G4double KW[24] =
642 {0, 0, 0, 0, 0, 0, 0, 93.2581, 0.0117, 6.7302, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
643 0, 0, 0, 0};
644
645 AddElement("K", 19, 24, *KN , *KA , *KS , *KW);
646
647 // Z = 20 --------------------------------------------------------------------
648 G4int CaN[24] =
649 {34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52,
650 53, 54, 55, 56, 57};
651
652 G4double CaA[24] =
653 {34.0141, 35.0048, 35.9931, 36.9859, 37.9763, 38.9707, 39.9626, 40.9623,
654 41.9586, 42.9588, 43.9555, 44.9562, 45.9537, 46.9545, 47.9525, 48.9557,
655 49.9575, 50.9615, 51.9651, 52.9701, 53.9747, 54.9806, 55.9858, 56.9924};
656
657 G4double CaS[24] =
658 {32, 70, 40, 24, 5, 19, 3, 4, 4, 5, 9, 10, 25, 25, 4, 4, 10, 10, 50, 54, 75,
659 75, 97, 107};
660
661 G4double CaW[24] =
662 {0, 0, 0, 0, 0, 0, 96.941, 0, 0.647, 0.135, 2.086, 0, 0.004, 0, 0.187, 0, 0,
663 0, 0, 0, 0, 0, 0, 0};
664
665 AddElement("Ca", 20, 24, *CaN , *CaA , *CaS , *CaW);
666
667 // Z = 21 --------------------------------------------------------------------
668 G4int ScN[24] =
669 {36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54,
670 55, 56, 57, 58, 59};
671
672 G4double ScA[24] =
673 {36.0149, 37.0031, 37.9947, 38.9848, 39.978, 40.9693, 41.9655, 42.9612,
674 43.9594, 44.9559, 45.9552, 46.9524, 47.9522, 48.95 , 49.9522, 50.9536,
675 51.9567, 52.9592, 53.963 , 54.9674, 55.9727, 56.977, 57.9831, 58.988 };
676
677 G4double ScS[24] =
678 {54, 32, 32, 26, 4, 3, 4, 20, 19, 12, 12, 22, 6, 4, 17, 22, 25, 32, 50, 110,
679 75, 75, 86, 97};
680
681 G4double ScW[24] =
682 {0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
683
684 AddElement("Sc", 21, 24, *ScN , *ScA , *ScS , *ScW);
685
686 // Z = 22 --------------------------------------------------------------------
687 G4int TiN[24] =
688 {38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56,
689 57, 58, 59, 60, 61};
690
691 G4double TiA[24] =
692 {38.0098, 39.0013, 39.9905, 40.9831, 41.973 , 42.9685, 43.9597, 44.9581,
693 45.9526, 46.9518, 47.9479, 48.9479, 49.9448, 50.9466, 51.9469, 52.9497,
694 53.9509, 54.9551, 55.958 , 56.9629, 57.9661, 58.972 , 59.9756, 60.982 };
695
696 G4double TiS[24] =
697 {27, 11, 17, 40, 6, 7, 8, 13, 12, 10, 10, 10, 11, 14, 8, 11, 25, 26, 30, 100,
698 75, 75, 86, 97};
699
700 G4double TiW[24] =
701 {0, 0, 0, 0, 0, 0, 0, 0, 8.25, 7.44, 73.72, 5.41, 5.18, 0, 0, 0, 0, 0, 0, 0,
702 0, 0, 0, 0};
703
704 AddElement("Ti", 22, 24, *TiN , *TiA , *TiS , *TiW);
705
706 // Z = 23 --------------------------------------------------------------------
707 G4int VN[24] =
708 {40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58,
709 59, 60, 61, 62, 63};
710
711 G4double VA[24] =
712 {40.0111, 40.9997, 41.9912, 42.9806, 43.9744, 44.9658, 45.9602, 46.9549,
713 47.9523, 48.9485, 49.9472, 50.944 , 51.9448, 52.9443, 53.9464, 54.9472,
714 55.9504, 56.9524, 57.9567, 58.9593, 59.9645, 60.9674, 61.9731, 62.9768};
715
716 G4double VS[24] =
717 {54, 27, 21, 25, 90, 18, 16, 12, 28, 14, 14, 14, 14, 4, 16, 11, 26, 27, 28,
718 35, 60, 75, 75, 97};
719
720 G4double VW[24] =
721 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.25, 99.75, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
722 0};
723
724 AddElement("V", 23, 24, *VN , *VA , *VS , *VW);
725
726 // Z = 24 --------------------------------------------------------------------
727 G4int CrN[24] =
728 {42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60,
729 61, 62, 63, 64, 65};
730
731 G4double CrA[24] =
732 {42.0064, 42.9977, 43.9855, 44.9792, 45.9684, 46.9629, 47.954 , 48.9513,
733 49.946 , 50.9448, 51.9405, 52.9407, 53.9389, 54.9408, 55.9406, 56.9438,
734 57.9442, 58.9486, 59.9497, 60.9541, 61.9558, 62.9619, 63.9642, 64.9704};
735
736 G4double CrS[24] =
737 {32, 90, 14, 11, 22, 15, 8, 28, 14, 14, 15, 15, 15, 16, 10, 10, 26, 27, 28,
738 30, 40, 75, 75, 97};
739
740 G4double CrW[24] =
741 {0, 0, 0, 0, 0, 0, 0, 0, 4.345, 0, 83.789, 9.501, 2.365, 0, 0, 0, 0, 0, 0, 0,
742 0, 0, 0, 0};
743
744 AddElement("Cr", 24, 24, *CrN , *CrA , *CrS , *CrW);
745
746 // Z = 25 --------------------------------------------------------------------
747 G4int MnN[24] =
748 {44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62,
749 63, 64, 65, 66, 67};
750
751 G4double MnA[24] =
752 {44.0069, 44.9945, 45.9867, 46.9761, 47.9689, 48.9596, 49.9542, 50.9482,
753 51.9456, 52.9413, 53.9404, 54.938 , 55.9389, 56.9383, 57.94 , 58.9404,
754 59.9432, 60.9445, 61.948 , 62.9498, 63.9537, 64.9561, 65.9608, 66.9638};
755
756 G4double MnS[24] =
757 {54, 32, 12, 17, 80, 26, 15, 14, 25, 15, 18, 14, 15, 4, 30, 30, 29, 28, 28,
758 30, 35, 60, 75, 86};
759
760 G4double MnW[24] =
761 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
762
763 AddElement("Mn", 25, 24, *MnN , *MnA , *MnS , *MnW);
764
765 // Z = 26 --------------------------------------------------------------------
766 G4int FeN[25] =
767 {45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
768 64, 65, 66, 67, 68, 69};
769
770 G4double FeA[25] =
771 {45.0146, 46.0008, 46.9929, 47.9806, 48.9736, 49.963 , 50.9568, 51.9481,
772 52.9453, 53.9396, 54.9383, 55.9349, 56.9354, 57.9333, 58.9349, 59.9341,
773 60.9367, 61.9368, 62.9401, 63.9409, 64.9449, 65.946 , 66.95 , 67.9525,
774 68.9577};
775
776 G4double FeS[25] =
777 {43, 38, 28, 11, 17, 60, 16, 11, 23, 14, 14, 15, 15, 15, 15, 4, 22, 16, 20,
778 30, 30, 35, 50, 75, 86};
779
780 G4double FeW[25] =
781 {0, 0, 0, 0, 0, 0, 0, 0, 0, 5.845, 0, 91.754, 2.119, 0.282, 0, 0, 0, 0, 0, 0,
782 0, 0, 0, 0, 0};
783
784 AddElement("Fe", 26, 25, *FeN , *FeA , *FeS , *FeW);
785
786 // Z = 27 --------------------------------------------------------------------
787 G4int CoN[25] =
788 {48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66,
789 67, 68, 69, 70, 71, 72};
790
791 G4double CoA[25] =
792 {48.0018, 48.9897, 49.9815, 50.9707, 51.9636, 52.9542, 53.9485, 54.942,
793 55.9398, 56.9363, 57.9358, 58.9332, 59.9338, 60.9325, 61.9341, 62.9336,
794 63.9358, 64.9365, 65.9398, 66.9406, 67.9444, 68.9452, 69.9498, 70.9517,
795 71.9564};
796
797 G4double CoS[25] =
798 {43, 28, 18, 16, 70, 19, 14, 15, 26, 15, 19, 15, 15, 17, 22, 22, 22, 14, 29,
799 30, 35, 40, 75, 86, 86};
800
801 G4double CoW[25] =
802 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
803
804 AddElement("Co", 27, 25, *CoN , *CoA , *CoS , *CoW);
805
806 // Z = 28 --------------------------------------------------------------------
807 G4int NiN[29] =
808 {50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68,
809 69, 70, 71, 72, 73, 74, 75, 76, 77, 78};
810
811 G4double NiA[29] =
812 {49.9959, 50.9877, 51.9757, 52.9685, 53.9579, 54.9513, 55.9421, 56.9398,
813 57.9353, 58.9344, 59.9308, 60.9311, 61.9283, 62.9297, 63.928 , 64.9301,
814 65.9291, 66.9316, 67.9318, 68.9352, 69.9361, 70.94 , 71.9413, 72.9461,
815 73.9479, 74.953 , 75.9553, 76.9608, 77.9638};
816
817 G4double NiS[29] =
818 {28, 28, 90, 17, 50, 12, 12, 3, 15, 15, 15, 15, 15, 15, 16, 16, 17, 20, 18,
819 15, 35, 40, 50, 64, 75, 86, 97, 107, 118};
820
821 G4double NiW[29] =
822 {0, 0, 0, 0, 0, 0, 0, 0, 68.0769, 0, 26.2231, 1.1399, 3.6345, 0, 0.9256, 0,
823 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
824
825 AddElement("Ni", 28, 29, *NiN , *NiA , *NiS , *NiW);
826
827 // Z = 29 --------------------------------------------------------------------
828 G4int CuN[29] =
829 {52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70,
830 71, 72, 73, 74, 75, 76, 77, 78, 79, 80};
831
832 G4double CuA[29] =
833 {51.9972, 52.9856, 53.9767, 54.9661, 55.9586, 56.9492, 57.9445, 58.9395,
834 59.9374, 60.9335, 61.9326, 62.9296, 63.9298, 64.9278, 65.9289, 66.9278,
835 67.9296, 68.9294, 69.9324, 70.9326, 71.9355, 72.9365, 73.9402, 74.9417,
836 75.946 , 76.948 , 77.9528, 78.9553, 79.9619};
837
838 G4double CuS[29] =
839 {28, 28, 23, 32, 15, 17, 27, 18, 27, 19, 4, 15, 15, 19, 19, 9, 50, 9, 16, 40,
840 21, 32, 43, 54, 64, 75, 86, 97, 97};
841
842 G4double CuW[29] =
843 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 69.17, 0, 30.83, 0, 0, 0, 0, 0, 0, 0, 0, 0,
844 0, 0, 0, 0, 0, 0};
845
846 AddElement("Cu", 29, 29, *CuN , *CuA , *CuS , *CuW);
847
848 // Z = 30 --------------------------------------------------------------------
849 G4int ZnN[29] =
850 {54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72,
851 73, 74, 75, 76, 77, 78, 79, 80, 81, 82};
852
853 G4double ZnA[29] =
854 {53.993 , 54.984 , 55.9724, 56.9649, 57.9546, 58.9493, 59.9418, 60.9395,
855 61.9343, 62.9332, 63.9291, 64.9292, 65.926 , 66.9271, 67.9248, 68.9266,
856 69.9253, 70.9277, 71.9269, 72.9298, 73.9295, 74.9329, 75.9334, 76.9371,
857 77.9386, 78.9427, 79.9444, 80.9505, 81.9548};
858
859 G4double ZnS[29] =
860 {43, 27, 28, 15, 50, 40, 11, 18, 11, 23, 18, 18, 16, 17, 17, 18, 4, 11, 7,
861 40, 50, 80, 13, 14, 17, 29, 18, 43, 43};
862
863 G4double ZnW[29] =
864 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 48.63, 0, 27.9, 4.1, 18.75, 0, 0.62, 0, 0, 0,
865 0, 0, 0, 0, 0, 0, 0, 0, 0};
866
867 AddElement("Zn", 30, 29, *ZnN , *ZnA , *ZnS , *ZnW);
868
869 // Z = 31 --------------------------------------------------------------------
870 G4int GaN[29] =
871 {56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74,
872 75, 76, 77, 78, 79, 80, 81, 82, 83, 84};
873
874 G4double GaA[29] =
875 {55.9949, 56.9829, 57.9742, 58.9634, 59.9571, 60.9492, 61.9442, 62.9391,
876 63.9368, 64.9327, 65.9316, 66.9282, 67.928 , 68.9256, 69.926 , 70.9247,
877 71.9264, 72.9252, 73.9269, 74.9265, 75.9289, 76.9293, 77.9317, 78.9329,
878 79.9366, 80.9377, 81.9432, 82.9469, 83.9523};
879
880 G4double GaS[29] =
881 {28, 28, 23, 18, 12, 21, 30, 11, 4, 19, 4, 19, 22, 3, 3, 19, 22, 7, 80, 7, 10,
882 60, 90, 13, 13, 21, 32, 54, 64};
883
884 G4double GaW[29] =
885 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 60.108, 0, 39.892, 0, 0, 0, 0, 0, 0,
886 0, 0, 0, 0, 0, 0, 0};
887
888 AddElement("Ga", 31, 29, *GaN , *GaA , *GaS , *GaW);
889
890 // Z = 32 --------------------------------------------------------------------
891 G4int GeN[29] =
892 {58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76,
893 77, 78, 79, 80, 81, 82, 83, 84, 85, 86};
894
895 G4double GeA[29] =
896 {57.991 , 58.9817, 59.9702, 60.9638, 61.9547, 62.9496, 63.9416, 64.9394,
897 65.9339, 66.9327, 67.9281, 68.928 , 69.9243, 70.925 , 71.9221, 72.9235,
898 73.9212, 74.9229, 75.9214, 76.9235, 77.9229, 78.9254, 79.9254, 80.9288,
899 81.9296, 82.9345, 83.9373, 84.9427, 85.9463};
900
901 G4double GeS[29] =
902 {34, 30, 25, 32, 15, 21, 27, 11, 30, 5, 7, 3, 19, 19, 16, 16, 16, 16, 16, 20,
903 4, 10, 25, 13, 26, 32, 43, 54, 64};
904
905 G4double GeW[29] =
906 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 20.84, 0, 27.54, 7.73, 36.28, 0, 7.61,
907 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
908
909 AddElement("Ge", 32, 29, *GeN , *GeA , *GeS , *GeW);
910
911 // Z = 33 --------------------------------------------------------------------
912 G4int AsN[30] =
913 {60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78,
914 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89};
915
916 G4double AsA[30] =
917 {59.9931, 60.9806, 61.9732, 62.9637, 63.9576, 64.9495, 65.9444, 66.9392,
918 67.9368, 68.9323, 69.9309, 70.9271, 71.9268, 72.9238, 73.9239, 74.9216,
919 75.9224, 76.9206, 77.9218, 78.9209, 79.9226, 80.9221, 81.9245, 82.925 ,
920 83.9291, 84.9318, 85.9362, 86.9396, 87.9446, 88.9492};
921
922 G4double AsS[30] =
923 {64, 64, 32, 54, 38, 42, 22, 11, 11, 30, 50, 5, 5, 4, 24, 18, 18, 23, 11, 6,
924 23, 6, 21, 24, 32, 32, 43, 54, 64, 64};
925
926 G4double AsW[30] =
927 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0,
928 0, 0, 0, 0, 0, 0};
929
930 AddElement("As", 33, 30, *AsN , *AsA , *AsS , *AsW);
931
932 // Z = 34 --------------------------------------------------------------------
933 G4int SeN[28] =
934 {65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83,
935 84, 85, 86, 87, 88, 89, 90, 91, 92};
936
937 G4double SeA[28] =
938 {64.9647, 65.9552, 66.9501, 67.9419, 68.9396, 69.9335, 70.9323, 71.9271,
939 72.9268, 73.9225, 74.9225, 75.9192, 76.9199, 77.9173, 78.9185, 79.9165,
940 80.918 , 81.9167, 82.9191, 83.9185, 84.9222, 85.9243, 86.9285, 87.9314,
941 88.936 , 89.9394, 90.9454, 91.9493};
942
943 G4double SeS[28] =
944 {64, 32, 21, 32, 40, 22, 22, 13, 12, 16, 16, 16, 16, 16, 16, 20, 21, 22, 4,
945 16, 30, 17, 40, 50, 32, 43, 54, 64};
946
947 G4double SeW[28] =
948 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0.89, 0, 9.37, 7.63, 23.77, 0, 49.61, 0, 8.73,
949 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
950
951 AddElement("Se", 34, 28, *SeN , *SeA , *SeS , *SeW);
952
953 // Z = 35 --------------------------------------------------------------------
954 G4int BrN[28] =
955 {67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85,
956 86, 87, 88, 89, 90, 91, 92, 93, 94};
957
958 G4double BrA[28] =
959 {66.9648, 67.9583, 68.9502, 69.9446, 70.9393, 71.9365, 72.9318, 73.9299,
960 74.9258, 75.9245, 76.9214, 77.9211, 78.9183, 79.9185, 80.9163, 81.9168,
961 82.9152, 83.9165, 84.9156, 85.9188, 86.9207, 87.9241, 88.9264, 89.9306,
962 90.934 , 91.9393, 92.9431, 93.9487};
963
964 G4double BrS[28] =
965 {54, 58, 34, 39, 32, 28, 14, 16, 15, 10, 3, 4, 20, 20, 3, 3, 5, 27, 21, 12,
966 19, 40, 60, 80, 80, 50, 32, 43};
967
968 G4double BrW[28] =
969 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 50.69, 0, 49.31, 0, 0, 0, 0, 0, 0, 0,
970 0, 0, 0, 0, 0, 0};
971
972 AddElement("Br", 35, 28, *BrN , *BrA , *BrS , *BrW);
973
974 // Z = 36 --------------------------------------------------------------------
975 G4int KrN[29] =
976 {69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87,
977 88, 89, 90, 91, 92, 93, 94, 95, 96, 97};
978
979 G4double KrA[29] =
980 {68.9653, 69.956 , 70.9505, 71.9419, 72.9389, 73.9333, 74.931 , 75.9259,
981 76.9247, 77.9204, 78.9201, 79.9164, 80.9166, 81.9135, 82.9141, 83.9115,
982 84.9125, 85.9106, 86.9134, 87.9144, 88.9176, 89.9195, 90.9234, 91.9262,
983 92.9313, 93.9344, 94.9398, 95.9431, 96.9486};
984
985 G4double KrS[29] =
986 {54, 43, 32, 29, 15, 60, 17, 11, 9, 7, 4, 4, 3, 28, 3, 3, 3, 12, 14, 14, 60,
987 20, 60, 13, 11, 32, 43, 54, 54};
988
989 G4double KrW[29] =
990 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0.35, 0, 2.28, 0, 11.58, 11.49, 57, 0, 17.3, 0,
991 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
992
993 AddElement("Kr", 36, 29, *KrN , *KrA , *KrS , *KrW);
994
995 // Z = 37 --------------------------------------------------------------------
996 G4int RbN[32] =
997 {71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89,
998 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102};
999
1000 G4double RbA[32] =
1001 {70.9653, 71.9591, 72.9504, 73.9445, 74.9386, 75.9351, 76.9304, 77.9281,
1002 78.924 , 79.9225, 80.919 , 81.9182, 82.9151, 83.9144, 84.9118, 85.9112,
1003 86.9092, 87.9113, 88.9123, 89.9148, 90.9165, 91.9197, 92.922 , 93.9264,
1004 94.9293, 95.9343, 96.9373, 97.9417, 98.9454, 99.9499, 100.953, 101.959};
1005
1006 G4double RbS[32] =
1007 {54, 54, 52, 77, 8, 8, 8, 8, 7, 8, 7, 8, 7, 3, 25, 25, 27, 5, 6, 9, 9, 7, 8,
1008 9, 21, 27, 30, 40, 16, 32, 18, 54};
1009
1010 G4double RbW[32] =
1011 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 72.17, 0, 27.83, 0, 0, 0, 0, 0,
1012 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1013
1014 AddElement("Rb", 37, 32, *RbN , *RbA , *RbS , *RbW);
1015
1016 // Z = 38 --------------------------------------------------------------------
1017 G4int SrN[32] =
1018 {73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91,
1019 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104};
1020
1021 G4double SrA[32] =
1022 {72.966 , 73.9563, 74.9499, 75.9416, 76.9378, 77.9322, 78.9297, 79.9245,
1023 80.9232, 81.9184, 82.9176, 83.9134, 84.9129, 85.9093, 86.9089, 87.9056,
1024 88.9075, 89.9077, 90.9102, 91.911 , 92.914 , 93.9154, 94.9194, 95.9217,
1025 96.9261, 97.9285, 98.9333, 99.9353, 100.941, 101.943, 102.949, 103.952};
1026
1027 G4double SrS[32] =
1028 {64, 54, 32, 32, 16, 8, 9, 8, 8, 6, 9, 4, 4, 24, 24, 24, 24, 29, 6, 7, 8, 8,
1029 8, 26, 20, 27, 15, 14, 13, 12, 54, 75};
1030
1031 G4double SrW[32] =
1032 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.56, 0, 9.86, 7, 82.58, 0, 0, 0, 0, 0,
1033 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1034
1035 AddElement("Sr", 38, 32, *SrN , *SrA , *SrS , *SrW);
1036
1037 // Z = 39 --------------------------------------------------------------------
1038 G4int YN[30] =
1039 {77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
1040 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106};
1041
1042 G4double YA[30] =
1043 {76.9496, 77.9435, 78.9373, 79.9343, 80.9291, 81.9268, 82.9223, 83.9204,
1044 84.9164, 85.9149, 86.9109, 87.9095, 88.9058, 89.9072, 90.9073, 91.9089,
1045 92.9096, 93.9116, 94.9128, 95.9159, 96.9181, 97.9222, 98.9246, 99.9278,
1046 100.93 , 101.934, 102.937, 103.941, 104.945, 105.95};
1047
1048 G4double YS[30] =
1049 {32, 43, 48, 43, 70, 11, 50, 10, 27, 15, 28, 29, 25, 25, 3, 10, 11, 8, 8, 23,
1050 13, 26, 26, 80, 10, 90, 32, 43, 54, 75};
1051
1052 G4double YW[30] =
1053 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1054 0, 0, 0, 0, 0};
1055
1056 AddElement("Y", 39, 30, *YN , *YA , *YS , *YW);
1057
1058 // Z = 40 --------------------------------------------------------------------
1059 G4int ZrN[30] =
1060 {79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97,
1061 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108};
1062
1063 G4double ZrA[30] =
1064 {78.9492, 79.9406, 80.9368, 81.9311, 82.9287, 83.9232, 84.9215, 85.9165,
1065 86.9148, 87.9102, 88.9089, 89.9047, 90.9056, 91.905 , 92.9065, 93.9063,
1066 94.908 , 95.9083, 96.911 , 97.9127, 98.9165, 99.9178, 100.921, 101.923,
1067 102.927, 103.929, 104.933, 105.936, 106.941, 107.944};
1068
1069 G4double ZrS[30] =
1070 {43, 32, 32, 55, 10, 21, 11, 30, 9, 11, 4, 23, 23, 23, 23, 25, 25, 3, 3, 21,
1071 21, 40, 30, 50, 12, 43, 43, 54, 64, 75};
1072
1073 G4double ZrW[30] =
1074 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 51.45, 11.22, 17.15, 0, 17.38, 0, 2.8, 0,
1075 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1076
1077 AddElement("Zr", 40, 30, *ZrN , *ZrA , *ZrS , *ZrW);
1078
1079 // Z = 41 --------------------------------------------------------------------
1080 G4int NbN[30] =
1081 {81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99,
1082 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110};
1083
1084 G4double NbA[30] =
1085 {80.949 , 81.9431, 82.9367, 83.9336, 84.9279, 85.925 , 86.9204, 87.918 ,
1086 88.9135, 89.9113, 90.907 , 91.9072, 92.9064, 93.9073, 94.9068, 95.9081,
1087 96.9081, 97.9103, 98.9116, 99.9142, 100.915, 101.918, 102.919, 103.922,
1088 104.924, 105.928, 106.93 , 107.935, 108.938, 109.943};
1089
1090 G4double NbS[30] =
1091 {43, 32, 34, 32, 24, 90, 70, 22, 40, 5, 3, 29, 24, 24, 20, 4, 28, 6, 14, 28,
1092 20, 40, 70, 12, 11, 32, 43, 54, 54, 64};
1093
1094 G4double NbW[30] =
1095 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1096 0, 0, 0, 0, 0};
1097
1098 AddElement("Nb", 41, 30, *NbN , *NbA , *NbS , *NbW);
1099
1100 // Z = 42 --------------------------------------------------------------------
1101 G4int MoN[31] =
1102 {83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101,
1103 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113};
1104
1105 G4double MoA[31] =
1106 {82.9487, 83.9401, 84.9366, 85.9307, 86.9273, 87.922 , 88.9195, 89.9139,
1107 90.9118, 91.9068, 92.9068, 93.9051, 94.9058, 95.9047, 96.906 , 97.9054,
1108 98.9077, 99.9075, 100.91 , 101.91 , 102.913, 103.914, 104.917, 105.918,
1109 106.922, 107.924, 108.928, 109.93 , 110.935, 111.937, 112.942};
1110
1111 G4double MoS[31] =
1112 {54, 43, 43, 47, 24, 22, 17, 7, 12, 4, 4, 20, 20, 20, 20, 20, 20, 6, 6, 22,
1113 70, 70, 80, 23, 17, 21, 32, 43, 54, 64, 64};
1114
1115 G4double MoW[31] =
1116 {0, 0, 0, 0, 0, 0, 0, 0, 0, 14.84, 0, 9.25, 15.92, 16.68, 9.55, 24.13, 0,
1117 9.63, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1118
1119 AddElement("Mo", 42, 31, *MoN , *MoA , *MoS , *MoW);
1120
1121 // Z = 43 --------------------------------------------------------------------
1122 G4int TcN[31] =
1123 {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102,
1124 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115};
1125
1126 G4double TcA[31] =
1127 {84.9489, 85.9429, 86.9365, 87.9328, 88.9275, 89.9236, 90.9184, 91.9153,
1128 92.9102, 93.9097, 94.9077, 95.9079, 96.9064, 97.9072, 98.9063, 99.9077,
1129 100.907, 101.909, 102.909, 103.911, 104.912, 105.914, 106.915, 107.918,
1130 108.92 , 109.923, 110.925, 111.929, 112.931, 113.936, 114.938};
1131
1132 G4double TcS[31] =
1133 {54, 32, 32, 32, 23, 26, 22, 28, 4, 5, 6, 6, 5, 4, 21, 23, 26, 10, 11, 50,
1134 60, 15, 16, 14, 23, 43, 43, 54, 64, 64, 75};
1135
1136 G4double TcW[31] =
1137 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1138 0, 0, 0, 0, 0, 0, 0};
1139
1140 AddElement("Tc", 43, 31, *TcN , *TcA , *TcS , *TcW);
1141
1142 // Z = 44 --------------------------------------------------------------------
1143 G4int RuN[32] =
1144 {87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104,
1145 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118};
1146
1147 G4double RuA[32] =
1148 {86.9492, 87.9404, 88.9361, 89.9298, 90.9264, 91.9201, 92.9171, 93.9114,
1149 94.9104, 95.9076, 96.9076, 97.9053, 98.9059, 99.9042, 100.906, 101.904,
1150 102.906, 103.905, 104.908, 105.907, 106.91 , 107.91 , 108.913, 109.914,
1151 110.918, 111.919, 112.923, 113.924, 114.928, 115.93 , 116.935, 117.937};
1152
1153 G4double RuS[32] =
1154 {64, 54, 54, 43, 54, 32, 90, 14, 13, 8, 9, 7, 21, 22, 22, 22, 22, 4, 4, 8,
1155 13, 13, 70, 25, 32, 58, 54, 39, 64, 75, 86, 97};
1156
1157 G4double RuW[32] =
1158 {0, 0, 0, 0, 0, 0, 0, 0, 0, 5.54, 0, 1.87, 12.76, 12.6, 17.06, 31.55, 0,
1159 18.62, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1160
1161 AddElement("Ru", 44, 32, *RuN , *RuA , *RuS , *RuW);
1162
1163 // Z = 45 --------------------------------------------------------------------
1164 G4int RhN[33] =
1165 {89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105,
1166 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120,
1167 121};
1168
1169 G4double RhA[33] =
1170 {88.9494, 89.9429, 90.9365, 91.932 , 92.9257, 93.9217, 94.9159, 95.9145,
1171 96.9113, 97.9107, 98.9081, 99.9081, 100.906, 101.907, 102.906, 103.907,
1172 104.906, 105.907, 106.907, 107.909, 108.909, 109.911, 110.912, 111.915,
1173 112.915, 113.919, 114.92 , 115.924, 116.925, 117.929, 118.931, 119.936,
1174 120.938};
1175
1176 G4double RhS[33] =
1177 {54, 54, 43, 43, 43, 48, 16, 14, 40, 13, 8, 22, 19, 5, 3, 3, 5, 8, 13, 11,
1178 13, 24, 22, 54, 43, 32, 54, 54, 64, 75, 86, 86, 97};
1179 G4double RhW[33] =
1180 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1181 0, 0, 0, 0, 0, 0, 0, 0, 0};
1182
1183 AddElement("Rh", 45, 33, *RhN , *RhA , *RhS , *RhW);
1184
1185 // Z = 46 --------------------------------------------------------------------
1186 G4int PdN[33] =
1187 {91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107,
1188 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122,
1189 123};
1190
1191 G4double PdA[33] =
1192 {90.9495, 91.9404, 92.9359, 93.9288, 94.9247, 95.9182, 96.9165, 97.9127,
1193 98.9118, 99.9085, 100.908, 101.906, 102.906, 103.904, 104.905, 105.903,
1194 106.905, 107.904, 108.906, 109.905, 110.908, 111.907, 112.91 , 113.91 ,
1195 114.914, 115.914, 116.918, 117.919, 118.923, 119.924, 120.928, 121.93 ,
1196 122.934};
1197
1198 G4double PdS[33] =
1199 {64, 54, 43, 43, 43, 16, 32, 23, 16, 12, 19, 3, 3, 5, 5, 5, 7, 4, 4, 12, 40,
1200 19, 40, 26, 70, 60, 32, 23, 32, 43, 54, 54, 64};
1201
1202 G4double PdW[33] =
1203 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.02, 0, 11.14, 22.33, 27.33, 0, 26.46, 0,
1204 11.72, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1205
1206 AddElement("Pd", 46, 33, *PdN , *PdA , *PdS , *PdW);
1207
1208 // Z = 47 --------------------------------------------------------------------
1209 G4int AgN[34] =
1210 {94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109,
1211 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124,
1212 125, 126, 127};
1213
1214 G4double AgA[34] =
1215 {93.9428, 94.9355, 95.9307, 96.924 , 97.9218, 98.9176, 99.9161, 100.913,
1216 101.912, 102.909, 103.909, 104.907, 105.907, 106.905, 107.906, 108.905,
1217 109.906, 110.905, 111.907, 112.907, 113.909, 114.909, 115.911, 116.912,
1218 117.915, 118.916, 119.919, 120.92 , 121.923, 122.925, 123.929, 124.931,
1219 125.934, 126.937};
1220
1221 G4double AgS[34] =
1222 {54, 43, 43, 43, 16, 16, 80, 11, 80, 18, 7, 12, 6, 6, 6, 3, 3, 4, 18, 18, 28,
1223 40, 50, 50, 70, 10, 80, 16, 22, 32, 43, 43, 43, 54};
1224
1225 G4double AgW[34] =
1226 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 51.839, 0, 48.161, 0, 0, 0, 0, 0, 0,
1227 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1228
1229 AddElement("Ag", 47, 34, *AgN , *AgA , *AgS , *AgW);
1230
1231 // Z = 48 --------------------------------------------------------------------
1232 G4int CdN[35] =
1233 {96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
1234 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126,
1235 127, 128, 129, 130};
1236
1237 G4double CdA[35] =
1238 {95.9398, 96.9349, 97.9276, 98.925 , 99.9202, 100.919, 101.915, 102.913,
1239 103.91 , 104.909, 105.906, 106.907, 107.904, 108.905, 109.903, 110.904,
1240 111.903, 112.904, 113.903, 114.905, 115.905, 116.907, 117.907, 118.91 ,
1241 119.91 , 120.913, 121.913, 122.917, 123.918, 124.921, 125.922, 126.926,
1242 127.928, 128.932, 129.934};
1243
1244 G4double CdS[35] =
1245 {54, 43, 22, 22, 10, 16, 80, 17, 10, 12, 6, 7, 6, 4, 3, 3, 30, 30, 30, 3, 3,
1246 4, 22, 90, 20, 90, 22, 40, 70, 70, 60, 80, 32, 43, 43};
1247
1248 G4double CdW[35] =
1249 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.25, 0, 0.89, 0, 12.49, 12.8, 24.13, 12.22,
1250 28.73, 0, 7.49, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1251
1252 AddElement("Cd", 48, 35, *CdN , *CdA , *CdS , *CdW);
1253
1254 // Z = 49 --------------------------------------------------------------------
1255 G4int InN[37] =
1256 {98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112,
1257 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127,
1258 128, 129, 130, 131, 132, 133, 134};
1259
1260 G4double InA[37] =
1261 {97.9422, 98.9346, 99.9312, 100.927, 101.925, 102.92 , 103.918, 104.915,
1262 105.913, 106.91 , 107.91 , 108.907, 109.907, 110.905, 111.906, 112.904,
1263 113.905, 114.904, 115.905, 116.905, 117.906, 118.906, 119.908, 120.908,
1264 121.91 , 122.91 , 123.913, 124.914, 125.916, 126.917, 127.92 , 128.922,
1265 129.925, 130.927, 131.933, 132.938, 133.945};
1266
1267 G4double InS[37] =
1268 {54, 54, 41, 32, 41, 27, 15, 19, 15, 14, 40, 6, 13, 6, 6, 4, 3, 5, 5, 6, 9,
1269 8, 40, 29, 50, 26, 50, 30, 40, 40, 50, 14, 50, 80, 70, 43, 54};
1270
1271 G4double InW[37] =
1272 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4.29, 0, 95.71, 0, 0, 0, 0, 0,
1273 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1274
1275 AddElement("In", 49, 37, *InN , *InA , *InS , *InW);
1276
1277 // Z = 50 --------------------------------------------------------------------
1278 G4int SnN[38] =
1279 {100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114,
1280 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129,
1281 130, 131, 132, 133, 134, 135, 136, 137};
1282
1283 G4double SnA[38] =
1284 {99.939 , 100.936, 101.93 , 102.928, 103.923, 104.921, 105.917, 106.916,
1285 107.912, 108.911, 109.908, 110.908, 111.905, 112.905, 113.903, 114.903,
1286 115.902, 116.903, 117.902, 118.903, 119.902, 120.904, 121.903, 122.906,
1287 123.905, 124.908, 125.908, 126.91 , 127.911, 128.913, 129.914, 130.917,
1288 131.918, 132.924, 133.928, 134.935, 135.939, 136.946};
1289
1290 G4double SnS[38] =
1291 {46, 54, 43, 32, 16, 10, 50, 90, 50, 11, 17, 8, 5, 4, 3, 3, 3, 3, 3, 3, 27,
1292 27, 29, 29, 15, 16, 11, 27, 29, 13, 30, 80, 28, 90, 11, 43, 54, 64};
1293
1294 G4double SnW[38] =
1295 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.97, 0, 0.66, 0.34, 14.54, 7.68, 24.22,
1296 8.59, 32.58, 0, 4.63, 0, 5.79, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1297
1298 AddElement("Sn", 50, 38, *SnN , *SnA , *SnS , *SnW);
1299
1300 // Z = 51 --------------------------------------------------------------------
1301 G4int SbN[37] =
1302 {103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117,
1303 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132,
1304 133, 134, 135, 136, 137, 138, 139};
1305
1306 G4double SbA[37] =
1307 {102.94 , 103.936, 104.932, 105.929, 106.924, 107.922, 108.918, 109.917,
1308 110.913, 111.912, 112.909, 113.909, 114.907, 115.907, 116.905, 117.906,
1309 118.904, 119.905, 120.904, 121.905, 122.904, 123.906, 124.905, 125.907,
1310 126.907, 127.909, 128.909, 129.912, 130.912, 131.914, 132.915, 133.921,
1311 134.925, 135.931, 136.935, 137.941, 138.946};
1312
1313 G4double SbS[37] =
1314 {54, 39, 17, 34, 32, 22, 20, 22, 22, 25, 24, 22, 22, 6, 10, 4, 9, 8, 24, 24,
1315 22, 22, 3, 30, 6, 27, 23, 27, 80, 25, 80, 60, 11, 32, 43, 54, 64};
1316
1317 G4double SbW[37] =
1318 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 57.21, 0, 42.79, 0, 0,
1319 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1320
1321 AddElement("Sb", 51, 37, *SbN , *SbA , *SbS , *SbW);
1322
1323 // Z = 52 --------------------------------------------------------------------
1324 G4int TeN[37] =
1325 {106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120,
1326 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135,
1327 136, 137, 138, 139, 140, 141, 142};
1328
1329 G4double TeA[37] =
1330 {105.938, 106.935, 107.929, 108.927, 109.922, 110.921, 111.917, 112.916,
1331 113.912, 114.912, 115.908, 116.909, 117.906, 118.906, 119.904, 120.905,
1332 121.903, 122.904, 123.903, 124.904, 125.903, 126.905, 127.904, 128.907,
1333 129.906, 130.909, 131.909, 132.911, 133.912, 134.916, 135.92 , 136.925,
1334 137.929, 138.935, 139.939, 140.944, 141.948};
1335
1336 G4double TeS[37] =
1337 {43, 32, 16, 80, 60, 80, 18, 22, 22, 11, 10, 20, 17, 9, 11, 27, 20, 19, 16,
1338 20, 20, 4, 19, 3, 21, 22, 12, 80, 40, 10, 50, 13, 22, 43, 54, 54, 64};
1339
1340 G4double TeW[37] =
1341 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.09, 0, 2.55, 0.89, 4.74, 7.07,
1342 18.84, 0, 31.74, 0, 34.08, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1343
1344 AddElement("Te", 52, 37, *TeN , *TeA , *TeS , *TeW);
1345
1346 // Z = 53 --------------------------------------------------------------------
1347 G4int IN[37] =
1348 {108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122,
1349 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137,
1350 138, 139, 140, 141, 142, 143, 144};
1351
1352 G4double IA[37] =
1353 {107.943, 108.938, 109.935, 110.93 , 111.928, 112.924, 113.922, 114.918,
1354 115.917, 116.914, 117.913, 118.91 , 119.91 , 120.907, 121.908, 122.906,
1355 123.906, 124.905, 125.906, 126.904, 127.906, 128.905, 129.907, 130.906,
1356 131.908, 132.908, 133.91 , 134.91 , 135.915, 136.918, 137.922, 138.926,
1357 139.931, 140.935, 141.94, 142.944, 143.95};
1358
1359 G4double IS[37] =
1360 {39, 16, 33, 32, 23, 60, 32, 50, 15, 80, 80, 70, 20, 12, 6, 4, 26, 20, 4, 4,
1361 4, 4, 4, 12, 11, 28, 16, 25, 50, 30, 90, 30, 23, 32, 43, 43, 54};
1362
1363 G4double IW[37] =
1364 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0,
1365 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1366
1367 AddElement("I", 53, 37, *IN , *IA , *IS , *IW);
1368
1369 // Z = 54 --------------------------------------------------------------------
1370 G4int XeN[38] =
1371 {110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124,
1372 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139,
1373 140, 141, 142, 143, 144, 145, 146, 147};
1374
1375 G4double XeA[38] =
1376 {109.944, 110.942, 111.936, 112.933, 113.928, 114.927, 115.922, 116.921,
1377 117.917, 118.916, 119.912, 120.911, 121.909, 122.908, 123.906, 124.906,
1378 125.904, 126.905, 127.904, 128.905, 129.904, 130.905, 131.904, 132.906,
1379 133.905, 134.907, 135.907, 136.912, 137.914, 138.919, 139.922, 140.927,
1380 141.93 , 142.935, 143.938, 144.944, 145.947, 146.953};
1381
1382 G4double XeS[38] =
1383 {43, 33, 16, 10, 22, 26, 26, 19, 107, 13, 50, 26, 90, 17, 21, 21, 7, 4, 15,
1384 9, 10, 10, 12, 4, 9, 11, 8, 8, 40, 23, 70, 10, 11, 24, 34, 43, 43, 54};
1385
1386 G4double XeW[38] =
1387 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.09, 0, 0.09, 0, 1.92, 26.44,
1388 4.08, 21.18, 26.89, 0, 10.44, 0, 8.87, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1389
1390 AddElement("Xe", 54, 38, *XeN , *XeA , *XeS , *XeW);
1391
1392 // Z = 55 --------------------------------------------------------------------
1393 G4int CsN[40] =
1394 {112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126,
1395 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141,
1396 142, 143, 144, 145, 146, 147, 148, 149, 150, 151};
1397
1398 G4double CsA[40] =
1399 {111.95 , 112.945, 113.941, 114.936, 115.933, 116.929, 117.927, 118.922,
1400 119.921, 120.917, 121.916, 122.913, 123.912, 124.91 , 125.909, 126.907,
1401 127.908, 128.906, 129.907, 130.905, 131.906, 132.905, 133.907, 134.906,
1402 135.907, 136.907, 137.911, 138.913, 139.917, 140.92 , 141.924, 142.927,
1403 143.932, 144.935, 145.94 , 146.944, 147.949, 148.953, 149.958, 150.962};
1404
1405 G4double CsS[40] =
1406 {33, 16, 33, 46, 38, 60, 14, 15, 11, 15, 18, 13, 13, 8, 13, 9, 6, 5, 9, 6, 4,
1407 3, 3, 3, 4, 3, 10, 5, 9, 11, 11, 24, 30, 50, 90, 16, 63, 32, 54, 75};
1408
1409 G4double CsW[40] =
1410 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0,
1411 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1412
1413 AddElement("Cs", 55, 40, *CsN , *CsA , *CsS , *CsW);
1414
1415 // Z = 56 --------------------------------------------------------------------
1416 G4int BaN[40] =
1417 {114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128,
1418 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143,
1419 144, 145, 146, 147, 148, 149, 150, 151, 152, 153};
1420
1421 G4double BaA[40] =
1422 {113.951, 114.948, 115.942, 116.939, 117.933, 118.931, 119.926, 120.924,
1423 121.92 , 122.919, 123.915, 124.915, 125.911, 126.911, 127.908, 128.909,
1424 129.906, 130.907, 131.905, 132.906, 133.905, 134.906, 135.905, 136.906,
1425 137.905, 138.909, 139.911, 140.914, 141.916, 142.921, 143.923, 144.927,
1426 145.93 , 146.934, 147.938, 148.942, 149.946, 150.951, 151.954, 152.96};
1427
1428 G4double BaS[40] =
1429 {48, 64, 54, 70, 54, 109, 32, 33, 32, 32, 15, 27, 15, 11, 12, 12, 7, 7, 3, 3,
1430 3, 3, 3, 3, 3, 3, 9, 9, 7, 14, 15, 60, 80, 10, 15, 43, 54, 64, 75, 97};
1431
1432 G4double BaW[40] =
1433 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.106, 0, 0.101, 0, 2.417,
1434 6.592, 7.854, 11.232, 71.698, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1435
1436 AddElement("Ba", 56, 40, *BaN , *BaA , *BaS , *BaW);
1437
1438 // Z = 57 --------------------------------------------------------------------
1439 G4int LaN[39] =
1440 {117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131,
1441 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146,
1442 147, 148, 149, 150, 151, 152, 153, 154, 155};
1443
1444 G4double LaA[39] =
1445 {116.95 , 117.947, 118.941, 119.938, 120.933, 121.931, 122.926, 123.925,
1446 124.921, 125.919, 126.916, 127.915, 128.913, 129.912, 130.91 , 131.91 ,
1447 132.908, 133.908, 134.907, 135.908, 136.906, 137.907, 138.906, 139.909,
1448 140.911, 141.914, 142.916, 143.92 , 144.922, 145.926, 146.928, 147.932,
1449 148.934, 149.939, 150.942, 151.946, 152.949, 153.954, 154.958};
1450
1451 G4double LaS[39] =
1452 {96, 86, 75, 64, 54, 54, 43, 32, 32, 32, 24, 43, 60, 22, 11, 50, 21, 28, 11,
1453 80, 50, 4, 3, 3, 5, 6, 16, 60, 70, 80, 80, 14, 32, 43, 54, 64, 75, 86, 97};
1454
1455 G4double LaW[39] =
1456 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.09, 99.91,
1457 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1458
1459 AddElement("La", 57, 39, *LaN , *LaA , *LaS , *LaW);
1460
1461 // Z = 58 --------------------------------------------------------------------
1462 G4int CeN[39] =
1463 {119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133,
1464 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148,
1465 149, 150, 151, 152, 153, 154, 155, 156, 157};
1466
1467 G4double CeA[39] =
1468 {118.953, 119.947, 120.944, 121.938, 122.936, 123.931, 124.929, 125.924,
1469 126.923, 127.919, 128.918, 129.915, 130.914, 131.911, 132.912, 133.909,
1470 134.909, 135.907, 136.908, 137.906, 138.907, 139.905, 140.908, 141.909,
1471 142.912, 143.914, 144.917, 145.919, 146.923, 147.924, 148.928, 149.93 ,
1472 150.934, 151.936, 152.941, 153.943, 154.948, 155.951, 156.956};
1473
1474 G4double CeS[39] =
1475 {97, 86, 75, 64, 54, 54, 43, 43, 32, 32, 22, 66, 44, 21, 21, 22, 12, 50, 50,
1476 11, 8, 3, 3, 4, 4, 4, 40, 70, 60, 13, 80, 13, 32, 43, 54, 64, 75, 86, 97};
1477
1478 G4double CeW[39] =
1479 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.185, 0, 0.251, 0,
1480 88.45, 0, 11.114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1481
1482 AddElement("Ce", 58, 39, *CeN , *CeA , *CeS , *CeW);
1483
1484 // Z = 59 --------------------------------------------------------------------
1485 G4int PrN[39] =
1486 {121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135,
1487 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150,
1488 151, 152, 153, 154, 155, 156, 157, 158, 159};
1489
1490 G4double PrA[39] =
1491 {120.955, 121.952, 122.946, 123.943, 124.938, 125.935, 126.931, 127.929,
1492 128.925, 129.923, 130.92 , 131.919, 132.916, 133.916, 134.913, 135.913,
1493 136.911, 137.911, 138.909, 139.909, 140.908, 141.91 , 142.911, 143.913,
1494 144.915, 145.918, 146.919, 147.922, 148.924, 149.927, 150.928, 151.932,
1495 152.934, 153.937, 154.94 , 155.944, 156.947, 157.952, 158.955};
1496
1497 G4double PrS[39] =
1498 {86, 86, 75, 64, 54, 54, 43, 43, 32, 32, 47, 21, 21, 32, 16, 50, 50, 16, 9,
1499 7, 3, 3, 3, 4, 8, 60, 40, 10, 11, 90, 40, 32, 32, 43, 54, 64, 75, 86, 97};
1500
1501 G4double PrW[39] =
1502 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0,
1503 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1504
1505 AddElement("Pr", 59, 39, *PrN , *PrA , *PrS , *PrW);
1506
1507 // Z = 60 --------------------------------------------------------------------
1508 G4int NdN[36] =
1509 {126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140,
1510 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155,
1511 156, 157, 158, 159, 160, 161};
1512
1513 G4double NdA[36] =
1514 {125.943, 126.941, 127.935, 128.933, 129.929, 130.927, 131.923, 132.922,
1515 133.919, 134.918, 135.915, 136.915, 137.912, 138.912, 139.909, 140.91 ,
1516 141.908, 142.91 , 143.91 , 144.913, 145.913, 146.916, 147.917, 148.92 ,
1517 149.921, 150.924, 151.925, 152.928, 153.929, 154.933, 155.935, 156.939,
1518 157.942, 158.946, 159.949, 160.954};
1519
1520 G4double NdS[36] =
1521 {75, 64, 64, 39, 54, 50, 32, 32, 36, 22, 60, 80, 22, 50, 21, 4, 3, 3, 3, 3,
1522 3, 3, 3, 3, 4, 4, 30, 29, 12, 16, 43, 54, 64, 75, 86, 97};
1523
1524 G4double NdW[36] =
1525 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 27.2, 12.2, 23.8, 8.3, 17.2,
1526 0, 5.7, 0, 5.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1527
1528 AddElement("Nd", 60, 36, *NdN , *NdA , *NdS , *NdW);
1529
1530 // Z = 61 --------------------------------------------------------------------
1531 G4int PmN[36] =
1532 {128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142,
1533 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
1534 158, 159, 160, 161, 162, 163};
1535
1536 G4double PmA[36] =
1537 {127.948, 128.943, 129.94 , 130.936, 131.934, 132.93 , 133.928, 134.925,
1538 135.923, 136.921, 137.919, 138.917, 139.916, 140.914, 141.913, 142.911,
1539 143.913, 144.913, 145.915, 146.915, 147.917, 148.918, 149.921, 150.921,
1540 151.923, 152.924, 153.927, 154.928, 155.931, 156.933, 157.937, 158.939,
1541 159.943, 160.946, 161.95 , 162.954};
1542
1543 G4double PmS[36] =
1544 {97, 86, 75, 64, 54, 54, 42, 35, 22, 15, 34, 60, 30, 29, 50, 4, 4, 4, 5, 3,
1545 7, 5, 22, 6, 80, 12, 80, 30, 40, 32, 43, 54, 64, 75, 86, 97};
1546
1547 G4double PmW[36] =
1548 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0,
1549 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1550
1551 AddElement("Pm", 61, 36, *PmN , *PmA , *PmS , *PmW);
1552
1553 // Z = 62 --------------------------------------------------------------------
1554 G4int SmN[36] =
1555 {130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144,
1556 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159,
1557 160, 161, 162, 163, 164, 165};
1558
1559 G4double SmA[36] =
1560 {129.949, 130.946, 131.941, 132.939, 133.934, 134.932, 135.928, 136.927,
1561 137.924, 138.922, 139.919, 140.918, 141.915, 142.915, 143.912, 144.913,
1562 145.913, 146.915, 147.915, 148.917, 149.917, 150.92 , 151.92 , 152.922,
1563 153.922, 154.925, 155.926, 156.928, 157.93 , 158.933, 159.935, 160.939,
1564 161.941, 162.945, 163.948, 164.953};
1565
1566 G4double SmS[36] =
1567 {97, 97, 75, 64, 54, 54, 43, 12, 32, 16, 16, 13, 11, 4, 4, 4, 4, 3, 3, 3, 3,
1568 3, 3, 3, 3, 3, 10, 50, 80, 32, 43, 54, 64, 75, 86, 97};
1569
1570 G4double SmW[36] =
1571 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.07, 0, 0, 14.99, 11.24, 13.82,
1572 7.38, 0, 26.75, 0, 22.75, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1573
1574 AddElement("Sm", 62, 36, *SmN , *SmA , *SmS , *SmW);
1575
1576 // Z = 63 --------------------------------------------------------------------
1577 G4int EuN[36] =
1578 {132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146,
1579 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161,
1580 162, 163, 164, 165, 166, 167};
1581
1582 G4double EuA[36] =
1583 {131.954, 132.949, 133.946, 134.942, 135.94 , 136.935, 137.933, 138.93 ,
1584 139.928, 140.925, 141.923, 142.92 , 143.919, 144.916, 145.917, 146.917,
1585 147.918, 148.918, 149.92 , 150.92 , 151.922, 152.921, 153.923, 154.923,
1586 155.925, 156.925, 157.928, 158.929, 159.932, 160.934, 161.937, 162.939,
1587 163.943, 164.946, 165.95 , 166.953};
1588
1589 G4double EuS[36] =
1590 {97, 97, 75, 64, 54, 54, 43, 16, 60, 30, 30, 14, 19, 5, 8, 4, 19, 5, 8, 3,
1591 3, 3, 3, 3, 6, 7, 80, 9, 22, 32, 43, 54, 64, 75, 86, 97};
1592
1593 G4double EuW[36] =
1594 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 47.81, 0, 52.19,
1595 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1596
1597 AddElement("Eu", 63, 36, *EuN , *EuA , *EuS , *EuW);
1598
1599 // Z = 64 --------------------------------------------------------------------
1600 G4int GdN[34] =
1601 {136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150,
1602 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165,
1603 166, 167, 168, 169};
1604
1605 G4double GdA[34] =
1606 {135.947, 136.945, 137.94 , 138.938, 139.934, 140.932, 141.928, 142.927,
1607 143.923, 144.922, 145.918, 146.919, 147.918, 148.919, 149.919, 150.92 ,
1608 151.92 , 152.922, 153.921, 154.923, 155.922, 156.924, 157.924, 158.926,
1609 159.927, 160.93 , 161.931, 162.934, 163.936, 164.939, 165.942, 166.946,
1610 167.948, 168.953};
1611
1612 G4double GdS[34] =
1613 {75, 64, 54, 54, 43, 32, 32, 22, 22, 40, 6, 4, 4, 5, 7, 4, 3, 3, 3, 3, 3, 3,
1614 3, 3, 3, 3, 5, 32, 43, 54, 64, 64, 75, 86};
1615
1616 G4double GdW[34] =
1617 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2, 0, 2.18, 14.8, 20.47,
1618 15.65, 24.84, 0, 21.86, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1619
1620 AddElement("Gd", 64, 34, *GdN , *GdA , *GdS , *GdW);
1621
1622 // Z = 65 --------------------------------------------------------------------
1623 G4int TbN[34] =
1624 {138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152,
1625 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167,
1626 168, 169, 170, 171};
1627
1628 G4double TbA[34] =
1629 {137.953, 138.948, 139.946, 140.941, 141.939, 142.935, 143.933, 144.929,
1630 145.927, 146.924, 147.924, 148.923, 149.924, 150.923, 151.924, 152.923,
1631 153.925, 154.923, 155.925, 156.924, 157.925, 158.925, 159.927, 160.928,
1632 161.929, 162.931, 163.933, 164.935, 165.938, 166.94 , 167.944, 168.946,
1633 169.95 , 170.953};
1634
1635 G4double TbS[34] =
1636 {86, 75, 96, 64, 82, 43, 32, 24, 50, 13, 30, 5, 9, 5, 40, 5, 50, 13, 5, 3,
1637 3, 3, 3, 3, 40, 5, 11, 21, 32, 43, 54, 64, 75, 86};
1638
1639 G4double TbW[34] =
1640 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0,
1641 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1642
1643 AddElement("Tb", 65, 34, *TbN , *TbA , *TbS , *TbW);
1644
1645 // Z = 66 --------------------------------------------------------------------
1646 G4int DyN[34] =
1647 {140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154,
1648 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169,
1649 170, 171, 172, 173};
1650
1651 G4double DyA[34] =
1652 {139.954, 140.951, 141.946, 142.944, 143.939, 144.937, 145.933, 146.931,
1653 147.927, 148.927, 149.926, 150.926, 151.925, 152.926, 153.924, 154.926,
1654 155.924, 156.925, 157.924, 158.926, 159.925, 160.927, 161.927, 162.929,
1655 163.929, 164.932, 165.933, 166.936, 167.937, 168.94 , 169.943, 170.946,
1656 171.949, 172.953};
1657
1658 G4double DyS[34] =
1659 {97, 75, 85, 54, 43, 32, 12, 60, 30, 12, 6, 5, 6, 5, 9, 13, 7, 7, 4, 3, 3, 3,
1660 3, 3, 3, 3, 3, 60, 32, 32, 43, 54, 64, 75};
1661
1662 G4double DyW[34] =
1663 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.06, 0, 0.1, 0, 2.34,
1664 18.91, 25.51, 24.9, 28.18, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1665
1666 AddElement("Dy", 66, 34, *DyN , *DyA , *DyS , *DyW);
1667
1668 // Z = 67 --------------------------------------------------------------------
1669 G4int HoN[34] =
1670 {142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156,
1671 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171,
1672 172, 173, 174, 175};
1673
1674 G4double HoA[34] =
1675 {141.96 , 142.955, 143.952, 144.947, 145.944, 146.94 , 147.937, 148.934,
1676 149.933, 150.932, 151.932, 152.93 , 153.931, 154.929, 155.93 , 156.928,
1677 157.929, 158.928, 159.929, 160.928, 161.929, 162.929, 163.93 , 164.93 ,
1678 165.932, 166.933, 167.935, 168.937, 169.94 , 170.941, 171.945, 172.947,
1679 173.951, 174.954};
1680
1681 G4double HoS[34] =
1682 {107, 75, 64, 64, 54, 43, 29, 23, 11, 13, 30, 6, 10, 25, 22, 50, 30, 4, 16,
1683 4, 5, 3, 3, 3, 3, 6, 30, 22, 50, 64, 43, 43, 54, 64};
1684
1685 G4double HoW[34] =
1686 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100,
1687 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1688
1689 AddElement("Ho", 67, 34, *HoN , *HoA , *HoS , *HoW);
1690
1691 // Z = 68 --------------------------------------------------------------------
1692 G4int ErN[34] =
1693 {144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158,
1694 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173,
1695 174, 175, 176, 177};
1696
1697 G4double ErA[34] =
1698 {143.961, 144.957, 145.952, 146.949, 147.944, 148.942, 149.938, 150.937,
1699 151.935, 152.935, 153.933, 154.933, 155.931, 156.932, 157.93 , 158.931,
1700 159.929, 160.93 , 161.929, 162.93 , 163.929, 164.931, 165.93 , 166.932,
1701 167.932, 168.935, 169.935, 170.938, 171.939, 172.942, 173.944, 174.948,
1702 175.95 , 176.954};
1703
1704 G4double ErS[34] =
1705 {86, 75, 64, 54, 43, 51, 11, 32, 30, 12, 6, 50, 80, 90, 11, 5, 50, 10, 4, 6,
1706 4, 4, 3, 3, 3, 3, 3, 3, 5, 21, 32, 43, 43, 64};
1707
1708 G4double ErW[34] =
1709 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.14, 0, 1.61, 0,
1710 33.61, 22.93, 26.78, 0, 14.93, 0, 0, 0, 0, 0, 0, 0};
1711
1712 AddElement("Er", 68, 34, *ErN , *ErA , *ErS , *ErW);
1713
1714 // Z = 69 --------------------------------------------------------------------
1715 G4int TmN[34] =
1716 {146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160,
1717 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175,
1718 176, 177, 178, 179};
1719
1720 G4double TmA[34] =
1721 {145.966, 146.961, 147.958, 148.953, 149.95 , 150.945, 151.944, 152.942,
1722 153.941, 154.939, 155.939, 156.937, 157.937, 158.935, 159.935, 160.933,
1723 161.934, 162.933, 163.933, 164.932, 165.934, 166.933, 167.934, 168.934,
1724 169.936, 170.936, 171.938, 172.94 , 173.942, 174.944, 175.947, 176.949,
1725 177.953, 178.955};
1726
1727 G4double TmS[34] =
1728 {75, 64, 75, 64, 54, 15, 32, 23, 12, 14, 60, 11, 13, 70, 33, 10, 30, 7, 20,
1729 4, 12, 3, 4, 3, 3, 3, 7, 6, 50, 50, 11, 32, 43, 54};
1730
1731 G4double TmW[34] =
1732 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100,
1733 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1734
1735 AddElement("Tm", 69, 34, *TmN , *TmA , *TmS , *TmW);
1736
1737 // Z = 70 --------------------------------------------------------------------
1738 G4int YbN[34] =
1739 {148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162,
1740 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177,
1741 178, 179, 180, 181};
1742
1743 G4double YbA[34] =
1744 {147.967, 148.963, 149.958, 150.955, 151.95 , 152.949, 153.946, 154.946,
1745 155.943, 156.943, 157.94 , 158.94 , 159.938, 160.938, 161.936, 162.936,
1746 163.935, 164.935, 165.934, 166.935, 167.934, 168.935, 169.935, 170.936,
1747 171.936, 172.938, 173.939, 174.941, 175.943, 176.945, 177.947, 178.95,
1748 179.952, 180.956};
1749
1750 G4double YbS[34] =
1751 {86, 75, 64, 34, 38, 32, 11, 32, 40, 60, 11, 10, 22, 24, 22, 11, 11, 22, 9,
1752 5, 5, 5, 3, 3, 30, 30, 30, 30, 3, 3, 11, 32, 43, 43};
1753
1754 G4double YbW[34] =
1755 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.13, 0, 3.04,
1756 14.28, 21.83, 16.13, 31.83, 0, 12.76, 0, 0, 0, 0, 0};
1757
1758 AddElement("Yb", 70, 34, *YbN , *YbA , *YbS , *YbW);
1759
1760 // Z = 71 --------------------------------------------------------------------
1761 G4int LuN[35] =
1762 {150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164,
1763 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179,
1764 180, 181, 182, 183, 184};
1765
1766 G4double LuA[35] =
1767 {149.973, 150.967, 151.964, 152.959, 153.957, 154.954, 155.953, 156.95 ,
1768 157.949, 158.947, 159.946, 160.944, 161.943, 162.941, 163.941, 164.94 ,
1769 165.94 , 166.938, 167.939, 168.938, 169.938, 170.938, 171.939, 172.939,
1770 173.94 , 174.941, 175.943, 176.944, 177.946, 178.947, 179.95 , 180.952,
1771 181.955, 182.958, 183.961};
1772
1773 G4double LuS[35] =
1774 {75, 65, 75, 64, 54, 14, 32, 23, 13, 50, 25, 26, 24, 24, 13, 90, 17, 11, 90,
1775 6, 20, 3, 4, 3, 30, 28, 28, 28, 3, 6, 80, 32, 32, 32, 43};
1776
1777 G4double LuW[35] =
1778 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1779 97.41, 2.59, 0, 0, 0, 0, 0, 0, 0, 0};
1780
1781 AddElement("Lu", 71, 35, *LuN , *LuA , *LuS , *LuW);
1782
1783 // Z = 72 --------------------------------------------------------------------
1784 G4int HfN[33] =
1785 {154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168,
1786 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183,
1787 184, 185, 186};
1788
1789 G4double HfA[33] =
1790 {153.964, 154.963, 155.959, 156.958, 157.955, 158.954, 159.951, 160.95 ,
1791 161.947, 162.947, 163.944, 164.945, 165.942, 166.943, 167.941, 168.941,
1792 169.94 , 170.94 , 171.939, 172.941, 173.94 , 174.942, 175.941, 176.943,
1793 177.944, 178.946, 179.947, 180.949, 181.951, 182.954, 183.955, 184.959,
1794 185.961};
1795
1796 G4double HfS[33] =
1797 {75, 64, 38, 32, 11, 32, 40, 80, 12, 34, 21, 40, 32, 22, 11, 90, 21, 21, 50,
1798 11, 3, 3, 29, 27, 27, 27, 27, 28, 7, 30, 40, 32, 32};
1799
1800 G4double HfW[33] =
1801 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.16, 0, 5.26,
1802 18.6, 27.28, 13.62, 35.08, 0, 0, 0, 0, 0, 0};
1803
1804 AddElement("Hf", 72, 33, *HfN , *HfA , *HfS , *HfW);
1805
1806 // Z = 73 --------------------------------------------------------------------
1807 G4int TaN[33] =
1808 {156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170,
1809 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185,
1810 186, 187, 188};
1811
1812 G4double TaA[33] =
1813 {155.972, 156.968, 157.966, 158.963, 159.961, 160.958, 161.957, 162.954,
1814 163.954, 164.951, 165.95 , 166.948, 167.948, 168.946, 169.946, 170.944,
1815 171.945, 172.944, 173.944, 174.944, 175.945, 176.944, 177.946, 178.946,
1816 179.947, 180.948, 181.95 , 182.951, 183.954, 184.956, 185.959, 186.96 ,
1817 187.964};
1818
1819 G4double TaS[33] =
1820 {64, 64, 54, 13, 33, 60, 14, 80, 43, 24, 32, 46, 39, 22, 21, 22, 20, 24, 90,
1821 11, 11, 4, 11, 6, 3, 3, 3, 3, 28, 15, 60, 32, 32};
1822
1823 G4double TaW[33] =
1824 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1825 0.012, 99.988, 0, 0, 0, 0, 0, 0, 0};
1826
1827 AddElement("Ta", 73, 33, *TaN , *TaA , *TaS , *TaW);
1828
1829 // Z = 74 --------------------------------------------------------------------
1830 G4int WN[33] =
1831 {158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172,
1832 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187,
1833 188, 189, 190};
1834
1835 G4double WA[33] =
1836 {157.974, 158.972, 159.968, 160.967, 161.963, 162.963, 163.959, 164.958,
1837 165.955, 166.955, 167.952, 168.952, 169.949, 170.949, 171.947, 172.948,
1838 173.946, 174.947, 175.946, 176.947, 177.946, 178.947, 179.947, 180.948,
1839 181.948, 182.95 , 183.951, 184.953, 185.954, 186.957, 187.958, 188.962,
1840 189.963};
1841
1842 G4double WS[33] =
1843 {75, 64, 38, 33, 11, 33, 40, 90, 13, 33, 21, 34, 51, 30, 29, 40, 32, 21, 21,
1844 32, 11, 17, 5, 6, 3, 29, 29, 30, 3, 3, 4, 21, 24};
1845
1846 G4double WW[33] =
1847 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.12, 0,
1848 26.5, 14.31, 30.64, 0, 28.43, 0, 0, 0, 0};
1849
1850 AddElement("W", 74, 33, *WN , *WA , *WS , *WW);
1851
1852 // Z = 75 --------------------------------------------------------------------
1853 G4int ReN[33] =
1854 {160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174,
1855 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189,
1856 190, 191, 192};
1857
1858 G4double ReA[33] =
1859 {159.981, 160.978, 161.976, 162.972, 163.97 , 164.967, 165.966, 166.963,
1860 167.962, 168.959, 169.958, 170.956, 171.955, 172.953, 173.953, 174.951,
1861 175.952, 176.95 , 177.951, 178.95 , 179.951, 180.95 , 181.951, 182.951,
1862 183.953, 184.953, 185.955, 186.956, 187.958, 188.959, 189.962, 190.963,
1863 191.966};
1864
1865 G4double ReS[33] =
1866 {64, 64, 55, 12, 33, 80, 15, 14, 43, 22, 43, 37, 33, 48, 44, 48, 21, 21, 22,
1867 60, 40, 15, 11, 9, 6, 30, 3, 30, 30, 9, 23, 11, 21};
1868
1869 G4double ReW[33] =
1870 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1871 37.4, 0, 62.6, 0, 0, 0, 0, 0};
1872
1873 AddElement("Re", 75, 33, *ReN , *ReA , *ReS , *ReW);
1874
1875 // Z = 76 --------------------------------------------------------------------
1876 G4int OsN[35] =
1877 {162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176,
1878 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191,
1879 192, 193, 194, 195, 196};
1880
1881 G4double OsA[35] =
1882 {161.984, 162.982, 163.978, 164.976, 165.973, 166.972, 167.968, 168.967,
1883 169.964, 170.963, 171.96 , 172.96 , 173.957, 174.957, 175.955, 176.955,
1884 177.953, 178.954, 179.952, 180.953, 181.952, 182.953, 183.952, 184.954,
1885 185.954, 186.956, 187.956, 188.958, 189.958, 190.961, 191.961, 192.964,
1886 193.965, 194.968, 195.97};
1887
1888 G4double OsS[35] =
1889 {75, 64, 38, 33, 11, 33, 40, 90, 14, 33, 21, 33, 50, 32, 22, 30, 22, 25, 20,
1890 22, 27, 11, 3, 3, 3, 30, 30, 30, 3, 3, 4, 4, 4, 54, 40};
1891
1892 G4double OsW[35] =
1893 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.02, 0,
1894 1.59, 1.96, 13.24, 16.15, 26.26, 0, 40.78, 0, 0, 0, 0};
1895
1896 AddElement("Os", 76, 35, *OsN , *OsA , *OsS , *OsW);
1897
1898 // Z = 77 --------------------------------------------------------------------
1899 G4int IrN[35] =
1900 {165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179,
1901 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194,
1902 195, 196, 197, 198, 199};
1903
1904 G4double IrA[35] =
1905 {164.988, 165.986, 166.982, 167.98 , 168.976, 169.975, 170.972, 171.971,
1906 172.968, 173.967, 174.964, 175.964, 176.961, 177.961, 178.959, 179.959,
1907 180.958, 181.958, 182.957, 183.957, 184.957, 185.958, 186.957, 187.959,
1908 188.959, 189.961, 190.961, 191.963, 192.963, 193.965, 194.966, 195.968,
1909 196.97 , 197.972, 198.974};
1910
1911 G4double IrS[35] =
1912 {43, 55, 11, 35, 10, 16, 14, 43, 25, 43, 37, 33, 48, 39, 43, 21, 23, 15, 15,
1913 29, 21, 22, 7, 8, 14, 21, 3, 3, 3, 3, 3, 40, 22, 21, 40};
1914
1915 G4double IrW[35] =
1916 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1917 37.3, 0, 62.7, 0, 0, 0, 0, 0, 0};
1918
1919 AddElement("Ir", 77, 35, *IrN , *IrA , *IrS , *IrW);
1920
1921 // Z = 78 --------------------------------------------------------------------
1922 G4int PtN[35] =
1923 {168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182,
1924 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197,
1925 198, 199, 200, 201, 202};
1926
1927 G4double PtA[35] =
1928 {167.988, 168.986, 169.982, 170.981, 171.977, 172.976, 173.973, 174.972,
1929 175.969, 176.968, 177.966, 178.965, 179.963, 180.963, 181.961, 182.962,
1930 183.96 , 184.961, 185.959, 186.961, 187.959, 188.961, 189.96 , 190.962,
1931 191.961, 192.963, 193.963, 194.965, 195.965, 196.967, 197.968, 198.971,
1932 199.971, 200.975, 201.976};
1933
1934 G4double PtS[35] =
1935 {38, 34, 11, 34, 40, 11, 14, 33, 21, 33, 50, 32, 22, 30, 22, 25, 20, 22, 30,
1936 20, 6, 12, 7, 5, 4, 3, 3, 3, 3, 3, 4, 5, 22, 50, 32};
1937
1938 G4double PtW[35] =
1939 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.014, 0,
1940 0.782, 0, 32.967, 33.832, 25.242, 0, 7.163, 0, 0, 0, 0};
1941
1942 AddElement("Pt", 78, 35, *PtN , *PtA , *PtS , *PtW);
1943
1944 // Z = 79 --------------------------------------------------------------------
1945 G4int AuN[35] =
1946 {171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185,
1947 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200,
1948 201, 202, 203, 204, 205};
1949
1950 G4double AuA[35] =
1951 {170.992, 171.99 , 172.986, 173.985, 174.982, 175.98 , 176.977, 177.976,
1952 178.973, 179.972, 180.97 , 181.97 , 182.968, 183.967, 184.966, 185.966,
1953 186.965, 187.965, 188.964, 189.965, 190.964, 191.965, 192.964, 193.965,
1954 194.965, 195.967, 196.967, 197.968, 198.969, 199.971, 200.972, 201.974,
1955 202.975, 203.978, 204.98 };
1956
1957 G4double AuS[35] =
1958 {27, 36, 11, 16, 26, 43, 25, 43, 37, 32, 48, 39, 43, 21, 23, 15, 16, 11, 22,
1959 17, 50, 17, 10, 12, 3, 4, 3, 3, 3, 60, 5, 18, 5, 22, 32};
1960
1961 G4double AuW[35] =
1962 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1963 100, 0, 0, 0, 0, 0, 0, 0, 0};
1964
1965 AddElement("Au", 79, 35, *AuN , *AuA , *AuS , *AuW);
1966
1967 // Z = 80 --------------------------------------------------------------------
1968 G4int HgN[34] =
1969 {175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189,
1970 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204,
1971 205, 206, 207, 208};
1972
1973 G4double HgA[34] =
1974 {174.991, 175.987, 176.986, 177.982, 178.982, 179.978, 180.978, 181.975,
1975 182.975, 183.972, 184.972, 185.969, 186.97 , 187.968, 188.968, 189.966,
1976 190.967, 191.966, 192.967, 193.965, 194.967, 195.966, 196.967, 197.967,
1977 198.968, 199.968, 200.97 , 201.971, 202.973, 203.973, 204.976, 205.977,
1978 206.983, 207.986};
1979
1980 G4double HgS[34] =
1981 {34, 40, 12, 16, 33, 21, 33, 50, 32, 22, 30, 22, 26, 19, 30, 16, 90, 30, 21,
1982 25, 50, 4, 5, 3, 3, 3, 3, 3, 3, 3, 5, 22, 16, 32};
1983
1984 G4double HgW[34] =
1985 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.15, 0, 9.97,
1986 16.87, 23.1, 13.18, 29.86, 0, 6.87, 0, 0, 0, 0};
1987
1988 AddElement("Hg", 80, 34, *HgN , *HgA , *HgS , *HgW);
1989
1990 // Z = 81 --------------------------------------------------------------------
1991 G4int TlN[34] =
1992 {177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191,
1993 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206,
1994 207, 208, 209, 210};
1995
1996 G4double TlA[34] =
1997 {176.997, 177.995, 178.991, 179.99 , 180.987, 181.986, 182.983, 183.982,
1998 184.979, 185.979, 186.976, 187.976, 188.974, 189.974, 190.972, 191.972,
1999 192.971, 193.971, 194.97 , 195.971, 196.97 , 197.97 , 198.97 , 199.971,
2000 200.971, 201.972, 202.972, 203.974, 204.974, 205.976, 206.977, 207.982,
2001 208.985, 209.99 };
2002
2003 G4double TlS[34] =
2004 {24, 23, 15, 48, 41, 43, 42, 32, 43, 39, 43, 24, 37, 46, 23, 22, 27, 22, 14,
2005 15, 30, 90, 11, 7, 16, 16, 3, 3, 3, 3, 6, 3, 10, 12};
2006
2007 G4double TlW[34] =
2008 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2009 29.524, 0, 70.476, 0, 0, 0, 0, 0};
2010
2011 AddElement("Tl", 81, 34, *TlN , *TlA , *TlS , *TlW);
2012
2013 // Z = 82 --------------------------------------------------------------------
2014 G4int PbN[34] =
2015 {181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195,
2016 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210,
2017 211, 212, 213, 214};
2018
2019 G4double PbA[34] =
2020 {180.997, 181.993, 182.992, 183.988, 184.988, 185.984, 186.984, 187.981,
2021 188.981, 189.978, 190.978, 191.976, 192.976, 193.974, 194.974, 195.973,
2022 196.973, 197.972, 198.973, 199.972, 200.973, 201.972, 202.973, 203.973,
2023 204.974, 205.974, 206.976, 207.977, 208.981, 209.984, 210.989, 211.992,
2024 212.996, 214};
2025
2026 G4double PbS[34] =
2027 {17, 18, 33, 21, 33, 50, 32, 22, 29, 22, 23, 19, 20, 16, 44, 15, 11, 10, 70,
2028 14, 30, 11, 7, 3, 3, 3, 3, 3, 3, 3, 3, 29, 11, 27};
2029
2030 G4double PbW[34] =
2031 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.4, 0,
2032 24.1, 22.1, 52.4, 0, 0, 0, 0, 0, 0};
2033
2034 AddElement("Pb", 82, 34, *PbN , *PbA , *PbS , *PbW);
2035
2036 // Z = 83 --------------------------------------------------------------------
2037 G4int BiN[32] =
2038 {185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199,
2039 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214,
2040 215, 216};
2041
2042 G4double BiA[32] =
2043 {184.998, 185.996, 186.993, 187.992, 188.99 , 189.989, 190.986, 191.985,
2044 192.983, 193.983, 194.981, 195.981, 196.979, 197.979, 198.978, 199.978,
2045 200.977, 201.978, 202.977, 203.978, 204.977, 205.978, 206.978, 207.98 ,
2046 208.98 , 209.984, 210.987, 211.991, 212.994, 213.999, 215.002, 216.006};
2047
2048 G4double BiS[32] =
2049 {24, 48, 41, 32, 43, 39, 43, 24, 37, 46, 23, 22, 26, 19, 13, 10, 30, 60, 23,
2050 28, 8, 9, 4, 4, 3, 3, 6, 3, 8, 12, 10, 11};
2051
2052 G4double BiW[32] =
2053 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100,
2054 0, 0, 0, 0, 0, 0, 0};
2055
2056 AddElement("Bi", 83, 32, *BiN , *BiA , *BiS , *BiW);
2057
2058 // Z = 84 --------------------------------------------------------------------
2059 G4int PoN[29] =
2060 {190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204,
2061 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218};
2062
2063 G4double PoA[29] =
2064 {189.995, 190.995, 191.992, 192.991, 193.988, 194.988, 195.986, 196.986,
2065 197.983, 198.984, 199.982, 200.982, 201.981, 202.981, 203.98 , 204.981,
2066 205.98 , 206.982, 207.981, 208.982, 209.983, 210.987, 211.989, 212.993,
2067 213.995, 214.999, 216.002, 217.006, 218.009};
2068
2069 G4double PoS[29] =
2070 {51, 32, 22, 30, 22, 24, 19, 21, 16, 44, 15, 11, 10, 70, 14, 30, 11, 8, 3, 3,
2071 3, 3, 3, 4, 3, 3, 29, 11, 27};
2072
2073 G4double PoW[29] =
2074 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0,
2075 0, 0, 0, 0};
2076
2077 AddElement("Po", 84, 29, *PoN , *PoA , *PoS , *PoW);
2078
2079 // Z = 85 --------------------------------------------------------------------
2080 G4int AtN[31] =
2081 {193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207,
2082 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222,
2083 223};
2084
2085 G4double AtA[31] =
2086 {193.000, 193.999, 194.997, 195.996, 196.993, 197.993, 198.991, 199.99,
2087 200.988, 201.988, 202.987, 203.987, 204.986, 205.987, 206.986, 207.987,
2088 208.986, 209.987, 210.987, 211.991, 212.993, 213.996, 214.999, 216.002,
2089 217.005, 218.009, 219.011, 220.015, 221.018, 222.022, 223.025};
2090
2091 G4double AtS[31] =
2092 {43, 43, 43, 25, 38, 46, 24, 22, 26, 19, 13, 10, 30, 60, 23, 28, 9, 9, 4, 4,
2093 6, 5, 8, 5, 8, 13, 90, 12, 32, 32, 43};
2094
2095 G4double AtW[31] =
2096 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0,
2097 0, 0, 0, 0, 0, 0};
2098
2099 AddElement("At", 85, 31, *AtN , *AtA , *AtS , *AtW);
2100
2101 // Z = 86 --------------------------------------------------------------------
2102 G4int RnN[33] =
2103 {196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210,
2104 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225,
2105 226, 227, 228};
2106
2107 G4double RnA[33] =
2108 {196.002, 197.002, 197.999, 198.998, 199.996, 200.996, 201.993, 202.993,
2109 203.991, 204.992, 205.99 , 206.991, 207.99 , 208.99 , 209.99 , 210.991,
2110 211.991, 212.994, 213.995, 214.999, 216.000, 217.004, 218.006, 219.009,
2111 220.011, 221.015, 222.018, 223.022, 224.024, 225.028, 226.031, 227.035,
2112 228.038};
2113
2114 G4double RnS[33] =
2115 {22, 30, 22, 24, 19, 21, 16, 44, 15, 12, 10, 70, 14, 30, 11, 8, 4, 8, 10, 9,
2116 8, 5, 4, 3, 29, 11, 27, 32, 32, 32, 43, 45, 50};
2117
2118 G4double RnW[33] =
2119 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
2120 100, 0, 0, 0, 0, 0, 0};
2121
2122 AddElement("Rn", 86, 33, *RnN , *RnA , *RnS , *RnW);
2123
2124 // Z = 87 --------------------------------------------------------------------
2125 G4int FrN[33] =
2126 {200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214,
2127 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229,
2128 230, 231, 232};
2129
2130 G4double FrA[33] =
2131 {200.006, 201.004, 202.003, 203.001, 204.001, 204.999, 205.998, 206.997,
2132 207.997, 208.996, 209.996, 210.996, 211.996, 212.996, 213.999, 215.000,
2133 216.003, 217.005, 218.008, 219.009, 220.012, 221.014, 222.018, 223.02 ,
2134 224.023, 225.026, 226.029, 227.032, 228.036, 229.038, 230.043, 231.045,
2135 232.05 };
2136
2137 G4double FrS[33] =
2138 {25, 38, 46, 25, 22, 25, 19, 12, 80, 30, 23, 22, 28, 9, 10, 8, 14, 8, 6, 8, 5,
2139 8, 23, 29, 50, 11, 10, 10, 22, 39, 48, 56, 69};
2140
2141 G4double FrW[33] =
2142 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0,
2143 0, 0, 0, 0, 0, 0, 0, 0};
2144
2145 AddElement("Fr", 87, 33, *FrN , *FrA , *FrS , *FrW);
2146
2147 // Z = 88 --------------------------------------------------------------------
2148 G4int RaN[32] =
2149 {203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217,
2150 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232,
2151 233, 234};
2152
2153 G4double RaA[32] =
2154 {203.009, 204.006, 205.006, 206.004, 207.004, 208.002, 209.002, 210 ,
2155 211.001, 212 , 213 , 214 , 215.003, 216.004, 217.006, 218.007,
2156 219.01 , 220.011, 221.014, 222.015, 223.018, 224.02 , 225.024, 226.025,
2157 227.029, 228.031, 229.035, 230.037, 231.041, 232.044, 233.048, 234.051};
2158
2159 G4double RaS[32] =
2160 {25, 19, 22, 16, 45, 15, 13, 10, 70, 15, 30, 12, 9, 10, 10, 12, 9, 11, 8, 5,
2161 3, 29, 3, 27, 27, 27, 70, 40, 32, 39, 50, 58};
2162
2163 G4double RaW[32] =
2164 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0,
2165 0, 0, 0, 0, 0, 0, 0};
2166
2167 AddElement("Ra", 88, 32, *RaN , *RaA , *RaS , *RaW);
2168
2169 // Z = 89 --------------------------------------------------------------------
2170 G4int AcN[30] =
2171 {207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
2172 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236};
2173
2174 G4double AcA[30] =
2175 {207.012, 208.011, 209.01 , 210.009, 211.008, 212.008, 213.007, 214.007,
2176 215.006, 216.009, 217.009, 218.012, 219.012, 220.015, 221.016, 222.018,
2177 223.019, 224.022, 225.023, 226.026, 227.028, 228.031, 229.033, 230.036,
2178 231.039, 232.042, 233.045, 234.048, 235.051, 236.055};
2179
2180 G4double AcS[30] =
2181 {25, 22, 26, 20, 14, 10, 60, 60, 60, 29, 14, 50, 50, 60, 50, 6, 8, 5, 8, 4,
2182 29, 28, 50, 11, 11, 11, 32, 43, 45, 54};
2183
2184 G4double AcW[30] =
2185 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0,
2186 0, 0, 0, 0, 0};
2187
2188 AddElement("Ac", 89, 30, *AcN , *AcA , *AcS , *AcW);
2189
2190 // Z = 90 --------------------------------------------------------------------
2191 G4int ThN[29] =
2192 {210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224,
2193 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238};
2194
2195 G4double ThA[29] =
2196 {210.015, 211.015, 212.013, 213.013, 214.011, 215.012, 216.011, 217.013,
2197 218.013, 219.016, 220.016, 221.018, 222.018, 223.021, 224.021, 225.024,
2198 226.025, 227.028, 228.029, 229.032, 230.033, 231.036, 232.038, 233.042,
2199 234.044, 235.048, 236.05 , 237.054, 238.056};
2200
2201 G4double ThS[29] =
2202 {17, 45, 15, 14, 10, 70, 17, 30, 15, 50, 24, 11, 14, 10, 13, 8, 5, 3, 29, 3,
2203 22, 22, 22, 22, 4, 50, 32, 39, 39};
2204
2205 G4double ThW[29] =
2206 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0,
2207 0, 0, 0, 0};
2208
2209 AddElement("Th", 90, 29, *ThN , *ThA , *ThS , *ThW);
2210
2211 // Z = 91 --------------------------------------------------------------------
2212 G4int PaN[28] =
2213 {213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227,
2214 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240};
2215
2216 G4double PaA[28] =
2217 {213.021, 214.021, 215.019, 216.019, 217.018, 218.02 , 219.02 , 220.022,
2218 221.022, 222.024, 223.024, 224.026, 225.026, 226.028, 227.029, 228.031,
2219 229.032, 230.035, 231.036, 232.039, 233.04 , 234.043, 235.045, 236.049,
2220 237.051, 238.054, 239.057, 240.061};
2221
2222 G4double PaS[28] =
2223 {27, 21, 15, 11, 80, 80, 80, 60, 60, 80, 80, 60, 80, 13, 8, 5, 10, 4, 28, 8,
2224 24, 5, 50, 21, 11, 60, 32, 32};
2225
2226 G4double PaW[28] =
2227 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0,
2228 0, 0, 0};
2229
2230 AddElement("Pa", 91, 28, *PaN , *PaA , *PaS , *PaW);
2231
2232 // Z = 92 --------------------------------------------------------------------
2233 G4int UN[25] =
2234 {218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232,
2235 233, 234, 235, 236, 237, 238, 239, 240, 241, 242};
2236
2237 G4double UA[25] =
2238 {218.023, 219.025, 220.025, 221.026, 222.026, 223.028, 224.028, 225.029,
2239 226.029, 227.031, 228.031, 229.033, 230.034, 231.036, 232.037, 233.04 ,
2240 234.041, 235.044, 236.046, 237.049, 238.051, 239.054, 240.057, 241.06 ,
2241 242.063};
2242
2243 G4double US[25] =
2244 {10, 90, 22, 11, 11, 80, 27, 50, 20, 18, 17, 9, 5, 4, 29, 3, 21, 21, 21, 21,
2245 21, 21, 6, 32, 22};
2246
2247 G4double UW[25] =
2248 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0055, 0.72, 0, 0, 99.2745,
2249 0, 0, 0, 0};
2250
2251 AddElement("U", 92, 25, *UN , *UA , *US , *UW);
2252
2253 // Z = 93 --------------------------------------------------------------------
2254 G4int NpN[20] =
2255 {225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239,
2256 240, 241, 242, 243, 244};
2257
2258 G4double NpA[20] =
2259 {225.034, 226.035, 227.035, 228.036, 229.036, 230.038, 231.038, 232.04 ,
2260 233.041, 234.043, 235.044, 236.047, 237.048, 238.051, 239.053, 240.056,
2261 241.058, 242.062, 243.064, 244.068};
2262
2263 G4double NpS[20] =
2264 {80, 10, 80, 21, 90, 60, 50, 11, 60, 9, 23, 50, 21, 21, 23, 16, 80, 22, 30,
2265 32};
2266
2267 G4double NpW[20] =
2268 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0};
2269
2270 AddElement("Np", 93, 20, *NpN , *NpA , *NpS , *NpW);
2271
2272 // Z = 94 --------------------------------------------------------------------
2273 G4int PuN[20] =
2274 {228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242,
2275 243, 244, 245, 246, 247};
2276
2277 G4double PuA[20] =
2278 {228.039, 229.04 , 230.04 , 231.041, 232.041, 233.043, 234.043, 235.045,
2279 236.046, 237.048, 238.05 , 239.052, 240.054, 241.057, 242.059, 243.062,
2280 244.064, 245.068, 246.07 , 247.074};
2281
2282 G4double PuS[20] =
2283 {30, 80, 26, 11, 20, 50, 8, 22, 29, 25, 21, 21, 21, 21, 21, 3, 5, 15, 16, 32};
2284
2285 G4double PuW[20] =
2286 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
2287
2288 AddElement("Pu", 94, 20, *PuN , *PuA , *PuS , *PuW);
2289
2290 // Z = 95 --------------------------------------------------------------------
2291 G4int AmN[19] =
2292 {231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245,
2293 246, 247, 248, 249};
2294
2295 G4double AmA[19] =
2296 {231.046, 232.047, 233.046, 234.048, 235.048, 236.05 , 237.05 , 238.052,
2297 239.053, 240.055, 241.057, 242.06 , 243.061, 244.064, 245.066, 246.07 ,
2298 247.072, 248.076, 249.078};
2299
2300 G4double AmS[19] =
2301 {32, 32, 23, 22, 22, 11, 60, 50, 3, 15, 21, 21, 23, 23, 4, 20, 11, 22, 32};
2302
2303 G4double AmW[19] =
2304 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0};
2305
2306 AddElement("Am", 95, 19, *AmN , *AmA , *AmS , *AmW);
2307
2308 // Z = 96 --------------------------------------------------------------------
2309 G4int CmN[20] =
2310 {233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247,
2311 248, 249, 250, 251, 252};
2312
2313 G4double CmA[20] =
2314 {233.051, 234.05 , 235.052, 236.051, 237.053, 238.053, 239.055, 240.056,
2315 241.058, 242.059, 243.061, 244.063, 245.065, 246.067, 247.07 , 248.072,
2316 249.076, 250.078, 251.082, 252.085};
2317
2318 G4double CmS[20] =
2319 {43, 32, 24, 22, 22, 40, 11, 29, 24, 21, 24, 21, 29, 24, 5, 5, 5, 12, 24, 32};
2320
2321 G4double CmW[20] =
2322 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
2323
2324 AddElement("Cm", 96, 20, *CmN , *CmA , *CmS , *CmW);
2325
2326 // Z = 97 --------------------------------------------------------------------
2327 G4int BkN[20] =
2328 {235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249,
2329 250, 251, 252, 253, 254};
2330
2331 G4double BkA[20] =
2332 {235.057, 236.057, 237.057, 238.058, 239.058, 240.06 , 241.06 , 242.062,
2333 243.063, 244.065, 245.066, 246.069, 247.07 , 248.073, 249.075, 250.078,
2334 251.081, 252.084, 253.087, 254.091};
2335
2336 G4double BkS[20] =
2337 {43, 43, 32, 31, 31, 16, 22, 22, 5, 16, 26, 60, 6, 80, 3, 4, 12, 22, 39, 32};
2338
2339 G4double BkW[20] =
2340 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0};
2341
2342 AddElement("Bk", 97, 20, *BkN , *BkA , *BkS , *BkW);
2343
2344 // Z = 98 --------------------------------------------------------------------
2345 G4int CfN[20] =
2346 {237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251,
2347 252, 253, 254, 255, 256};
2348
2349 G4double CfA[20] =
2350 {237.062, 238.061, 239.063, 240.062, 241.064, 242.064, 243.065, 244.066,
2351 245.068, 246.069, 247.071, 248.072, 249.075, 250.076, 251.08 , 252.082,
2352 253.085, 254.087, 255.091, 256.093};
2353
2354 G4double CfS[20] =
2355 {54, 43, 25, 22, 27, 40, 15, 4, 11, 24, 9, 6, 3, 24, 5, 5, 7, 13, 22, 32};
2356
2357 G4double CfW[20] =
2358 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
2359
2360 AddElement("Cf", 98, 20, *CfN , *CfA , *CfS , *CfW);
2361
2362 // Z = 99 --------------------------------------------------------------------
2363 G4int EsN[18] =
2364 {240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254,
2365 255, 256, 257};
2366
2367 G4double EsA[18] =
2368 {240.069, 241.069, 242.07 , 243.07 , 244.071, 245.071, 246.073, 247.074,
2369 248.075, 249.076, 250.079, 251.08 , 252.083, 253.085, 254.088, 255.09 ,
2370 256.094, 257.096};
2371
2372 G4double EsS[18] =
2373 {43, 32, 35, 31, 20, 22, 24, 30, 60, 30, 11, 7, 50, 3, 5, 12, 11, 44};
2374
2375 G4double EsW[18] =
2376 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
2377
2378 AddElement("Es", 99, 18, *EsN , *EsA , *EsS , *EsW);
2379
2380 // Z = 100 -------------------------------------------------------------------
2381 G4int FmN[18] =
2382 {242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256,
2383 257, 258, 259};
2384
2385 G4double FmA[18] =
2386 {242.073, 243.075, 244.074, 245.075, 246.075, 247.077, 248.077, 249.079,
2387 250.08 , 251.082, 252.082, 253.085, 254.087, 255.09 , 256.092, 257.095,
2388 258.097, 259.101};
2389
2390 G4double FmS[18] =
2391 {43, 25, 31, 30, 40, 16, 13, 15, 13, 9, 6, 5, 3, 5, 8, 7, 22, 30};
2392
2393 G4double FmW[18] =
2394 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0};
2395
2396 AddElement("Fm", 100, 18, *FmN , *FmA , *FmS , *FmW);
2397
2398 // Z = 101 -------------------------------------------------------------------
2399 G4int MdN[16] =
2400 {245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259,
2401 260};
2402
2403 G4double MdA[16] =
2404 {245.081, 246.082, 247.082, 248.083, 249.083, 250.084, 251.085, 252.087,
2405 253.087, 254.09 , 255.091, 256.094, 257.096, 258.098, 259.101, 260.104};
2406
2407 G4double MdS[16] =
2408 {40, 41, 39, 26, 24, 32, 22, 21, 22, 11, 8, 60, 4, 5, 22, 34};
2409
2410 G4double MdW[16] =
2411 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0};
2412
2413 AddElement("Md", 101, 16, *MdN , *MdA , *MdS , *MdW);
2414
2415 // Z = 102 -------------------------------------------------------------------
2416 G4int NoN[14] =
2417 {249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262};
2418
2419 G4double NoA[14] =
2420 {249.088, 250.087, 251.089, 252.089, 253.091, 254.091, 255.093, 256.094,
2421 257.097, 258.098, 259.101, 260.103, 261.106, 262.108};
2422
2423 G4double NoS[14] =
2424 {37, 22, 19, 14, 26, 19, 13, 9, 30, 22, 11, 22, 32, 58};
2425
2426 G4double NoW[14] =
2427 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
2428
2429 AddElement("No", 102, 14, *NoN , *NoA , *NoS , *NoW);
2430
2431 // Z = 103 -------------------------------------------------------------------
2432 G4int LrN[13] =
2433 {251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263};
2434
2435 G4double LrA[13] =
2436 {251.094, 252.095, 253.095, 254.097, 255.097, 256.099, 257.1 , 258.102,
2437 259.103, 260.106, 261.107, 262.11 , 263.111};
2438
2439 G4double LrS[13] =
2440 {32, 32, 24, 36, 22, 24, 22, 11, 80, 12, 22, 32, 39};
2441
2442 G4double LrW[13] =
2443 {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0};
2444
2445 AddElement("Lr", 103, 13, *LrN , *LrA , *LrS , *LrW);
2446
2447 // Z = 104 -------------------------------------------------------------------
2448 G4int RfN[12] =
2449 {253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264};
2450
2451 G4double RfA[12] =
2452 {253.101, 254.1 , 255.101, 256.101, 257.103, 258.104, 259.106, 260.106,
2453 261.109, 262.11 , 263.113, 264.114};
2454
2455 G4double RfS[12] =
2456 {49, 31, 22, 29, 29, 22, 80, 22, 11, 30, 20, 48};
2457
2458 G4double RfW[12] =
2459 {0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
2460
2461 AddElement("Rf", 104, 12, *RfN , *RfA , *RfS , *RfW);
2462
2463 // Z = 105 -------------------------------------------------------------------
2464 G4int DbN[11] =
2465 {255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265};
2466
2467 G4double DbA[11] =
2468 {255.107, 256.108, 257.108, 258.109, 259.11 , 260.111, 261.112, 262.114,
2469 263.115, 264.117, 265.119};
2470
2471 G4double DbS[11] =
2472 {45, 39, 25, 37, 31, 25, 25, 20, 18, 25, 30};
2473
2474 G4double DbW[11] =
2475 {0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
2476
2477 AddElement("Db", 105, 11, *DbN , *DbA , *DbS , *DbW);
2478
2479 // Z = 106 -------------------------------------------------------------------
2480 G4int SgN[9] =
2481 {258, 259, 260, 261, 262, 263, 264, 265, 266};
2482
2483 G4double SgA[9] =
2484 {258.113, 259.115, 260.114, 261.116, 262.116, 263.118, 264.119, 265.121,
2485 266.122};
2486
2487 G4double SgS[9] =
2488 {45, 23, 40, 30, 30, 13, 30, 15, 31};
2489
2490 G4double SgW[9] =
2491 {0, 0, 0, 0, 0, 0, 0, 0, 100};
2492
2493 AddElement("Sg", 106, 9, *SgN , *SgA , *SgS , *SgW);
2494
2495 // Z = 107 -------------------------------------------------------------------
2496 G4int BhN[8] =
2497 {260, 261, 262, 263, 264, 265, 266, 267};
2498
2499 G4double BhA[8] =
2500 {260.122, 261.122, 262.123, 263.123, 264.125, 265.125, 266.127, 267.128};
2501
2502 G4double BhS[8] =
2503 {66, 26, 40, 45, 30, 41, 38, 37};
2504
2505 G4double BhW[8] =
2506 {0, 0, 0, 0, 100, 0, 0, 0};
2507
2508 AddElement("Bh", 107, 8, *BhN , *BhA , *BhS , *BhW);
2509
2510 for(G4int i=0; i<maxNumElements; ++i) {elmNames.push_back(elmSymbol[i]);}
2511
2512 if(0<verbose) {
2513 G4cout << "G4NistElementBuilder: " << maxNumElements-1 << " Elements "
2514 << index << " Isotopes"
2515 << G4endl;
2516 }
2517}
2518
2519//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
imax
static const G4int imax
Definition: G4ElectroNuclearCrossSection.cc:98
G4ElementTable
std::vector< G4Element * > G4ElementTable
Definition: G4ElementTable.hh:42
maxNumElements
const G4int maxNumElements
Definition: G4NistElementBuilder.hh:65
maxAbundance
const G4int maxAbundance
Definition: G4NistElementBuilder.hh:66
HA
static const G4double HA[nHA]
Definition: G4PhotoNuclearCrossSection.cc:71
mole
static constexpr double mole
Definition: G4SIunits.hh:279
eV
static constexpr double eV
Definition: G4SIunits.hh:201
g
static constexpr double g
Definition: G4SIunits.hh:168
GeV
static constexpr double GeV
Definition: G4SIunits.hh:203
G4MUTEX_INITIALIZER
#define G4MUTEX_INITIALIZER
Definition: G4Threading.hh:85
G4MUTEXLOCK
#define G4MUTEXLOCK(mutex)
Definition: G4Threading.hh:251
G4MUTEXUNLOCK
#define G4MUTEXUNLOCK(mutex)
Definition: G4Threading.hh:254
G4Mutex
std::mutex G4Mutex
Definition: G4Threading.hh:81
G4double
double G4double
Definition: G4Types.hh:83
G4bool
bool G4bool
Definition: G4Types.hh:86
G4int
int G4int
Definition: G4Types.hh:85
Z
const G4int Z[17]
Definition: G4WaterStopping.cc:51
A
const G4double A[17]
Definition: G4WaterStopping.cc:53
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:397
G4Element::AddIsotope
void AddIsotope(G4Isotope *isotope, G4double RelativeAbundance)
Definition: G4Element.cc:151
G4Element::GetIndex
size_t GetIndex() const
Definition: G4Element.hh:182
G4Element::SetNaturalAbundanceFlag
void SetNaturalAbundanceFlag(G4bool)
Definition: G4Element.hh:268
G4NistElementBuilder::elmIndex
G4int elmIndex[maxNumElements]
Definition: G4NistElementBuilder.hh:145
G4NistElementBuilder::atomicMass
G4double atomicMass[maxNumElements]
Definition: G4NistElementBuilder.hh:139
G4NistElementBuilder::GetAtomicMass
G4double GetAtomicMass(G4int Z, G4int N) const
Definition: G4NistElementBuilder.hh:182
G4NistElementBuilder::BuildElement
G4Element * BuildElement(G4int Z)
Definition: G4NistElementBuilder.cc:156
G4NistElementBuilder::GetAtomicMassAmu
G4double GetAtomicMassAmu(const G4String &symb) const
Definition: G4NistElementBuilder.cc:95
G4NistElementBuilder::sigMass
G4double sigMass[maxAbundance]
Definition: G4NistElementBuilder.hh:148
G4NistElementBuilder::FindOrBuildElement
G4Element * FindOrBuildElement(G4int Z, G4bool buildIsotopes=true)
Definition: G4NistElementBuilder.cc:105
G4NistElementBuilder::relAbundance
G4double relAbundance[maxAbundance]
Definition: G4NistElementBuilder.hh:149
G4NistElementBuilder::elmSymbol
G4String elmSymbol[maxNumElements]
Definition: G4NistElementBuilder.hh:138
G4NistElementBuilder::GetZ
G4int GetZ(const G4String &symb) const
Definition: G4NistElementBuilder.cc:85
G4NistElementBuilder::FindElement
G4Element * FindElement(G4int Z) const
Definition: G4NistElementBuilder.hh:255
G4NistElementBuilder::elmNames
std::vector< G4String > elmNames
Definition: G4NistElementBuilder.hh:154
G4NistElementBuilder::idxIsotopes
G4int idxIsotopes[maxNumElements]
Definition: G4NistElementBuilder.hh:143
G4NistElementBuilder::nIsotopes
G4int nIsotopes[maxNumElements]
Definition: G4NistElementBuilder.hh:141
G4NistElementBuilder::AddElement
void AddElement(const G4String &symbol, G4int Z, G4int NumberOfIsotopes, const G4int &N, const G4double &A, const G4double &sigmaA, const G4double &W)
Definition: G4NistElementBuilder.cc:237
G4NistElementBuilder::bindingEnergy
G4double bindingEnergy[maxNumElements]
Definition: G4NistElementBuilder.hh:140
G4NistElementBuilder::massIsotopes
G4double massIsotopes[maxAbundance]
Definition: G4NistElementBuilder.hh:147
G4NistElementBuilder::PrintElement
void PrintElement(G4int Z) const
Definition: G4NistElementBuilder.cc:203
G4NistElementBuilder::nFirstIsotope
G4int nFirstIsotope[maxNumElements]
Definition: G4NistElementBuilder.hh:142
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition: G4InuclParticleNames.hh:76
anonymous_namespace{paraMaker.cc}::CA
static const G4double CA[5]
Definition: paraMaker.cc:45
source.hepunit.electron_mass_c2
float electron_mass_c2
Definition: hepunit.py:273
source.hepunit.proton_mass_c2
float proton_mass_c2
Definition: hepunit.py:274
source.hepunit.amu_c2
float amu_c2
Definition: hepunit.py:276
NS
#define NS(x)
Definition: xmltok.cc:1624
Generated by doxygen 1.9.3