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