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