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