| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 25
26 // ------------------------------------------- 26 // -------------------------------------------------------------------
27 // 27 //
28 // GEANT4 Class file 28 // GEANT4 Class file
29 // 29 //
30 // File name: G4NistElementBuilder 30 // File name: G4NistElementBuilder
31 // 31 //
32 // Author: Vladimir Ivanchenko 32 // Author: Vladimir Ivanchenko
33 // 33 //
34 // Creation date: 23.12.2004 34 // Creation date: 23.12.2004
35 // 35 //
36 // Modifications: 36 // Modifications:
37 // 02.05.2006 Subtract mass of atomic electron 37 // 02.05.2006 Subtract mass of atomic electrons from NIST mass (VI)
38 // 11.05.2006 Do not subtract mass of atomic e 38 // 11.05.2006 Do not subtract mass of atomic electrons from NIST mass (VI)
39 // 17.10.2006 Add natiral abundances flag to e 39 // 17.10.2006 Add natiral abundances flag to element and
40 // use G4 units for isotope mass ve 40 // use G4 units for isotope mass vector (VI)
41 // 10.05.2007 Add protection agains Z>101 (VI) 41 // 10.05.2007 Add protection agains Z>101 (VI)
42 // 26.07.2007 Create one and only one Nist ele 42 // 26.07.2007 Create one and only one Nist element with given Z and
43 // allow users to create there own 43 // allow users to create there own elements with the same Z (VI)
44 // 44 //
45 // ------------------------------------------- 45 // -------------------------------------------------------------------
46 // 46 //
47 // Class Description: 47 // Class Description:
48 // 48 //
49 // Element data from the NIST DB on Atomic Wei 49 // Element data from the NIST DB on Atomic Weights and Isotope Compositions
50 // http://physics.nist.gov/PhysRefData/Composi 50 // http://physics.nist.gov/PhysRefData/Compositions/index.html
51 51
52 #include "G4NistElementBuilder.hh" 52 #include "G4NistElementBuilder.hh"
53 53
54 #include "G4AutoLock.hh" 54 #include "G4AutoLock.hh"
55 #include "G4PhysicalConstants.hh" 55 #include "G4PhysicalConstants.hh"
56 #include "G4SystemOfUnits.hh" 56 #include "G4SystemOfUnits.hh"
57 57
58 #include <sstream> 58 #include <sstream>
59 59
60 namespace 60 namespace
61 { 61 {
62 G4Mutex nistElementMutex = G4MUTEX_INITIALIZER 62 G4Mutex nistElementMutex = G4MUTEX_INITIALIZER;
63 } 63 }
64 64
65 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 66
67 G4NistElementBuilder::G4NistElementBuilder(G4i 67 G4NistElementBuilder::G4NistElementBuilder(G4int vb) : verbose(vb)
68 { 68 {
69 atomicMass[0] = 0.0; <<
70 nFirstIsotope[0] = 0; 69 nFirstIsotope[0] = 0;
71 nIsotopes[0] = 0; 70 nIsotopes[0] = 0;
72 idxIsotopes[0] = 0; << 71 relAbundance[0] = 0.0;
73 elmSymbol.push_back(" "); <<
74 Initialise(); 72 Initialise();
75 for (G4int i=0; i<maxNumElements; ++i) { 73 for (G4int i=0; i<maxNumElements; ++i) {
76 elmIndex[i] = -1; 74 elmIndex[i] = -1;
77 } 75 }
78 } 76 }
79 77
80 //....oooOO0OOooo........oooOO0OOooo........oo 78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
81 79
82 G4int G4NistElementBuilder::GetZ(const G4Strin 80 G4int G4NistElementBuilder::GetZ(const G4String& name) const
83 { 81 {
84 G4int Z = maxNumElements; 82 G4int Z = maxNumElements;
85 // Loop checking, 07-Aug-2015, Vladimir Ivan 83 // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
86 do { 84 do {
87 --Z; 85 --Z;
88 } while (Z > 0 && elmSymbol[Z] != name); 86 } while (Z > 0 && elmSymbol[Z] != name);
89 return Z; 87 return Z;
90 } 88 }
91 89
92 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
93 91
94 G4double G4NistElementBuilder::GetAtomicMassAm 92 G4double G4NistElementBuilder::GetAtomicMassAmu(const G4String& name) const
95 { 93 {
96 G4int Z = maxNumElements; 94 G4int Z = maxNumElements;
97 // Loop checking, 07-Aug-2015, Vladimir Ivan 95 // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
98 do { 96 do {
99 --Z; 97 --Z;
100 } while (Z > 0 && elmSymbol[Z] != name); 98 } while (Z > 0 && elmSymbol[Z] != name);
101 return GetAtomicMassAmu(Z); 99 return GetAtomicMassAmu(Z);
102 } 100 }
103 101
104 //....oooOO0OOooo........oooOO0OOooo........oo 102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
105 103
106 G4Element* G4NistElementBuilder::FindOrBuildEl 104 G4Element* G4NistElementBuilder::FindOrBuildElement(G4int Z, G4bool)
107 { 105 {
108 G4Element* anElement = FindElement(Z); 106 G4Element* anElement = FindElement(Z);
109 if (anElement == nullptr && Z > 0 && Z < max 107 if (anElement == nullptr && Z > 0 && Z < maxNumElements) {
110 anElement = BuildElement(Z); 108 anElement = BuildElement(Z);
111 } 109 }
112 return anElement; 110 return anElement;
113 } 111 }
114 112
115 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
116 114
117 G4Element* G4NistElementBuilder::FindOrBuildEl 115 G4Element* G4NistElementBuilder::FindOrBuildElement(const G4String& symb, G4bool)
118 { 116 {
119 G4Element* elm = nullptr; 117 G4Element* elm = nullptr;
120 const G4ElementTable* theElementTable = G4El 118 const G4ElementTable* theElementTable = G4Element::GetElementTable();
121 std::size_t nelm = theElementTable->size(); << 119 size_t nelm = theElementTable->size();
122 for (std::size_t i = 0; i < nelm; ++i) { << 120 for (size_t i = 0; i < nelm; ++i) {
123 if (symb == ((*theElementTable)[i])->GetSy 121 if (symb == ((*theElementTable)[i])->GetSymbol()) {
124 elm = (*theElementTable)[i]; 122 elm = (*theElementTable)[i];
125 break; 123 break;
126 } 124 }
127 } 125 }
128 if (nullptr == elm) { 126 if (nullptr == elm) {
129 for (G4int Z = 1; Z < maxNumElements; ++Z) 127 for (G4int Z = 1; Z < maxNumElements; ++Z) {
130 if (symb == elmSymbol[Z]) { 128 if (symb == elmSymbol[Z]) {
131 elm = BuildElement(Z); 129 elm = BuildElement(Z);
132 break; 130 break;
133 } 131 }
134 } 132 }
135 } 133 }
136 return elm; 134 return elm;
137 } 135 }
138 136
139 //....oooOO0OOooo........oooOO0OOooo........oo 137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
140 138
141 G4Element* G4NistElementBuilder::BuildElement( 139 G4Element* G4NistElementBuilder::BuildElement(G4int Z)
142 { 140 {
143 G4Element* theElement = nullptr; 141 G4Element* theElement = nullptr;
144 if (Z < 1 || Z >= maxNumElements) { 142 if (Z < 1 || Z >= maxNumElements) {
145 return theElement; 143 return theElement;
146 } 144 }
147 const G4ElementTable* theTable = G4Element:: <<
148 // Nist element was already built <<
149 if (0 <= elmIndex[Z]) { <<
150 return (*theTable)[elmIndex[Z]]; <<
151 } <<
152 145
153 G4AutoLock l(&nistElementMutex); 146 G4AutoLock l(&nistElementMutex);
154 if (0 <= elmIndex[Z]) { 147 if (0 <= elmIndex[Z]) {
>> 148 // Nist element was already built
>> 149 const G4ElementTable* theTable = G4Element::GetElementTable();
155 theElement = (*theTable)[elmIndex[Z]]; 150 theElement = (*theTable)[elmIndex[Z]];
156 } 151 }
157 else { 152 else {
158 G4double Aeff = atomicMass[Z]; 153 G4double Aeff = atomicMass[Z];
159 if (verbose > 1) { 154 if (verbose > 1) {
160 G4cout << "G4NistElementBuilder: Build E 155 G4cout << "G4NistElementBuilder: Build Element <" << elmSymbol[Z] << "> Z= " << Z
161 << " Aeff= " << Aeff; 156 << " Aeff= " << Aeff;
162 G4cout << " with natural isotope compos 157 G4cout << " with natural isotope composition" << G4endl;
163 } 158 }
164 159
165 // build Element with its Isotopes 160 // build Element with its Isotopes
166 // 161 //
167 G4int nc = nIsotopes[Z]; 162 G4int nc = nIsotopes[Z];
168 G4int n0 = nFirstIsotope[Z]; 163 G4int n0 = nFirstIsotope[Z];
169 G4int idx = idxIsotopes[Z]; 164 G4int idx = idxIsotopes[Z];
170 std::vector<G4Isotope*> iso; 165 std::vector<G4Isotope*> iso;
>> 166 G4Isotope* ist;
171 for (G4int i = 0; i < nc; ++i) { 167 for (G4int i = 0; i < nc; ++i) {
172 if (relAbundance[idx + i] > 0.0) { 168 if (relAbundance[idx + i] > 0.0) {
173 G4int n1 = n0 + i; <<
174 std::ostringstream os; 169 std::ostringstream os;
175 os << n1; << 170 os << elmSymbol[Z] << n0 + i;
176 const G4String& ss = elmSymbol[Z] + os << 171 ist = new G4Isotope(os.str(), Z, n0 + i, GetAtomicMass(Z, n0 + i) * g / (mole * amu_c2));
177 iso.push_back(new G4Isotope(ss, Z, n1, <<
178 CLHEP::g / (CLHEP::mole * CLHEP::a <<
179 /* 172 /*
180 G4cout << " Z= " << Z << " N= " << n 173 G4cout << " Z= " << Z << " N= " << n0 + i
181 << " miso(amu)= " << GetIsotopeMass 174 << " miso(amu)= " << GetIsotopeMass(Z, n0 + i)/amu_c2
182 << " matom(amu)= " << GetAtomicMass( 175 << " matom(amu)= " << GetAtomicMass(Z, n0 + i)/amu_c2 << G4endl;
183 */ 176 */
>> 177 iso.push_back(ist);
184 } 178 }
185 } 179 }
186 auto ni = (G4int)iso.size(); 180 auto ni = (G4int)iso.size();
>> 181 G4double w;
187 theElement = new G4Element(elmSymbol[Z], e 182 theElement = new G4Element(elmSymbol[Z], elmSymbol[Z], ni);
188 for (G4int j = 0; j < ni; ++j) { 183 for (G4int j = 0; j < ni; ++j) {
189 G4double w = relAbundance[idx + (iso[j]) << 184 w = relAbundance[idx + (iso[j])->GetN() - n0];
190 theElement->AddIsotope(iso[j], w); << 185 ist = iso[j];
>> 186 theElement->AddIsotope(ist, w);
191 } 187 }
192 theElement->SetNaturalAbundanceFlag(true); 188 theElement->SetNaturalAbundanceFlag(true);
193 elmIndex[Z] = (G4int)theElement->GetIndex( 189 elmIndex[Z] = (G4int)theElement->GetIndex();
194 } 190 }
195 l.unlock(); 191 l.unlock();
196 return theElement; 192 return theElement;
197 } 193 }
198 194
199 //....oooOO0OOooo........oooOO0OOooo........oo 195 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
200 196
201 void G4NistElementBuilder::PrintElement(G4int 197 void G4NistElementBuilder::PrintElement(G4int Z) const
202 { 198 {
203 G4int imin = Z; 199 G4int imin = Z;
204 G4int imax = Z + 1; 200 G4int imax = Z + 1;
205 if (Z == 0) { 201 if (Z == 0) {
206 imin = 1; 202 imin = 1;
207 imax = maxNumElements; 203 imax = maxNumElements;
208 } 204 }
209 if (imax > maxNumElements) { 205 if (imax > maxNumElements) {
210 imax = maxNumElements; 206 imax = maxNumElements;
211 } 207 }
212 208
213 for (G4int i = imin; i < imax; ++i) { 209 for (G4int i = imin; i < imax; ++i) {
214 G4int nc = nIsotopes[i]; 210 G4int nc = nIsotopes[i];
215 G4cout << "Nist Element: <" << elmSymbol[i 211 G4cout << "Nist Element: <" << elmSymbol[i] << "> Z= " << i << " Aeff(amu)= " << atomicMass[i]
216 << " " << nc << " isotopes:" << G4 212 << " " << nc << " isotopes:" << G4endl;
217 G4int j; 213 G4int j;
218 G4int idx = idxIsotopes[i]; 214 G4int idx = idxIsotopes[i];
219 G4int n0 = nFirstIsotope[i]; 215 G4int n0 = nFirstIsotope[i];
220 G4cout << " N: "; 216 G4cout << " N: ";
221 for (j = 0; j < nc; ++j) { 217 for (j = 0; j < nc; ++j) {
222 G4cout << n0 + j << " "; 218 G4cout << n0 + j << " ";
223 } 219 }
224 G4cout << G4endl; 220 G4cout << G4endl;
225 G4cout << " mass(amu): "; 221 G4cout << " mass(amu): ";
226 for (j = 0; j < nc; ++j) { 222 for (j = 0; j < nc; ++j) {
227 G4cout << GetAtomicMass(i, n0 + j) << " 223 G4cout << GetAtomicMass(i, n0 + j) << " ";
228 } 224 }
229 G4cout << G4endl; 225 G4cout << G4endl;
230 G4cout << " abundance: "; 226 G4cout << " abundance: ";
231 for (j = 0; j < nc; ++j) { 227 for (j = 0; j < nc; ++j) {
232 G4cout << relAbundance[idx + j] << " "; 228 G4cout << relAbundance[idx + j] << " ";
233 } 229 }
234 G4cout << G4endl; 230 G4cout << G4endl;
235 } 231 }
236 } 232 }
237 233
238 //....oooOO0OOooo........oooOO0OOooo........oo 234 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
239 235
240 void G4NistElementBuilder::AddElement(const G4 236 void G4NistElementBuilder::AddElement(const G4String& name, G4int Z, G4int nc, const G4int& N,
241 const G4double& A, const G4double& sA, const 237 const G4double& A, const G4double& sA, const G4double& W)
242 { 238 {
243 if (verbose > 1) { 239 if (verbose > 1) {
244 G4cout << "AddElement " << name << " Z= " 240 G4cout << "AddElement " << name << " Z= " << Z << " nc= " << nc << G4endl;
245 } 241 }
246 if (Z >= maxNumElements) { 242 if (Z >= maxNumElements) {
247 G4cout << "G4NistElementBuilder::AddElemen 243 G4cout << "G4NistElementBuilder::AddElement: Warning! Z= " << Z << " is too big" << G4endl;
248 return; 244 return;
249 } 245 }
250 246
251 if (index + nc >= maxAbundance) { 247 if (index + nc >= maxAbundance) {
252 G4cout << "G4NistElementBuilder::AddElemen 248 G4cout << "G4NistElementBuilder::AddElement: ERROR index= " << index << " + nc= " << nc
253 << " is above array size " << maxAb 249 << " is above array size " << maxAbundance << G4endl;
254 return; 250 return;
255 } 251 }
256 252
257 elmSymbol.push_back(name); << 253 elmSymbol[Z] = name;
258 atomicMass[Z] = 0.0; 254 atomicMass[Z] = 0.0;
259 idxIsotopes[Z] = index; 255 idxIsotopes[Z] = index;
260 nIsotopes[Z] = nc; 256 nIsotopes[Z] = nc;
261 257
262 nFirstIsotope[Z] = (&N)[0]; 258 nFirstIsotope[Z] = (&N)[0];
263 G4double ww = 0.0; 259 G4double ww = 0.0;
264 G4double www; 260 G4double www;
265 261
266 for (G4int i = 0; i < nc; ++i) { 262 for (G4int i = 0; i < nc; ++i) {
267 www = 0.01 * (&W)[i]; 263 www = 0.01 * (&W)[i];
268 // mass of the isotope in G4 units 264 // mass of the isotope in G4 units
269 massIsotopes[index] = (&A)[i] * amu_c2 - Z 265 massIsotopes[index] = (&A)[i] * amu_c2 - Z * electron_mass_c2 + bindingEnergy[Z];
270 sigMass[index] = (&sA)[i] * amu_c2; 266 sigMass[index] = (&sA)[i] * amu_c2;
271 relAbundance[index] = www; 267 relAbundance[index] = www;
272 268
273 // computation of mean atomic mass of the 269 // computation of mean atomic mass of the element in atomic units
274 atomicMass[Z] += www * (&A)[i]; 270 atomicMass[Z] += www * (&A)[i];
275 ww += www; 271 ww += www;
276 ++index; 272 ++index;
277 } 273 }
278 274
279 if (ww > 0.0) { << 275 if (ww != 1.0) {
280 G4int idx = idxIsotopes[Z]; 276 G4int idx = idxIsotopes[Z];
281 atomicMass[Z] /= ww; 277 atomicMass[Z] /= ww;
282 for (G4int j = 0; j < nc; ++j) { 278 for (G4int j = 0; j < nc; ++j) {
283 relAbundance[idx + j] /= ww; 279 relAbundance[idx + j] /= ww;
284 } 280 }
285 } 281 }
286 282
287 if (verbose > 1) { 283 if (verbose > 1) {
288 PrintElement(Z); 284 PrintElement(Z);
289 } 285 }
290 else if (1 == verbose) { 286 else if (1 == verbose) {
291 G4cout << "Nist Element " << elmSymbol[Z] 287 G4cout << "Nist Element " << elmSymbol[Z] << " Z= " << Z << " Aeff(amu)= " << atomicMass[Z]
292 << " " << nIsotopes[Z] << " isotop 288 << " " << nIsotopes[Z] << " isotopes:" << G4endl;
293 } 289 }
294 } 290 }
295 291
296 //....oooOO0OOooo........oooOO0OOooo........oo 292 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
297 293
298 void G4NistElementBuilder::Initialise() 294 void G4NistElementBuilder::Initialise()
299 { 295 {
300 // Parameterisation from D.Lunney,J.M.Pearso 296 // Parameterisation from D.Lunney,J.M.Pearson,C.Thibault,
301 // Rev.Mod.Phys. 75 (2003) 1021 297 // Rev.Mod.Phys. 75 (2003) 1021
302 bindingEnergy[0] = 0.0; 298 bindingEnergy[0] = 0.0;
303 for (G4int i = 1; i < maxNumElements; ++i) { 299 for (G4int i = 1; i < maxNumElements; ++i) {
304 auto Z = G4double(i); 300 auto Z = G4double(i);
305 bindingEnergy[i] = (14.4381 * std::pow(Z, 301 bindingEnergy[i] = (14.4381 * std::pow(Z, 2.39) + 1.55468e-6 * std::pow(Z, 5.35)) * eV;
306 } 302 }
307 303
308 // NIST data 304 // NIST data
309 index = 0; 305 index = 0;
310 306
311 // Z = 1 ----------------------------------- 307 // Z = 1 ---------------------------------------------------------------------
312 G4int HN[6] = {1, 2, 3, 4, 5, 6}; 308 G4int HN[6] = {1, 2, 3, 4, 5, 6};
313 309
314 G4double HA[6] = {1.00783, 2.0141, 3.01605, 310 G4double HA[6] = {1.00783, 2.0141, 3.01605, 4.02783, 5.03954, 6.04494};
315 311
316 // Garantee consistence with G4 masses 312 // Garantee consistence with G4 masses
317 HA[0] = (proton_mass_c2 + electron_mass_c2 - 313 HA[0] = (proton_mass_c2 + electron_mass_c2 - bindingEnergy[1]) / amu_c2;
318 HA[1] = (1.875613 * GeV + electron_mass_c2 - 314 HA[1] = (1.875613 * GeV + electron_mass_c2 - bindingEnergy[1]) / amu_c2;
319 HA[2] = (2.80925 * GeV + electron_mass_c2 - 315 HA[2] = (2.80925 * GeV + electron_mass_c2 - bindingEnergy[1]) / amu_c2;
320 316
321 G4double HS[6] = {4, 4, 11, 12, 102, 28}; 317 G4double HS[6] = {4, 4, 11, 12, 102, 28};
322 318
323 G4double HW[6] = {99.9885, 0.0115, 0, 0, 0, 319 G4double HW[6] = {99.9885, 0.0115, 0, 0, 0, 0};
324 320
325 AddElement("H", 1, 6, *HN, *HA, *HS, *HW); 321 AddElement("H", 1, 6, *HN, *HA, *HS, *HW);
326 322
327 // Z = 2 ----------------------------------- 323 // Z = 2 ---------------------------------------------------------------------
328 G4int HeN[8] = {3, 4, 5, 6, 7, 8, 9, 10}; 324 G4int HeN[8] = {3, 4, 5, 6, 7, 8, 9, 10};
329 325
330 G4double HeA[8] = {3.01603, 4.0026, 5.01222, 326 G4double HeA[8] = {3.01603, 4.0026, 5.01222, 6.01889, 7.02803, 8.03392, 9.04382, 10.0524};
331 327
332 // Garantee consistence with G4 masses 328 // Garantee consistence with G4 masses
333 HeA[0] = (2.80923 * GeV + 2.0 * electron_mas 329 HeA[0] = (2.80923 * GeV + 2.0 * electron_mass_c2 - bindingEnergy[2]) / amu_c2;
334 HeA[1] = (3.727417 * GeV + 2.0 * electron_ma 330 HeA[1] = (3.727417 * GeV + 2.0 * electron_mass_c2 - bindingEnergy[2]) / amu_c2;
335 331
336 G4double HeS[8] = {9, 10, 50, 11, 30, 8, 70, 332 G4double HeS[8] = {9, 10, 50, 11, 30, 8, 70, 80};
337 333
338 G4double HeW[8] = {0.000137, 99.9999, 0, 0, 334 G4double HeW[8] = {0.000137, 99.9999, 0, 0, 0, 0, 0, 0};
339 335
340 AddElement("He", 2, 8, *HeN, *HeA, *HeS, *He 336 AddElement("He", 2, 8, *HeN, *HeA, *HeS, *HeW);
341 337
342 // Z = 3 ----------------------------------- 338 // Z = 3 ---------------------------------------------------------------------
343 G4int LiN[9] = {4, 5, 6, 7, 8, 9, 10, 11, 12 339 G4int LiN[9] = {4, 5, 6, 7, 8, 9, 10, 11, 12};
344 G4double LiA[9] = {4.02718, 5.01254, 6.01512 340 G4double LiA[9] = {4.02718, 5.01254, 6.01512, 7.016, 8.02249, 9.02679, 10.0355, 11.0438, 12.0538};
345 341
346 G4double LiS[9] = {23, 50, 5, 5, 5, 21, 16, 342 G4double LiS[9] = {23, 50, 5, 5, 5, 21, 16, 29, 107};
347 343
348 G4double LiW[9] = {0, 0, 7.59, 92.41, 0, 0, 344 G4double LiW[9] = {0, 0, 7.59, 92.41, 0, 0, 0, 0, 0};
349 345
350 AddElement("Li", 3, 9, *LiN, *LiA, *LiS, *Li 346 AddElement("Li", 3, 9, *LiN, *LiA, *LiS, *LiW);
351 347
352 // Z = 4 ----------------------------------- 348 // Z = 4 ---------------------------------------------------------------------
353 G4int BeN[10] = {5, 6, 7, 8, 9, 10, 11, 12, 349 G4int BeN[10] = {5, 6, 7, 8, 9, 10, 11, 12, 13, 14};
354 350
355 G4double BeA[10] = { 351 G4double BeA[10] = {
356 5.04079, 6.01973, 7.01693, 8.00531, 9.0121 352 5.04079, 6.01973, 7.01693, 8.00531, 9.01218, 10.0135, 11.0217, 12.0269, 13.0361, 14.0428};
357 353
358 G4double BeS[10] = {429, 6, 5, 4, 4, 4, 7, 1 354 G4double BeS[10] = {429, 6, 5, 4, 4, 4, 7, 16, 54, 12};
359 355
360 G4double BeW[10] = {0, 0, 0, 0, 100, 0, 0, 0 356 G4double BeW[10] = {0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
361 357
362 AddElement("Be", 4, 10, *BeN, *BeA, *BeS, *B 358 AddElement("Be", 4, 10, *BeN, *BeA, *BeS, *BeW);
363 359
364 // Z = 5 ----------------------------------- 360 // Z = 5 ---------------------------------------------------------------------
365 G4int BN[13] = {7, 8, 9, 10, 11, 12, 13, 14, 361 G4int BN[13] = {7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19};
366 362
367 G4double BA[13] = {7.02992, 8.02461, 9.01333 363 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.063 364 15.0311, 16.0398, 17.0469, 18.0562, 19.0637};
369 365
370 G4double BS[13] = {80, 12, 11, 4, 5, 15, 12, 366 G4double BS[13] = {80, 12, 11, 4, 5, 15, 12, 23, 24, 60, 15, 86, 43};
371 367
372 G4double BW[13] = {0, 0, 0, 19.9, 80.1, 0, 0 368 G4double BW[13] = {0, 0, 0, 19.9, 80.1, 0, 0, 0, 0, 0, 0, 0, 0};
373 369
374 AddElement("B", 5, 13, *BN, *BA, *BS, *BW); 370 AddElement("B", 5, 13, *BN, *BA, *BS, *BW);
375 371
376 // Z = 6 ----------------------------------- 372 // Z = 6 ---------------------------------------------------------------------
377 G4int CN[15] = {8, 9, 10, 11, 12, 13, 14, 15 373 G4int CN[15] = {8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22};
378 374
379 G4double CA[15] = {8.03768, 9.03104, 10.0169 375 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.049 376 17.0226, 18.0268, 19.0353, 20.0403, 21.0493, 22.0565};
381 377
382 G4double CS[15] = {25, 23, 4, 10, 0, 10, 4, 378 G4double CS[15] = {25, 23, 4, 10, 0, 10, 4, 9, 4, 19, 30, 12, 22, 54, 97};
383 379
384 G4double CW[15] = {0, 0, 0, 0, 98.93, 1.07, 380 G4double CW[15] = {0, 0, 0, 0, 98.93, 1.07, 0, 0, 0, 0, 0, 0, 0, 0, 0};
385 381
386 AddElement("C", 6, 15, *CN, *CA, *CS, *CW); 382 AddElement("C", 6, 15, *CN, *CA, *CS, *CW);
387 383
388 // Z = 7 ----------------------------------- 384 // Z = 7 ---------------------------------------------------------------------
389 G4int NN[15] = {10, 11, 12, 13, 14, 15, 16, 385 G4int NN[15] = {10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};
390 386
391 G4double NA[15] = {10.0426, 11.0268, 12.0186 387 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 388 18.0141, 19.017, 20.0234, 21.0271, 22.0344, 23.0405, 24.0505};
393 389
394 G4double NS[15] = {43, 19, 11, 29, 9, 9, 28, 390 G4double NS[15] = {43, 19, 11, 29, 9, 9, 28, 16, 21, 18, 60, 10, 21, 76, 54};
395 391
396 G4double NW[15] = {0, 0, 0, 0, 99.632, 0.368 392 G4double NW[15] = {0, 0, 0, 0, 99.632, 0.368, 0, 0, 0, 0, 0, 0, 0, 0, 0};
397 393
398 AddElement("N", 7, 15, *NN, *NA, *NS, *NW); 394 AddElement("N", 7, 15, *NN, *NA, *NS, *NW);
399 395
400 // Z = 8 ----------------------------------- 396 // Z = 8 ---------------------------------------------------------------------
401 G4int ON[15] = {12, 13, 14, 15, 16, 17, 18, 397 G4int ON[15] = {12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26};
402 398
403 G4double OA[15] = {12.0344, 13.0248, 14.0086 399 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, 400 20.0041, 21.0087, 22.01, 23.0157, 24.0204, 25.0291, 26.0377};
405 401
406 G4double Os[15] = {20, 10, 8, 5, 15, 22, 9, 402 G4double Os[15] = {20, 10, 8, 5, 15, 22, 9, 3, 13, 13, 60, 11, 33, 40, 46};
407 403
408 G4double OW[15] = {0, 0, 0, 0, 99.757, 0.038 404 G4double OW[15] = {0, 0, 0, 0, 99.757, 0.038, 0.205, 0, 0, 0, 0, 0, 0, 0, 0};
409 405
410 AddElement("O", 8, 15, *ON, *OA, *Os, *OW); 406 AddElement("O", 8, 15, *ON, *OA, *Os, *OW);
411 407
412 // Z = 9 ----------------------------------- 408 // Z = 9 ---------------------------------------------------------------------
413 G4int FN[16] = {14, 15, 16, 17, 18, 19, 20, 409 G4int FN[16] = {14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29};
414 410
415 G4double FA[16] = {14.0361, 15.018, 16.0115, 411 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.026 412 23.0036, 24.0081, 25.0121, 26.0196, 27.0269, 28.0357, 29.0433};
417 413
418 G4double FS[16] = {43, 14, 9, 27, 6, 7, 9, 1 414 G4double FS[16] = {43, 14, 9, 27, 6, 7, 9, 19, 13, 90, 70, 80, 13, 45, 55, 62};
419 415
420 G4double FW[16] = {0, 0, 0, 0, 0, 100, 0, 0, 416 G4double FW[16] = {0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
421 417
422 AddElement("F", 9, 16, *FN, *FA, *FS, *FW); 418 AddElement("F", 9, 16, *FN, *FA, *FS, *FW);
423 419
424 // Z = 10 ---------------------------------- 420 // Z = 10 --------------------------------------------------------------------
425 G4int NeN[17] = {16, 17, 18, 19, 20, 21, 22, 421 G4int NeN[17] = {16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32};
426 422
427 G4double NeA[17] = {16.0258, 17.0177, 18.005 423 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.012 424 23.9936, 24.9978, 26.0005, 27.0076, 28.0121, 29.0193, 30.0239, 31.0331, 32.0399};
429 425
430 G4double NeS[17] = {22, 50, 16, 6, 20, 4, 23 426 G4double NeS[17] = {22, 50, 16, 6, 20, 4, 23, 26, 11, 50, 60, 10, 12, 32, 88, 97, 94};
431 427
432 G4double NeW[17] = {0, 0, 0, 0, 90.48, 0.27, 428 G4double NeW[17] = {0, 0, 0, 0, 90.48, 0.27, 9.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
433 429
434 AddElement("Ne", 10, 17, *NeN, *NeA, *NeS, * 430 AddElement("Ne", 10, 17, *NeN, *NeA, *NeS, *NeW);
435 431
436 // Z = 11 ---------------------------------- 432 // Z = 11 --------------------------------------------------------------------
437 G4int NaN[18] = {18, 19, 20, 21, 22, 23, 24, 433 G4int NaN[18] = {18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35};
438 434
439 G4double NaA[18] = {18.0272, 19.0139, 20.007 435 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 436 26.994, 27.9989, 29.0028, 30.0092, 31.0136, 32.0196, 33.0274, 34.0349, 35.0442};
441 437
442 G4double NaS[18] = {43, 13, 7, 8, 5, 23, 23, 438 G4double NaS[18] = {43, 13, 7, 8, 5, 23, 23, 13, 15, 40, 80, 10, 10, 18, 52, 160, 113, 166};
443 439
444 G4double NaW[18] = {0, 0, 0, 0, 0, 100, 0, 0 440 G4double NaW[18] = {0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
445 441
446 AddElement("Na", 11, 18, *NaN, *NaA, *NaS, * 442 AddElement("Na", 11, 18, *NaN, *NaA, *NaS, *NaW);
447 443
448 // Z = 12 ---------------------------------- 444 // Z = 12 --------------------------------------------------------------------
449 G4int MgN[18] = {20, 21, 22, 23, 24, 25, 26, 445 G4int MgN[18] = {20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37};
450 446
451 G4double MgA[18] = {20.0189, 21.0117, 21.999 447 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.999 448 27.9839, 28.9886, 29.9905, 30.9965, 31.9992, 33.0056, 34.0091, 35.0175, 36.0224, 37.0312};
453 449
454 G4double MgS[18] = {29, 18, 15, 13, 20, 20, 450 G4double MgS[18] = {29, 18, 15, 13, 20, 20, 21, 21, 22, 30, 70, 80, 10, 16, 28, 47, 97, 97};
455 451
456 G4double MgW[18] = {0, 0, 0, 0, 78.99, 10, 1 452 G4double MgW[18] = {0, 0, 0, 0, 78.99, 10, 11.01, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
457 453
458 AddElement("Mg", 12, 18, *MgN, *MgA, *MgS, * 454 AddElement("Mg", 12, 18, *MgN, *MgA, *MgS, *MgW);
459 455
460 // Z = 13 ---------------------------------- 456 // Z = 13 --------------------------------------------------------------------
461 G4int AlN[19] = {21, 22, 23, 24, 25, 26, 27, 457 G4int AlN[19] = {21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39};
462 458
463 G4double AlA[19] = {21.028, 22.0195, 23.0073 459 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 460 28.9804, 29.983, 30.9839, 31.9881, 32.9909, 33.9969, 34.9999, 36.0063, 37.0103, 38.0169,
465 39.0219}; 461 39.0219};
466 462
467 G4double AlS[19] = {32, 10, 27, 4, 7, 21, 14 463 G4double AlS[19] = {32, 10, 27, 4, 7, 21, 14, 15, 13, 15, 22, 90, 70, 10, 15, 29, 58, 60, 64};
468 464
469 G4double AlW[19] = {0, 0, 0, 0, 0, 0, 100, 0 465 G4double AlW[19] = {0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
470 466
471 AddElement("Al", 13, 19, *AlN, *AlA, *AlS, * 467 AddElement("Al", 13, 19, *AlN, *AlA, *AlS, *AlW);
472 468
473 // Z = 14 ---------------------------------- 469 // Z = 14 --------------------------------------------------------------------
474 G4int SiN[21] = { 470 G4int SiN[21] = {
475 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32 471 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42};
476 472
477 G4double SiA[21] = {22.0345, 23.0255, 24.011 473 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 474 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}; 475 41.0127, 42.0161};
480 476
481 G4double SiS[21] = { 477 G4double SiS[21] = {
482 22, 21, 21, 11, 3, 17, 20, 3, 5, 7, 23, 17 478 22, 21, 21, 11, 3, 17, 20, 3, 5, 7, 23, 17, 15, 40, 11, 13, 29, 43, 54, 64, 75};
483 479
484 G4double SiW[21] = { 480 G4double SiW[21] = {
485 0, 0, 0, 0, 0, 0, 92.2297, 4.6832, 3.0872, 481 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 482
487 AddElement("Si", 14, 21, *SiN, *SiA, *SiS, * 483 AddElement("Si", 14, 21, *SiN, *SiA, *SiS, *SiW);
488 484
489 // Z = 15 ---------------------------------- 485 // Z = 15 --------------------------------------------------------------------
490 G4int PN[23] = { 486 G4int PN[23] = {
491 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34 487 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46};
492 488
493 G4double PA[23] = {24.0343, 25.0203, 26.0118 489 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.978 490 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.023 491 42.0001, 43.0033, 44.0099, 45.0151, 46.0238};
496 492
497 G4double PS[23] = { 493 G4double PS[23] = {
498 54, 21, 21, 40, 4, 8, 4, 20, 20, 12, 5, 20 494 54, 21, 21, 40, 4, 8, 4, 20, 20, 12, 5, 20, 14, 40, 15, 16, 21, 50, 54, 54, 75, 86, 97};
499 495
500 G4double PW[23] = {0, 0, 0, 0, 0, 0, 0, 100, 496 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 497
502 AddElement("P", 15, 23, *PN, *PA, *PS, *PW); 498 AddElement("P", 15, 23, *PN, *PA, *PS, *PW);
503 499
504 // Z = 16 ---------------------------------- 500 // Z = 16 --------------------------------------------------------------------
505 G4int SN[24] = { 501 G4int SN[24] = {
506 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36 502 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 503
508 G4double SA[24] = {26.0279, 27.0188, 28.0044 504 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 505 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} 506 44.9948, 45.9996, 47.0076, 48.013, 49.022};
511 507
512 G4double SS[24] = { 508 G4double SS[24] = {
513 32, 22, 17, 50, 3, 16, 12, 12, 11, 10, 25, 509 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 510
515 G4double SW[24] = { 511 G4double SW[24] = {
516 0, 0, 0, 0, 0, 0, 94.93, 0.76, 4.29, 0, 0. 512 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 513
518 AddElement("S", 16, 24, *SN, *SA, *SS, *SW); 514 AddElement("S", 16, 24, *SN, *SA, *SS, *SW);
519 515
520 // Z = 17 ---------------------------------- 516 // Z = 17 --------------------------------------------------------------------
521 G4int ClN[24] = { 517 G4int ClN[24] = {
522 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38 518 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 519
524 G4double ClA[24] = {28.0285, 29.0141, 30.004 520 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, 521 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.013 522 46.9879, 47.9948, 48.9999, 50.0077, 51.0135};
527 523
528 G4double ClS[24] = { 524 G4double ClS[24] = {
529 54, 21, 21, 50, 7, 6, 13, 4, 8, 5, 12, 19, 525 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 526
531 G4double ClW[24] = { 527 G4double ClW[24] = {
532 0, 0, 0, 0, 0, 0, 0, 75.78, 0, 24.22, 0, 0 528 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 529
534 AddElement("Cl", 17, 24, *ClN, *ClA, *ClS, * 530 AddElement("Cl", 17, 24, *ClN, *ClA, *ClS, *ClW);
535 531
536 // Z = 18 ---------------------------------- 532 // Z = 18 --------------------------------------------------------------------
537 G4int ArN[24] = { 533 G4int ArN[24] = {
538 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40 534 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 535
540 G4double ArA[24] = {30.0216, 31.0121, 31.997 536 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.963 537 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.998 538 47.9751, 48.9822, 49.9859, 50.9932, 51.9982, 53.0062};
543 539
544 G4double ArS[24] = { 540 G4double ArS[24] = {
545 32, 22, 50, 30, 3, 8, 27, 3, 5, 5, 3, 7, 4 541 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 542
547 G4double ArW[24] = { 543 G4double ArW[24] = {
548 0, 0, 0, 0, 0, 0, 0.3365, 0, 0.0632, 0, 99 544 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 545
550 AddElement("Ar", 18, 24, *ArN, *ArA, *ArS, * 546 AddElement("Ar", 18, 24, *ArN, *ArA, *ArS, *ArW);
551 547
552 // Z = 19 ---------------------------------- 548 // Z = 19 --------------------------------------------------------------------
553 G4int KN[24] = { 549 G4int KN[24] = {
554 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42 550 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 551
556 G4double KA[24] = {32.0219, 33.0073, 33.9984 552 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.960 553 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 554 50.9764, 51.9826, 52.9871, 53.994, 54.9994};
559 555
560 G4double KS[24] = { 556 G4double KS[24] = {
561 54, 21, 32, 21, 8, 29, 8, 3, 29, 28, 3, 10 557 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 558
563 G4double KW[24] = { 559 G4double KW[24] = {
564 0, 0, 0, 0, 0, 0, 0, 93.2581, 0.0117, 6.73 560 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 561
566 AddElement("K", 19, 24, *KN, *KA, *KS, *KW); 562 AddElement("K", 19, 24, *KN, *KA, *KS, *KW);
567 563
568 // Z = 20 ---------------------------------- 564 // Z = 20 --------------------------------------------------------------------
569 G4int CaN[24] = { 565 G4int CaN[24] = {
570 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44 566 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 567
572 G4double CaA[24] = {34.0141, 35.0048, 35.993 568 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.953 569 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.985 570 51.9651, 52.9701, 53.9747, 54.9806, 55.9858, 56.9924};
575 571
576 G4double CaS[24] = { 572 G4double CaS[24] = {
577 32, 70, 40, 24, 5, 19, 3, 4, 4, 5, 9, 10, 573 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 574
579 G4double CaW[24] = {0, 0, 0, 0, 0, 0, 96.941 575 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}; 576 0, 0, 0, 0, 0, 0};
581 577
582 AddElement("Ca", 20, 24, *CaN, *CaA, *CaS, * 578 AddElement("Ca", 20, 24, *CaN, *CaA, *CaS, *CaW);
583 579
584 // Z = 21 ---------------------------------- 580 // Z = 21 --------------------------------------------------------------------
585 G4int ScN[24] = { 581 G4int ScN[24] = {
586 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46 582 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 583
588 G4double ScA[24] = {36.0149, 37.0031, 37.994 584 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.952 585 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} 586 54.9674, 55.9727, 56.977, 57.9831, 58.988};
591 587
592 G4double ScS[24] = { 588 G4double ScS[24] = {
593 54, 32, 32, 26, 4, 3, 4, 20, 19, 12, 12, 2 589 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 590
595 G4double ScW[24] = {0, 0, 0, 0, 0, 0, 0, 0, 591 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 592
597 AddElement("Sc", 21, 24, *ScN, *ScA, *ScS, * 593 AddElement("Sc", 21, 24, *ScN, *ScA, *ScS, *ScW);
598 594
599 // Z = 22 ---------------------------------- 595 // Z = 22 --------------------------------------------------------------------
600 G4int TiN[24] = { 596 G4int TiN[24] = {
601 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48 597 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 598
603 G4double TiA[24] = {38.0098, 39.0013, 39.990 599 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.944 600 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, 601 55.958, 56.9629, 57.9661, 58.972, 59.9756, 60.982};
606 602
607 G4double TiS[24] = { 603 G4double TiS[24] = {
608 27, 11, 17, 40, 6, 7, 8, 13, 12, 10, 10, 1 604 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 605
610 G4double TiW[24] = { 606 G4double TiW[24] = {
611 0, 0, 0, 0, 0, 0, 0, 0, 8.25, 7.44, 73.72, 607 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 608
613 AddElement("Ti", 22, 24, *TiN, *TiA, *TiS, * 609 AddElement("Ti", 22, 24, *TiN, *TiA, *TiS, *TiW);
614 610
615 // Z = 23 ---------------------------------- 611 // Z = 23 --------------------------------------------------------------------
616 G4int VN[24] = { 612 G4int VN[24] = {
617 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50 613 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 614
619 G4double VA[24] = {40.0111, 40.9997, 41.9912 615 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 616 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.973 617 57.9567, 58.9593, 59.9645, 60.9674, 61.9731, 62.9768};
622 618
623 G4double VS[24] = { 619 G4double VS[24] = {
624 54, 27, 21, 25, 90, 18, 16, 12, 28, 14, 14 620 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 621
626 G4double VW[24] = {0, 0, 0, 0, 0, 0, 0, 0, 0 622 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 623
628 AddElement("V", 23, 24, *VN, *VA, *VS, *VW); 624 AddElement("V", 23, 24, *VN, *VA, *VS, *VW);
629 625
630 // Z = 24 ---------------------------------- 626 // Z = 24 --------------------------------------------------------------------
631 G4int CrN[24] = { 627 G4int CrN[24] = {
632 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52 628 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 629
634 G4double CrA[24] = {42.0064, 42.9977, 43.985 630 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.940 631 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.970 632 60.9541, 61.9558, 62.9619, 63.9642, 64.9704};
637 633
638 G4double CrS[24] = { 634 G4double CrS[24] = {
639 32, 90, 14, 11, 22, 15, 8, 28, 14, 14, 15, 635 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 636
641 G4double CrW[24] = { 637 G4double CrW[24] = {
642 0, 0, 0, 0, 0, 0, 0, 0, 4.345, 0, 83.789, 638 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 639
644 AddElement("Cr", 24, 24, *CrN, *CrA, *CrS, * 640 AddElement("Cr", 24, 24, *CrN, *CrA, *CrS, *CrW);
645 641
646 // Z = 25 ---------------------------------- 642 // Z = 25 --------------------------------------------------------------------
647 G4int MnN[24] = { 643 G4int MnN[24] = {
648 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54 644 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 645
650 G4double MnA[24] = {44.0069, 44.9945, 45.986 646 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 647 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.963 648 62.9498, 63.9537, 64.9561, 65.9608, 66.9638};
653 649
654 G4double MnS[24] = { 650 G4double MnS[24] = {
655 54, 32, 12, 17, 80, 26, 15, 14, 25, 15, 18 651 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 652
657 G4double MnW[24] = {0, 0, 0, 0, 0, 0, 0, 0, 653 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 654
659 AddElement("Mn", 25, 24, *MnN, *MnA, *MnS, * 655 AddElement("Mn", 25, 24, *MnN, *MnA, *MnS, *MnW);
660 656
661 // Z = 26 ---------------------------------- 657 // Z = 26 --------------------------------------------------------------------
662 G4int FeN[25] = {45, 46, 47, 48, 49, 50, 51, 658 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}; 659 65, 66, 67, 68, 69};
664 660
665 G4double FeA[25] = {45.0146, 46.0008, 46.992 661 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.935 662 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, 663 62.9401, 63.9409, 64.9449, 65.946, 66.95, 67.9525, 68.9577};
668 664
669 G4double FeS[25] = {43, 38, 28, 11, 17, 60, 665 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}; 666 30, 35, 50, 75, 86};
671 667
672 G4double FeW[25] = { 668 G4double FeW[25] = {
673 0, 0, 0, 0, 0, 0, 0, 0, 0, 5.845, 0, 91.75 669 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 670
675 AddElement("Fe", 26, 25, *FeN, *FeA, *FeS, * 671 AddElement("Fe", 26, 25, *FeN, *FeA, *FeS, *FeW);
676 672
677 // Z = 27 ---------------------------------- 673 // Z = 27 --------------------------------------------------------------------
678 G4int CoN[25] = {48, 49, 50, 51, 52, 53, 54, 674 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}; 675 68, 69, 70, 71, 72};
680 676
681 G4double CoA[25] = {48.0018, 48.9897, 49.981 677 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.933 678 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.949 679 65.9398, 66.9406, 67.9444, 68.9452, 69.9498, 70.9517, 71.9564};
684 680
685 G4double CoS[25] = {43, 28, 18, 16, 70, 19, 681 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}; 682 30, 35, 40, 75, 86, 86};
687 683
688 G4double CoW[25] = {0, 0, 0, 0, 0, 0, 0, 0, 684 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 685
690 AddElement("Co", 27, 25, *CoN, *CoA, *CoS, * 686 AddElement("Co", 27, 25, *CoN, *CoA, *CoS, *CoW);
691 687
692 // Z = 28 ---------------------------------- 688 // Z = 28 --------------------------------------------------------------------
693 G4int NiN[29] = {50, 51, 52, 53, 54, 55, 56, 689 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}; 690 70, 71, 72, 73, 74, 75, 76, 77, 78};
695 691
696 G4double NiA[29] = {49.9959, 50.9877, 51.975 692 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.928 693 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, 694 67.9318, 68.9352, 69.9361, 70.94, 71.9413, 72.9461, 73.9479, 74.953, 75.9553, 76.9608, 77.9638};
699 695
700 G4double NiS[29] = {28, 28, 90, 17, 50, 12, 696 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}; 697 35, 40, 50, 64, 75, 86, 97, 107, 118};
702 698
703 G4double NiW[29] = {0, 0, 0, 0, 0, 0, 0, 0, 699 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}; 700 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
705 701
706 AddElement("Ni", 28, 29, *NiN, *NiA, *NiS, * 702 AddElement("Ni", 28, 29, *NiN, *NiA, *NiS, *NiW);
707 703
708 // Z = 29 ---------------------------------- 704 // Z = 29 --------------------------------------------------------------------
709 G4int CuN[29] = {52, 53, 54, 55, 56, 57, 58, 705 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}; 706 72, 73, 74, 75, 76, 77, 78, 79, 80};
711 707
712 G4double CuA[29] = {51.9972, 52.9856, 53.976 708 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.929 709 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.940 710 69.9324, 70.9326, 71.9355, 72.9365, 73.9402, 74.9417, 75.946, 76.948, 77.9528, 78.9553,
715 79.9619}; 711 79.9619};
716 712
717 G4double CuS[29] = {28, 28, 23, 32, 15, 17, 713 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}; 714 21, 32, 43, 54, 64, 75, 86, 97, 97};
719 715
720 G4double CuW[29] = { 716 G4double CuW[29] = {
721 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 69.17, 0, 717 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 718
723 AddElement("Cu", 29, 29, *CuN, *CuA, *CuS, * 719 AddElement("Cu", 29, 29, *CuN, *CuA, *CuS, *CuW);
724 720
725 // Z = 30 ---------------------------------- 721 // Z = 30 --------------------------------------------------------------------
726 G4int ZnN[29] = {54, 55, 56, 57, 58, 59, 60, 722 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}; 723 74, 75, 76, 77, 78, 79, 80, 81, 82};
728 724
729 G4double ZnA[29] = {53.993, 54.984, 55.9724, 725 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 726 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.937 727 72.9298, 73.9295, 74.9329, 75.9334, 76.9371, 77.9386, 78.9427, 79.9444, 80.9505, 81.9548};
732 728
733 G4double ZnS[29] = {43, 27, 28, 15, 50, 40, 729 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}; 730 50, 80, 13, 14, 17, 29, 18, 43, 43};
735 731
736 G4double ZnW[29] = {0, 0, 0, 0, 0, 0, 0, 0, 732 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}; 733 0, 0, 0, 0, 0, 0, 0, 0};
738 734
739 AddElement("Zn", 30, 29, *ZnN, *ZnA, *ZnS, * 735 AddElement("Zn", 30, 29, *ZnN, *ZnA, *ZnS, *ZnW);
740 736
741 // Z = 31 ---------------------------------- 737 // Z = 31 --------------------------------------------------------------------
742 G4int GaN[29] = {56, 57, 58, 59, 60, 61, 62, 738 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}; 739 76, 77, 78, 79, 80, 81, 82, 83, 84};
744 740
745 G4double GaA[29] = {55.9949, 56.9829, 57.974 741 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 742 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.932 743 74.9265, 75.9289, 76.9293, 77.9317, 78.9329, 79.9366, 80.9377, 81.9432, 82.9469, 83.9523};
748 744
749 G4double GaS[29] = {28, 28, 23, 18, 12, 21, 745 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}; 746 60, 90, 13, 13, 21, 32, 54, 64};
751 747
752 G4double GaW[29] = {0, 0, 0, 0, 0, 0, 0, 0, 748 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}; 749 0, 0, 0, 0, 0, 0};
754 750
755 AddElement("Ga", 31, 29, *GaN, *GaA, *GaS, * 751 AddElement("Ga", 31, 29, *GaN, *GaA, *GaS, *GaW);
756 752
757 // Z = 32 ---------------------------------- 753 // Z = 32 --------------------------------------------------------------------
758 G4int GeN[29] = {58, 59, 60, 61, 62, 63, 64, 754 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}; 755 78, 79, 80, 81, 82, 83, 84, 85, 86};
760 756
761 G4double GeA[29] = {57.991, 58.9817, 59.9702 757 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 758 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.928 759 76.9235, 77.9229, 78.9254, 79.9254, 80.9288, 81.9296, 82.9345, 83.9373, 84.9427, 85.9463};
764 760
765 G4double GeS[29] = {34, 30, 25, 32, 15, 21, 761 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}; 762 4, 10, 25, 13, 26, 32, 43, 54, 64};
767 763
768 G4double GeW[29] = {0, 0, 0, 0, 0, 0, 0, 0, 764 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}; 765 0, 0, 0, 0, 0, 0, 0, 0, 0};
770 766
771 AddElement("Ge", 32, 29, *GeN, *GeA, *GeS, * 767 AddElement("Ge", 32, 29, *GeN, *GeA, *GeS, *GeW);
772 768
773 // Z = 33 ---------------------------------- 769 // Z = 33 --------------------------------------------------------------------
774 G4int AsN[30] = {60, 61, 62, 63, 64, 65, 66, 770 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}; 771 80, 81, 82, 83, 84, 85, 86, 87, 88, 89};
776 772
777 G4double AsA[30] = {59.9931, 60.9806, 61.973 773 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.926 774 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.924 775 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}; 776 87.9446, 88.9492};
781 777
782 G4double AsS[30] = {64, 64, 32, 54, 38, 42, 778 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}; 779 23, 6, 21, 24, 32, 32, 43, 54, 64, 64};
784 780
785 G4double AsW[30] = { 781 G4double AsW[30] = {
786 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 782 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 783
788 AddElement("As", 33, 30, *AsN, *AsA, *AsS, * 784 AddElement("As", 33, 30, *AsN, *AsA, *AsS, *AsW);
789 785
790 // Z = 34 ---------------------------------- 786 // Z = 34 --------------------------------------------------------------------
791 G4int SeN[28] = {65, 66, 67, 68, 69, 70, 71, 787 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}; 788 85, 86, 87, 88, 89, 90, 91, 92};
793 789
794 G4double SeA[28] = {64.9647, 65.9552, 66.950 790 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.919 791 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.928 792 82.9191, 83.9185, 84.9222, 85.9243, 86.9285, 87.9314, 88.936, 89.9394, 90.9454, 91.9493};
797 793
798 G4double SeS[28] = {64, 32, 21, 32, 40, 22, 794 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}; 795 30, 17, 40, 50, 32, 43, 54, 64};
800 796
801 G4double SeW[28] = {0, 0, 0, 0, 0, 0, 0, 0, 797 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}; 798 0, 0, 0, 0, 0, 0, 0, 0, 0};
803 799
804 AddElement("Se", 34, 28, *SeN, *SeA, *SeS, * 800 AddElement("Se", 34, 28, *SeN, *SeA, *SeS, *SeW);
805 801
806 // Z = 35 ---------------------------------- 802 // Z = 35 --------------------------------------------------------------------
807 G4int BrN[28] = {67, 68, 69, 70, 71, 72, 73, 803 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}; 804 87, 88, 89, 90, 91, 92, 93, 94};
809 805
810 G4double BrA[28] = {66.9648, 67.9583, 68.950 806 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.918 807 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.926 808 84.9156, 85.9188, 86.9207, 87.9241, 88.9264, 89.9306, 90.934, 91.9393, 92.9431, 93.9487};
813 809
814 G4double BrS[28] = {54, 58, 34, 39, 32, 28, 810 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}; 811 40, 60, 80, 80, 50, 32, 43};
816 812
817 G4double BrW[28] = { 813 G4double BrW[28] = {
818 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 50.69, 814 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 815
820 AddElement("Br", 35, 28, *BrN, *BrA, *BrS, * 816 AddElement("Br", 35, 28, *BrN, *BrA, *BrS, *BrW);
821 817
822 // Z = 36 ---------------------------------- 818 // Z = 36 --------------------------------------------------------------------
823 G4int KrN[29] = {69, 70, 71, 72, 73, 74, 75, 819 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}; 820 89, 90, 91, 92, 93, 94, 95, 96, 97};
825 821
826 G4double KrA[29] = {68.9653, 69.956, 70.9505 822 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.913 823 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.926 824 87.9144, 88.9176, 89.9195, 90.9234, 91.9262, 92.9313, 93.9344, 94.9398, 95.9431, 96.9486};
829 825
830 G4double KrS[29] = {54, 43, 32, 29, 15, 60, 826 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}; 827 20, 60, 13, 11, 32, 43, 54, 54};
832 828
833 G4double KrW[29] = {0, 0, 0, 0, 0, 0, 0, 0, 829 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}; 830 0, 0, 0, 0, 0, 0, 0, 0, 0};
835 831
836 AddElement("Kr", 36, 29, *KrN, *KrA, *KrS, * 832 AddElement("Kr", 36, 29, *KrN, *KrA, *KrS, *KrW);
837 833
838 // Z = 37 ---------------------------------- 834 // Z = 37 --------------------------------------------------------------------
839 G4int RbN[32] = {71, 72, 73, 74, 75, 76, 77, 835 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, 1 836 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102};
841 837
842 G4double RbA[32] = {70.9653, 71.9591, 72.950 838 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, 839 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 840 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}; 841 99.9499, 100.953, 101.959};
846 842
847 G4double RbS[32] = {54, 54, 52, 77, 8, 8, 8, 843 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}; 844 21, 27, 30, 40, 16, 32, 18, 54};
849 845
850 G4double RbW[32] = {0, 0, 0, 0, 0, 0, 0, 0, 846 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}; 847 0, 0, 0, 0, 0, 0, 0, 0, 0};
852 848
853 AddElement("Rb", 37, 32, *RbN, *RbA, *RbS, * 849 AddElement("Rb", 37, 32, *RbN, *RbA, *RbS, *RbW);
854 850
855 // Z = 38 ---------------------------------- 851 // Z = 38 --------------------------------------------------------------------
856 G4int SrN[32] = {73, 74, 75, 76, 77, 78, 79, 852 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, 853 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104};
858 854
859 G4double SrA[32] = {72.966, 73.9563, 74.9499 855 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.912 856 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, 857 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}; 858 101.943, 102.949, 103.952};
863 859
864 G4double SrS[32] = {64, 54, 32, 32, 16, 8, 9 860 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}; 861 26, 20, 27, 15, 14, 13, 12, 54, 75};
866 862
867 G4double SrW[32] = {0, 0, 0, 0, 0, 0, 0, 0, 863 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}; 864 0, 0, 0, 0, 0, 0, 0, 0, 0};
869 865
870 AddElement("Sr", 38, 32, *SrN, *SrA, *SrS, * 866 AddElement("Sr", 38, 32, *SrN, *SrA, *SrS, *SrW);
871 867
872 // Z = 39 ---------------------------------- 868 // Z = 39 --------------------------------------------------------------------
873 G4int YN[30] = {77, 78, 79, 80, 81, 82, 83, 869 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, 870 97, 98, 99, 100, 101, 102, 103, 104, 105, 106};
875 871
876 G4double YA[30] = {76.9496, 77.9435, 78.9373 872 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.905 873 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.924 874 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}; 875 104.945, 105.95};
880 876
881 G4double YS[30] = {32, 43, 48, 43, 70, 11, 5 877 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}; 878 13, 26, 26, 80, 10, 90, 32, 43, 54, 75};
883 879
884 G4double YW[30] = { 880 G4double YW[30] = {
885 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0 881 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 882
887 AddElement("Y", 39, 30, *YN, *YA, *YS, *YW); 883 AddElement("Y", 39, 30, *YN, *YA, *YS, *YW);
888 884
889 // Z = 40 ---------------------------------- 885 // Z = 40 --------------------------------------------------------------------
890 G4int ZrN[30] = {79, 80, 81, 82, 83, 84, 85, 886 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 887 99, 100, 101, 102, 103, 104, 105, 106, 107, 108};
892 888
893 G4double ZrA[30] = {78.9492, 79.9406, 80.936 889 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.905 890 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.92 891 97.9127, 98.9165, 99.9178, 100.921, 101.923, 102.927, 103.929, 104.933, 105.936, 106.941,
896 107.944}; 892 107.944};
897 893
898 G4double ZrS[30] = {43, 32, 32, 55, 10, 21, 894 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}; 895 21, 40, 30, 50, 12, 43, 43, 54, 64, 75};
900 896
901 G4double ZrW[30] = {0, 0, 0, 0, 0, 0, 0, 0, 897 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}; 898 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
903 899
904 AddElement("Zr", 40, 30, *ZrN, *ZrA, *ZrS, * 900 AddElement("Zr", 40, 30, *ZrN, *ZrA, *ZrS, *ZrW);
905 901
906 // Z = 41 ---------------------------------- 902 // Z = 41 --------------------------------------------------------------------
907 G4int NbN[30] = {81, 82, 83, 84, 85, 86, 87, 903 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, 10 904 101, 102, 103, 104, 105, 106, 107, 108, 109, 110};
909 905
910 G4double NbA[30] = {80.949, 81.9431, 82.9367 906 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 907 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.92 908 99.9142, 100.915, 101.918, 102.919, 103.922, 104.924, 105.928, 106.93, 107.935, 108.938,
913 109.943}; 909 109.943};
914 910
915 G4double NbS[30] = {43, 32, 34, 32, 24, 90, 911 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}; 912 20, 40, 70, 12, 11, 32, 43, 54, 54, 64};
917 913
918 G4double NbW[30] = { 914 G4double NbW[30] = {
919 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0 915 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 916
921 AddElement("Nb", 41, 30, *NbN, *NbA, *NbS, * 917 AddElement("Nb", 41, 30, *NbN, *NbA, *NbS, *NbW);
922 918
923 // Z = 42 ---------------------------------- 919 // Z = 42 --------------------------------------------------------------------
924 G4int MoN[31] = {83, 84, 85, 86, 87, 88, 89, 920 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, 11 921 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113};
926 922
927 G4double MoA[31] = {82.9487, 83.9401, 84.936 923 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.905 924 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 925 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}; 926 112.942};
931 927
932 G4double MoS[31] = {54, 43, 43, 47, 24, 22, 928 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}; 929 70, 80, 23, 17, 21, 32, 43, 54, 64, 64};
934 930
935 G4double MoW[31] = {0, 0, 0, 0, 0, 0, 0, 0, 931 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}; 932 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
937 933
938 AddElement("Mo", 42, 31, *MoN, *MoA, *MoS, * 934 AddElement("Mo", 42, 31, *MoN, *MoA, *MoS, *MoW);
939 935
940 // Z = 43 ---------------------------------- 936 // Z = 43 --------------------------------------------------------------------
941 G4int TcN[31] = {85, 86, 87, 88, 89, 90, 91, 937 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, 11 938 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115};
943 939
944 G4double TcA[31] = {84.9489, 85.9429, 86.936 940 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.906 941 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.91 942 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}; 943 112.931, 113.936, 114.938};
948 944
949 G4double TcS[31] = {54, 32, 32, 32, 23, 26, 945 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}; 946 15, 16, 14, 23, 43, 43, 54, 64, 64, 75};
951 947
952 G4double TcW[31] = { 948 G4double TcW[31] = {
953 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100 949 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 950
955 AddElement("Tc", 43, 31, *TcN, *TcA, *TcS, * 951 AddElement("Tc", 43, 31, *TcN, *TcA, *TcS, *TcW);
956 952
957 // Z = 44 ---------------------------------- 953 // Z = 44 --------------------------------------------------------------------
958 G4int RuN[32] = {87, 88, 89, 90, 91, 92, 93, 954 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, 11 955 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118};
960 956
961 G4double RuA[32] = {86.9492, 87.9404, 88.936 957 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.905 958 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, 959 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}; 960 115.93, 116.935, 117.937};
965 961
966 G4double RuS[32] = {64, 54, 54, 43, 54, 32, 962 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 963 13, 70, 25, 32, 58, 54, 39, 64, 75, 86, 97};
968 964
969 G4double RuW[32] = {0, 0, 0, 0, 0, 0, 0, 0, 965 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}; 966 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
971 967
972 AddElement("Ru", 44, 32, *RuN, *RuA, *RuS, * 968 AddElement("Ru", 44, 32, *RuN, *RuA, *RuS, *RuW);
973 969
974 // Z = 45 ---------------------------------- 970 // Z = 45 --------------------------------------------------------------------
975 G4int RhN[33] = {89, 90, 91, 92, 93, 94, 95, 971 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, 11 972 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121};
977 973
978 G4double RhA[33] = {88.9494, 89.9429, 90.936 974 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.90 975 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.91 976 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.93 977 116.925, 117.929, 118.931, 119.936, 120.938};
982 978
983 G4double RhS[33] = {54, 54, 43, 43, 43, 48, 979 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 980 24, 22, 54, 43, 32, 54, 54, 64, 75, 86, 86, 97};
985 G4double RhW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 981 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}; 982 0, 0, 0, 0, 0, 0, 0, 0};
987 983
988 AddElement("Rh", 45, 33, *RhN, *RhA, *RhS, * 984 AddElement("Rh", 45, 33, *RhN, *RhA, *RhS, *RhW);
989 985
990 // Z = 46 ---------------------------------- 986 // Z = 46 --------------------------------------------------------------------
991 G4int PdN[33] = {91, 92, 93, 94, 95, 96, 97, 987 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, 11 988 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123};
993 989
994 G4double PdA[33] = {90.9495, 91.9404, 92.935 990 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.90 991 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 992 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}; 993 119.924, 120.928, 121.93, 122.934};
998 994
999 G4double PdS[33] = {64, 54, 43, 43, 43, 16, 995 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, 5 996 19, 40, 26, 70, 60, 32, 23, 32, 43, 54, 54, 64};
1001 997
1002 G4double PdW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 998 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 999 11.72, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1004 1000
1005 AddElement("Pd", 46, 33, *PdN, *PdA, *PdS, 1001 AddElement("Pd", 46, 33, *PdN, *PdA, *PdS, *PdW);
1006 1002
1007 // Z = 47 --------------------------------- 1003 // Z = 47 --------------------------------------------------------------------
1008 G4int AgN[34] = {94, 95, 96, 97, 98, 99, 10 1004 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, 1 1005 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127};
1010 1006
1011 G4double AgA[34] = {93.9428, 94.9355, 95.93 1007 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.9 1008 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.9 1009 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.9 1010 121.923, 122.925, 123.929, 124.931, 125.934, 126.937};
1015 1011
1016 G4double AgS[34] = {54, 43, 43, 43, 16, 16, 1012 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, 4 1013 40, 50, 50, 70, 10, 80, 16, 22, 32, 43, 43, 43, 54};
1018 1014
1019 G4double AgW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1015 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}; 1016 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1021 1017
1022 AddElement("Ag", 47, 34, *AgN, *AgA, *AgS, 1018 AddElement("Ag", 47, 34, *AgN, *AgA, *AgS, *AgW);
1023 1019
1024 // Z = 48 --------------------------------- 1020 // Z = 48 --------------------------------------------------------------------
1025 G4int CdN[35] = {96, 97, 98, 99, 100, 101, 1021 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, 1 1022 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130};
1027 1023
1028 G4double CdA[35] = {95.9398, 96.9349, 97.92 1024 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.9 1025 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.9 1026 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.9 1027 125.922, 126.926, 127.928, 128.932, 129.934};
1032 1028
1033 G4double CdS[35] = {54, 43, 22, 22, 10, 16, 1029 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, 3 1030 22, 90, 20, 90, 22, 40, 70, 70, 60, 80, 32, 43, 43};
1035 1031
1036 G4double CdW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1032 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 1033 28.73, 0, 7.49, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1038 1034
1039 AddElement("Cd", 48, 35, *CdN, *CdA, *CdS, 1035 AddElement("Cd", 48, 35, *CdN, *CdA, *CdS, *CdW);
1040 1036
1041 // Z = 49 --------------------------------- 1037 // Z = 49 --------------------------------------------------------------------
1042 G4int InN[37] = {98, 99, 100, 101, 102, 103 1038 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, 1 1039 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132,
1044 133, 134}; 1040 133, 134};
1045 1041
1046 G4double InA[37] = {97.9422, 98.9346, 99.93 1042 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 1043 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.9 1044 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.93 1045 127.92, 128.922, 129.925, 130.927, 131.933, 132.938, 133.945};
1050 1046
1051 G4double InS[37] = {54, 54, 41, 32, 41, 27, 1047 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, 5 1048 40, 29, 50, 26, 50, 30, 40, 40, 50, 14, 50, 80, 70, 43, 54};
1053 1049
1054 G4double InW[37] = {0, 0, 0, 0, 0, 0, 0, 0, 1050 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}; 1051 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1056 1052
1057 AddElement("In", 49, 37, *InN, *InA, *InS, 1053 AddElement("In", 49, 37, *InN, *InA, *InS, *InW);
1058 1054
1059 // Z = 50 --------------------------------- 1055 // Z = 50 --------------------------------------------------------------------
1060 G4int SnN[38] = {100, 101, 102, 103, 104, 1 1056 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, 1 1057 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134,
1062 135, 136, 137}; 1058 135, 136, 137};
1063 1059
1064 G4double SnA[38] = {99.939, 100.936, 101.93 1060 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.9 1061 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.9 1062 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.9 1063 128.913, 129.914, 130.917, 131.918, 132.924, 133.928, 134.935, 135.939, 136.946};
1068 1064
1069 G4double SnS[38] = {46, 54, 43, 32, 16, 10, 1065 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, 2 1066 29, 29, 15, 16, 11, 27, 29, 13, 30, 80, 28, 90, 11, 43, 54, 64};
1071 1067
1072 G4double SnW[38] = {0, 0, 0, 0, 0, 0, 0, 0, 1068 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 1069 8.59, 32.58, 0, 4.63, 0, 5.79, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1074 1070
1075 AddElement("Sn", 50, 38, *SnN, *SnA, *SnS, 1071 AddElement("Sn", 50, 38, *SnN, *SnA, *SnS, *SnW);
1076 1072
1077 // Z = 51 --------------------------------- 1073 // Z = 51 --------------------------------------------------------------------
1078 G4int SbN[37] = {103, 104, 105, 106, 107, 1 1074 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, 1 1075 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137,
1080 138, 139}; 1076 138, 139};
1081 1077
1082 G4double SbA[37] = {102.94, 103.936, 104.93 1078 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.9 1079 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.9 1080 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.9 1081 130.912, 131.914, 132.915, 133.921, 134.925, 135.931, 136.935, 137.941, 138.946};
1086 1082
1087 G4double SbS[37] = {54, 39, 17, 34, 32, 22, 1083 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, 1084 22, 22, 3, 30, 6, 27, 23, 27, 80, 25, 80, 60, 11, 32, 43, 54, 64};
1089 1085
1090 G4double SbW[37] = {0, 0, 0, 0, 0, 0, 0, 0, 1086 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} 1087 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1092 1088
1093 AddElement("Sb", 51, 37, *SbN, *SbA, *SbS, 1089 AddElement("Sb", 51, 37, *SbN, *SbA, *SbS, *SbW);
1094 1090
1095 // Z = 52 --------------------------------- 1091 // Z = 52 --------------------------------------------------------------------
1096 G4int TeN[37] = {106, 107, 108, 109, 110, 1 1092 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, 1 1093 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140,
1098 141, 142}; 1094 141, 142};
1099 1095
1100 G4double TeA[37] = {105.938, 106.935, 107.9 1096 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.9 1097 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.9 1098 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.92 1099 133.912, 134.916, 135.92, 136.925, 137.929, 138.935, 139.939, 140.944, 141.948};
1104 1100
1105 G4double TeS[37] = {43, 32, 16, 80, 60, 80, 1101 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, 1102 20, 4, 19, 3, 21, 22, 12, 80, 40, 10, 50, 13, 22, 43, 54, 54, 64};
1107 1103
1108 G4double TeW[37] = {0, 0, 0, 0, 0, 0, 0, 0, 1104 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, 1105 18.84, 0, 31.74, 0, 34.08, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1110 1106
1111 AddElement("Te", 52, 37, *TeN, *TeA, *TeS, 1107 AddElement("Te", 52, 37, *TeN, *TeA, *TeS, *TeW);
1112 1108
1113 // Z = 53 --------------------------------- 1109 // Z = 53 --------------------------------------------------------------------
1114 G4int IN[37] = {108, 109, 110, 111, 112, 11 1110 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, 1 1111 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142,
1116 143, 144}; 1112 143, 144};
1117 1113
1118 G4double IA[37] = {107.943, 108.938, 109.93 1114 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 1115 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.9 1116 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.94 1117 138.926, 139.931, 140.935, 141.94, 142.944, 143.95};
1122 1118
1123 G4double IS[37] = {39, 16, 33, 32, 23, 60, 1119 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, 1120 4, 4, 12, 11, 28, 16, 25, 50, 30, 90, 30, 23, 32, 43, 43, 54};
1125 1121
1126 G4double IW[37] = {0, 0, 0, 0, 0, 0, 0, 0, 1122 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}; 1123 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1128 1124
1129 AddElement("I", 53, 37, *IN, *IA, *IS, *IW) 1125 AddElement("I", 53, 37, *IN, *IA, *IS, *IW);
1130 1126
1131 // Z = 54 --------------------------------- 1127 // Z = 54 --------------------------------------------------------------------
1132 G4int XeN[38] = {110, 111, 112, 113, 114, 1 1128 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, 1 1129 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144,
1134 145, 146, 147}; 1130 145, 146, 147};
1135 1131
1136 G4double XeA[38] = {109.944, 110.942, 111.9 1132 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.9 1133 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.9 1134 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.9 1135 137.914, 138.919, 139.922, 140.927, 141.93, 142.935, 143.938, 144.944, 145.947, 146.953};
1140 1136
1141 G4double XeS[38] = {43, 33, 16, 10, 22, 26, 1137 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, 1 1138 10, 10, 12, 4, 9, 11, 8, 8, 40, 23, 70, 10, 11, 24, 34, 43, 43, 54};
1143 1139
1144 G4double XeW[38] = {0, 0, 0, 0, 0, 0, 0, 0, 1140 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, 1141 21.18, 26.89, 0, 10.44, 0, 8.87, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1146 1142
1147 AddElement("Xe", 54, 38, *XeN, *XeA, *XeS, 1143 AddElement("Xe", 54, 38, *XeN, *XeA, *XeS, *XeW);
1148 1144
1149 // Z = 55 --------------------------------- 1145 // Z = 55 --------------------------------------------------------------------
1150 G4int CsN[40] = {112, 113, 114, 115, 116, 1 1146 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, 1 1147 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}; 1148 147, 148, 149, 150, 151};
1153 1149
1154 G4double CsA[40] = {111.95, 112.945, 113.94 1150 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.9 1151 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.9 1152 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.93 1153 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}; 1154 150.962};
1159 1155
1160 G4double CsS[40] = {33, 16, 33, 46, 38, 60, 1156 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, 1157 3, 3, 3, 4, 3, 10, 5, 9, 11, 11, 24, 30, 50, 90, 16, 63, 32, 54, 75};
1162 1158
1163 G4double CsW[40] = {0, 0, 0, 0, 0, 0, 0, 0, 1159 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, 1160 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1165 1161
1166 AddElement("Cs", 55, 40, *CsN, *CsA, *CsS, 1162 AddElement("Cs", 55, 40, *CsN, *CsA, *CsS, *CsW);
1167 1163
1168 // Z = 56 --------------------------------- 1164 // Z = 56 --------------------------------------------------------------------
1169 G4int BaN[40] = {114, 115, 116, 117, 118, 1 1165 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, 1 1166 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}; 1167 149, 150, 151, 152, 153};
1172 1168
1173 G4double BaA[40] = {113.951, 114.948, 115.9 1169 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.91 1170 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.9 1171 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.9 1172 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}; 1173 151.954, 152.96};
1178 1174
1179 G4double BaS[40] = {48, 64, 54, 70, 54, 109 1175 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 1176 3, 3, 3, 3, 3, 3, 9, 9, 7, 14, 15, 60, 80, 10, 15, 43, 54, 64, 75, 97};
1181 1177
1182 G4double BaW[40] = {0, 0, 0, 0, 0, 0, 0, 0, 1178 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, 1179 6.592, 7.854, 11.232, 71.698, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1184 1180
1185 AddElement("Ba", 56, 40, *BaN, *BaA, *BaS, 1181 AddElement("Ba", 56, 40, *BaN, *BaA, *BaS, *BaW);
1186 1182
1187 // Z = 57 --------------------------------- 1183 // Z = 57 --------------------------------------------------------------------
1188 G4int LaN[39] = {117, 118, 119, 120, 121, 1 1184 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, 1 1185 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151,
1190 152, 153, 154, 155}; 1186 152, 153, 154, 155};
1191 1187
1192 G4double LaA[39] = {116.95, 117.947, 118.94 1188 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.9 1189 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.9 1190 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.9 1191 145.926, 146.928, 147.932, 148.934, 149.939, 150.942, 151.946, 152.949, 153.954, 154.958};
1196 1192
1197 G4double LaS[39] = {96, 86, 75, 64, 54, 54, 1193 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 1194 80, 50, 4, 3, 3, 5, 6, 16, 60, 70, 80, 80, 14, 32, 43, 54, 64, 75, 86, 97};
1199 1195
1200 G4double LaW[39] = {0, 0, 0, 0, 0, 0, 0, 0, 1196 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, 1197 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1202 1198
1203 AddElement("La", 57, 39, *LaN, *LaA, *LaS, 1199 AddElement("La", 57, 39, *LaN, *LaA, *LaS, *LaW);
1204 1200
1205 // Z = 58 --------------------------------- 1201 // Z = 58 --------------------------------------------------------------------
1206 G4int CeN[39] = {119, 120, 121, 122, 123, 1 1202 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, 1 1203 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153,
1208 154, 155, 156, 157}; 1204 154, 155, 156, 157};
1209 1205
1210 G4double CeA[39] = {118.953, 119.947, 120.9 1206 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.9 1207 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.9 1208 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.93 1209 146.923, 147.924, 148.928, 149.93, 150.934, 151.936, 152.941, 153.943, 154.948, 155.951,
1214 156.956}; 1210 156.956};
1215 1211
1216 G4double CeS[39] = {97, 86, 75, 64, 54, 54, 1212 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, 1213 11, 8, 3, 3, 4, 4, 4, 40, 70, 60, 13, 80, 13, 32, 43, 54, 64, 75, 86, 97};
1218 1214
1219 G4double CeW[39] = {0, 0, 0, 0, 0, 0, 0, 0, 1215 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, 1216 0, 11.114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1221 1217
1222 AddElement("Ce", 58, 39, *CeN, *CeA, *CeS, 1218 AddElement("Ce", 58, 39, *CeN, *CeA, *CeS, *CeW);
1223 1219
1224 // Z = 59 --------------------------------- 1220 // Z = 59 --------------------------------------------------------------------
1225 G4int PrN[39] = {121, 122, 123, 124, 125, 1 1221 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, 1 1222 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155,
1227 156, 157, 158, 159}; 1223 156, 157, 158, 159};
1228 1224
1229 G4double PrA[39] = {120.955, 121.952, 122.9 1225 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.91 1226 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.91 1227 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.9 1228 148.924, 149.927, 150.928, 151.932, 152.934, 153.937, 154.94, 155.944, 156.947, 157.952,
1233 158.955}; 1229 158.955};
1234 1230
1235 G4double PrS[39] = {86, 86, 75, 64, 54, 54, 1231 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 1232 3, 3, 3, 4, 8, 60, 40, 10, 11, 90, 40, 32, 32, 43, 54, 64, 75, 86, 97};
1237 1233
1238 G4double PrW[39] = {0, 0, 0, 0, 0, 0, 0, 0, 1234 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} 1235 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1240 1236
1241 AddElement("Pr", 59, 39, *PrN, *PrA, *PrS, 1237 AddElement("Pr", 59, 39, *PrN, *PrA, *PrS, *PrW);
1242 1238
1243 // Z = 60 --------------------------------- 1239 // Z = 60 --------------------------------------------------------------------
1244 G4int NdN[36] = {126, 127, 128, 129, 130, 1 1240 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, 1 1241 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160,
1246 161}; 1242 161};
1247 1243
1248 G4double NdA[36] = {125.943, 126.941, 127.9 1244 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.9 1245 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.9 1246 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.9 1247 154.933, 155.935, 156.939, 157.942, 158.946, 159.949, 160.954};
1252 1248
1253 G4double NdS[36] = {75, 64, 64, 39, 54, 50, 1249 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 1250 3, 3, 3, 4, 4, 30, 29, 12, 16, 43, 54, 64, 75, 86, 97};
1255 1251
1256 G4double NdW[36] = {0, 0, 0, 0, 0, 0, 0, 0, 1252 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 1253 0, 5.7, 0, 5.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1258 1254
1259 AddElement("Nd", 60, 36, *NdN, *NdA, *NdS, 1255 AddElement("Nd", 60, 36, *NdN, *NdA, *NdS, *NdW);
1260 1256
1261 // Z = 61 --------------------------------- 1257 // Z = 61 --------------------------------------------------------------------
1262 G4int PmN[36] = {128, 129, 130, 131, 132, 1 1258 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, 1 1259 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162,
1264 163}; 1260 163};
1265 1261
1266 G4double PmA[36] = {127.948, 128.943, 129.9 1262 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.9 1263 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.9 1264 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.9 1265 156.933, 157.937, 158.939, 159.943, 160.946, 161.95, 162.954};
1270 1266
1271 G4double PmS[36] = {97, 86, 75, 64, 54, 54, 1267 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, 1268 5, 22, 6, 80, 12, 80, 30, 40, 32, 43, 54, 64, 75, 86, 97};
1273 1269
1274 G4double PmW[36] = {0, 0, 0, 0, 0, 0, 0, 0, 1270 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}; 1271 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1276 1272
1277 AddElement("Pm", 61, 36, *PmN, *PmA, *PmS, 1273 AddElement("Pm", 61, 36, *PmN, *PmA, *PmS, *PmW);
1278 1274
1279 // Z = 62 --------------------------------- 1275 // Z = 62 --------------------------------------------------------------------
1280 G4int SmN[36] = {130, 131, 132, 133, 134, 1 1276 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, 1 1277 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164,
1282 165}; 1278 165};
1283 1279
1284 G4double SmA[36] = {129.949, 130.946, 131.9 1280 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.9 1281 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 1282 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.9 1283 158.933, 159.935, 160.939, 161.941, 162.945, 163.948, 164.953};
1288 1284
1289 G4double SmS[36] = {97, 97, 75, 64, 54, 54, 1285 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, 7 1286 3, 3, 3, 3, 10, 50, 80, 32, 43, 54, 64, 75, 86, 97};
1291 1287
1292 G4double SmW[36] = {0, 0, 0, 0, 0, 0, 0, 0, 1288 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, 1289 7.38, 0, 26.75, 0, 22.75, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1294 1290
1295 AddElement("Sm", 62, 36, *SmN, *SmA, *SmS, 1291 AddElement("Sm", 62, 36, *SmN, *SmA, *SmS, *SmW);
1296 1292
1297 // Z = 63 --------------------------------- 1293 // Z = 63 --------------------------------------------------------------------
1298 G4int EuN[36] = {132, 133, 134, 135, 136, 1 1294 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, 1 1295 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166,
1300 167}; 1296 167};
1301 1297
1302 G4double EuA[36] = {131.954, 132.949, 133.9 1298 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.91 1299 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.9 1300 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.9 1301 161.937, 162.939, 163.943, 164.946, 165.95, 166.953};
1306 1302
1307 G4double EuS[36] = {97, 97, 75, 64, 54, 54, 1303 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, 1304 3, 3, 3, 6, 7, 80, 9, 22, 32, 43, 54, 64, 75, 86, 97};
1309 1305
1310 G4double EuW[36] = {0, 0, 0, 0, 0, 0, 0, 0, 1306 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}; 1307 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1312 1308
1313 AddElement("Eu", 63, 36, *EuN, *EuA, *EuS, 1309 AddElement("Eu", 63, 36, *EuN, *EuA, *EuS, *EuW);
1314 1310
1315 // Z = 64 --------------------------------- 1311 // Z = 64 --------------------------------------------------------------------
1316 G4int GdN[34] = {136, 137, 138, 139, 140, 1 1312 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, 1 1313 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169};
1318 1314
1319 G4double GdA[34] = {135.947, 136.945, 137.9 1315 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.9 1316 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.9 1317 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.9 1318 164.939, 165.942, 166.946, 167.948, 168.953};
1323 1319
1324 G4double GdS[34] = {75, 64, 54, 54, 43, 32, 1320 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}; 1321 3, 3, 3, 5, 32, 43, 54, 64, 64, 75, 86};
1326 1322
1327 G4double GdW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1323 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, 1324 15.65, 24.84, 0, 21.86, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1329 1325
1330 AddElement("Gd", 64, 34, *GdN, *GdA, *GdS, 1326 AddElement("Gd", 64, 34, *GdN, *GdA, *GdS, *GdW);
1331 1327
1332 // Z = 65 --------------------------------- 1328 // Z = 65 --------------------------------------------------------------------
1333 G4int TbN[34] = {138, 139, 140, 141, 142, 1 1329 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, 1 1330 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171};
1335 1331
1336 G4double TbA[34] = {137.953, 138.948, 139.9 1332 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.9 1333 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.9 1334 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 1335 165.938, 166.94, 167.944, 168.946, 169.95, 170.953};
1340 1336
1341 G4double TbS[34] = {86, 75, 96, 64, 82, 43, 1337 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, 7 1338 3, 3, 3, 40, 5, 11, 21, 32, 43, 54, 64, 75, 86};
1343 1339
1344 G4double TbW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1340 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}; 1341 0, 0, 0, 0, 0, 0, 0, 0, 0};
1346 1342
1347 AddElement("Tb", 65, 34, *TbN, *TbA, *TbS, 1343 AddElement("Tb", 65, 34, *TbN, *TbA, *TbS, *TbW);
1348 1344
1349 // Z = 66 --------------------------------- 1345 // Z = 66 --------------------------------------------------------------------
1350 G4int DyN[34] = {140, 141, 142, 143, 144, 1 1346 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, 1 1347 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173};
1352 1348
1353 G4double DyA[34] = {139.954, 140.951, 141.9 1349 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.9 1350 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.9 1351 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.94 1352 167.937, 168.94, 169.943, 170.946, 171.949, 172.953};
1357 1353
1358 G4double DyS[34] = {97, 75, 85, 54, 43, 32, 1354 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 1355 3, 3, 3, 3, 3, 60, 32, 32, 43, 54, 64, 75};
1360 1356
1361 G4double DyW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1357 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, 1358 25.51, 24.9, 28.18, 0, 0, 0, 0, 0, 0, 0, 0, 0};
1363 1359
1364 AddElement("Dy", 66, 34, *DyN, *DyA, *DyS, 1360 AddElement("Dy", 66, 34, *DyN, *DyA, *DyS, *DyW);
1365 1361
1366 // Z = 67 --------------------------------- 1362 // Z = 67 --------------------------------------------------------------------
1367 G4int HoN[34] = {142, 143, 144, 145, 146, 1 1363 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, 1 1364 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175};
1369 1365
1370 G4double HoA[34] = {141.96, 142.955, 143.95 1366 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.92 1367 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 1368 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}; 1369 172.947, 173.951, 174.954};
1374 1370
1375 G4double HoS[34] = {107, 75, 64, 64, 54, 43 1371 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, 1372 5, 3, 3, 3, 3, 6, 30, 22, 50, 64, 43, 43, 54, 64};
1377 1373
1378 G4double HoW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1374 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}; 1375 0, 0, 0, 0, 0, 0, 0, 0, 0};
1380 1376
1381 AddElement("Ho", 67, 34, *HoN, *HoA, *HoS, 1377 AddElement("Ho", 67, 34, *HoN, *HoA, *HoS, *HoW);
1382 1378
1383 // Z = 68 --------------------------------- 1379 // Z = 68 --------------------------------------------------------------------
1384 G4int ErN[34] = {144, 145, 146, 147, 148, 1 1380 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, 1 1381 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177};
1386 1382
1387 G4double ErA[34] = {143.961, 144.957, 145.9 1383 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.9 1384 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 1385 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}; 1386 173.944, 174.948, 175.95, 176.954};
1391 1387
1392 G4double ErS[34] = {86, 75, 64, 54, 43, 51, 1388 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, 6 1389 4, 3, 3, 3, 3, 3, 3, 5, 21, 32, 43, 43, 64};
1394 1390
1395 G4double ErW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1391 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, 1392 22.93, 26.78, 0, 14.93, 0, 0, 0, 0, 0, 0, 0};
1397 1393
1398 AddElement("Er", 68, 34, *ErN, *ErA, *ErS, 1394 AddElement("Er", 68, 34, *ErN, *ErA, *ErS, *ErW);
1399 1395
1400 // Z = 69 --------------------------------- 1396 // Z = 69 --------------------------------------------------------------------
1401 G4int TmN[34] = {146, 147, 148, 149, 150, 1 1397 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, 1 1398 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179};
1403 1399
1404 G4double TmA[34] = {145.966, 146.961, 147.9 1400 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.9 1401 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.9 1402 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.9 1403 173.942, 174.944, 175.947, 176.949, 177.953, 178.955};
1408 1404
1409 G4double TmS[34] = {75, 64, 75, 64, 54, 15, 1405 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, 1406 12, 3, 4, 3, 3, 3, 7, 6, 50, 50, 11, 32, 43, 54};
1411 1407
1412 G4double TmW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1408 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}; 1409 0, 0, 0, 0, 0, 0, 0, 0, 0};
1414 1410
1415 AddElement("Tm", 69, 34, *TmN, *TmA, *TmS, 1411 AddElement("Tm", 69, 34, *TmN, *TmA, *TmS, *TmW);
1416 1412
1417 // Z = 70 --------------------------------- 1413 // Z = 70 --------------------------------------------------------------------
1418 G4int YbN[34] = {148, 149, 150, 151, 152, 1 1414 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, 1 1415 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181};
1420 1416
1421 G4double YbA[34] = {147.967, 148.963, 149.9 1417 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 1418 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.9 1419 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.95 1420 176.945, 177.947, 178.95, 179.952, 180.956};
1425 1421
1426 G4double YbS[34] = {86, 75, 64, 34, 38, 32, 1422 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, 1423 5, 5, 3, 3, 30, 30, 30, 30, 3, 3, 11, 32, 43, 43};
1428 1424
1429 G4double YbW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1425 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, 1426 14.28, 21.83, 16.13, 31.83, 0, 12.76, 0, 0, 0, 0, 0};
1431 1427
1432 AddElement("Yb", 70, 34, *YbN, *YbA, *YbS, 1428 AddElement("Yb", 70, 34, *YbN, *YbA, *YbS, *YbW);
1433 1429
1434 // Z = 71 --------------------------------- 1430 // Z = 71 --------------------------------------------------------------------
1435 G4int LuN[35] = {150, 151, 152, 153, 154, 1 1431 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, 1 1432 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184};
1437 1433
1438 G4double LuA[35] = {149.973, 150.967, 151.9 1434 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.9 1435 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.9 1436 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.95 1437 178.947, 179.95, 180.952, 181.955, 182.958, 183.961};
1442 1438
1443 G4double LuS[35] = {75, 65, 75, 64, 54, 14, 1439 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 1440 20, 3, 4, 3, 30, 28, 28, 28, 3, 6, 80, 32, 32, 32, 43};
1445 1441
1446 G4double LuW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1442 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}; 1443 97.41, 2.59, 0, 0, 0, 0, 0, 0, 0, 0};
1448 1444
1449 AddElement("Lu", 71, 35, *LuN, *LuA, *LuS, 1445 AddElement("Lu", 71, 35, *LuN, *LuA, *LuS, *LuW);
1450 1446
1451 // Z = 72 --------------------------------- 1447 // Z = 72 --------------------------------------------------------------------
1452 G4int HfN[33] = {154, 155, 156, 157, 158, 1 1448 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, 1 1449 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186};
1454 1450
1455 G4double HfA[33] = {153.964, 154.963, 155.9 1451 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.9 1452 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.94 1453 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}; 1454 182.954, 183.955, 184.959, 185.961};
1459 1455
1460 G4double HfS[33] = {75, 64, 38, 32, 11, 32, 1456 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, 1457 11, 3, 3, 29, 27, 27, 27, 27, 28, 7, 30, 40, 32, 32};
1462 1458
1463 G4double HfW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1459 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, 1460 18.6, 27.28, 13.62, 35.08, 0, 0, 0, 0, 0, 0};
1465 1461
1466 AddElement("Hf", 72, 33, *HfN, *HfA, *HfS, 1462 AddElement("Hf", 72, 33, *HfN, *HfA, *HfS, *HfW);
1467 1463
1468 // Z = 73 --------------------------------- 1464 // Z = 73 --------------------------------------------------------------------
1469 G4int TaN[33] = {156, 157, 158, 159, 160, 1 1465 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, 1 1466 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188};
1471 1467
1472 G4double TaA[33] = {155.972, 156.968, 157.9 1468 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.94 1469 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.9 1470 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.96 1471 183.954, 184.956, 185.959, 186.96, 187.964};
1476 1472
1477 G4double TaS[33] = {64, 64, 54, 13, 33, 60, 1473 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, 1474 11, 11, 4, 11, 6, 3, 3, 3, 3, 28, 15, 60, 32, 32};
1479 1475
1480 G4double TaW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1476 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}; 1477 99.988, 0, 0, 0, 0, 0, 0, 0};
1482 1478
1483 AddElement("Ta", 73, 33, *TaN, *TaA, *TaS, 1479 AddElement("Ta", 73, 33, *TaN, *TaA, *TaS, *TaW);
1484 1480
1485 // Z = 74 --------------------------------- 1481 // Z = 74 --------------------------------------------------------------------
1486 G4int WN[33] = {158, 159, 160, 161, 162, 16 1482 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, 1 1483 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190};
1488 1484
1489 G4double WA[33] = {157.974, 158.972, 159.96 1485 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.9 1486 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.9 1487 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.9 1488 185.954, 186.957, 187.958, 188.962, 189.963};
1493 1489
1494 G4double WS[33] = {75, 64, 38, 33, 11, 33, 1490 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, 1491 11, 17, 5, 6, 3, 29, 29, 30, 3, 3, 4, 21, 24};
1496 1492
1497 G4double WW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1493 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} 1494 26.5, 14.31, 30.64, 0, 28.43, 0, 0, 0, 0};
1499 1495
1500 AddElement("W", 74, 33, *WN, *WA, *WS, *WW) 1496 AddElement("W", 74, 33, *WN, *WA, *WS, *WW);
1501 1497
1502 // Z = 75 --------------------------------- 1498 // Z = 75 --------------------------------------------------------------------
1503 G4int ReN[33] = {160, 161, 162, 163, 164, 1 1499 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, 1 1500 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192};
1505 1501
1506 G4double ReA[33] = {159.981, 160.978, 161.9 1502 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.9 1503 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 1504 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}; 1505 188.959, 189.962, 190.963, 191.966};
1510 1506
1511 G4double ReS[33] = {64, 64, 55, 12, 33, 80, 1507 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, 2 1508 60, 40, 15, 11, 9, 6, 30, 3, 30, 30, 9, 23, 11, 21};
1513 1509
1514 G4double ReW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1510 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}; 1511 37.4, 0, 62.6, 0, 0, 0, 0, 0};
1516 1512
1517 AddElement("Re", 75, 33, *ReN, *ReA, *ReS, 1513 AddElement("Re", 75, 33, *ReN, *ReA, *ReS, *ReW);
1518 1514
1519 // Z = 76 --------------------------------- 1515 // Z = 76 --------------------------------------------------------------------
1520 G4int OsN[35] = {162, 163, 164, 165, 166, 1 1516 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, 1 1517 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196};
1522 1518
1523 G4double OsA[35] = {161.984, 162.982, 163.9 1519 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 1520 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.9 1521 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.9 1522 190.961, 191.961, 192.964, 193.965, 194.968, 195.97};
1527 1523
1528 G4double OsS[35] = {75, 64, 38, 33, 11, 33, 1524 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, 1525 22, 27, 11, 3, 3, 3, 30, 30, 30, 3, 3, 4, 4, 4, 54, 40};
1530 1526
1531 G4double OsW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1527 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 1528 1.59, 1.96, 13.24, 16.15, 26.26, 0, 40.78, 0, 0, 0, 0};
1533 1529
1534 AddElement("Os", 76, 35, *OsN, *OsA, *OsS, 1530 AddElement("Os", 76, 35, *OsN, *OsA, *OsS, *OsW);
1535 1531
1536 // Z = 77 --------------------------------- 1532 // Z = 77 --------------------------------------------------------------------
1537 G4int IrN[35] = {165, 166, 167, 168, 169, 1 1533 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, 1 1534 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199};
1539 1535
1540 G4double IrA[35] = {164.988, 165.986, 166.9 1536 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.9 1537 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.9 1538 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.9 1539 192.963, 193.965, 194.966, 195.968, 196.97, 197.972, 198.974};
1544 1540
1545 G4double IrS[35] = {43, 55, 11, 35, 10, 16, 1541 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, 1542 29, 21, 22, 7, 8, 14, 21, 3, 3, 3, 3, 3, 40, 22, 21, 40};
1547 1543
1548 G4double IrW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1544 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}; 1545 37.3, 0, 62.7, 0, 0, 0, 0, 0, 0};
1550 1546
1551 AddElement("Ir", 77, 35, *IrN, *IrA, *IrS, 1547 AddElement("Ir", 77, 35, *IrN, *IrA, *IrS, *IrW);
1552 1548
1553 // Z = 78 --------------------------------- 1549 // Z = 78 --------------------------------------------------------------------
1554 G4int PtN[35] = {168, 169, 170, 171, 172, 1 1550 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, 1 1551 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202};
1556 1552
1557 G4double PtA[35] = {167.988, 168.986, 169.9 1553 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.9 1554 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.9 1555 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.9 1556 195.965, 196.967, 197.968, 198.971, 199.971, 200.975, 201.976};
1561 1557
1562 G4double PtS[35] = {38, 34, 11, 34, 40, 11, 1558 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, 1559 20, 6, 12, 7, 5, 4, 3, 3, 3, 3, 3, 4, 5, 22, 50, 32};
1564 1560
1565 G4double PtW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1561 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.16 1562 0.782, 0, 32.967, 33.832, 25.242, 0, 7.163, 0, 0, 0, 0};
1567 1563
1568 AddElement("Pt", 78, 35, *PtN, *PtA, *PtS, 1564 AddElement("Pt", 78, 35, *PtN, *PtA, *PtS, *PtW);
1569 1565
1570 // Z = 79 --------------------------------- 1566 // Z = 79 --------------------------------------------------------------------
1571 G4int AuN[35] = {171, 172, 173, 174, 175, 1 1567 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, 1 1568 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205};
1573 1569
1574 G4double AuA[35] = {170.992, 171.99, 172.98 1570 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 1571 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.9 1572 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.9 1573 200.972, 201.974, 202.975, 203.978, 204.98};
1578 1574
1579 G4double AuS[35] = {27, 36, 11, 16, 26, 43, 1575 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, 1576 17, 50, 17, 10, 12, 3, 4, 3, 3, 3, 60, 5, 18, 5, 22, 32};
1581 1577
1582 G4double AuW[35] = {0, 0, 0, 0, 0, 0, 0, 0, 1578 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}; 1579 100, 0, 0, 0, 0, 0, 0, 0, 0};
1584 1580
1585 AddElement("Au", 79, 35, *AuN, *AuA, *AuS, 1581 AddElement("Au", 79, 35, *AuN, *AuA, *AuS, *AuW);
1586 1582
1587 // Z = 80 --------------------------------- 1583 // Z = 80 --------------------------------------------------------------------
1588 G4int HgN[34] = {175, 176, 177, 178, 179, 1 1584 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, 1 1585 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208};
1590 1586
1591 G4double HgA[34] = {174.991, 175.987, 176.9 1587 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.9 1588 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.9 1589 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.9 1590 202.973, 203.973, 204.976, 205.977, 206.983, 207.986};
1595 1591
1596 G4double HgS[34] = {34, 40, 12, 16, 33, 21, 1592 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, 1593 25, 50, 4, 5, 3, 3, 3, 3, 3, 3, 3, 5, 22, 16, 32};
1598 1594
1599 G4double HgW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1595 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, 1596 16.87, 23.1, 13.18, 29.86, 0, 6.87, 0, 0, 0, 0};
1601 1597
1602 AddElement("Hg", 80, 34, *HgN, *HgA, *HgS, 1598 AddElement("Hg", 80, 34, *HgN, *HgA, *HgS, *HgW);
1603 1599
1604 // Z = 81 --------------------------------- 1600 // Z = 81 --------------------------------------------------------------------
1605 G4int TlN[34] = {177, 178, 179, 180, 181, 1 1601 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, 2 1602 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210};
1607 1603
1608 G4double TlA[34] = {176.997, 177.995, 178.9 1604 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.9 1605 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, 1606 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.9 1607 205.976, 206.977, 207.982, 208.985, 209.99};
1612 1608
1613 G4double TlS[34] = {24, 23, 15, 48, 41, 43, 1609 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, 1610 15, 30, 90, 11, 7, 16, 16, 3, 3, 3, 3, 6, 3, 10, 12};
1615 1611
1616 G4double TlW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1612 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}; 1613 29.524, 0, 70.476, 0, 0, 0, 0, 0};
1618 1614
1619 AddElement("Tl", 81, 34, *TlN, *TlA, *TlS, 1615 AddElement("Tl", 81, 34, *TlN, *TlA, *TlS, *TlW);
1620 1616
1621 // Z = 82 --------------------------------- 1617 // Z = 82 --------------------------------------------------------------------
1622 G4int PbN[34] = {181, 182, 183, 184, 185, 1 1618 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, 2 1619 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214};
1624 1620
1625 G4double PbA[34] = {180.997, 181.993, 182.9 1621 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.9 1622 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.9 1623 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.9 1624 208.981, 209.984, 210.989, 211.992, 212.996, 214};
1629 1625
1630 G4double PbS[34] = {17, 18, 33, 21, 33, 50, 1626 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 1627 14, 30, 11, 7, 3, 3, 3, 3, 3, 3, 3, 3, 29, 11, 27};
1632 1628
1633 G4double PbW[34] = {0, 0, 0, 0, 0, 0, 0, 0, 1629 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}; 1630 24.1, 22.1, 52.4, 0, 0, 0, 0, 0, 0};
1635 1631
1636 AddElement("Pb", 82, 34, *PbN, *PbA, *PbS, 1632 AddElement("Pb", 82, 34, *PbN, *PbA, *PbS, *PbW);
1637 1633
1638 // Z = 83 --------------------------------- 1634 // Z = 83 --------------------------------------------------------------------
1639 G4int BiN[32] = {185, 186, 187, 188, 189, 1 1635 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, 2 1636 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216};
1641 1637
1642 G4double BiA[32] = {184.998, 185.996, 186.9 1638 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.9 1639 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.9 1640 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}; 1641 213.999, 215.002, 216.006};
1646 1642
1647 G4double BiS[32] = {24, 48, 41, 32, 43, 39, 1643 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 1644 28, 8, 9, 4, 4, 3, 3, 6, 3, 8, 12, 10, 11};
1649 1645
1650 G4double BiW[32] = {0, 0, 0, 0, 0, 0, 0, 0, 1646 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}; 1647 0, 0, 0, 0, 0, 0, 0};
1652 1648
1653 AddElement("Bi", 83, 32, *BiN, *BiA, *BiS, 1649 AddElement("Bi", 83, 32, *BiN, *BiA, *BiS, *BiW);
1654 1650
1655 // Z = 84 --------------------------------- 1651 // Z = 84 --------------------------------------------------------------------
1656 G4int PoN[29] = {190, 191, 192, 193, 194, 1 1652 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, 2 1653 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218};
1658 1654
1659 G4double PoA[29] = {189.995, 190.995, 191.9 1655 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.9 1656 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.9 1657 208.982, 209.983, 210.987, 211.989, 212.993, 213.995, 214.999, 216.002, 217.006, 218.009};
1662 1658
1663 G4double PoS[29] = {51, 32, 22, 30, 22, 24, 1659 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}; 1660 3, 3, 3, 4, 3, 3, 29, 11, 27};
1665 1661
1666 G4double PoW[29] = { 1662 G4double PoW[29] = {
1667 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1663 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 1664
1669 AddElement("Po", 84, 29, *PoN, *PoA, *PoS, 1665 AddElement("Po", 84, 29, *PoN, *PoA, *PoS, *PoW);
1670 1666
1671 // Z = 85 --------------------------------- 1667 // Z = 85 --------------------------------------------------------------------
1672 G4int AtN[31] = {193, 194, 195, 196, 197, 1 1668 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, 2 1669 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223};
1674 1670
1675 G4double AtA[31] = {193.000, 193.999, 194.9 1671 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.9 1672 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.9 1673 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}; 1674 221.018, 222.022, 223.025};
1679 1675
1680 G4double AtS[31] = {43, 43, 43, 25, 38, 46, 1676 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}; 1677 5, 8, 5, 8, 13, 90, 12, 32, 32, 43};
1682 1678
1683 G4double AtW[31] = { 1679 G4double AtW[31] = {
1684 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1680 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 1681
1686 AddElement("At", 85, 31, *AtN, *AtA, *AtS, 1682 AddElement("At", 85, 31, *AtN, *AtA, *AtS, *AtW);
1687 1683
1688 // Z = 86 --------------------------------- 1684 // Z = 86 --------------------------------------------------------------------
1689 G4int RnN[33] = {196, 197, 198, 199, 200, 2 1685 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, 2 1686 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228};
1691 1687
1692 G4double RnA[33] = {196.002, 197.002, 197.9 1688 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 1689 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.0 1690 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}; 1691 225.028, 226.031, 227.035, 228.038};
1696 1692
1697 G4double RnS[33] = {22, 30, 22, 24, 19, 21, 1693 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, 1694 5, 4, 3, 29, 11, 27, 32, 32, 32, 43, 45, 50};
1699 1695
1700 G4double RnW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1696 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}; 1697 100, 0, 0, 0, 0, 0, 0};
1702 1698
1703 AddElement("Rn", 86, 33, *RnN, *RnA, *RnS, 1699 AddElement("Rn", 86, 33, *RnN, *RnA, *RnS, *RnW);
1704 1700
1705 // Z = 87 --------------------------------- 1701 // Z = 87 --------------------------------------------------------------------
1706 G4int FrN[33] = {200, 201, 202, 203, 204, 2 1702 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, 2 1703 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232};
1708 1704
1709 G4double FrA[33] = {200.006, 201.004, 202.0 1705 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.9 1706 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.0 1707 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.0 1708 228.036, 229.038, 230.043, 231.045, 232.05};
1713 1709
1714 G4double FrS[33] = {25, 38, 46, 25, 22, 25, 1710 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 1711 8, 23, 29, 50, 11, 10, 10, 22, 39, 48, 56, 69};
1716 1712
1717 G4double FrW[33] = {0, 0, 0, 0, 0, 0, 0, 0, 1713 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}; 1714 0, 0, 0, 0, 0, 0, 0, 0};
1719 1715
1720 AddElement("Fr", 87, 33, *FrN, *FrA, *FrS, 1716 AddElement("Fr", 87, 33, *FrN, *FrA, *FrS, *FrW);
1721 1717
1722 // Z = 88 --------------------------------- 1718 // Z = 88 --------------------------------------------------------------------
1723 G4int RaN[32] = {203, 204, 205, 206, 207, 2 1719 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, 2 1720 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234};
1725 1721
1726 G4double RaA[32] = {203.009, 204.006, 205.0 1722 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, 1723 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.03 1724 224.02, 225.024, 226.025, 227.029, 228.031, 229.035, 230.037, 231.041, 232.044, 233.048,
1729 234.051}; 1725 234.051};
1730 1726
1731 G4double RaS[32] = {25, 19, 22, 16, 45, 15, 1727 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 1728 29, 3, 27, 27, 27, 70, 40, 32, 39, 50, 58};
1733 1729
1734 G4double RaW[32] = {0, 0, 0, 0, 0, 0, 0, 0, 1730 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}; 1731 0, 0, 0, 0, 0, 0, 0};
1736 1732
1737 AddElement("Ra", 88, 32, *RaN, *RaA, *RaS, 1733 AddElement("Ra", 88, 32, *RaN, *RaA, *RaS, *RaW);
1738 1734
1739 // Z = 89 --------------------------------- 1735 // Z = 89 --------------------------------------------------------------------
1740 G4int AcN[30] = {207, 208, 209, 210, 211, 2 1736 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, 2 1737 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236};
1742 1738
1743 G4double AcA[30] = {207.012, 208.011, 209.0 1739 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.0 1740 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.0 1741 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}; 1742 235.051, 236.055};
1747 1743
1748 G4double AcS[30] = {25, 22, 26, 20, 14, 10, 1744 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}; 1745 28, 50, 11, 11, 11, 32, 43, 45, 54};
1750 1746
1751 G4double AcW[30] = { 1747 G4double AcW[30] = {
1752 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1748 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 1749
1754 AddElement("Ac", 89, 30, *AcN, *AcA, *AcS, 1750 AddElement("Ac", 89, 30, *AcN, *AcA, *AcS, *AcW);
1755 1751
1756 // Z = 90 --------------------------------- 1752 // Z = 90 --------------------------------------------------------------------
1757 G4int ThN[29] = {210, 211, 212, 213, 214, 2 1753 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, 2 1754 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238};
1759 1755
1760 G4double ThA[29] = {210.015, 211.015, 212.0 1756 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.0 1757 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.0 1758 228.029, 229.032, 230.033, 231.036, 232.038, 233.042, 234.044, 235.048, 236.05, 237.054,
1763 238.056}; 1759 238.056};
1764 1760
1765 G4double ThS[29] = {17, 45, 15, 14, 10, 70, 1761 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}; 1762 22, 22, 22, 22, 4, 50, 32, 39, 39};
1767 1763
1768 G4double ThW[29] = { 1764 G4double ThW[29] = {
1769 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1765 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 1766
1771 AddElement("Th", 90, 29, *ThN, *ThA, *ThS, 1767 AddElement("Th", 90, 29, *ThN, *ThA, *ThS, *ThW);
1772 1768
1773 // Z = 91 --------------------------------- 1769 // Z = 91 --------------------------------------------------------------------
1774 G4int PaN[28] = {213, 214, 215, 216, 217, 2 1770 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, 2 1771 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240};
1776 1772
1777 G4double PaA[28] = {213.021, 214.021, 215.0 1773 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.0 1774 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.04 1775 232.039, 233.04, 234.043, 235.045, 236.049, 237.051, 238.054, 239.057, 240.061};
1780 1776
1781 G4double PaS[28] = {27, 21, 15, 11, 80, 80, 1777 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}; 1778 24, 5, 50, 21, 11, 60, 32, 32};
1783 1779
1784 G4double PaW[28] = { 1780 G4double PaW[28] = {
1785 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1781 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 1782
1787 AddElement("Pa", 91, 28, *PaN, *PaA, *PaS, 1783 AddElement("Pa", 91, 28, *PaN, *PaA, *PaS, *PaW);
1788 1784
1789 // Z = 92 --------------------------------- 1785 // Z = 92 --------------------------------------------------------------------
1790 G4int UN[25] = {218, 219, 220, 221, 222, 22 1786 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, 2 1787 234, 235, 236, 237, 238, 239, 240, 241, 242};
1792 1788
1793 G4double UA[25] = {218.023, 219.025, 220.02 1789 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.0 1790 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.0 1791 236.046, 237.049, 238.051, 239.054, 240.057, 241.06, 242.063};
1796 1792
1797 G4double US[25] = { 1793 G4double US[25] = {
1798 10, 90, 22, 11, 11, 80, 27, 50, 20, 18, 1 1794 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 1795
1800 G4double UW[25] = { 1796 G4double UW[25] = {
1801 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1797 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 1798
1803 AddElement("U", 92, 25, *UN, *UA, *US, *UW) 1799 AddElement("U", 92, 25, *UN, *UA, *US, *UW);
1804 1800
1805 // Z = 93 --------------------------------- 1801 // Z = 93 --------------------------------------------------------------------
1806 G4int NpN[20] = {225, 226, 227, 228, 229, 2 1802 G4int NpN[20] = {225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240,
1807 241, 242, 243, 244}; 1803 241, 242, 243, 244};
1808 1804
1809 G4double NpA[20] = {225.034, 226.035, 227.0 1805 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.0 1806 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}; 1807 243.064, 244.068};
1812 1808
1813 G4double NpS[20] = { 1809 G4double NpS[20] = {
1814 80, 10, 80, 21, 90, 60, 50, 11, 60, 9, 23 1810 80, 10, 80, 21, 90, 60, 50, 11, 60, 9, 23, 50, 21, 21, 23, 16, 80, 22, 30, 32};
1815 1811
1816 G4double NpW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 1812 G4double NpW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0};
1817 1813
1818 AddElement("Np", 93, 20, *NpN, *NpA, *NpS, 1814 AddElement("Np", 93, 20, *NpN, *NpA, *NpS, *NpW);
1819 1815
1820 // Z = 94 --------------------------------- 1816 // Z = 94 --------------------------------------------------------------------
1821 G4int PuN[20] = {228, 229, 230, 231, 232, 2 1817 G4int PuN[20] = {228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243,
1822 244, 245, 246, 247}; 1818 244, 245, 246, 247};
1823 1819
1824 G4double PuA[20] = {228.039, 229.04, 230.04 1820 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.05 1821 237.048, 238.05, 239.052, 240.054, 241.057, 242.059, 243.062, 244.064, 245.068, 246.07,
1826 247.074}; 1822 247.074};
1827 1823
1828 G4double PuS[20] = {30, 80, 26, 11, 20, 50, 1824 G4double PuS[20] = {30, 80, 26, 11, 20, 50, 8, 22, 29, 25, 21, 21, 21, 21, 21, 3, 5, 15, 16, 32};
1829 1825
1830 G4double PuW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 1826 G4double PuW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
1831 1827
1832 AddElement("Pu", 94, 20, *PuN, *PuA, *PuS, 1828 AddElement("Pu", 94, 20, *PuN, *PuA, *PuS, *PuW);
1833 1829
1834 // Z = 95 --------------------------------- 1830 // Z = 95 --------------------------------------------------------------------
1835 G4int AmN[19] = { 1831 G4int AmN[19] = {
1836 231, 232, 233, 234, 235, 236, 237, 238, 2 1832 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249};
1837 1833
1838 G4double AmA[19] = {231.046, 232.047, 233.0 1834 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.06 1835 240.055, 241.057, 242.06, 243.061, 244.064, 245.066, 246.07, 247.072, 248.076, 249.078};
1840 1836
1841 G4double AmS[19] = {32, 32, 23, 22, 22, 11, 1837 G4double AmS[19] = {32, 32, 23, 22, 22, 11, 60, 50, 3, 15, 21, 21, 23, 23, 4, 20, 11, 22, 32};
1842 1838
1843 G4double AmW[19] = {0, 0, 0, 0, 0, 0, 0, 0, 1839 G4double AmW[19] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0};
1844 1840
1845 AddElement("Am", 95, 19, *AmN, *AmA, *AmS, 1841 AddElement("Am", 95, 19, *AmN, *AmA, *AmS, *AmW);
1846 1842
1847 // Z = 96 --------------------------------- 1843 // Z = 96 --------------------------------------------------------------------
1848 G4int CmN[20] = {233, 234, 235, 236, 237, 2 1844 G4int CmN[20] = {233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248,
1849 249, 250, 251, 252}; 1845 249, 250, 251, 252};
1850 1846
1851 G4double CmA[20] = {233.051, 234.05, 235.05 1847 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.0 1848 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}; 1849 251.082, 252.085};
1854 1850
1855 G4double CmS[20] = {43, 32, 24, 22, 22, 40, 1851 G4double CmS[20] = {43, 32, 24, 22, 22, 40, 11, 29, 24, 21, 24, 21, 29, 24, 5, 5, 5, 12, 24, 32};
1856 1852
1857 G4double CmW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 1853 G4double CmW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
1858 1854
1859 AddElement("Cm", 96, 20, *CmN, *CmA, *CmS, 1855 AddElement("Cm", 96, 20, *CmN, *CmA, *CmS, *CmW);
1860 1856
1861 // Z = 97 --------------------------------- 1857 // Z = 97 --------------------------------------------------------------------
1862 G4int BkN[20] = {235, 236, 237, 238, 239, 2 1858 G4int BkN[20] = {235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250,
1863 251, 252, 253, 254}; 1859 251, 252, 253, 254};
1864 1860
1865 G4double BkA[20] = {235.057, 236.057, 237.0 1861 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.07 1862 244.065, 245.066, 246.069, 247.07, 248.073, 249.075, 250.078, 251.081, 252.084, 253.087,
1867 254.091}; 1863 254.091};
1868 1864
1869 G4double BkS[20] = {43, 43, 32, 31, 31, 16, 1865 G4double BkS[20] = {43, 43, 32, 31, 31, 16, 22, 22, 5, 16, 26, 60, 6, 80, 3, 4, 12, 22, 39, 32};
1870 1866
1871 G4double BkW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 1867 G4double BkW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0, 0, 0};
1872 1868
1873 AddElement("Bk", 97, 20, *BkN, *BkA, *BkS, 1869 AddElement("Bk", 97, 20, *BkN, *BkA, *BkS, *BkW);
1874 1870
1875 // Z = 98 --------------------------------- 1871 // Z = 98 --------------------------------------------------------------------
1876 G4int CfN[20] = {237, 238, 239, 240, 241, 2 1872 G4int CfN[20] = {237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252,
1877 253, 254, 255, 256}; 1873 253, 254, 255, 256};
1878 1874
1879 G4double CfA[20] = {237.062, 238.061, 239.0 1875 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.0 1876 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}; 1877 255.091, 256.093};
1882 1878
1883 G4double CfS[20] = {54, 43, 25, 22, 27, 40, 1879 G4double CfS[20] = {54, 43, 25, 22, 27, 40, 15, 4, 11, 24, 9, 6, 3, 24, 5, 5, 7, 13, 22, 32};
1884 1880
1885 G4double CfW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 1881 G4double CfW[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
1886 1882
1887 AddElement("Cf", 98, 20, *CfN, *CfA, *CfS, 1883 AddElement("Cf", 98, 20, *CfN, *CfA, *CfS, *CfW);
1888 1884
1889 // Z = 99 --------------------------------- 1885 // Z = 99 --------------------------------------------------------------------
1890 G4int EsN[18] = { 1886 G4int EsN[18] = {
1891 240, 241, 242, 243, 244, 245, 246, 247, 2 1887 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257};
1892 1888
1893 G4double EsA[18] = {240.069, 241.069, 242.0 1889 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.08 1890 249.076, 250.079, 251.08, 252.083, 253.085, 254.088, 255.09, 256.094, 257.096};
1895 1891
1896 G4double EsS[18] = {43, 32, 35, 31, 20, 22, 1892 G4double EsS[18] = {43, 32, 35, 31, 20, 22, 24, 30, 60, 30, 11, 7, 50, 3, 5, 12, 11, 44};
1897 1893
1898 G4double EsW[18] = {0, 0, 0, 0, 0, 0, 0, 0, 1894 G4double EsW[18] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0, 0, 0};
1899 1895
1900 AddElement("Es", 99, 18, *EsN, *EsA, *EsS, 1896 AddElement("Es", 99, 18, *EsN, *EsA, *EsS, *EsW);
1901 1897
1902 // Z = 100 -------------------------------- 1898 // Z = 100 -------------------------------------------------------------------
1903 G4int FmN[18] = { 1899 G4int FmN[18] = {
1904 242, 243, 244, 245, 246, 247, 248, 249, 2 1900 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259};
1905 1901
1906 G4double FmA[18] = {242.073, 243.075, 244.0 1902 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.08 1903 250.08, 251.082, 252.082, 253.085, 254.087, 255.09, 256.092, 257.095, 258.097, 259.101};
1908 1904
1909 G4double FmS[18] = {43, 25, 31, 30, 40, 16, 1905 G4double FmS[18] = {43, 25, 31, 30, 40, 16, 13, 15, 13, 9, 6, 5, 3, 5, 8, 7, 22, 30};
1910 1906
1911 G4double FmW[18] = {0, 0, 0, 0, 0, 0, 0, 0, 1907 G4double FmW[18] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0};
1912 1908
1913 AddElement("Fm", 100, 18, *FmN, *FmA, *FmS, 1909 AddElement("Fm", 100, 18, *FmN, *FmA, *FmS, *FmW);
1914 1910
1915 // Z = 101 -------------------------------- 1911 // Z = 101 -------------------------------------------------------------------
1916 G4int MdN[16] = {245, 246, 247, 248, 249, 2 1912 G4int MdN[16] = {245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260};
1917 1913
1918 G4double MdA[16] = {245.081, 246.082, 247.0 1914 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.09 1915 253.087, 254.09, 255.091, 256.094, 257.096, 258.098, 259.101, 260.104};
1920 1916
1921 G4double MdS[16] = {40, 41, 39, 26, 24, 32, 1917 G4double MdS[16] = {40, 41, 39, 26, 24, 32, 22, 21, 22, 11, 8, 60, 4, 5, 22, 34};
1922 1918
1923 G4double MdW[16] = {0, 0, 0, 0, 0, 0, 0, 0, 1919 G4double MdW[16] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0};
1924 1920
1925 AddElement("Md", 101, 16, *MdN, *MdA, *MdS, 1921 AddElement("Md", 101, 16, *MdN, *MdA, *MdS, *MdW);
1926 1922
1927 // Z = 102 -------------------------------- 1923 // Z = 102 -------------------------------------------------------------------
1928 G4int NoN[14] = {249, 250, 251, 252, 253, 2 1924 G4int NoN[14] = {249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262};
1929 1925
1930 G4double NoA[14] = {249.088, 250.087, 251.0 1926 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.1 1927 257.097, 258.098, 259.101, 260.103, 261.106, 262.108};
1932 1928
1933 G4double NoS[14] = {37, 22, 19, 14, 26, 19, 1929 G4double NoS[14] = {37, 22, 19, 14, 26, 19, 13, 9, 30, 22, 11, 22, 32, 58};
1934 1930
1935 G4double NoW[14] = {0, 0, 0, 0, 0, 0, 0, 0, 1931 G4double NoW[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
1936 1932
1937 AddElement("No", 102, 14, *NoN, *NoA, *NoS, 1933 AddElement("No", 102, 14, *NoN, *NoA, *NoS, *NoW);
1938 1934
1939 // Z = 103 -------------------------------- 1935 // Z = 103 -------------------------------------------------------------------
1940 G4int LrN[13] = {251, 252, 253, 254, 255, 2 1936 G4int LrN[13] = {251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263};
1941 1937
1942 G4double LrA[13] = {251.094, 252.095, 253.0 1938 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}; 1939 260.106, 261.107, 262.11, 263.111};
1944 1940
1945 G4double LrS[13] = {32, 32, 24, 36, 22, 24, 1941 G4double LrS[13] = {32, 32, 24, 36, 22, 24, 22, 11, 80, 12, 22, 32, 39};
1946 1942
1947 G4double LrW[13] = {0, 0, 0, 0, 0, 0, 0, 0, 1943 G4double LrW[13] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100, 0};
1948 1944
1949 AddElement("Lr", 103, 13, *LrN, *LrA, *LrS, 1945 AddElement("Lr", 103, 13, *LrN, *LrA, *LrS, *LrW);
1950 1946
1951 // Z = 104 -------------------------------- 1947 // Z = 104 -------------------------------------------------------------------
1952 G4int RfN[12] = {253, 254, 255, 256, 257, 2 1948 G4int RfN[12] = {253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264};
1953 1949
1954 G4double RfA[12] = {253.101, 254.1, 255.101 1950 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}; 1951 262.11, 263.113, 264.114};
1956 1952
1957 G4double RfS[12] = {49, 31, 22, 29, 29, 22, 1953 G4double RfS[12] = {49, 31, 22, 29, 29, 22, 80, 22, 11, 30, 20, 48};
1958 1954
1959 G4double RfW[12] = {0, 0, 0, 0, 0, 0, 0, 0, 1955 G4double RfW[12] = {0, 0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
1960 1956
1961 AddElement("Rf", 104, 12, *RfN, *RfA, *RfS, 1957 AddElement("Rf", 104, 12, *RfN, *RfA, *RfS, *RfW);
1962 1958
1963 // Z = 105 -------------------------------- 1959 // Z = 105 -------------------------------------------------------------------
1964 G4int DbN[11] = {255, 256, 257, 258, 259, 2 1960 G4int DbN[11] = {255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265};
1965 1961
1966 G4double DbA[11] = {255.107, 256.108, 257.1 1962 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}; 1963 263.115, 264.117, 265.119};
1968 1964
1969 G4double DbS[11] = {45, 39, 25, 37, 31, 25, 1965 G4double DbS[11] = {45, 39, 25, 37, 31, 25, 25, 20, 18, 25, 30};
1970 1966
1971 G4double DbW[11] = {0, 0, 0, 0, 0, 0, 0, 10 1967 G4double DbW[11] = {0, 0, 0, 0, 0, 0, 0, 100, 0, 0, 0};
1972 1968
1973 AddElement("Db", 105, 11, *DbN, *DbA, *DbS, 1969 AddElement("Db", 105, 11, *DbN, *DbA, *DbS, *DbW);
1974 1970
1975 // Z = 106 -------------------------------- 1971 // Z = 106 -------------------------------------------------------------------
1976 G4int SgN[9] = {258, 259, 260, 261, 262, 26 1972 G4int SgN[9] = {258, 259, 260, 261, 262, 263, 264, 265, 266};
1977 1973
1978 G4double SgA[9] = { 1974 G4double SgA[9] = {
1979 258.113, 259.115, 260.114, 261.116, 262.1 1975 258.113, 259.115, 260.114, 261.116, 262.116, 263.118, 264.119, 265.121, 266.122};
1980 1976
1981 G4double SgS[9] = {45, 23, 40, 30, 30, 13, 1977 G4double SgS[9] = {45, 23, 40, 30, 30, 13, 30, 15, 31};
1982 1978
1983 G4double SgW[9] = {0, 0, 0, 0, 0, 0, 0, 0, 1979 G4double SgW[9] = {0, 0, 0, 0, 0, 0, 0, 0, 100};
1984 1980
1985 AddElement("Sg", 106, 9, *SgN, *SgA, *SgS, 1981 AddElement("Sg", 106, 9, *SgN, *SgA, *SgS, *SgW);
1986 1982
1987 // Z = 107 -------------------------------- 1983 // Z = 107 -------------------------------------------------------------------
1988 G4int BhN[8] = {260, 261, 262, 263, 264, 26 1984 G4int BhN[8] = {260, 261, 262, 263, 264, 265, 266, 267};
1989 1985
1990 G4double BhA[8] = {260.122, 261.122, 262.12 1986 G4double BhA[8] = {260.122, 261.122, 262.123, 263.123, 264.125, 265.125, 266.127, 267.128};
1991 1987
1992 G4double BhS[8] = {66, 26, 40, 45, 30, 41, 1988 G4double BhS[8] = {66, 26, 40, 45, 30, 41, 38, 37};
1993 1989
1994 G4double BhW[8] = {0, 0, 0, 0, 100, 0, 0, 0 1990 G4double BhW[8] = {0, 0, 0, 0, 100, 0, 0, 0};
1995 1991
1996 AddElement("Bh", 107, 8, *BhN, *BhA, *BhS, 1992 AddElement("Bh", 107, 8, *BhN, *BhA, *BhS, *BhW);
>> 1993
>> 1994 for (auto& i : elmSymbol) {
>> 1995 elmNames.push_back(i);
>> 1996 }
1997 1997
1998 if (0 < verbose) { 1998 if (0 < verbose) {
1999 G4cout << "G4NistElementBuilder: " << max 1999 G4cout << "G4NistElementBuilder: " << maxNumElements - 1 << " Elements " << index
2000 << " Isotopes" << G4endl; 2000 << " Isotopes" << G4endl;
2001 } 2001 }
2002 } 2002 }
2003 2003