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