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