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