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