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