Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
G4QNeutronCaptureRatio.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// $Id$
28//
29//
30// G4 Physics class: G4QNeutronCaptureRatio for N+A Diffraction Interactions
31// Created: M.V. Kossov, CERN/ITEP(Moscow), 10-OCT-01
32// The last update: M.V. Kossov, CERN/ITEP (Moscow) 15-Oct-06
33//
34// ----------------------------------------------------------------------
35// Short description: (n,gamma) capture is a part of the incoherent
36// (inelastic) interaction. This part is calculated in the class.
37// ----------------------------------------------------------------------
38
39//#define debug
40//#define pdebug
41//#define fdebug
42//#define nandebug
43
44#include "G4QNeutronCaptureRatio.hh"
45#include "G4SystemOfUnits.hh"
46
47// Returns Pointer to the class G4QNeutronCaptureRatio
48G4QNeutronCaptureRatio* G4QNeutronCaptureRatio::GetPointer()
49{
50 static G4QNeutronCaptureRatio theRatios; // *** Static body of the NeutCaptureRatio ***
51 return &theRatios;
52}
53
54// Calculation of netronCapture/inelastic ratio
55G4double G4QNeutronCaptureRatio::GetRatio(G4double pIU, G4int tgZ, G4int tgN)
56{
57 // Table parameters
58 static const G4int npp=100; // Number of steps in the R(s) LinTable
59 static const G4int mpp=npp+1; // Number of elements in the R(s) LinTable
60 static const G4double pma=6.; // The first LinTabEl(mom=0)=1., mom>pma -> logTab
61 static const G4double dp=pma/npp; // Step of the linear Table
62 static const G4int nls=150; // Number of steps in the R(lns) logTable
63 static const G4int mls=nls+1; // Number of elements in the R(lns) logTable
64 static const G4double lpi=1.79; // The min ln(p) logTabEl(p=5.99 < pma=6.)
65 static const G4double lpa=8.; // The max ln(p) logTabEl(p=5.99 - 2981 GeV)
66 static const G4double mi=std::exp(lpi);// The min mom of logTabEl(~ 5.99 GeV)
67 static const G4double max_s=std::exp(lpa);// The max mom of logTabEl(~ 2981 GeV)
68 static const G4double dl=(lpa-lpi)/nls;// Step of the logarithmic Table
69 static const G4double edl=std::exp(dl);// Multiplication step of the logarithmic Table
70 static const G4double toler=.0001; // Tolarence (GeV) defining the same momentum
71 static G4double lastP=0.; // Last mometum value for which R was calculated
72 static G4double lastR=0.; // Last ratio R which was calculated
73 // Local Associative Data Base:
74 static std::vector<G4int> vZ; // Vector of calculated Z (target)
75 static std::vector<G4int> vN; // Vector of calculated N (target)
76 static std::vector<G4double> vH; // Vector of max mom initialized in the LinTable
77 static std::vector<G4int> vJ; // Vector of topBin number initialized in LinTable
78 static std::vector<G4double> vM; // Vector of relMax ln(p) initialized in LogTable
79 static std::vector<G4int> vK; // Vector of topBin number initialized in LogTable
80 static std::vector<G4double*> vT; // Vector of pointers to LinTable in C++ heap
81 static std::vector<G4double*> vL; // Vector of pointers to LogTable in C++ heap
82 // Last values of the Associative Data Base:
83 static G4int lastZ=0; // theLast of calculated A
84 static G4int lastN=0; // theLast of calculated A
85 static G4double lastH=0.; // theLast of max mom initialized in the LinTable
86 static G4int lastJ=0; // theLast of topBinNumber initialized in LinTable
87 static G4double lastM=0.; // theLast of relMax ln(p) initialized in LogTab.
88 static G4int lastK=0; // theLast of topBinNumber initialized in LogTable
89 static G4double* lastT=0; // theLast of pointer to LinTable in the C++ heap
90 static G4double* lastL=0; // theLast of pointer to LogTable in the C++ heap
91 // LogTable is created only if necessary. R(p>2981GeV) calcul by formula for any nuclei
92 G4int A=tgN+tgZ;
93 if(pIU > 50) return 0.;
94 if(pIU > 30 && ((tgN==1 && tgZ==1) || (tgN==8 && tgZ==7))) return 0.;
95 if(pIU > 20 && tgN==2 && tgZ==1) return 0.;
96 if(pIU > 15 && ((tgN==1 && tgZ==2) || (tgN==8 && tgZ==8))) return 0.;
97 if(pIU<toler || A<1) return 1.; // Fake use of toler as non zero number
98 if(A>247)
99 {
100 G4cout<<"-*-Warning-*-G4NeutronCaptureRatio::GetRatio:A="<<A<<">247, return 0"<<G4endl;
101 return 0.;
102 }
103 G4int nDB=vZ.size(); // A number of nuclei already initialized in AMDB
104 if(nDB && lastZ==tgZ && lastN==tgN && std::fabs(pIU-lastP)<toler) return lastR;
105 if(pIU>max_s)
106 {
107 lastR=CalcCap2In_Ratio(s,tgZ,tgN); //@@ Probably user ought to be notified about bigP
108 return lastR;
109 }
110 G4bool found=false;
111 G4int i=-1;
112 if(nDB) for (i=0; i<nDB; i++) if(tgZ==vZ[i] && tgN==vN[i]) // Sirch for this Z,N in AMDB
113 {
114 found=true; // The (Z,N) is found
115 break;
116 }
117 if(!nDB || !found) // Create new line in the AMDB
118 {
119 lastZ = tgZ;
120 lastN = tgN;
121 lastT = new G4double[mpp]; // Create the linear Table
122 lastJ = static_cast<int>(pIU/dp)+1; // MaxBin to be initialized
123 if(lastJ>npp)
124 {
125 lastJ=npp;
126 lastH=pma;
127 }
128 else lastH = lastJ*dp; // Calculate max initialized s for LinTab
129 G4double pv=0;
130 lastT[0]=1.;
131 for(G4int j=1; j<=lastJ; j++) // Calculate LinTab values
132 {
133 pv+=dp;
134 lastT[j]=CalcCap2In_Ratio(pv,tgZ,tgN); // ??
135 }
136 lastL=new G4double[mls]; // Create the logarithmic Table
137 G4double ls=std::log(s);
138 lastK = static_cast<int>((ls-lpi)/dl)+1; // MaxBin to be initialized in LogTaB
139 if(lastK>nls)
140 {
141 lastK=nls;
142 lastM=lpa-lpi;
143 }
144 else lastM = lastK*dl; // Calculate max initialized ln(s)-lpi for LogTab
145 pv=mi;
146 for(G4int j=0; j<=lastK; j++) // Calculate LogTab values
147 {
148 lastL[j]=CalcCap2In_Ratio(pv,tgZ,tgN);
149 if(j!=lastK) pv*=edl;
150 }
151 i++; // Make a new record to AMDB and position on it
152 vZ.push_back(lastZ);
153 vN.push_back(lastN);
154 vH.push_back(lastH);
155 vJ.push_back(lastJ);
156 vM.push_back(lastM);
157 vK.push_back(lastK);
158 vT.push_back(lastT);
159 vL.push_back(lastL);
160 }
161 else // The A value was found in AMDB
162 {
163 lastZ=vZ[i];
164 lastN=vN[i];
165 lastH=vH[i];
166 lastJ=vJ[i];
167 lastM=vM[i];
168 lastK=vK[i];
169 lastT=vT[i];
170 lastL=vL[i];
171 if(s>lastH) // At least LinTab must be updated
172 {
173 G4int nextN=lastJ+1; // The next bin to be initialized
174 if(lastJ<npp)
175 {
176 lastJ = static_cast<int>(pIU/dp)+1;// MaxBin to be initialized
177 G4double pv=lastH;
178 if(lastJ>npp)
179 {
180 lastJ=npp;
181 lastH=pma;
182 }
183 else lastH = lastJ*dp; // Calculate max initialized s for LinTab
184 for(G4int j=nextN; j<=lastJ; j++)// Calculate LogTab values
185 {
186 pv+=dp;
187 lastT[j]=CalcCap2In_Ratio(pv,tgZ,tgN);
188 }
189 } // End of LinTab update
190 if(lastJ>=nextN)
191 {
192 vH[i]=lastH;
193 vJ[i]=lastJ;
194 }
195 G4int nextK=lastK+1;
196 if(pIU>pma && lastK<nls) // LogTab must be updated
197 {
198 G4double pv=std::exp(lastM+lpi); // Define starting poit (lastM will be changed)
199 G4double ls=std::log(s);
200 lastK = static_cast<int>((ls-lpi)/dl)+1; // MaxBin to be initialized in LogTaB
201 if(lastK>nls)
202 {
203 lastK=nls;
204 lastM=lpa-lpi;
205 }
206 else lastM = lastK*dl; // Calculate max initialized ln(p)-lpi for LogTab
207 for(G4int j=nextK; j<=lastK; j++)// Calculate LogTab values
208 {
209 pv*=edl;
210 lastL[j]=CalcCap2In_Ratio(pv,tgZ,tgN);
211 }
212 } // End of LogTab update
213 if(lastK>=nextK)
214 {
215 vM[i]=lastM;
216 vK[i]=lastK;
217 }
218 }
219 }
220 // Now one can use tabeles to calculate the value
221 if(pIU<pma) // Use linear table
222 {
223 G4int n=static_cast<int>(pIU/dp); // Low edge number of the bin
224 G4double d=s-n*dp; // Linear shift
225 G4double v=lastT[n]; // Base
226 lastR=v+d*(lastT[n+1]-v)/dp; // Result
227 }
228 else // Use log table
229 {
230 G4double ls=std::log(pIU)-lpi; // ln(p)-l_min
231 G4int n=static_cast<int>(ls/dl); // Low edge number of the bin
232 G4double d=ls-n*dl; // Log shift
233 G4double v=lastL[n]; // Base
234 lastR=v+d*(lastL[n+1]-v)/dl; // Result
235 }
236 if(lastR<0.) lastR=0.;
237 if(lastR>1.) lastR=1.;
238 return lastR;
239} // End of GetRatio
240
241// Calculate Capture/Inelastic Ratio as a function of total momentum (in GeV/c)
242G4double G4QNeutronCaptureRatio::CalcCap2In_Ratio(G4double p, G4int Z, G4int N)
243{
244 //==> n (Z=0)
245 static const G4int N0=1;
246 static const G4double pZ0N1[5]={.0001, 40., 0., 0., 1.}; // *** No Capture ?
247 static const std::pair<G4int, const G4double*> Z0N1(1,pZ0N1);
248 static const std::pair<G4int, const G4double*> Z0[N0]={Z0N1};
249 //==> H (Z=1) *** no protons, which are treated separately ***
250 static const G4int N1=2;
251 static const G4double pZ1N1[5]={.07, 40., .0006, .0001, .01};
252 static const std::pair<G4int, const G4double*> Z1N1(1,pZ1N1);
253 static const G4double pZ1N2[5]={.0001, 40., 0., 0., 1.}; // *** No Capture ?
254 static const std::pair<G4int, const G4double*> Z1N2(2,pZ1N2);
255 static const std::pair<G4int, const G4double*> Z1[N1]={Z1N1, Z1N2};
256 //==> He(Z=2)
257 static const G4int N2=2;
258 static const G4double pZ2N1[5]={.1, 40., 0., 0., 1.}; // *** Unknown threshold
259 static const std::pair<G4int, const G4double*> Z2N1(1,pZ2N1);
260 static const G4double pZ2N2[5]={.0001, 40., 0., 0., 1.}; // *** No Capture ?
261 static const std::pair<G4int, const G4double*> Z2N2(2,pZ2N2);
262 static const std::pair<G4int, const G4double*> Z2[N2]={Z2N1, Z2N2};
263 //==> Li(Z=3)
264 static const G4int N3=2;
265 static const G4double pZ3N3[5]={.001, 40., 3.E-5, .0001, .1};
266 static const std::pair<G4int, const G4double*> Z3N1(3,pZ3N3);
267 static const G4double pZ3N4[5]={.022, 19., 3.E-5, .0001, .04};
268 static const std::pair<G4int, const G4double*> Z3N2(4,pZ3N4);
269 static const std::pair<G4int, const G4double*> Z3[N3]={Z3N1, Z3N2};
270 //==> Be(Z=4)
271 static const G4int N4=1;
272 static const G4double pZ4N5[5]={.0004, 40., 0., 0., 1.};
273 static const std::pair<G4int, const G4double*> Z4N5(5,pZ4N5);
274 static const std::pair<G4int, const G4double*> Z4[N4]={Z4N5};
275 //==> B (Z=5)
276 static const G4int N5=2;
277 static const G4double pZ5N5[5]={.011, 9., .0002, .0001, .002};
278 static const std::pair<G4int, const G4double*> Z5N5(5,pZ5N5);
279 static const G4double pZ5N6[5]={.027, 9., 0., 0., 1.};
280 static const std::pair<G4int, const G4double*> Z5N6(6,pZ5N6);
281 static const std::pair<G4int, const G4double*> Z5[N5]={Z5N5, Z5N6};
282 //==> C (Z=6)
283 static const G4int N6=2;
284 static const G4double pZ6N6[5]={.08, 40., .0003, .0001, .07}; // *** Only Nat Mix ***
285 static const std::pair<G4int, const G4double*> Z6N6(6,pZ6N6);
286 static const G4double pZ6N7[5]={.08, 40., .0003, .0001, .07}; // *** Only Nat Mix ***
287 static const std::pair<G4int, const G4double*> Z6N7(7,pZ6N7);
288 static const std::pair<G4int, const G4double*> Z6[N6]={Z6N6, Z6N7};
289 //==> N (Z=7)
290 static const G4int N7=2;
291 static const G4double pZ7N7[5]={.005, 3., 0., 0., 1.};
292 static const std::pair<G4int, const G4double*> Z7N7(7,pZ7N7);
293 static const G4double pZ7N8[5]={.084, 40., .0001, .0001, .015};
294 static const std::pair<G4int, const G4double*> Z7N8(8,pZ7N8);
295 static const std::pair<G4int, const G4double*> Z7[N7]={Z7N7, Z7N8};
296 //==> O (Z=8)
297 static const G4int N8=3;
298 static const G4double pZ8N8[5]={.08, 40., 0., 0., 1.};
299 static const std::pair<G4int, const G4double*> Z8N8(8,pZ8N8);
300 static const G4double pZ8N9[5]={.0065, 5., .0013, .0001, .02};
301 static const std::pair<G4int, const G4double*> Z8N9(9,pZ8N9);
302 static const G4double pZ8N10[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
303 static const std::pair<G4int, const G4double*> Z8N10(10,pZ8N10);
304 static const std::pair<G4int, const G4double*> Z8[N8]={Z8N8, Z8N9, Z8N10};
305 //==> F (Z=9)
306 static const G4int N9=1;
307 static const G4double pZ9N10[5]={.013, 27., .0001, .0001, .02};
308 static const std::pair<G4int, const G4double*> Z9N10(10,pZ9N10);
309 static const std::pair<G4int, const G4double*> Z9[N9]={Z9N10};
310 //==> Ne(Z=10)
311 static const G4int N10=3;
312 static const G4double pZ10N10[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
313 static const std::pair<G4int, const G4double*> Z10N10(10,pZ10N10);
314 static const G4double pZ10N11[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
315 static const std::pair<G4int, const G4double*> Z10N11(11,pZ10N11);
316 static const G4double pZ10N12[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
317 static const std::pair<G4int, const G4double*> Z10N12(12,pZ10N12);
318 static const std::pair<G4int, const G4double*> Z10[N10]={Z10N10, Z10N11, Z10N12};
319 //==> Na(Z=11)
320 static const G4int N11=1;
321 static const G4double pZ11N12[5]={.024, 17., .0005, .0001, .03};
322 static const std::pair<G4int, const G4double*> Z11N12(12,pZ11N12);
323 static const std::pair<G4int, const G4double*> Z11[N11]={Z11N12};
324 //==> Mg(Z=12)
325 static const G4int N12=3;
326 static const G4double pZ12N12[5]={.045, 40., .0003, .0001, .02};
327 static const std::pair<G4int, const G4double*> Z12N12(12,pZ12N12);
328 static const G4double pZ12N13[5]={.019, 7., .0002, .0001, .01};
329 static const std::pair<G4int, const G4double*> Z12N13(13,pZ12N13);
330 static const G4double pZ12N14[5]={.053, 40., .0006, .0001, .007};
331 static const std::pair<G4int, const G4double*> Z12N14(14,pZ12N14);
332 static const std::pair<G4int, const G4double*> Z12[N12]={Z12N12, Z12N13, Z12N14};
333 //==> Al(Z=13)
334 static const G4int N13=1;
335 static const G4double pZ13N14[5]={.035, 17., .001, .03, 1.};
336 static const std::pair<G4int, const G4double*> Z13N14(14,pZ13N14);
337 static const std::pair<G4int, const G4double*> Z13[N13]={Z13N14};
338 //==> Si(Z=14)
339 static const G4int N14=3;
340 static const G4double pZ14N14[5]={.052, 40., .002, .0001, .008};
341 static const std::pair<G4int, const G4double*> Z14N14(14,pZ14N14);
342 static const G4double pZ14N15[5]={.048, 40., .0001, .0004, .02};
343 static const std::pair<G4int, const G4double*> Z14N15(15,pZ14N15);
344 static const G4double pZ14N16[5]={.06, 40., .0015, .0001, .01};
345 static const std::pair<G4int, const G4double*> Z14N16(16,pZ14N16);
346 static const std::pair<G4int, const G4double*> Z14[N14]={Z14N14, Z14N15, Z14N16};
347 //==> P (Z=15)
348 static const G4int N15=1;
349 static const G4double pZ15N16[5]={.024, 7., .0008, .0001, .03};
350 static const std::pair<G4int, const G4double*> Z15N16(16,pZ15N16);
351 static const std::pair<G4int, const G4double*> Z15[N15]={Z15N16};
352 //==> S (Z=16)
353 static const G4int N16=4;
354 static const G4double pZ16N16[5]={.036, 12., .0003, .03, .004};
355 static const std::pair<G4int, const G4double*> Z16N16(16,pZ16N16);
356 static const G4double pZ16N17[5]={.018, 40., .0033, .0001, .002};
357 static const std::pair<G4int, const G4double*> Z16N17(17,pZ16N17);
358 static const G4double pZ16N18[5]={.053, 25., .002, .0001, .0043};
359 static const std::pair<G4int, const G4double*> Z16N18(18,pZ16N18);
360 static const G4double pZ16N20[5]={.065, 25., .002, .0001, .0043};
361 static const std::pair<G4int, const G4double*> Z16N20(20,pZ16N20);
362 static const std::pair<G4int, const G4double*> Z16[N16]={Z16N16, Z16N17, Z16N18, Z16N20};
363 //==> Cl(Z=17)
364 static const G4int N17=2;
365 static const G4double pZ17N18[5]={.014, 4., .0004, .175, .002};
366 static const std::pair<G4int, const G4double*> Z17N18(18,pZ17N18);
367 static const G4double pZ17N20[5]={.035, 8., .008, 18., .0005};
368 static const std::pair<G4int, const G4double*> Z17N20(20,pZ17N20);
369 static const std::pair<G4int, const G4double*> Z17[N17]={Z17N18, Z17N20};
370 //==> Ar(Z=18)
371 static const G4int N18=3;
372 static const G4double pZ18N18[5]={.036, 8., .0005, .1, .01};
373 static const std::pair<G4int, const G4double*> Z18N18(18,pZ18N18);
374 static const G4double pZ18N20[5]={.025, 6., .0027, .19, .0003};
375 static const std::pair<G4int, const G4double*> Z18N20(20,pZ18N20);
376 static const G4double pZ18N22[5]={.028, 6., .001, .19, .0003};
377 static const std::pair<G4int, const G4double*> Z18N22(22,pZ18N22);
378 static const std::pair<G4int, const G4double*> Z18[N18]={Z18N18, Z18N20, Z18N22};
379 //==> K (Z=19)
380 static const G4int N19=3;
381 static const G4double pZ19N20[5]={.04, 8., 0., 0., 1.};
382 static const std::pair<G4int, const G4double*> Z19N20(20,pZ19N20);
383 static const G4double pZ19N21[5]={.049, 5., 0., 0., 1.};
384 static const std::pair<G4int, const G4double*> Z19N21(21,pZ19N21);
385 static const G4double pZ19N22[5]={.04, 11., .005, .0001, .005};
386 static const std::pair<G4int, const G4double*> Z19N22(22,pZ19N22);
387 static const std::pair<G4int, const G4double*> Z19[N19]={Z19N20, Z19N21, Z19N22};
388 //==> Ca(Z=20)
389 static const G4int N20=6;
390 static const G4double pZ20N20[5]={.05, 14., .0006, .09, .009};
391 static const std::pair<G4int, const G4double*> Z20N20(20,pZ20N20);
392 static const G4double pZ20N22[5]={.047, 30., .003, .0001, .014};
393 static const std::pair<G4int, const G4double*> Z20N22(22,pZ20N22);
394 static const G4double pZ20N23[5]={.01, 3.5, .0015, .16, .002};
395 static const std::pair<G4int, const G4double*> Z20N23(23,pZ20N23);
396 static const G4double pZ20N24[5]={.04, 30., .002, .0001, .008};
397 static const std::pair<G4int, const G4double*> Z20N24(24,pZ20N24);
398 static const G4double pZ20N26[5]={.044, 40., .001, .0001, .01};
399 static const std::pair<G4int, const G4double*> Z20N26(26,pZ20N26);
400 static const G4double pZ20N28[5]={.055, 14., .001, .18, .001};
401 static const std::pair<G4int, const G4double*> Z20N28(28,pZ20N28);
402 static const std::pair<G4int, const G4double*> Z20[N20]={Z20N20, Z20N22, Z20N23,
403 Z20N24, Z20N26, Z20N28};
404 //==> Sc(Z=21)
405 static const G4int N21=1;
406 static const G4double pZ21N24[5]={.014, 4., 0., 0., 1.};
407 static const std::pair<G4int, const G4double*> Z21N24(24,pZ21N24);
408 static const std::pair<G4int, const G4double*> Z21[N21]={Z21N24};
409 //==> Ti(Z=22)
410 static const G4int N22=5;
411 static const G4double pZ22N24[5]={.036, 27., .007, .0001, .005};
412 static const std::pair<G4int, const G4double*> Z22N24(24,pZ22N24);
413 static const G4double pZ22N25[5]={.013, 9., .017, .0001, .005};
414 static const std::pair<G4int, const G4double*> Z22N25(25,pZ22N25);
415 static const G4double pZ22N26[5]={.043, 40., .002, .0001, .01};
416 static const std::pair<G4int, const G4double*> Z22N26(26,pZ22N26);
417 static const G4double pZ22N27[5]={.047, 30., .007, .0001, .01};
418 static const std::pair<G4int, const G4double*> Z22N27(27,pZ22N27);
419 static const G4double pZ22N28[5]={.052, 40., .0005, .0001, .01};
420 static const std::pair<G4int, const G4double*> Z22N28(28,pZ22N28);
421 static const std::pair<G4int, const G4double*> Z22[N22]={Z22N24, Z22N25, Z22N26,
422 Z22N27, Z22N28};
423 //==> V (Z=23)
424 static const G4int N23=2;
425 static const G4double pZ23N27[5]={.023, 30., .01, .0001, .003}; // *** Only Nat mix ***
426 static const std::pair<G4int, const G4double*> Z23N27(27,pZ23N27);
427 static const G4double pZ23N28[5]={.023, 30., .01, .0001, .003}; // *** Only Nat mix ***
428 static const std::pair<G4int, const G4double*> Z23N28(28,pZ23N28);
429 static const std::pair<G4int, const G4double*> Z23[N23]={Z23N27, Z23N28};
430 //==> Cr(Z=24)
431 static const G4int N24=4;
432 static const G4double pZ24N26[5]={.035, 27., .004, .0001, .01};
433 static const std::pair<G4int, const G4double*> Z24N26(26,pZ24N26);
434 static const G4double pZ24N28[5]={.049, 40., .001, .0001, .016};
435 static const std::pair<G4int, const G4double*> Z24N28(28,pZ24N28);
436 static const G4double pZ24N29[5]={.032, 30., .005, .0001, .005};
437 static const std::pair<G4int, const G4double*> Z24N29(29,pZ24N29);
438 static const G4double pZ24N30[5]={.034, 30., .002, .0001, .008};
439 static const std::pair<G4int, const G4double*> Z24N30(30,pZ24N30);
440 static const std::pair<G4int, const G4double*> Z24[N24]={Z24N26, Z24N28, Z24N29, Z24N30};
441 //==> Mn(Z=25)
442 static const G4int N25=1;
443 static const G4double pZ25N30[5]={.0145, 10., .01, .0001, .003};
444 static const std::pair<G4int, const G4double*> Z25N30(30,pZ25N30);
445 static const std::pair<G4int, const G4double*> Z25[N25]={Z25N30};
446 //==> Fe(Z=26)
447 static const G4int N26=4;
448 static const G4double pZ26N28[5]={.048, 27., .0016, .0001, .016};
449 static const std::pair<G4int, const G4double*> Z26N28(28,pZ26N28);
450 static const G4double pZ26N30[5]={.036, 27., .004, .0001, .006};
451 static const std::pair<G4int, const G4double*> Z26N30(30,pZ26N30);
452 static const G4double pZ26N31[5]={.036, 27., .005, .0001, .005};
453 static const std::pair<G4int, const G4double*> Z26N31(31,pZ26N31);
454 static const G4double pZ26N32[5]={.036, 27., .005, .0001, .007};
455 static const std::pair<G4int, const G4double*> Z26N32(32,pZ26N32);
456 static const std::pair<G4int, const G4double*> Z26[N26]={Z26N28, Z26N30, Z26N31, Z26N32};
457 //==> Co(Z=27)
458 static const G4int N27=1;
459 static const G4double pZ27N32[5]={.044, 22., .002, .0001, .016};
460 static const std::pair<G4int, const G4double*> Z27N32(32,pZ27N32);
461 static const std::pair<G4int, const G4double*> Z27[N27]={Z27N32};
462 //==> Ni(Z=28)
463 static const G4int N28=5;
464 static const G4double pZ28N30[5]={.045, 20., .003, .0001, .01};
465 static const std::pair<G4int, const G4double*> Z28N30(30,pZ28N30);
466 static const G4double pZ28N32[5]={.046, 20., .016, .0001, .01};
467 static const std::pair<G4int, const G4double*> Z28N32(32,pZ28N32);
468 static const G4double pZ28N33[5]={.046, 20., .016, .0001, .01};
469 static const std::pair<G4int, const G4double*> Z28N33(33,pZ28N33);
470 static const G4double pZ28N34[5]={.045, 20., .005, .0001, .007};
471 static const std::pair<G4int, const G4double*> Z28N34(34,pZ28N34);
472 static const G4double pZ28N36[5]={.045, 20., .005, .0001, .007};
473 static const std::pair<G4int, const G4double*> Z28N36(36,pZ28N36);
474 static const std::pair<G4int, const G4double*> Z28[N28]={Z28N30, Z28N32, Z28N33,
475 Z28N34, Z28N36};
476 //==> Cu(Z=29)
477 static const G4int N29=2;
478 static const G4double pZ29N34[5]={.035, 15., .008, .0001, .015};
479 static const std::pair<G4int, const G4double*> Z29N34(34,pZ29N34);
480 static const G4double pZ29N36[5]={.036, 15., .003, .0001, .013};
481 static const std::pair<G4int, const G4double*> Z29N36(36,pZ29N36);
482 static const std::pair<G4int, const G4double*> Z29[N29]={Z29N34, Z29N36};
483 //==> Zn(Z=30)
484 static const G4int N30=5;
485 static const G4double pZ30N34[5]={.041, 20., .008, .0001, .02}; // *** only NAT mix ***
486 static const std::pair<G4int, const G4double*> Z30N34(34,pZ30N34);
487 static const G4double pZ30N36[5]={.041, 20., .008, .0001, .02}; // *** only NAT mix ***
488 static const std::pair<G4int, const G4double*> Z30N36(36,pZ30N36);
489 static const G4double pZ30N37[5]={.041, 20., .008, .0001, .02}; // *** only NAT mix ***
490 static const std::pair<G4int, const G4double*> Z30N37(37,pZ30N37);
491 static const G4double pZ30N38[5]={.041, 20., .008, .0001, .02}; // *** only NAT mix ***
492 static const std::pair<G4int, const G4double*> Z30N38(38,pZ30N38);
493 static const G4double pZ30N40[5]={.041, 20., .008, .0001, .02}; // *** only NAT mix ***
494 static const std::pair<G4int, const G4double*> Z30N40(40,pZ30N40);
495 static const std::pair<G4int, const G4double*> Z30[N30]={Z30N34, Z30N36, Z30N37,
496 Z30N38, Z30N40};
497 //==> Ga(Z=31)
498 static const G4int N31=2;
499 static const G4double pZ31N38[5]={.024, 7., .03, .01, .003};
500 static const std::pair<G4int, const G4double*> Z31N38(38,pZ31N38);
501 static const G4double pZ31N40[5]={.026, 9., .015, .01, .003};
502 static const std::pair<G4int, const G4double*> Z31N40(40,pZ31N40);
503 static const std::pair<G4int, const G4double*> Z31[N31]={Z31N38, Z31N40};
504 //==> Ge(Z=32)
505 static const G4int N32=5;
506 static const G4double pZ32N38[5]={.037, 12., .15, .025, .003};
507 static const std::pair<G4int, const G4double*> Z32N38(38,pZ32N38);
508 static const G4double pZ32N40[5]={.035, 20., .015, .01, .0035};
509 static const std::pair<G4int, const G4double*> Z32N40(40,pZ32N40);
510 static const G4double pZ32N41[5]={.009, 3., .02, .03, .0001};
511 static const std::pair<G4int, const G4double*> Z32N41(41,pZ32N41);
512 static const G4double pZ32N42[5]={.027, 12., .003, .0001, .01};
513 static const std::pair<G4int, const G4double*> Z32N42(42,pZ32N42);
514 static const G4double pZ32N44[5]={.031, 20., .025, .0005, .0045};
515 static const std::pair<G4int, const G4double*> Z32N44(44,pZ32N44);
516 static const std::pair<G4int, const G4double*> Z32[N32]={Z32N38, Z32N40, Z32N41,
517 Z32N42, Z32N44};
518 //==> As(Z=33)
519 static const G4int N33=1;
520 static const G4double pZ33N42[5]={.017, 5., .004, .05, .006};
521 static const std::pair<G4int, const G4double*> Z33N42(42,pZ33N42);
522 static const std::pair<G4int, const G4double*> Z33[N33]={Z33N42};
523 //==> Se(Z=34)
524 static const G4int N34=6;
525 static const G4double pZ34N40[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
526 static const std::pair<G4int, const G4double*> Z34N40(40,pZ34N40);
527 static const G4double pZ34N42[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
528 static const std::pair<G4int, const G4double*> Z34N42(42,pZ34N42);
529 static const G4double pZ34N43[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
530 static const std::pair<G4int, const G4double*> Z34N43(43,pZ34N43);
531 static const G4double pZ34N44[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
532 static const std::pair<G4int, const G4double*> Z34N44(44,pZ34N44);
533 static const G4double pZ34N46[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
534 static const std::pair<G4int, const G4double*> Z34N46(46,pZ34N46);
535 static const G4double pZ34N48[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
536 static const std::pair<G4int, const G4double*> Z34N48(48,pZ34N48);
537 static const std::pair<G4int, const G4double*> Z34[N34]={Z34N40, Z34N42, Z34N43,
538 Z34N44, Z34N46, Z34N48};
539 //==> Br(Z=35)
540 static const G4int N35=2;
541 static const G4double pZ35N44[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
542 static const std::pair<G4int, const G4double*> Z35N44(44,pZ35N44);
543 static const G4double pZ35N46[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
544 static const std::pair<G4int, const G4double*> Z35N46(46,pZ35N46);
545 static const std::pair<G4int, const G4double*> Z35[N35]={Z35N44, Z35N46};
546 //==> Kr(Z=36)
547 static const G4int N36=6;
548 static const G4double pZ36N42[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
549 static const std::pair<G4int, const G4double*> Z36N42(42,pZ36N42);
550 static const G4double pZ36N44[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
551 static const std::pair<G4int, const G4double*> Z36N44(44,pZ36N44);
552 static const G4double pZ36N46[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
553 static const std::pair<G4int, const G4double*> Z36N46(46,pZ36N46);
554 static const G4double pZ36N47[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
555 static const std::pair<G4int, const G4double*> Z36N47(47,pZ36N47);
556 static const G4double pZ36N48[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
557 static const std::pair<G4int, const G4double*> Z36N48(48,pZ36N48);
558 static const G4double pZ36N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
559 static const std::pair<G4int, const G4double*> Z36N50(50,pZ36N50);
560 static const std::pair<G4int, const G4double*> Z36[N36]={Z36N42, Z36N44, Z36N46,
561 Z36N47, Z36N48, Z36N50};
562 //==> Rb(Z=37)
563 static const G4int N37=2;
564 static const G4double pZ37N48[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
565 static const std::pair<G4int, const G4double*> Z37N48(48,pZ37N48);
566 static const G4double pZ37N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
567 static const std::pair<G4int, const G4double*> Z37N50(50,pZ37N50);
568 static const std::pair<G4int, const G4double*> Z37[N37]={Z37N48, Z37N50};
569 //==> Sr(Z=38)
570 static const G4int N38=4;
571 static const G4double pZ38N46[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
572 static const std::pair<G4int, const G4double*> Z38N46(46,pZ38N46);
573 static const G4double pZ38N48[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
574 static const std::pair<G4int, const G4double*> Z38N48(48,pZ38N48);
575 static const G4double pZ38N49[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
576 static const std::pair<G4int, const G4double*> Z38N49(49,pZ38N49);
577 static const G4double pZ38N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
578 static const std::pair<G4int, const G4double*> Z38N50(50,pZ38N50);
579 static const std::pair<G4int, const G4double*> Z38[N38]={Z38N46, Z38N48, Z38N49, Z38N50};
580 //==> Y (Z=39)
581 static const G4int N39=1;
582 static const G4double pZ39N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
583 static const std::pair<G4int, const G4double*> Z39N50(50,pZ39N50);
584 static const std::pair<G4int, const G4double*> Z39[N39]={Z39N50};
585 //==> Zr(Z=40)
586 static const G4int N40=5;
587 static const G4double pZ40N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
588 static const std::pair<G4int, const G4double*> Z40N50(50,pZ40N50);
589 static const G4double pZ40N51[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
590 static const std::pair<G4int, const G4double*> Z40N51(51,pZ40N51);
591 static const G4double pZ40N52[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
592 static const std::pair<G4int, const G4double*> Z40N52(52,pZ40N52);
593 static const G4double pZ40N54[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
594 static const std::pair<G4int, const G4double*> Z40N54(54,pZ40N54);
595 static const G4double pZ40N56[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
596 static const std::pair<G4int, const G4double*> Z40N56(56,pZ40N56);
597 static const std::pair<G4int, const G4double*> Z40[N40]={Z40N50, Z40N51, Z40N52,
598 Z40N54, Z40N56};
599 //==> Nb(Z=41)
600 static const G4int N41=1;
601 static const G4double pZ41N52[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
602 static const std::pair<G4int, const G4double*> Z41N52(52,pZ41N52);
603 static const std::pair<G4int, const G4double*> Z41[N41]={Z41N52};
604 //==> Mo(Z=42)
605 static const G4int N42=7;
606 static const G4double pZ42N50[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
607 static const std::pair<G4int, const G4double*> Z42N50(50,pZ42N50);
608 static const G4double pZ42N52[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
609 static const std::pair<G4int, const G4double*> Z42N52(52,pZ42N52);
610 static const G4double pZ42N53[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
611 static const std::pair<G4int, const G4double*> Z42N53(53,pZ42N53);
612 static const G4double pZ42N54[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
613 static const std::pair<G4int, const G4double*> Z42N54(54,pZ42N54);
614 static const G4double pZ42N55[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
615 static const std::pair<G4int, const G4double*> Z42N55(55,pZ42N55);
616 static const G4double pZ42N56[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
617 static const std::pair<G4int, const G4double*> Z42N56(56,pZ42N56);
618 static const G4double pZ42N58[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
619 static const std::pair<G4int, const G4double*> Z42N58(58,pZ42N58);
620 static const std::pair<G4int, const G4double*> Z42[N42]={Z42N50, Z42N52, Z42N53, Z42N54,
621 Z42N55, Z42N56, Z42N58};
622 //==> Mo(Z=43)
623 static const G4int N43=1;
624 static const G4double pZ43N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
625 static const std::pair<G4int, const G4double*> Z43N0(0,pZ43N0);
626 static const std::pair<G4int, const G4double*> Z43[N43]={Z43N0};
627 //==> Ru(Z=44)
628 static const G4int N44=7;
629 static const G4double pZ44N52[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
630 static const std::pair<G4int, const G4double*> Z44N52(52,pZ44N52);
631 static const G4double pZ44N54[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
632 static const std::pair<G4int, const G4double*> Z44N54(54,pZ44N54);
633 static const G4double pZ44N55[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
634 static const std::pair<G4int, const G4double*> Z44N55(55,pZ44N55);
635 static const G4double pZ44N56[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
636 static const std::pair<G4int, const G4double*> Z44N56(56,pZ44N56);
637 static const G4double pZ44N57[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
638 static const std::pair<G4int, const G4double*> Z44N57(57,pZ44N57);
639 static const G4double pZ44N58[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
640 static const std::pair<G4int, const G4double*> Z44N58(58,pZ44N58);
641 static const G4double pZ44N60[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
642 static const std::pair<G4int, const G4double*> Z44N60(60,pZ44N60);
643 static const std::pair<G4int, const G4double*> Z44[N44]={Z44N52, Z44N54, Z44N55, Z44N56,
644 Z44N57, Z44N58, Z44N60};
645 //==> Rh(Z=45)
646 static const G4int N45=1;
647 static const G4double pZ45N58[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
648 static const std::pair<G4int, const G4double*> Z45N58(58,pZ45N58);
649 static const std::pair<G4int, const G4double*> Z45[N45]={Z45N58};
650 //==> Pd(Z=46)
651 static const G4int N46=6;
652 static const G4double pZ46N56[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
653 static const std::pair<G4int, const G4double*> Z46N56(56,pZ46N56);
654 static const G4double pZ46N58[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
655 static const std::pair<G4int, const G4double*> Z46N58(58,pZ46N58);
656 static const G4double pZ46N59[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
657 static const std::pair<G4int, const G4double*> Z46N59(59,pZ46N59);
658 static const G4double pZ46N60[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
659 static const std::pair<G4int, const G4double*> Z46N60(60,pZ46N60);
660 static const G4double pZ46N62[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
661 static const std::pair<G4int, const G4double*> Z46N62(62,pZ46N62);
662 static const G4double pZ46N64[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
663 static const std::pair<G4int, const G4double*> Z46N64(64,pZ46N64);
664 static const std::pair<G4int, const G4double*> Z46[N46]={Z46N56, Z46N58, Z46N59,
665 Z46N60, Z46N62, Z46N64};
666 //==> Ag(Z=47)
667 static const G4int N47=2;
668 static const G4double pZ47N60[5]={.018, 5., .1, .004, .003};
669 static const std::pair<G4int, const G4double*> Z47N60(60,pZ47N60);
670 static const G4double pZ47N62[5]={.018, 4., .015, .06, .0008};
671 static const std::pair<G4int, const G4double*> Z47N62(62,pZ47N62);
672 static const std::pair<G4int, const G4double*> Z47[N47]={Z47N60, Z47N62};
673 //==> Cd(Z=48)
674 static const G4int N48=8;
675 static const G4double pZ48N58[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
676 static const std::pair<G4int, const G4double*> Z48N58(58,pZ48N58);
677 static const G4double pZ48N60[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
678 static const std::pair<G4int, const G4double*> Z48N60(60,pZ48N60);
679 static const G4double pZ48N62[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
680 static const std::pair<G4int, const G4double*> Z48N62(62,pZ48N62);
681 static const G4double pZ48N63[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
682 static const std::pair<G4int, const G4double*> Z48N63(63,pZ48N63);
683 static const G4double pZ48N64[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
684 static const std::pair<G4int, const G4double*> Z48N64(64,pZ48N64);
685 static const G4double pZ48N65[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
686 static const std::pair<G4int, const G4double*> Z48N65(65,pZ48N65);
687 static const G4double pZ48N66[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
688 static const std::pair<G4int, const G4double*> Z48N66(66,pZ48N66);
689 static const G4double pZ48N68[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
690 static const std::pair<G4int, const G4double*> Z48N68(68,pZ48N68);
691 static const std::pair<G4int, const G4double*> Z48[N48]={Z48N58, Z48N60, Z48N62, Z48N63,
692 Z48N64, Z48N65, Z48N66, Z48N68};
693 //==> In(Z=49)
694 static const G4int N49=2;
695 static const G4double pZ49N64[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
696 static const std::pair<G4int, const G4double*> Z49N64(64,pZ49N64);
697 static const G4double pZ49N66[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
698 static const std::pair<G4int, const G4double*> Z49N66(66,pZ49N66);
699 static const std::pair<G4int, const G4double*> Z49[N49]={Z49N64, Z49N66};
700 //==> Sn(Z=50)
701 static const G4int N50=10;
702 static const G4double pZ50N62[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
703 static const std::pair<G4int, const G4double*> Z50N62(62,pZ50N62);
704 static const G4double pZ50N64[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
705 static const std::pair<G4int, const G4double*> Z50N64(64,pZ50N64);
706 static const G4double pZ50N65[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
707 static const std::pair<G4int, const G4double*> Z50N65(65,pZ50N65);
708 static const G4double pZ50N66[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
709 static const std::pair<G4int, const G4double*> Z50N66(66,pZ50N66);
710 static const G4double pZ50N67[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
711 static const std::pair<G4int, const G4double*> Z50N67(67,pZ50N67);
712 static const G4double pZ50N68[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
713 static const std::pair<G4int, const G4double*> Z50N68(68,pZ50N68);
714 static const G4double pZ50N69[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
715 static const std::pair<G4int, const G4double*> Z50N69(69,pZ50N69);
716 static const G4double pZ50N70[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
717 static const std::pair<G4int, const G4double*> Z50N70(70,pZ50N70);
718 static const G4double pZ50N72[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
719 static const std::pair<G4int, const G4double*> Z50N72(72,pZ50N72);
720 static const G4double pZ50N74[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
721 static const std::pair<G4int, const G4double*> Z50N74(74,pZ50N74);
722 static const std::pair<G4int, const G4double*> Z50[N50]={Z50N62, Z50N64, Z50N65, Z50N66,
723 Z50N67, Z50N68, Z50N69, Z50N70,
724 Z50N72, Z50N74};
725 //==> Sb(Z=51)
726 static const G4int N51=2;
727 static const G4double pZ51N70[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
728 static const std::pair<G4int, const G4double*> Z51N70(70,pZ51N70);
729 static const G4double pZ51N72[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
730 static const std::pair<G4int, const G4double*> Z51N72(72,pZ51N72);
731 static const std::pair<G4int, const G4double*> Z51[N51]={Z51N70, Z51N72};
732 //==> Te(Z=52)
733 static const G4int N52=8;
734 static const G4double pZ52N68[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
735 static const std::pair<G4int, const G4double*> Z52N68(68,pZ52N68);
736 static const G4double pZ52N70[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
737 static const std::pair<G4int, const G4double*> Z52N70(70,pZ52N70);
738 static const G4double pZ52N71[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
739 static const std::pair<G4int, const G4double*> Z52N71(71,pZ52N71);
740 static const G4double pZ52N72[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
741 static const std::pair<G4int, const G4double*> Z52N72(72,pZ52N72);
742 static const G4double pZ52N73[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
743 static const std::pair<G4int, const G4double*> Z52N73(73,pZ52N73);
744 static const G4double pZ52N74[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
745 static const std::pair<G4int, const G4double*> Z52N74(74,pZ52N74);
746 static const G4double pZ52N76[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
747 static const std::pair<G4int, const G4double*> Z52N76(76,pZ52N76);
748 static const G4double pZ52N78[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
749 static const std::pair<G4int, const G4double*> Z52N78(78,pZ52N78);
750 static const std::pair<G4int, const G4double*> Z52[N52]={Z52N68, Z52N70, Z52N71, Z52N72,
751 Z52N73, Z52N74, Z52N76, Z52N78};
752 //==> I (Z=53)
753 static const G4int N53=1;
754 static const G4double pZ53N74[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
755 static const std::pair<G4int, const G4double*> Z53N74(74,pZ53N74);
756 static const std::pair<G4int, const G4double*> Z53[N53]={Z53N74};
757 //==> Xe(Z=54)
758 static const G4int N54=9;
759 static const G4double pZ54N70[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
760 static const std::pair<G4int, const G4double*> Z54N70(70,pZ54N70);
761 static const G4double pZ54N72[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
762 static const std::pair<G4int, const G4double*> Z54N72(72,pZ54N72);
763 static const G4double pZ54N74[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
764 static const std::pair<G4int, const G4double*> Z54N74(74,pZ54N74);
765 static const G4double pZ54N75[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
766 static const std::pair<G4int, const G4double*> Z54N75(75,pZ54N75);
767 static const G4double pZ54N76[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
768 static const std::pair<G4int, const G4double*> Z54N76(76,pZ54N76);
769 static const G4double pZ54N77[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
770 static const std::pair<G4int, const G4double*> Z54N77(77,pZ54N77);
771 static const G4double pZ54N78[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
772 static const std::pair<G4int, const G4double*> Z54N78(78,pZ54N78);
773 static const G4double pZ54N80[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
774 static const std::pair<G4int, const G4double*> Z54N80(80,pZ54N80);
775 static const G4double pZ54N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
776 static const std::pair<G4int, const G4double*> Z54N82(82,pZ54N82);
777 static const std::pair<G4int, const G4double*> Z54[N54]={Z54N70, Z54N72, Z54N74,
778 Z54N75, Z54N76, Z54N77,
779 Z54N78, Z54N80, Z54N82};
780 //==> Cs(Z=55)
781 static const G4int N55=1;
782 static const G4double pZ55N78[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
783 static const std::pair<G4int, const G4double*> Z55N78(78,pZ55N78);
784 static const std::pair<G4int, const G4double*> Z55[N55]={Z55N78};
785 //==> Ba(Z=56)
786 static const G4int N56=7;
787 static const G4double pZ56N74[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
788 static const std::pair<G4int, const G4double*> Z56N74(70,pZ56N74);
789 static const G4double pZ56N76[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
790 static const std::pair<G4int, const G4double*> Z56N76(71,pZ56N76);
791 static const G4double pZ56N78[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
792 static const std::pair<G4int, const G4double*> Z56N78(72,pZ56N78);
793 static const G4double pZ56N79[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
794 static const std::pair<G4int, const G4double*> Z56N79(73,pZ56N79);
795 static const G4double pZ56N80[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
796 static const std::pair<G4int, const G4double*> Z56N80(74,pZ56N80);
797 static const G4double pZ56N81[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
798 static const std::pair<G4int, const G4double*> Z56N81(76,pZ56N81);
799 static const G4double pZ56N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
800 static const std::pair<G4int, const G4double*> Z56N82(78,pZ56N82);
801 static const std::pair<G4int, const G4double*> Z56[N56]={Z56N74, Z56N76, Z56N78, Z56N79,
802 Z56N80, Z56N81, Z56N82};
803 //==> La(Z=57)
804 static const G4int N57=2;
805 static const G4double pZ57N81[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
806 static const std::pair<G4int, const G4double*> Z57N81(81,pZ57N81);
807 static const G4double pZ57N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
808 static const std::pair<G4int, const G4double*> Z57N82(82,pZ57N82);
809 static const std::pair<G4int, const G4double*> Z57[N57]={Z57N81, Z57N82};
810 //==> Ce(Z=58)
811 static const G4int N58=4;
812 static const G4double pZ58N78[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
813 static const std::pair<G4int, const G4double*> Z58N78(78,pZ58N78);
814 static const G4double pZ58N80[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
815 static const std::pair<G4int, const G4double*> Z58N80(80,pZ58N80);
816 static const G4double pZ58N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
817 static const std::pair<G4int, const G4double*> Z58N82(82,pZ58N82);
818 static const G4double pZ58N84[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
819 static const std::pair<G4int, const G4double*> Z58N84(84,pZ58N84);
820 static const std::pair<G4int, const G4double*> Z58[N58]={Z58N78, Z58N80, Z58N82, Z58N84};
821 //==> Pr(Z=59)
822 static const G4int N59=1;
823 static const G4double pZ59N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
824 static const std::pair<G4int, const G4double*> Z59N82(82,pZ59N82);
825 static const std::pair<G4int, const G4double*> Z59[N59]={Z59N82};
826 //==> Nd(Z=60)
827 static const G4int N60=7;
828 static const G4double pZ60N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
829 static const std::pair<G4int, const G4double*> Z60N82(82,pZ60N82);
830 static const G4double pZ60N83[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
831 static const std::pair<G4int, const G4double*> Z60N83(83,pZ60N83);
832 static const G4double pZ60N84[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
833 static const std::pair<G4int, const G4double*> Z60N84(84,pZ60N84);
834 static const G4double pZ60N85[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
835 static const std::pair<G4int, const G4double*> Z60N85(85,pZ60N85);
836 static const G4double pZ60N86[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
837 static const std::pair<G4int, const G4double*> Z60N86(86,pZ60N86);
838 static const G4double pZ60N88[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
839 static const std::pair<G4int, const G4double*> Z60N88(88,pZ60N88);
840 static const G4double pZ60N90[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
841 static const std::pair<G4int, const G4double*> Z60N90(90,pZ60N90);
842 static const std::pair<G4int, const G4double*> Z60[N60]={Z60N82, Z60N83, Z60N84, Z60N85,
843 Z60N86, Z60N88, Z60N90};
844 //==> Mo(Z=61)
845 static const G4int N61=1;
846 static const G4double pZ61N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
847 static const std::pair<G4int, const G4double*> Z61N0(0,pZ61N0);
848 static const std::pair<G4int, const G4double*> Z61[N61]={Z61N0};
849 //==> Sm(Z=62)
850 static const G4int N62=7;
851 static const G4double pZ62N82[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
852 static const std::pair<G4int, const G4double*> Z62N82(82,pZ62N82);
853 static const G4double pZ62N85[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
854 static const std::pair<G4int, const G4double*> Z62N85(85,pZ62N85);
855 static const G4double pZ62N86[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
856 static const std::pair<G4int, const G4double*> Z62N86(86,pZ62N86);
857 static const G4double pZ62N87[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
858 static const std::pair<G4int, const G4double*> Z62N87(87,pZ62N87);
859 static const G4double pZ62N88[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
860 static const std::pair<G4int, const G4double*> Z62N88(88,pZ62N88);
861 static const G4double pZ62N90[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
862 static const std::pair<G4int, const G4double*> Z62N90(90,pZ62N90);
863 static const G4double pZ62N92[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
864 static const std::pair<G4int, const G4double*> Z62N92(92,pZ62N92);
865 static const std::pair<G4int, const G4double*> Z62[N62]={Z62N82, Z62N85, Z62N86, Z62N87,
866 Z62N88, Z62N90, Z62N92};
867 //==> Eu(Z=63)
868 static const G4int N63=2;
869 static const G4double pZ63N88[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
870 static const std::pair<G4int, const G4double*> Z63N88(88,pZ63N88);
871 static const G4double pZ63N90[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
872 static const std::pair<G4int, const G4double*> Z63N90(90,pZ63N90);
873 static const std::pair<G4int, const G4double*> Z63[N63]={Z63N88, Z63N90};
874 //==> Gd(Z=64)
875 static const G4int N64=7;
876 static const G4double pZ64N88[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
877 static const std::pair<G4int, const G4double*> Z64N88(88,pZ64N88);
878 static const G4double pZ64N90[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
879 static const std::pair<G4int, const G4double*> Z64N90(90,pZ64N90);
880 static const G4double pZ64N91[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
881 static const std::pair<G4int, const G4double*> Z64N91(91,pZ64N91);
882 static const G4double pZ64N92[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
883 static const std::pair<G4int, const G4double*> Z64N92(92,pZ64N92);
884 static const G4double pZ64N93[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
885 static const std::pair<G4int, const G4double*> Z64N93(93,pZ64N93);
886 static const G4double pZ64N94[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
887 static const std::pair<G4int, const G4double*> Z64N94(94,pZ64N94);
888 static const G4double pZ64N96[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
889 static const std::pair<G4int, const G4double*> Z64N96(96,pZ64N96);
890 static const std::pair<G4int, const G4double*> Z64[N64]={Z64N88, Z64N90, Z64N91, Z64N92,
891 Z64N93, Z64N94, Z64N96};
892 //==> Tb(Z=65)
893 static const G4int N65=1;
894 static const G4double pZ65N94[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
895 static const std::pair<G4int, const G4double*> Z65N94(82,pZ65N94);
896 static const std::pair<G4int, const G4double*> Z65[N65]={Z65N94};
897 //==> Dy(Z=66)
898 static const G4int N66=7;
899 static const G4double pZ66N90[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
900 static const std::pair<G4int, const G4double*> Z66N90(90,pZ66N90);
901 static const G4double pZ66N92[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
902 static const std::pair<G4int, const G4double*> Z66N92(92,pZ66N92);
903 static const G4double pZ66N94[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
904 static const std::pair<G4int, const G4double*> Z66N94(94,pZ66N94);
905 static const G4double pZ66N95[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
906 static const std::pair<G4int, const G4double*> Z66N95(95,pZ66N95);
907 static const G4double pZ66N96[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
908 static const std::pair<G4int, const G4double*> Z66N96(96,pZ66N96);
909 static const G4double pZ66N97[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
910 static const std::pair<G4int, const G4double*> Z66N97(97,pZ66N97);
911 static const G4double pZ66N98[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
912 static const std::pair<G4int, const G4double*> Z66N98(98,pZ66N98);
913 static const std::pair<G4int, const G4double*> Z66[N66]={Z66N90, Z66N92, Z66N94, Z66N95,
914 Z66N96, Z66N97, Z66N98};
915 //==> Ho(Z=67)
916 static const G4int N67=1;
917 static const G4double pZ67N98[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
918 static const std::pair<G4int, const G4double*> Z67N98(98,pZ67N98);
919 static const std::pair<G4int, const G4double*> Z67[N67]={Z67N98};
920 //==> Er(Z=68)
921 static const G4int N68=6;
922 static const G4double pZ68N94[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
923 static const std::pair<G4int, const G4double*> Z68N94(94,pZ68N94);
924 static const G4double pZ68N96[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
925 static const std::pair<G4int, const G4double*> Z68N96(96,pZ68N96);
926 static const G4double pZ68N98[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
927 static const std::pair<G4int, const G4double*> Z68N98(98,pZ68N98);
928 static const G4double pZ68N99[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
929 static const std::pair<G4int, const G4double*> Z68N99(99,pZ68N99);
930 static const G4double pZ68N100[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
931 static const std::pair<G4int, const G4double*> Z68N100(100,pZ68N100);
932 static const G4double pZ68N102[5]={.01, 27., 0., 0., 1.}; // *** NotImplemented ***
933 static const std::pair<G4int, const G4double*> Z68N102(102,pZ68N102);
934 static const std::pair<G4int, const G4double*> Z68[N68]={Z68N94, Z68N96, Z68N98,
935 Z68N99, Z68N100, Z68N102};
936 //==> Tm(Z=69)
937 static const G4int N69=1;
938 static const G4double pZ69N100[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
939 static const std::pair<G4int, const G4double*> Z69N100(100,pZ69N100);
940 static const std::pair<G4int, const G4double*> Z69[N69]={Z69N100};
941 //==> Yb(Z=70)
942 static const G4int N70=7;
943 static const G4double pZ70N98[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
944 static const std::pair<G4int, const G4double*> Z70N98(98,pZ70N98);
945 static const G4double pZ70N100[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
946 static const std::pair<G4int, const G4double*> Z70N100(100,pZ70N100);
947 static const G4double pZ70N101[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
948 static const std::pair<G4int, const G4double*> Z70N101(101,pZ70N101);
949 static const G4double pZ70N102[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
950 static const std::pair<G4int, const G4double*> Z70N102(102,pZ70N102);
951 static const G4double pZ70N103[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
952 static const std::pair<G4int, const G4double*> Z70N103(103,pZ70N103);
953 static const G4double pZ70N104[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
954 static const std::pair<G4int, const G4double*> Z70N104(104,pZ70N104);
955 static const G4double pZ70N106[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
956 static const std::pair<G4int, const G4double*> Z70N106(106,pZ70N106);
957 static const std::pair<G4int, const G4double*> Z70[N70]={Z70N98, Z70N100, Z70N101,
958 Z70N102, Z70N103, Z70N104,
959 Z70N106};
960 //==> Lu(Z=71)
961 static const G4int N71=2;
962 static const G4double pZ71N104[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
963 static const std::pair<G4int, const G4double*> Z71N104(104,pZ71N104);
964 static const G4double pZ71N105[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
965 static const std::pair<G4int, const G4double*> Z71N105(105,pZ71N105);
966 static const std::pair<G4int, const G4double*> Z71[N71]={Z71N104, Z71N105};
967 //==> Hf(Z=72)
968 static const G4int N72=6;
969 static const G4double pZ72N102[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
970 static const std::pair<G4int, const G4double*> Z72N102(102,pZ72N102);
971 static const G4double pZ72N104[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
972 static const std::pair<G4int, const G4double*> Z72N104(104,pZ72N104);
973 static const G4double pZ72N105[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
974 static const std::pair<G4int, const G4double*> Z72N105(105,pZ72N105);
975 static const G4double pZ72N106[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
976 static const std::pair<G4int, const G4double*> Z72N106(106,pZ72N106);
977 static const G4double pZ72N107[5]={.018, 5., .1, .03, .002}; // *** NotImplemented ***
978 static const std::pair<G4int, const G4double*> Z72N107(107,pZ72N107);
979 static const G4double pZ72N108[5]={.01, 27., 0., 0., 1.}; // *** NotImplemented ***
980 static const std::pair<G4int, const G4double*> Z72N108(108,pZ72N108);
981 static const std::pair<G4int, const G4double*> Z72[N72]={Z72N102, Z72N104, Z72N105,
982 Z72N106, Z72N107, Z72N108};
983 //==> Ta(Z=73)
984 static const G4int N73=1;
985 static const G4double pZ73N108[5]={.0065, 2., .15, .001, .0012};
986 static const std::pair<G4int, const G4double*> Z73N108(108,pZ73N108);
987 static const std::pair<G4int, const G4double*> Z73[N73]={Z73N108};
988 //==> W (Z=74)
989 static const G4int N74=5;
990 static const G4double pZ74N106[5]={.014, 5., .03, .04, .001}; // *** No DATA ***
991 static const std::pair<G4int, const G4double*> Z74N106(106,pZ74N106);
992 static const G4double pZ74N108[5]={.013, 4.5, .35, .045, .001};
993 static const std::pair<G4int, const G4double*> Z74N108(108,pZ74N108);
994 static const G4double pZ74N109[5]={.012, 4., .02, .04, .001};
995 static const std::pair<G4int, const G4double*> Z74N109(109,pZ74N109);
996 static const G4double pZ74N110[5]={.014, 6., .02, .03, .0015};
997 static const std::pair<G4int, const G4double*> Z74N110(110,pZ74N110);
998 static const G4double pZ74N112[5]={.014, 6., .02, .03, .0015};
999 static const std::pair<G4int, const G4double*> Z74N112(112,pZ74N112);
1000 static const std::pair<G4int, const G4double*> Z74[N74]={Z74N106, Z74N108, Z74N109,
1001 Z74N110, Z74N112};
1002 //==> Re(Z=75)
1003 static const G4int N75=2;
1004 static const G4double pZ75N110[5]={.015, 4., .2, .0001, .0021};
1005 static const std::pair<G4int, const G4double*> Z75N110(110,pZ75N110);
1006 static const G4double pZ75N112[5]={.015, 4., .1, .0001, .002};
1007 static const std::pair<G4int, const G4double*> Z75N112(112,pZ75N112);
1008 static const std::pair<G4int, const G4double*> Z75[N75]={Z75N110, Z75N112};
1009 //==> Os(Z=76)
1010 static const G4int N76=7;
1011 static const G4double pZ76N108[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1012 static const std::pair<G4int, const G4double*> Z76N108(108,pZ76N108);
1013 static const G4double pZ76N110[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1014 static const std::pair<G4int, const G4double*> Z76N110(110,pZ76N110);
1015 static const G4double pZ76N111[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1016 static const std::pair<G4int, const G4double*> Z76N111(111,pZ76N111);
1017 static const G4double pZ76N112[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1018 static const std::pair<G4int, const G4double*> Z76N112(112,pZ76N112);
1019 static const G4double pZ76N113[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1020 static const std::pair<G4int, const G4double*> Z76N113(113,pZ76N113);
1021 static const G4double pZ76N114[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1022 static const std::pair<G4int, const G4double*> Z76N114(114,pZ76N114);
1023 static const G4double pZ76N116[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1024 static const std::pair<G4int, const G4double*> Z76N116(116,pZ76N116);
1025 static const std::pair<G4int, const G4double*> Z76[N76]={Z76N108, Z76N110, Z76N111,
1026 Z76N112, Z76N113, Z76N114,
1027 Z76N116};
1028 //==> Ir(Z=77)
1029 static const G4int N77=2;
1030 static const G4double pZ77N114[5]={.012, 3., .1, .0001, .003};
1031 static const std::pair<G4int, const G4double*> Z77N114(114,pZ77N114);
1032 static const G4double pZ77N116[5]={.012, 3., .08, .0001, .002};
1033 static const std::pair<G4int, const G4double*> Z77N116(116,pZ77N116);
1034 static const std::pair<G4int, const G4double*> Z77[N77]={Z77N114, Z77N116};
1035 //==> Pt(Z=78)
1036 static const G4int N78=6;
1037 static const G4double pZ78N112[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1038 static const std::pair<G4int, const G4double*> Z78N112(112,pZ78N112);
1039 static const G4double pZ78N114[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1040 static const std::pair<G4int, const G4double*> Z78N114(114,pZ78N114);
1041 static const G4double pZ78N116[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1042 static const std::pair<G4int, const G4double*> Z78N116(116,pZ78N116);
1043 static const G4double pZ78N117[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1044 static const std::pair<G4int, const G4double*> Z78N117(117,pZ78N117);
1045 static const G4double pZ78N118[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1046 static const std::pair<G4int, const G4double*> Z78N118(118,pZ78N118);
1047 static const G4double pZ78N120[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1048 static const std::pair<G4int, const G4double*> Z78N120(120,pZ78N120);
1049 static const std::pair<G4int, const G4double*> Z78[N78]={Z78N112, Z78N114, Z78N116,
1050 Z78N117, Z78N118, Z78N120};
1051 //==> Au(Z=79)
1052 static const G4int N79=1;
1053 static const G4double pZ79N118[5]={.012, 4., .2, .0001, .002};
1054 static const std::pair<G4int, const G4double*> Z79N118(118,pZ79N118);
1055 static const std::pair<G4int, const G4double*> Z79[N79]={Z79N118};
1056 //==> Hg(Z=80)
1057 static const G4int N80=7;
1058 static const G4double pZ80N116[5]={.022, 11., .006, .044, .001};
1059 static const std::pair<G4int, const G4double*> Z80N116(116,pZ80N116);
1060 static const G4double pZ80N118[5]={.024, 10., .04, .018, .0016};
1061 static const std::pair<G4int, const G4double*> Z80N118(118,pZ80N118);
1062 static const G4double pZ80N119[5]={.017, 8., .02, .02, .0016};
1063 static const std::pair<G4int, const G4double*> Z80N119(119,pZ80N119);
1064 static const G4double pZ80N120[5]={.021, 9., .03, .02, .0016};
1065 static const std::pair<G4int, const G4double*> Z80N120(120,pZ80N120);
1066 static const G4double pZ80N121[5]={.01, 7., .025, .025, .0016};
1067 static const std::pair<G4int, const G4double*> Z80N121(121,pZ80N121);
1068 static const G4double pZ80N122[5]={.023, 9., .008, .045, .0012};
1069 static const std::pair<G4int, const G4double*> Z80N122(122,pZ80N122);
1070 static const G4double pZ80N124[5]={.018, 8., .007, .04, .0013};
1071 static const std::pair<G4int, const G4double*> Z80N124(124,pZ80N124);
1072 static const std::pair<G4int, const G4double*> Z80[N80]={Z80N116, Z80N118, Z80N119,
1073 Z80N120, Z80N121, Z80N122,
1074 Z80N124};
1075 //==> Tl(Z=81)
1076 static const G4int N81=2;
1077 static const G4double pZ81N122[5]={.018, 5., .1, .03, .002}; // *** No DATA ***
1078 static const std::pair<G4int, const G4double*> Z81N122(122,pZ81N122);
1079 static const G4double pZ81N124[5]={.01, 27., 0., 0., 1.}; // *** No DATA ***
1080 static const std::pair<G4int, const G4double*> Z81N124(124,pZ81N124);
1081 static const std::pair<G4int, const G4double*> Z81[N81]={Z81N122, Z81N124};
1082 //==> Pb(Z=82)
1083 static const G4int N82=4;
1084 static const G4double pZ82N122[5]={.032, 9., .001, .13, .004};
1085 static const std::pair<G4int, const G4double*> Z82N122(122,pZ82N122);
1086 static const G4double pZ82N124[5]={.02, 7., .005, .14, .002};
1087 static const std::pair<G4int, const G4double*> Z82N124(124,pZ82N124);
1088 static const G4double pZ82N125[5]={.021, 9., .0012, .05, .01};
1089 static const std::pair<G4int, const G4double*> Z82N125(125,pZ82N125);
1090 static const G4double pZ82N126[5]={.049, 14., .0007, .145, .001};
1091 static const std::pair<G4int, const G4double*> Z82N126(126,pZ82N126);
1092 static const std::pair<G4int, const G4double*> Z82[N82]={Z82N122, Z82N124, Z82N125,
1093 Z82N126};
1094 //==> Bi(Z=83)
1095 static const G4int N83=1;
1096 static const G4double pZ83N126[5]={.033, 10., .001, .13, .006};
1097 static const std::pair<G4int, const G4double*> Z83N126(126,pZ83N126);
1098 static const std::pair<G4int, const G4double*> Z83[N83]={Z83N126};
1099 //==> Po(Z=84)
1100 static const G4int N84=1;
1101 static const G4double pZ84N0[5]={.01, 27., 0., 0., 1.}; // *** NoStableIsotopes ***
1102 static const std::pair<G4int, const G4double*> Z84N0(0,pZ84N0);
1103 static const std::pair<G4int, const G4double*> Z84[N84]={Z84N0};
1104 //==> At(Z=85)
1105 static const G4int N85=1;
1106 static const G4double pZ85N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
1107 static const std::pair<G4int, const G4double*> Z85N0(0,pZ85N0);
1108 static const std::pair<G4int, const G4double*> Z85[N85]={Z85N0};
1109 //==> Rn(Z=86)
1110 static const G4int N86=1;
1111 static const G4double pZ86N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
1112 static const std::pair<G4int, const G4double*> Z86N0(0,pZ86N0);
1113 static const std::pair<G4int, const G4double*> Z86[N86]={Z86N0};
1114 //==> Fr(Z=87)
1115 static const G4int N87=1;
1116 static const G4double pZ87N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
1117 static const std::pair<G4int, const G4double*> Z87N0(0,pZ87N0);
1118 static const std::pair<G4int, const G4double*> Z87[N87]={Z87N0};
1119 //==> Ra(Z=88)
1120 static const G4int N88=1;
1121 static const G4double pZ88N138[5]={.012, 4.4, .068, .051, .0008};
1122 static const std::pair<G4int, const G4double*> Z88N138(138,pZ88N138);
1123 static const std::pair<G4int, const G4double*> Z88[N88]={Z88N138};
1124 //==> Ac(Z=89)
1125 static const G4int N89=1;
1126 static const G4double pZ89N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
1127 static const std::pair<G4int, const G4double*> Z89N0(0,pZ89N0);
1128 static const std::pair<G4int, const G4double*> Z89[N89]={Z89N0};
1129 //==> Th(Z=90)
1130 static const G4int N90=1;
1131 static const G4double pZ90N142[5]={.01, 3.6, .07, .029, .0009};
1132 static const std::pair<G4int, const G4double*> Z90N142(142,pZ90N142);
1133 static const std::pair<G4int, const G4double*> Z90[N90]={Z90N142};
1134 //==> Pa(Z=91)
1135 static const G4int N91=1;
1136 static const G4double pZ91N0[5]={.01, 27., 0., 0., 1.}; // *** NoStableIsotopes ***
1137 static const std::pair<G4int, const G4double*> Z91N0(0,pZ91N0);
1138 static const std::pair<G4int, const G4double*> Z91[N91]={Z91N0};
1139 //==> U (Z=92)
1140 static const G4int N92=2;
1141 static const G4double pZ92N143[5]={.005, 2.5, .06, .008, .002};
1142 static const std::pair<G4int, const G4double*> Z92N143(143,pZ92N143);
1143 static const G4double pZ92N146[5]={.009, 3.1, .04, .025, .001};
1144 static const std::pair<G4int, const G4double*> Z92N146(146,pZ92N146);
1145 static const std::pair<G4int, const G4double*> Z92[N92]={Z92N143, Z92N146};
1146 //==> Np(Z=93)
1147 static const G4int N93=1;
1148 static const G4double pZ93N144[5]={.009, 2.35, .35, .003, .0007};
1149 static const std::pair<G4int, const G4double*> Z93N144(144,pZ93N144);
1150 static const std::pair<G4int, const G4double*> Z93[N93]={Z93N144};
1151 //==> Pu(Z=94)
1152 static const G4int N94=3;
1153 static const G4double pZ94N145[5]={.005, 2.6, .04, .023, .0002}; // *** Artificial ***
1154 static const std::pair<G4int, const G4double*> Z94N145(145,pZ94N145);
1155 static const G4double pZ94N148[5]={.009, 3.1, .06, .02, .001}; // *** Artificial ***
1156 static const std::pair<G4int, const G4double*> Z94N148(148,pZ94N148);
1157 static const G4double pZ94N150[5]={.01, 27., 0., 0., 1.};
1158 static const std::pair<G4int, const G4double*> Z94N150(150,pZ94N150);
1159 static const std::pair<G4int, const G4double*> Z94[N94]={Z94N145, Z94N148, Z94N150};
1160 //==> Am(Z=95)
1161 static const G4int N95=1;
1162 static const G4double pZ95N0[5]={.018, 5., .1, .03, .002}; // *** NoStableIsotopes ***
1163 static const std::pair<G4int, const G4double*> Z95N0(0,pZ95N0);
1164 static const std::pair<G4int, const G4double*> Z95[N95]={Z95N0};
1165 //==> Cm(Z=96)
1166 static const G4int N96=1;
1167 static const G4double pZ96N151[5]={.005, 2.5, .07, .027, .0009};
1168 static const std::pair<G4int, const G4double*> Z96N151(151,pZ96N151);
1169 static const std::pair<G4int, const G4double*> Z96[N96]={Z96N151};
1170
1171 static const G4int NZ=97; // #of Elements covered by CHIPS
1172 static const std::pair<G4int, const G4double*>* Pars[NZ]={Z0,Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8,Z9,
1173 Z10,Z11,Z12,Z13,Z14,Z15,Z16,Z17,Z18,Z19,Z20,Z21,Z22,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z30,
1174 Z31,Z32,Z33,Z34,Z35,Z36,Z37,Z38,Z39,Z40,Z41,Z42,Z43,Z44,Z45,Z46,Z47,Z48,Z49,Z50,Z51,
1175 Z52,Z53,Z54,Z55,Z56,Z57,Z58,Z59,Z60,Z61,Z62,Z63,Z64,Z65,Z66,Z67,Z68,Z69,Z70,Z71,Z72,
1176 Z73,Z74,Z75,Z76,Z77,Z78,Z79,Z80,Z81,Z82,Z83,Z84,Z85,Z86,Z87,Z88,Z89,Z90,Z91,Z92,Z93,
1177 Z94,Z95,Z96};
1178 static const G4int NIso[NZ]={N0,N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15,N16,
1179 N17,N18,N19,N20,N21,N22,N23,N24,N25,N26,N27,N28,N29,N30,N31,N32,N33,N34,N35,N36,N37,
1180 N38,N39,N40,N41,N42,N43,N44,N45,N46,N47,N48,N49,N50,N51,N52,N53,N54,N55,N56,N57,N58,
1181 N59,N60,N61,N62,N63,N64,N65,N66,N67,N68,N69,N70,N71,N72,N73,N74,N75,N76,N77,N78,N79,
1182 N80,N81,N82,N83,N84,N85,N86,N87,N88,N89,N90,N91,N92,N93,N94,N95,N96};
1183 static G4int mZ=0;
1184 static G4int mN=0;
1185 static G4double P=.1; // s=SQRT(M_N^2+M_h^2+2*E_h*M_N)
1186 static G4double R=0.; // Prototype of the result
1187 static G4double p0=0.;
1188 static G4double p1=0.;
1189 static G4double p2=0.;
1190 static G4double p3=0.;
1191 static G4double p4=0.;
1192 G4int A=Z+N;
1193 if(A<=1) return 0.;
1194 if(p<=0.001) return 1.;
1195 if(Z!=mZ || N!=mN) // Recalculate the parameters for different isotope
1196 {
1197 mZ=Z;
1198 mN=N;
1199 P=p;
1200 if(Z<97 && N<248) // General solution (*** Z/A limits ***)
1201 {
1202 p0=.01; // Default guess if there is no data
1203 p1=0.;
1204 p2=27.;
1205 p3=0.;
1206 p4=1.;
1207 G4int nn=NIso[Z];
1208 G4bool nfound=true;
1209 if(nn) for (G4int in=0; in<nn; in++)
1210 {
1211 std::pair<G4int, const G4double*> curIs=Pars[Z][in];
1212 if(curIs.first == N)
1213 {
1214 const G4double* curT=curIs.second;
1215 p0 = curT[0];
1216 if(p < p0)
1217 {
1218 R=1.;
1219 return R;
1220 }
1221 p1 = curT[1];
1222 p2 = curT[2];
1223 p3 = curT[3];
1224 p4 = curT[4];
1225 nfound = false;
1226 break;
1227 }
1228 }
1229 if(nfound) G4cout<<"-Warning-G4QNeutronCaptureRatio::CSLin: Z="<<Z<<", N="
1230 <<N<<" isotope is not implemented in CHIPS"<<G4endl;
1231 R=std::pow(p0/p,p1);
1232 if(p2>0.)
1233 {
1234 G4double dp=p-p3;
1235 R+=p2*std::exp(-dp*dp/p4);
1236 }
1237 if(R>1.) R=1.;
1238 }
1239 else G4cerr<<"-Warning-G4QNeutronCaptureRatio::CalcR:*Bad A* Z="<<Z<<", N="<<N<<G4endl;
1240 }
1241 else if(std::fabs(p-P)/P<.0001) return R; // Normally used at high energies (direct) only
1242 return R;
1243} // End of CalcCap2IN_Ratio
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4bool
bool G4bool
Definition: G4Types.hh:67
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cerr
G4DLLIMPORT std::ostream G4cerr
G4cout
G4DLLIMPORT std::ostream G4cout
G4QNeutronCaptureRatio::GetRatio
G4double GetRatio(G4double pIU, G4int tgZ, G4int tgN)
Definition: G4QNeutronCaptureRatio.cc:55
G4QNeutronCaptureRatio::GetPointer
static G4QNeutronCaptureRatio * GetPointer()
Definition: G4QNeutronCaptureRatio.cc:48