Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/cross_sections/src/G4ElectroNuclearCrossSection.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/cross_sections/src/G4ElectroNuclearCrossSection.cc (Version 11.3.0) and /processes/hadronic/cross_sections/src/G4ElectroNuclearCrossSection.cc (Version 3.0)


  1 //                                                  1 
  2 // *******************************************    
  3 // * License and Disclaimer                       
  4 // *                                              
  5 // * The  Geant4 software  is  copyright of th    
  6 // * the Geant4 Collaboration.  It is provided    
  7 // * conditions of the Geant4 Software License    
  8 // * LICENSE and available at  http://cern.ch/    
  9 // * include a list of copyright holders.         
 10 // *                                              
 11 // * Neither the authors of this software syst    
 12 // * institutes,nor the agencies providing fin    
 13 // * work  make  any representation or  warran    
 14 // * regarding  this  software system or assum    
 15 // * use.  Please see the license in the file     
 16 // * for the full disclaimer and the limitatio    
 17 // *                                              
 18 // * This  code  implementation is the result     
 19 // * technical work of the GEANT4 collaboratio    
 20 // * By using,  copying,  modifying or  distri    
 21 // * any work based  on the software)  you  ag    
 22 // * use  in  resulting  scientific  publicati    
 23 // * acceptance of all terms of the Geant4 Sof    
 24 // *******************************************    
 25 //                                                
 26 //                                                
 27 // G4 Physics class: G4ElectroNuclearCrossSect    
 28 // Created: M.V. Kossov, CERN/ITEP(Moscow), 10    
 29 // The last update: M.V. Kossov, CERN/ITEP (Mo    
 30                                                   
 31 // A. Dotti 2-May-2013: I keep exactly the sam    
 32 //                      general assumptions ab    
 33 //                      and "lastSeen" reuse),    
 34 //                      mechanism and redue th    
 35 //                      On simple setups this     
 36 //                      faster                    
 37 // W. Pokorski 29-08-2013: Major clean up, res    
 38 //                      of the cross section i    
 39 //                      but per element (Z). A    
 40                                                   
 41 #include <iostream>                               
 42                                                   
 43 #include "G4SystemOfUnits.hh"                     
 44 #include "G4PhysicalConstants.hh"                 
 45 #include "G4HadTmpUtil.hh"                        
 46 #include "G4ElectroNuclearCrossSection.hh"        
 47                                                   
 48 // factory                                        
 49 #include "G4CrossSectionFactory.hh"               
 50 //                                                
 51 G4_DECLARE_XS_FACTORY(G4ElectroNuclearCrossSec    
 52                                                   
 53 //                                                
 54 //A. Dotti 2-May-2013                             
 55 //Optimizing code                                 
 56                                                   
 57 // W. Pokorski 22-August-2013                     
 58 // Moved cache element to header file             
 59 // all static variables moved to private membe    
 60                                                   
 61 // removing all the static consts and putting     
 62                                                   
 63 // Parametrization of the PhotoNucCS              
 64 static const G4double sLowEnergyLimit = 100.;     
 65 static const G4double shd=1.0734;                 
 66 static const G4double poc=0.0375;                 
 67 static const G4double pos=16.5;                   
 68 static const G4double reg=.11;                    
 69 static const G4double mel=0.5109989;              
 70 static const G4double mel2=mel*mel;     // Squ    
 71 static const G4double lmel=G4Log(mel);      //    
 72 //                                                
 73 static const G4int nE=336;                        
 74 static const G4int mLL=nE-1;                      
 75 static const G4double EMi=2.0612; // Minimum t    
 76 static const G4double EMa=50000.; // Maximum t    
 77 static const G4double EMa2=EMa*EMa; // Maximum    
 78 static const G4double lEMi=G4Log(EMi);  // Min    
 79 static const G4double lEMa=G4Log(EMa);  // Max    
 80 static const G4double lEMa2=lEMa*lEMa;  // Max    
 81 static const G4double dlnE=(lEMa-lEMi)/mLL; //    
 82 static const G4double alop=1./137.036/3.141592    
 83 static const G4double le1=(lEMa-1.)*EMa;  // (    
 84 static const G4double leh=(lEMa-.5)*EMa2; // (    
 85 //                                                
 86 static const G4double ha=poc*.5;        // a/2    
 87 static const G4double hab=ha*pos;    // a*b/2     
 88 static const G4double ab=poc*pos;      // a*b     
 89 static const G4double d1 = 1 - reg;               
 90 static const G4double d2 = 2 - reg;               
 91 static const G4double cd=shd/reg;    // c/d       
 92 static const G4double cd1=shd/d1;    // c/d1      
 93 static const G4double cd2=shd/d2;    // c/d2      
 94 static const G4double ele=G4Exp(-reg*lEMa); //    
 95 static const G4double ele1=G4Exp(d1*lEMa); //     
 96 static const G4double ele2=G4Exp(d2*lEMa); //     
 97 //                                                
 98 static const G4double phte=poc*(lEMa-pos)+shd*    
 99 static const G4int    imax=27;   // Not more t    
100 static const G4double eps=0.001; // Accuracy w    
101 //                                                
102 static const G4double infEn = 9.e27;              
103 //                                                
104 static const G4double dM=938.27+939.57; // Mea    
105 static const G4double Q0=843.;          // Coe    
106 static const G4double Q02=Q0*Q0;        // Squ    
107 static const G4double blK0=G4Log(185.);   // C    
108 static const G4double bp=0.85;          // Pow    
109 static const G4double clK0=G4Log(1390.);  // C    
110 static const G4double cp=3.;            // Pow    
111 //                                                
112 static const G4int nN=14;                         
113 //                                                
114 static const G4int LL[nN]={138, 2, 32, 75, 26,    
115 // !! @@ Change it from ln(A) to A approximati    
116 static const G4double A[nN]={1.,2.,3.,4.,6.,7.    
117 //                                                
118 static const G4double P00[nE]={                   
119     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
120     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
121     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
122     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
123     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
124     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
125     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
126     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
127     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
128     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
129     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
130     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
131     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
132     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
133     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
134     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
135     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
136     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
137     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
138     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
139     5.658793e-09,1.254527e-08,2.393031e-08,4.5    
140     8.030425e-07,1.694925e-06,3.434914e-06,6.4    
141     3.581584e-05,4.713652e-05,6.004450e-05,7.4    
142     1.556326e-04,1.828616e-04,2.132878e-04,2.4    
143     3.999890e-04,4.378523e-04,4.735003e-04,5.0    
144     6.088412e-04,6.283321e-04,6.460063e-04,6.6    
145     7.170748e-04,7.295161e-04,7.418666e-04,7.5    
146     8.096392e-04,8.258506e-04,8.432877e-04,8.6    
147     9.410623e-04,9.595795e-04,9.769243e-04,9.9    
148     1.051075e-03,1.065714e-03,1.080833e-03,1.0    
149     1.151571e-03,1.163605e-03,1.175135e-03,1.1    
150     1.229133e-03,1.239639e-03,1.250100e-03,1.2    
151     1.301668e-03,1.311824e-03,1.321929e-03,1.3    
152     1.371743e-03,1.381573e-03,1.391360e-03,1.4    
153     1.439674e-03,1.449216e-03,1.458718e-03,1.4    
154     1.505649e-03,1.514923e-03,1.524160e-03,1.5    
155     1.569826e-03,1.578859e-03,1.587860e-03,1.5    
156     1.632405e-03,1.641227e-03,1.650021e-03,1.6    
157     1.693598e-03,1.702238e-03,1.710855e-03,1.7    
158     1.753604e-03,1.762091e-03,1.770557e-03,1.7    
159     1.812608e-03,1.820965e-03,1.829305e-03,1.8    
160     1.870774e-03,1.879024e-03,1.887260e-03,1.8    
161     1.928252e-03,1.936415e-03,1.944567e-03,1.9    
162     1.985179e-03,1.993274e-03,2.001361e-03,2.0    
163     2.041684e-03,2.049728e-03,2.057767e-03,2.0    
164     2.097884e-03,2.105894e-03,2.113901e-03,2.1    
165     2.153892e-03,2.161884e-03,2.169874e-03,2.1    
166     2.209815e-03,2.217803e-03,2.225791e-03,2.2    
167 static const G4double P10[nE]={                   
168     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
169     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
170     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
171     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
172     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
173     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
174     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
175     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
176     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
177     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
178     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
179     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
180     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
181     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
182     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
183     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
184     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
185     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
186     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
187     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
188     7.765990e-07,1.754155e-06,3.420430e-06,6.6    
189     1.331907e-04,2.897632e-04,6.045196e-04,1.1    
190     7.118943e-03,9.570755e-03,1.245202e-02,1.5    
191     3.600931e-02,4.329446e-02,5.168438e-02,6.1    
192     1.082176e-01,1.207266e-01,1.328639e-01,1.4    
193     1.830362e-01,1.909891e-01,1.984219e-01,2.0    
194     2.310550e-01,2.373266e-01,2.437434e-01,2.5    
195     2.824425e-01,2.925383e-01,3.037299e-01,3.1    
196     3.725249e-01,3.867626e-01,4.005068e-01,4.1    
197     4.647653e-01,4.786698e-01,4.934691e-01,5.0    
198     5.691075e-01,5.832202e-01,5.971567e-01,6.1    
199     6.686149e-01,6.838333e-01,6.994497e-01,7.1    
200     7.837715e-01,8.019397e-01,8.205697e-01,8.3    
201     9.211652e-01,9.428799e-01,9.651625e-01,9.8    
202     1.085653e+00,1.111684e+00,1.138400e+00,1.1    
203     1.282939e+00,1.314184e+00,1.346258e+00,1.3    
204     1.519941e+00,1.557522e+00,1.596116e+00,1.6    
205     1.805343e+00,1.850668e+00,1.897234e+00,1.9    
206     2.149995e+00,2.204816e+00,2.261160e+00,2.3    
207     2.567376e+00,2.633867e+00,2.702233e+00,2.7    
208     3.074206e+00,3.155064e+00,3.238230e+00,3.3    
209     3.691226e+00,3.789797e+00,3.891216e+00,3.9    
210     4.444197e+00,4.564641e+00,4.688605e+00,4.8    
211     5.365163e+00,5.512657e+00,5.664510e+00,5.8    
212     6.494036e+00,6.675037e+00,6.861440e+00,7.0    
213     7.880607e+00,8.103173e+00,8.332447e+00,8.5    
214     9.587066e+00,9.861271e+00,1.014381e+01,1.0    
215     1.169119e+01,1.202964e+01,1.237847e+01,1.2    
216 static const G4double P20[nE]={                   
217     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
218     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
219     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
220     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
221     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
222     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
223     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
224     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
225     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
226     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
227     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
228     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
229     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
230     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
231     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
232     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
233     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
234     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
235     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
236     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
237     1.065910e-04,2.453685e-04,4.892572e-04,9.7    
238     2.212593e-02,4.961468e-02,1.065569e-01,2.1    
239     1.419220e+00,1.950269e+00,2.593461e+00,3.3    
240     8.409767e+00,1.035907e+01,1.267273e+01,1.5    
241     2.982262e+01,3.395557e+01,3.808836e+01,4.2    
242     5.672268e+01,5.996799e+01,6.309401e+01,6.6    
243     7.810711e+01,8.126885e+01,8.460302e+01,8.8    
244     1.067429e+02,1.130306e+02,1.202141e+02,1.2    
245     1.687086e+02,1.796566e+02,1.905484e+02,2.0    
246     2.463419e+02,2.595499e+02,2.740378e+02,2.8    
247     3.550669e+02,3.716193e+02,3.884660e+02,4.0    
248     4.832092e+02,5.052552e+02,5.285702e+02,5.5    
249     6.667118e+02,6.992147e+02,7.335640e+02,7.6    
250     9.370917e+02,9.850666e+02,1.035803e+03,1.0    
251     1.336863e+03,1.407884e+03,1.483006e+03,1.5    
252     1.929001e+03,2.034272e+03,2.145649e+03,2.2    
253     2.807480e+03,2.963847e+03,3.129341e+03,3.3    
254     4.113928e+03,4.346817e+03,4.593401e+03,4.8    
255     6.062289e+03,6.410140e+03,6.778595e+03,7.1    
256     8.976188e+03,9.497193e+03,1.004927e+04,1.0    
257     1.334591e+04,1.412831e+04,1.495767e+04,1.5    
258     1.991550e+04,2.109332e+04,2.234227e+04,2.3    
259     2.981612e+04,3.159335e+04,3.347850e+04,3.5    
260     4.477059e+04,4.745818e+04,5.030985e+04,5.3    
261     6.740741e+04,7.148022e+04,7.580293e+04,8.0    
262     1.017437e+05,1.079281e+05,1.144938e+05,1.2    
263     1.539289e+05,1.633378e+05,1.733295e+05,1.8    
264     2.333914e+05,2.477325e+05,2.629658e+05,2.7    
265 static const G4double P01[nE]={                   
266     0.000000e+00,0.000000e+00,0.000000e+00,2.9    
267     3.199604e-04,4.032996e-04,4.911973e-04,5.8    
268     1.004812e-03,1.124386e-03,1.250223e-03,1.3    
269     1.980284e-03,2.147940e-03,2.320852e-03,2.4    
270     3.231140e-03,3.423010e-03,3.617256e-03,3.8    
271     4.547357e-03,4.723746e-03,4.895195e-03,5.0    
272     5.670813e-03,5.809570e-03,5.942803e-03,6.0    
273     6.530413e-03,6.633483e-03,6.732384e-03,6.8    
274     7.166336e-03,7.242378e-03,7.315204e-03,7.3    
275     7.635427e-03,7.691607e-03,7.745512e-03,7.7    
276     7.983384e-03,8.025299e-03,8.065520e-03,8.1    
277     8.243782e-03,8.275321e-03,8.305622e-03,8.3    
278     8.440441e-03,8.464389e-03,8.487429e-03,8.5    
279     8.590305e-03,8.608648e-03,8.626317e-03,8.6    
280     8.705479e-03,8.719642e-03,8.733300e-03,8.7    
281     8.794677e-03,8.805692e-03,8.816325e-03,8.8    
282     8.864237e-03,8.872859e-03,8.881188e-03,8.8    
283     8.918810e-03,8.925597e-03,8.932158e-03,8.9    
284     8.961852e-03,8.967219e-03,8.972411e-03,8.9    
285     8.995958e-03,9.000224e-03,9.004354e-03,9.0    
286     9.023195e-03,9.026706e-03,9.030244e-03,9.0    
287     9.059416e-03,9.071006e-03,9.084867e-03,9.1    
288     9.181864e-03,9.207136e-03,9.234874e-03,9.2    
289     9.422441e-03,9.471977e-03,9.525862e-03,9.5    
290     9.839376e-03,9.902903e-03,9.963776e-03,1.0    
291     1.021107e-02,1.024941e-02,1.028482e-02,1.0    
292     1.043213e-02,1.045828e-02,1.048413e-02,1.0    
293     1.061960e-02,1.064973e-02,1.068119e-02,1.0    
294     1.084929e-02,1.088177e-02,1.091287e-02,1.0    
295     1.104790e-02,1.107250e-02,1.109694e-02,1.1    
296     1.121276e-02,1.123411e-02,1.125504e-02,1.1    
297     1.135576e-02,1.137542e-02,1.139497e-02,1.1    
298     1.149159e-02,1.151073e-02,1.152981e-02,1.1    
299     1.162429e-02,1.164299e-02,1.166162e-02,1.1    
300     1.175368e-02,1.177187e-02,1.178998e-02,1.1    
301     1.187939e-02,1.189704e-02,1.191463e-02,1.1    
302     1.200146e-02,1.201862e-02,1.203572e-02,1.2    
303     1.212025e-02,1.213697e-02,1.215364e-02,1.2    
304     1.223616e-02,1.225251e-02,1.226881e-02,1.2    
305     1.234962e-02,1.236566e-02,1.238165e-02,1.2    
306     1.246103e-02,1.247679e-02,1.249253e-02,1.2    
307     1.257072e-02,1.258627e-02,1.260179e-02,1.2    
308     1.267901e-02,1.269438e-02,1.270973e-02,1.2    
309     1.278618e-02,1.280141e-02,1.281663e-02,1.2    
310     1.289249e-02,1.290763e-02,1.292275e-02,1.2    
311     1.299819e-02,1.301326e-02,1.302831e-02,1.3    
312     1.310350e-02,1.311853e-02,1.313355e-02,1.3    
313     1.320863e-02,1.322365e-02,1.323866e-02,1.3    
314 static const G4double P11[nE]={                   
315     0.000000e+00,0.000000e+00,0.000000e+00,5.1    
316     6.042723e-04,7.780058e-04,9.681932e-04,1.1    
317     2.214927e-03,2.536789e-03,2.887800e-03,3.2    
318     5.161429e-03,5.738401e-03,6.354262e-03,7.0    
319     9.954245e-03,1.079211e-02,1.166922e-02,1.2    
320     1.631002e-02,1.728180e-02,1.825789e-02,1.9    
321     2.312152e-02,2.408180e-02,2.503409e-02,2.5    
322     2.965126e-02,3.054421e-02,3.142876e-02,3.2    
323     3.568998e-02,3.651243e-02,3.732532e-02,3.8    
324     4.124621e-02,4.200316e-02,4.275250e-02,4.3    
325     4.637734e-02,4.707965e-02,4.777482e-02,4.8    
326     5.115054e-02,5.180682e-02,5.245710e-02,5.3    
327     5.562592e-02,5.624410e-02,5.685739e-02,5.7    
328     5.985567e-02,6.044255e-02,6.102544e-02,6.1    
329     6.388395e-02,6.444520e-02,6.500318e-02,6.5    
330     6.774747e-02,6.828779e-02,6.882546e-02,6.9    
331     7.147659e-02,7.199987e-02,7.252099e-02,7.3    
332     7.509630e-02,7.560571e-02,7.611336e-02,7.6    
333     7.862701e-02,7.912515e-02,7.962186e-02,8.0    
334     8.208616e-02,8.257558e-02,8.306395e-02,8.3    
335     8.550172e-02,8.599961e-02,8.651664e-02,8.7    
336     9.125803e-02,9.328999e-02,9.579442e-02,9.8    
337     1.151048e-01,1.205775e-01,1.267685e-01,1.3    
338     1.729166e-01,1.861694e-01,2.010268e-01,2.1    
339     2.959020e-01,3.168904e-01,3.376171e-01,3.5    
340     4.294320e-01,4.450769e-01,4.599676e-01,4.7    
341     5.276464e-01,5.408286e-01,5.542584e-01,5.6    
342     6.315125e-01,6.502781e-01,6.704690e-01,6.9    
343     7.887092e-01,8.136834e-01,8.383295e-01,8.6    
344     9.552702e-01,9.786340e-01,1.002560e+00,1.0    
345     1.126566e+00,1.151612e+00,1.176902e+00,1.2    
346     1.310225e+00,1.338704e+00,1.367897e+00,1.3    
347     1.525899e+00,1.560130e+00,1.595301e+00,1.6    
348     1.786109e+00,1.827424e+00,1.869844e+00,1.9    
349     2.099449e+00,2.149066e+00,2.199987e+00,2.2    
350     2.475340e+00,2.534822e+00,2.595869e+00,2.6    
351     2.926136e+00,2.997533e+00,3.070832e+00,3.1    
352     3.467838e+00,3.553766e+00,3.642020e+00,3.7    
353     4.120669e+00,4.224403e+00,4.330993e+00,4.4    
354     4.909864e+00,5.035479e+00,5.164607e+00,5.2    
355     5.866778e+00,6.019333e+00,6.176218e+00,6.3    
356     7.030342e+00,7.216120e+00,7.407242e+00,7.6    
357     8.448929e+00,8.675743e+00,8.909163e+00,9.1    
358     1.018273e+01,1.046031e+01,1.074606e+01,1.1    
359     1.230674e+01,1.264721e+01,1.299783e+01,1.3    
360     1.491454e+01,1.533306e+01,1.576417e+01,1.6    
361     1.812306e+01,1.863856e+01,1.916973e+01,1.9    
362     2.207856e+01,2.271477e+01,2.337048e+01,2.4    
363 static const G4double P21[nE]={                   
364     0.000000e+00,0.000000e+00,0.000000e+00,8.9    
365     1.143952e-03,1.506169e-03,1.917732e-03,2.3    
366     4.952288e-03,5.818752e-03,6.797967e-03,7.9    
367     1.389665e-02,1.588243e-02,1.807616e-02,2.0    
368     3.234740e-02,3.600656e-02,3.996753e-02,4.4    
369     6.317608e-02,6.853047e-02,7.408791e-02,7.9    
370     1.046534e-01,1.112996e-01,1.181068e-01,1.2    
371     1.544640e-01,1.622007e-01,1.701127e-01,1.7    
372     2.120436e-01,2.209398e-01,2.300138e-01,2.3    
373     2.781207e-01,2.883204e-01,2.987380e-01,3.0    
374     3.540870e-01,3.658557e-01,3.778718e-01,3.9    
375     4.419246e-01,4.555823e-01,4.695390e-01,4.8    
376     5.441662e-01,5.601251e-01,5.764518e-01,5.9    
377     6.640062e-01,6.827845e-01,7.020148e-01,7.2    
378     8.054363e-01,8.276787e-01,8.504769e-01,8.7    
379     9.734160e-01,9.999218e-01,1.027113e+00,1.0    
380     1.174094e+00,1.205855e+00,1.238460e+00,1.2    
381     1.415083e+00,1.453324e+00,1.492606e+00,1.5    
382     1.705796e+00,1.752032e+00,1.799555e+00,1.8    
383     2.057951e+00,2.114105e+00,2.171862e+00,2.2    
384     2.487871e+00,2.558479e+00,2.634058e+00,2.7    
385     3.406090e+00,3.762360e+00,4.214883e+00,4.7    
386     8.066375e+00,9.251616e+00,1.063353e+01,1.2    
387     2.200875e+01,2.555457e+01,2.965149e+01,3.4    
388     5.841585e+01,6.535063e+01,7.240850e+01,7.9    
389     1.065618e+02,1.129463e+02,1.192092e+02,1.2    
390     1.503616e+02,1.570072e+02,1.639851e+02,1.7    
391     2.081256e+02,2.198120e+02,2.327710e+02,2.4    
392     3.160962e+02,3.353012e+02,3.548335e+02,3.7    
393     4.563033e+02,4.784941e+02,5.019153e+02,5.2    
394     6.349401e+02,6.643178e+02,6.948910e+02,7.2    
395     8.716998e+02,9.129543e+02,9.565398e+02,1.0    
396     1.215413e+03,1.276652e+03,1.341501e+03,1.4    
397     1.727580e+03,1.818857e+03,1.915447e+03,2.0    
398     2.489147e+03,2.624520e+03,2.767701e+03,2.9    
399     3.617331e+03,3.817744e+03,4.029722e+03,4.2    
400     5.288205e+03,5.585276e+03,5.899591e+03,6.2    
401     7.767793e+03,8.209307e+03,8.676651e+03,9.1    
402     1.145821e+04,1.211642e+04,1.281345e+04,1.3    
403     1.696773e+04,1.795200e+04,1.899477e+04,2.0    
404     2.521774e+04,2.669389e+04,2.825838e+04,2.9    
405     3.760630e+04,3.982616e+04,4.217976e+04,4.4    
406     5.625866e+04,5.960545e+04,6.315512e+04,6.6    
407     8.441148e+04,8.946939e+04,9.483567e+04,1.0    
408     1.270030e+05,1.346642e+05,1.427950e+05,1.5    
409     1.915807e+05,2.032100e+05,2.155559e+05,2.2    
410     2.896999e+05,3.073887e+05,3.261726e+05,3.4    
411     4.390797e+05,4.660376e+05,4.946723e+05,5.2    
412 static const G4double P02[nE]={                   
413     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
414     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
415     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
416     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
417     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
418     2.266269e-04,3.153714e-04,4.017279e-04,4.8    
419     7.941235e-04,8.648089e-04,9.328938e-04,9.9    
420     1.237327e-03,1.291871e-03,1.344747e-03,1.3    
421     1.586647e-03,1.631791e-03,1.676264e-03,1.7    
422     1.895742e-03,1.940822e-03,1.987146e-03,2.0    
423     2.250027e-03,2.312541e-03,2.379824e-03,2.4    
424     2.813598e-03,2.923687e-03,3.041745e-03,3.1    
425     3.705488e-03,3.840963e-03,3.973777e-03,4.1    
426     4.572513e-03,4.676977e-03,4.776189e-03,4.8    
427     5.197723e-03,5.268357e-03,5.334907e-03,5.3    
428     5.613724e-03,5.660086e-03,5.703745e-03,5.7    
429     5.886946e-03,5.917549e-03,5.946428e-03,5.9    
430     6.068552e-03,6.089154e-03,6.108666e-03,6.1    
431     6.192190e-03,6.206496e-03,6.220122e-03,6.2    
432     6.278332e-03,6.288730e-03,6.299041e-03,6.3    
433     6.351039e-03,6.362103e-03,6.373708e-03,6.3    
434     6.448411e-03,6.468623e-03,6.491167e-03,6.5    
435     6.644722e-03,6.684861e-03,6.728714e-03,6.7    
436     7.015759e-03,7.088847e-03,7.167329e-03,7.2    
437     7.613307e-03,7.703728e-03,7.791106e-03,7.8    
438     8.159311e-03,8.219026e-03,8.274939e-03,8.3    
439     8.515509e-03,8.559288e-03,8.602682e-03,8.6    
440     8.826098e-03,8.873546e-03,8.921888e-03,8.9    
441     9.163124e-03,9.207992e-03,9.250957e-03,9.2    
442     9.437672e-03,9.470422e-03,9.502096e-03,9.5    
443     9.650264e-03,9.678770e-03,9.707115e-03,9.7    
444     9.848065e-03,9.876273e-03,9.904506e-03,9.9    
445     1.004587e-02,1.007412e-02,1.010233e-02,1.0    
446     1.024251e-02,1.027031e-02,1.029800e-02,1.0    
447     1.043486e-02,1.046188e-02,1.048878e-02,1.0    
448     1.062144e-02,1.064761e-02,1.067366e-02,1.0    
449     1.080215e-02,1.082752e-02,1.085277e-02,1.0    
450     1.097749e-02,1.100214e-02,1.102670e-02,1.1    
451     1.114815e-02,1.117219e-02,1.119615e-02,1.1    
452     1.131482e-02,1.133834e-02,1.136179e-02,1.1    
453     1.147813e-02,1.150122e-02,1.152426e-02,1.1    
454     1.163867e-02,1.166140e-02,1.168409e-02,1.1    
455     1.179693e-02,1.181938e-02,1.184180e-02,1.1    
456     1.195338e-02,1.197561e-02,1.199781e-02,1.2    
457     1.210843e-02,1.213049e-02,1.215253e-02,1.2    
458     1.226246e-02,1.228440e-02,1.230633e-02,1.2    
459     1.241582e-02,1.243769e-02,1.245956e-02,1.2    
460     1.256884e-02,1.259069e-02,1.261254e-02,1.2    
461 static const G4double P12[nE]={                   
462     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
463     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
464     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
465     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
466     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
467     1.179082e-03,1.668007e-03,2.159647e-03,2.6    
468     4.620522e-03,5.109707e-03,5.596353e-03,6.0    
469     7.989358e-03,8.461919e-03,8.934839e-03,9.4    
470     1.131287e-02,1.180118e-02,1.229762e-02,1.2    
471     1.499297e-02,1.560045e-02,1.624454e-02,1.6    
472     2.027454e-02,2.132233e-02,2.248557e-02,2.3    
473     3.076594e-02,3.305745e-02,3.559175e-02,3.8    
474     5.125391e-02,5.475126e-02,5.828676e-02,6.1    
475     7.573356e-02,7.907571e-02,8.234849e-02,8.5    
476     9.755076e-02,1.003496e-01,1.030683e-01,1.0    
477     1.155154e-01,1.177894e-01,1.199970e-01,1.2    
478     1.301187e-01,1.319759e-01,1.337825e-01,1.3    
479     1.421304e-01,1.436768e-01,1.451865e-01,1.4    
480     1.522486e-01,1.535763e-01,1.548799e-01,1.5    
481     1.609624e-01,1.621554e-01,1.633749e-01,1.6    
482     1.701192e-01,1.716882e-01,1.733847e-01,1.7    
483     1.854234e-01,1.889661e-01,1.930388e-01,1.9    
484     2.236126e-01,2.323048e-01,2.420924e-01,2.5    
485     3.126625e-01,3.322153e-01,3.538544e-01,3.7    
486     4.887811e-01,5.186554e-01,5.484082e-01,5.7    
487     6.852353e-01,7.096034e-01,7.331188e-01,7.5    
488     8.437111e-01,8.657805e-01,8.883258e-01,9.1    
489     1.015650e+00,1.045194e+00,1.076217e+00,1.1    
490     1.245702e+00,1.280202e+00,1.314250e+00,1.3    
491     1.475916e+00,1.507016e+00,1.538017e+00,1.5    
492     1.696751e+00,1.730185e+00,1.764449e+00,1.7    
493     1.951130e+00,1.991990e+00,2.034139e+00,2.0    
494     2.265373e+00,2.315904e+00,2.367919e+00,2.4    
495     2.651032e+00,2.712435e+00,2.775494e+00,2.8    
496     3.116791e+00,3.190508e+00,3.266141e+00,3.3    
497     3.674707e+00,3.762868e+00,3.853312e+00,3.9    
498     4.342009e+00,4.447529e+00,4.555814e+00,4.6    
499     5.141606e+00,5.268252e+00,5.398277e+00,5.5    
500     6.102709e+00,6.255223e+00,6.411885e+00,6.5    
501     7.261893e+00,7.446187e+00,7.635583e+00,7.8    
502     8.664658e+00,8.888076e+00,9.117781e+00,9.3    
503     1.036753e+01,1.063919e+01,1.091862e+01,1.1    
504     1.244074e+01,1.277199e+01,1.311284e+01,1.3    
505     1.497170e+01,1.537667e+01,1.579352e+01,1.6    
506     1.806931e+01,1.856562e+01,1.907666e+01,1.9    
507     2.186947e+01,2.247911e+01,2.310704e+01,2.3    
508     2.654196e+01,2.729244e+01,2.806568e+01,2.8    
509     3.229921e+01,3.322497e+01,3.417906e+01,3.5    
510 static const G4double P22[nE]={                   
511     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
512     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
513     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
514     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
515     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
516     6.137882e-03,8.831785e-03,1.163101e-02,1.4    
517     2.709822e-02,3.048398e-02,3.396264e-02,3.7    
518     5.281322e-02,5.690771e-02,6.113789e-02,6.5    
519     8.456445e-02,8.984669e-02,9.538880e-02,1.0    
520     1.285578e-01,1.367447e-01,1.457007e-01,1.5    
521     2.076059e-01,2.251690e-01,2.452818e-01,2.6    
522     4.036571e-01,4.513583e-01,5.057654e-01,5.6    
523     8.760602e-01,9.663525e-01,1.060475e+00,1.1    
524     1.569856e+00,1.676792e+00,1.784762e+00,1.8    
525     2.334084e+00,2.444997e+00,2.556071e+00,2.6    
526     3.112811e+00,3.224358e+00,3.335994e+00,3.4    
527     3.896283e+00,4.009001e+00,4.122034e+00,4.2    
528     4.693754e+00,4.809831e+00,4.926647e+00,5.0    
529     5.524919e+00,5.648145e+00,5.772858e+00,5.8    
530     6.409661e+00,6.546552e+00,6.690783e+00,6.8    
531     7.567228e+00,7.789765e+00,8.037769e+00,8.3    
532     9.981511e+00,1.060252e+01,1.133833e+01,1.2    
533     1.743708e+01,1.931953e+01,2.150424e+01,2.4    
534     3.888631e+01,4.411756e+01,5.008434e+01,5.6    
535     9.098150e+01,1.008525e+02,1.109841e+02,1.2    
536     1.619260e+02,1.718706e+02,1.817614e+02,1.9    
537     2.326990e+02,2.438250e+02,2.555393e+02,2.6    
538     3.282386e+02,3.466366e+02,3.665465e+02,3.8    
539     4.858444e+02,5.123742e+02,5.393579e+02,5.6    
540     6.796009e+02,7.091381e+02,7.394822e+02,7.7    
541     9.098582e+02,9.490768e+02,9.904994e+02,1.0    
542     1.238214e+03,1.297406e+03,1.360335e+03,1.4    
543     1.739279e+03,1.829682e+03,1.925585e+03,2.0    
544     2.498462e+03,2.634120e+03,2.777703e+03,2.9    
545     3.630477e+03,3.831604e+03,4.044271e+03,4.2    
546     5.304913e+03,5.601955e+03,5.916011e+03,6.2    
547     7.778144e+03,8.217189e+03,8.681530e+03,9.1    
548     1.143804e+04,1.208877e+04,1.277731e+04,1.3    
549     1.687087e+04,1.783860e+04,1.886306e+04,1.9    
550     2.496305e+04,2.640711e+04,2.793656e+04,2.9    
551     3.705657e+04,3.921839e+04,4.150907e+04,4.3    
552     5.518666e+04,5.843278e+04,6.187379e+04,6.5    
553     8.244579e+04,8.733369e+04,9.251705e+04,9.8    
554     1.235418e+05,1.309211e+05,1.387492e+05,1.4    
555     1.856555e+05,1.968232e+05,2.086742e+05,2.2    
556     2.797587e+05,2.966988e+05,3.146809e+05,3.3    
557     4.226461e+05,4.483978e+05,4.757419e+05,5.0    
558     6.400668e+05,6.792939e+05,7.209586e+05,7.6    
559 static const G4double P03[nE]={                   
560     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
561     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
562     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
563     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
564     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
565     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
566     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
567     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
568     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
569     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
570     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
571     3.120096e-04,5.006502e-04,7.064650e-04,9.2    
572     1.899130e-03,2.146133e-03,2.388721e-03,2.6    
573     3.483317e-03,3.673901e-03,3.854657e-03,4.0    
574     4.618564e-03,4.745669e-03,4.865110e-03,4.9    
575     5.361367e-03,5.443076e-03,5.519761e-03,5.5    
576     5.838252e-03,5.890835e-03,5.940264e-03,5.9    
577     6.146889e-03,6.181306e-03,6.213770e-03,6.2    
578     6.351124e-03,6.374368e-03,6.396429e-03,6.4    
579     6.491149e-03,6.507570e-03,6.523412e-03,6.5    
580     6.597557e-03,6.612423e-03,6.627837e-03,6.6    
581     6.726353e-03,6.753276e-03,6.783463e-03,6.8    
582     6.991787e-03,7.046648e-03,7.106649e-03,7.1    
583     7.498466e-03,7.597591e-03,7.703710e-03,7.8    
584     8.304010e-03,8.426239e-03,8.544896e-03,8.6    
585     9.053136e-03,9.136896e-03,9.215579e-03,9.2    
586     9.554447e-03,9.615520e-03,9.675636e-03,9.7    
587     9.976259e-03,1.003820e-02,1.010084e-02,1.0    
588     1.041285e-02,1.047195e-02,1.052910e-02,1.0    
589     1.078376e-02,1.082922e-02,1.087331e-02,1.0    
590     1.107975e-02,1.111933e-02,1.115860e-02,1.1    
591     1.135262e-02,1.139120e-02,1.142974e-02,1.1    
592     1.162169e-02,1.165987e-02,1.169795e-02,1.1    
593     1.188649e-02,1.192375e-02,1.196086e-02,1.1    
594     1.214372e-02,1.217975e-02,1.221559e-02,1.2    
595     1.239202e-02,1.242676e-02,1.246133e-02,1.2    
596     1.263159e-02,1.266515e-02,1.269855e-02,1.2    
597     1.286333e-02,1.289586e-02,1.292826e-02,1.2    
598     1.308833e-02,1.311999e-02,1.315153e-02,1.3    
599     1.330763e-02,1.333855e-02,1.336937e-02,1.3    
600     1.352215e-02,1.355246e-02,1.358269e-02,1.3    
601     1.373274e-02,1.376254e-02,1.379228e-02,1.3    
602     1.394010e-02,1.396950e-02,1.399885e-02,1.4    
603     1.414489e-02,1.417397e-02,1.420301e-02,1.4    
604     1.434768e-02,1.437652e-02,1.440533e-02,1.4    
605     1.454901e-02,1.457768e-02,1.460633e-02,1.4    
606     1.474934e-02,1.477791e-02,1.480646e-02,1.4    
607     1.494914e-02,1.497766e-02,1.500618e-02,1.5    
608 static const G4double P13[nE]={                   
609     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
610     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
611     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
612     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
613     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
614     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
615     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
616     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
617     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
618     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
619     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
620     6.214464e-03,1.014119e-02,1.455950e-02,1.9    
621     4.271502e-02,4.909153e-02,5.554923e-02,6.2    
622     8.744455e-02,9.354199e-02,9.950466e-02,1.0    
623     1.270507e-01,1.320872e-01,1.369666e-01,1.4    
624     1.591164e-01,1.631241e-01,1.670017e-01,1.7    
625     1.845939e-01,1.877850e-01,1.908772e-01,1.9    
626     2.049976e-01,2.075809e-01,2.100927e-01,2.1    
627     2.217034e-01,2.238606e-01,2.259710e-01,2.2    
628     2.358717e-01,2.377555e-01,2.396290e-01,2.4    
629     2.492319e-01,2.513402e-01,2.535933e-01,2.5    
630     2.694712e-01,2.741902e-01,2.796437e-01,2.8    
631     3.211310e-01,3.330113e-01,3.464030e-01,3.6    
632     4.427220e-01,4.692400e-01,4.984989e-01,5.3    
633     6.801126e-01,7.204962e-01,7.608999e-01,8.0    
634     9.498326e-01,9.840127e-01,1.017105e+00,1.0    
635     1.172873e+00,1.203659e+00,1.234892e+00,1.2    
636     1.406103e+00,1.444672e+00,1.484870e+00,1.5    
637     1.704122e+00,1.749563e+00,1.794852e+00,1.8    
638     2.015447e+00,2.058627e+00,2.101780e+00,2.1    
639     2.322931e+00,2.369348e+00,2.416823e+00,2.4    
640     2.673759e+00,2.729644e+00,2.787177e+00,2.8    
641     3.101119e+00,3.169415e+00,3.239621e+00,3.3    
642     3.620375e+00,3.702702e+00,3.787172e+00,3.8    
643     4.243205e+00,4.341489e+00,4.442258e+00,4.5    
644     4.985593e+00,5.102639e+00,5.222653e+00,5.3    
645     5.870177e+00,6.009805e+00,6.153034e+00,6.2    
646     6.926941e+00,7.094076e+00,7.265614e+00,7.4    
647     8.194045e+00,8.394879e+00,8.601118e+00,8.8    
648     9.719220e+00,9.961464e+00,1.021036e+01,1.0    
649     1.156181e+01,1.185504e+01,1.215646e+01,1.2    
650     1.379550e+01,1.415161e+01,1.451783e+01,1.4    
651     1.651185e+01,1.694562e+01,1.739189e+01,1.7    
652     1.982475e+01,2.035460e+01,2.089992e+01,2.1    
653     2.387613e+01,2.452501e+01,2.519307e+01,2.5    
654     2.884305e+01,2.963961e+01,3.046001e+01,3.1    
655     3.494673e+01,3.592682e+01,3.693656e+01,3.7    
656     4.246392e+01,4.367240e+01,4.491780e+01,4.6    
657 static const G4double P23[nE]={                   
658     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
659     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
660     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
661     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
662     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
663     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
664     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
665     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
666     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
667     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
668     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
669     1.238128e-01,2.055575e-01,3.004149e-01,4.0    
670     9.663775e-01,1.131003e+00,1.302922e+00,1.4    
671     2.234130e+00,2.429223e+00,2.625932e+00,2.8    
672     3.621155e+00,3.820738e+00,4.020088e+00,4.2    
673     5.010628e+00,5.207217e+00,5.403298e+00,5.5    
674     6.376896e+00,6.570570e+00,6.764031e+00,6.9    
675     7.730848e+00,7.924762e+00,8.119112e+00,8.3    
676     9.102464e+00,9.302681e+00,9.504585e+00,9.7    
677     1.054144e+01,1.075759e+01,1.097915e+01,1.1    
678     1.222528e+01,1.252429e+01,1.285366e+01,1.3    
679     1.541751e+01,1.624472e+01,1.723002e+01,1.8    
680     2.550749e+01,2.808038e+01,3.106957e+01,3.4    
681     5.479148e+01,6.188616e+01,6.995395e+01,7.9    
682     1.250017e+02,1.383453e+02,1.521040e+02,1.6    
683     2.224704e+02,2.364195e+02,2.503384e+02,2.6    
684     3.220770e+02,3.375972e+02,3.538252e+02,3.7    
685     4.515193e+02,4.755366e+02,5.013340e+02,5.2    
686     6.556933e+02,6.906377e+02,7.265303e+02,7.6    
687     9.179815e+02,9.589917e+02,1.001230e+03,1.0    
688     1.238592e+03,1.293040e+03,1.350433e+03,1.4    
689     1.691327e+03,1.772284e+03,1.858181e+03,1.9    
690     2.372614e+03,2.494796e+03,2.624240e+03,2.7    
691     3.394638e+03,3.576525e+03,3.768859e+03,3.9    
692     4.908260e+03,5.176412e+03,5.459761e+03,5.7    
693     7.136172e+03,7.530536e+03,7.947272e+03,8.3    
694     1.041452e+04,1.099549e+04,1.160967e+04,1.2    
695     1.525129e+04,1.611006e+04,1.701843e+04,1.7    
696     2.241358e+04,2.368790e+04,2.503657e+04,2.6    
697     3.306037e+04,3.495851e+04,3.696842e+04,3.9    
698     4.894528e+04,5.178262e+04,5.478849e+04,5.7    
699     7.272641e+04,7.698154e+04,8.149143e+04,8.6    
700     1.084410e+05,1.148416e+05,1.216281e+05,1.2    
701     1.622330e+05,1.718876e+05,1.821283e+05,1.9    
702     2.434706e+05,2.580712e+05,2.735637e+05,2.9    
703     3.664650e+05,3.885990e+05,4.120927e+05,4.3    
704     5.531172e+05,5.867475e+05,6.224550e+05,6.6    
705     8.369994e+05,8.882063e+05,9.425923e+05,1.0    
706 static const G4double P04[nE]={                   
707     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
708     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
709     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
710     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
711     1.519363e-05,2.703219e-05,4.014211e-05,5.5    
712     1.367177e-04,1.629335e-04,1.925910e-04,2.2    
713     3.991538e-04,4.544632e-04,5.151945e-04,5.8    
714     9.055657e-04,1.001169e-03,1.101697e-03,1.2    
715     1.665438e-03,1.784771e-03,1.904473e-03,2.0    
716     2.491853e-03,2.605055e-03,2.717349e-03,2.8    
717     3.290343e-03,3.417114e-03,3.552511e-03,3.6    
718     4.435314e-03,4.658182e-03,4.892305e-03,5.1    
719     6.043927e-03,6.241290e-03,6.422941e-03,6.5    
720     7.110833e-03,7.211089e-03,7.301407e-03,7.3    
721     7.635730e-03,7.684395e-03,7.728406e-03,7.7    
722     7.893992e-03,7.918663e-03,7.941175e-03,7.9    
723     8.028333e-03,8.041823e-03,8.054305e-03,8.0    
724     8.104959e-03,8.113293e-03,8.121181e-03,8.1    
725     8.155681e-03,8.161895e-03,8.167971e-03,8.1    
726     8.196173e-03,8.201721e-03,8.207490e-03,8.2    
727     8.243208e-03,8.253043e-03,8.264530e-03,8.2    
728     8.366621e-03,8.400357e-03,8.439779e-03,8.4    
729     8.733284e-03,8.813450e-03,8.901722e-03,8.9    
730     9.483175e-03,9.630255e-03,9.787493e-03,9.9    
731     1.067675e-02,1.085922e-02,1.103733e-02,1.1    
732     1.181358e-02,1.194350e-02,1.206587e-02,1.2    
733     1.259247e-02,1.268636e-02,1.277816e-02,1.2    
734     1.322486e-02,1.331430e-02,1.340408e-02,1.3    
735     1.384946e-02,1.393508e-02,1.401867e-02,1.4    
736     1.440114e-02,1.447087e-02,1.453873e-02,1.4    
737     1.485754e-02,1.491855e-02,1.497899e-02,1.5    
738     1.527579e-02,1.533444e-02,1.539291e-02,1.5    
739     1.568263e-02,1.573999e-02,1.579713e-02,1.5    
740     1.607900e-02,1.613455e-02,1.618979e-02,1.6    
741     1.646145e-02,1.651485e-02,1.656795e-02,1.6    
742     1.682885e-02,1.688015e-02,1.693115e-02,1.6    
743     1.718208e-02,1.723147e-02,1.728062e-02,1.7    
744     1.752281e-02,1.757058e-02,1.761813e-02,1.7    
745     1.785289e-02,1.789928e-02,1.794549e-02,1.7    
746     1.817402e-02,1.821925e-02,1.826434e-02,1.8    
747     1.848769e-02,1.853197e-02,1.857613e-02,1.8    
748     1.879521e-02,1.883871e-02,1.888210e-02,1.8    
749     1.909771e-02,1.914058e-02,1.918336e-02,1.9    
750     1.939619e-02,1.943856e-02,1.948086e-02,1.9    
751     1.969155e-02,1.973354e-02,1.977548e-02,1.9    
752     1.998458e-02,2.002629e-02,2.006798e-02,2.0    
753     2.027601e-02,2.031756e-02,2.035908e-02,2.0    
754     2.056653e-02,2.060799e-02,2.064945e-02,2.0    
755 static const G4double P14[nE]={                   
756     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
757     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
758     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
759     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
760     6.371064e-05,1.154186e-04,1.746380e-04,2.4    
761     6.607827e-04,8.052700e-04,9.741837e-04,1.1    
762     2.281454e-03,2.664366e-03,3.098601e-03,3.5    
763     6.190306e-03,7.018768e-03,7.918090e-03,8.8    
764     1.348004e-02,1.477087e-02,1.610715e-02,1.7    
765     2.331419e-02,2.483956e-02,2.640072e-02,2.8    
766     3.516681e-02,3.729155e-02,3.963237e-02,4.2    
767     5.648894e-02,6.112768e-02,6.615316e-02,7.1    
768     9.323627e-02,9.833048e-02,1.031653e-01,1.0    
769     1.231347e-01,1.263418e-01,1.293208e-01,1.3    
770     1.413425e-01,1.432706e-01,1.450684e-01,1.4    
771     1.524436e-01,1.536536e-01,1.547918e-01,1.5    
772     1.595992e-01,1.604178e-01,1.611987e-01,1.6    
773     1.646587e-01,1.652843e-01,1.658945e-01,1.6    
774     1.688138e-01,1.693906e-01,1.699718e-01,1.7    
775     1.729235e-01,1.735601e-01,1.742424e-01,1.7    
776     1.788979e-01,1.802932e-01,1.819727e-01,1.8    
777     1.985139e-01,2.044278e-01,2.115504e-01,2.2    
778     2.700636e-01,2.874240e-01,3.071261e-01,3.2    
779     4.500720e-01,4.894190e-01,5.327720e-01,5.8    
780     8.018245e-01,8.621139e-01,9.227611e-01,9.8    
781     1.211429e+00,1.264444e+00,1.315911e+00,1.3    
782     1.557947e+00,1.605276e+00,1.652968e+00,1.7    
783     1.907216e+00,1.962904e+00,2.020515e+00,2.0    
784     2.333525e+00,2.399370e+00,2.465608e+00,2.5    
785     2.797093e+00,2.863310e+00,2.929727e+00,2.9    
786     3.271270e+00,3.342824e+00,3.415884e+00,3.4    
787     3.808870e+00,3.893828e+00,3.981122e+00,4.0    
788     4.454920e+00,4.557525e+00,4.662858e+00,4.7    
789     5.232077e+00,5.354782e+00,5.480568e+00,5.6    
790     6.157992e+00,6.303681e+00,6.452959e+00,6.6    
791     7.256411e+00,7.429221e+00,7.606326e+00,7.7    
792     8.560577e+00,8.766096e+00,8.976834e+00,9.1    
793     1.011428e+01,1.035969e+01,1.061148e+01,1.0    
794     1.197310e+01,1.226741e+01,1.256956e+01,1.2    
795     1.420644e+01,1.456086e+01,1.492492e+01,1.5    
796     1.690054e+01,1.732896e+01,1.776928e+01,1.8    
797     2.016239e+01,2.068210e+01,2.121649e+01,2.1    
798     2.412492e+01,2.475737e+01,2.540797e+01,2.6    
799     2.895345e+01,2.972536e+01,3.051974e+01,3.1    
800     3.485392e+01,3.579859e+01,3.677113e+01,3.7    
801     4.208324e+01,4.324227e+01,4.443590e+01,4.5    
802     5.096240e+01,5.238778e+01,5.385619e+01,5.5    
803     6.189289e+01,6.364971e+01,6.546010e+01,6.7    
804 static const G4double P24[nE]={                   
805     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
806     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
807     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
808     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
809     2.672072e-04,4.930762e-04,7.606037e-04,1.0    
810     3.213192e-03,4.009598e-03,4.971730e-03,6.1    
811     1.326287e-02,1.591403e-02,1.901913e-02,2.2    
812     4.355667e-02,5.073638e-02,5.878240e-02,6.7    
813     1.137711e-01,1.277352e-01,1.426538e-01,1.5    
814     2.312639e-01,2.518194e-01,2.735251e-01,2.9    
815     4.079082e-01,4.435227e-01,4.839951e-01,5.3    
816     8.064902e-01,9.030479e-01,1.010929e+00,1.1    
817     1.649089e+00,1.780588e+00,1.909280e+00,2.0    
818     2.490116e+00,2.592716e+00,2.690981e+00,2.7    
819     3.124097e+00,3.200497e+00,3.273938e+00,3.3    
820     3.603063e+00,3.662412e+00,3.719967e+00,3.7    
821     3.985641e+00,4.035323e+00,4.084172e+00,4.1    
822     4.320973e+00,4.367929e+00,4.415151e+00,4.4    
823     4.662638e+00,4.716179e+00,4.771788e+00,4.8    
824     5.081309e+00,5.154369e+00,5.235071e+00,5.3    
825     5.843034e+00,6.040974e+00,6.286575e+00,6.5    
826     8.971439e+00,1.000824e+01,1.129521e+01,1.2    
827     2.298186e+01,2.674168e+01,3.113945e+01,3.6    
828     6.634717e+01,7.687414e+01,8.882816e+01,1.0    
829     1.703848e+02,1.903058e+02,2.109585e+02,2.3    
830     3.185090e+02,3.401451e+02,3.617925e+02,3.8    
831     4.732483e+02,4.971083e+02,5.218876e+02,5.4    
832     6.668734e+02,7.015497e+02,7.385227e+02,7.7    
833     9.589179e+02,1.009554e+03,1.062052e+03,1.1    
834     1.349889e+03,1.412779e+03,1.477789e+03,1.5    
835     1.844372e+03,1.928304e+03,2.016627e+03,2.1    
836     2.537958e+03,2.661032e+03,2.791360e+03,2.9    
837     3.567658e+03,3.751220e+03,3.945429e+03,4.1    
838     5.097080e+03,5.368173e+03,5.654580e+03,5.9    
839     7.347046e+03,7.744537e+03,8.164282e+03,8.6    
840     1.064318e+04,1.122542e+04,1.184040e+04,1.2    
841     1.547628e+04,1.633140e+04,1.723507e+04,1.8    
842     2.258721e+04,2.384817e+04,2.518152e+04,2.6    
843     3.309379e+04,3.496121e+04,3.693706e+04,3.9    
844     4.868365e+04,5.146069e+04,5.440064e+04,5.7    
845     7.190880e+04,7.605431e+04,8.044530e+04,8.5    
846     1.066357e+05,1.128456e+05,1.194264e+05,1.2    
847     1.587342e+05,1.680664e+05,1.779602e+05,1.8    
848     2.371343e+05,2.511997e+05,2.661175e+05,2.8    
849     3.554481e+05,3.767048e+05,3.992580e+05,4.2    
850     5.344641e+05,5.666699e+05,6.008519e+05,6.3    
851     8.059888e+05,8.548986e+05,9.068263e+05,9.6    
852     1.218769e+06,1.293211e+06,1.372270e+06,1.4    
853 static const G4double P05[nE]={                   
854     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
855     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
856     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
857     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
858     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
859     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
860     5.861537e-07,3.061765e-06,5.963955e-06,9.3    
861     3.261650e-05,4.273714e-05,5.639991e-05,7.3    
862     2.302289e-04,3.053102e-04,4.003118e-04,5.1    
863     1.049782e-03,1.203029e-03,1.370490e-03,1.5    
864     2.457122e-03,2.712762e-03,2.971183e-03,3.2    
865     4.197027e-03,4.422883e-03,4.643939e-03,4.8    
866     5.727123e-03,5.946622e-03,6.168285e-03,6.3    
867     7.295793e-03,7.519457e-03,7.740099e-03,7.9    
868     8.765403e-03,8.950129e-03,9.127152e-03,9.2    
869     9.895275e-03,1.002632e-02,1.015036e-02,1.0    
870     1.067556e-02,1.076370e-02,1.084690e-02,1.0    
871     1.119809e-02,1.125710e-02,1.131288e-02,1.1    
872     1.155017e-02,1.159052e-02,1.162888e-02,1.1    
873     1.179510e-02,1.182421e-02,1.185232e-02,1.1    
874     1.198450e-02,1.201129e-02,1.203924e-02,1.2    
875     1.222166e-02,1.227215e-02,1.232884e-02,1.2    
876     1.271890e-02,1.282083e-02,1.293183e-02,1.3    
877     1.364395e-02,1.382089e-02,1.400920e-02,1.4    
878     1.506624e-02,1.528290e-02,1.549464e-02,1.5    
879     1.642228e-02,1.657827e-02,1.672530e-02,1.6    
880     1.735724e-02,1.746933e-02,1.757862e-02,1.7    
881     1.810471e-02,1.820886e-02,1.831308e-02,1.8    
882     1.882813e-02,1.892735e-02,1.902442e-02,1.9    
883     1.947187e-02,1.955394e-02,1.963391e-02,1.9    
884     2.001002e-02,2.008193e-02,2.015313e-02,2.0    
885     2.050193e-02,2.057069e-02,2.063918e-02,2.0    
886     2.097784e-02,2.104476e-02,2.111138e-02,2.1    
887     2.143959e-02,2.150419e-02,2.156841e-02,2.1    
888     2.188391e-02,2.194588e-02,2.200747e-02,2.2    
889     2.230994e-02,2.236936e-02,2.242845e-02,2.2    
890     2.271896e-02,2.277613e-02,2.283300e-02,2.2    
891     2.311312e-02,2.316835e-02,2.322333e-02,2.3    
892     2.349466e-02,2.354825e-02,2.360164e-02,2.3    
893     2.386561e-02,2.391784e-02,2.396991e-02,2.4    
894     2.422776e-02,2.427888e-02,2.432984e-02,2.4    
895     2.458267e-02,2.463286e-02,2.468294e-02,2.4    
896     2.493168e-02,2.498112e-02,2.503047e-02,2.5    
897     2.527595e-02,2.532480e-02,2.537359e-02,2.5    
898     2.561653e-02,2.566494e-02,2.571329e-02,2.5    
899     2.595436e-02,2.600245e-02,2.605050e-02,2.6    
900     2.629030e-02,2.633818e-02,2.638604e-02,2.6    
901     2.662513e-02,2.667291e-02,2.672069e-02,2.6    
902 static const G4double P15[nE]={                   
903     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
904     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
905     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
906     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
907     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
908     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
909     3.974782e-06,2.111703e-05,4.187126e-05,6.6    
910     2.547110e-04,3.424807e-04,4.647771e-04,6.2    
911     2.209563e-03,3.022272e-03,4.083389e-03,5.3    
912     1.208961e-02,1.415500e-02,1.648372e-02,1.9    
913     3.316793e-02,3.745194e-02,4.191903e-02,4.6    
914     6.515418e-02,6.985421e-02,7.459850e-02,7.9    
915     1.001215e-01,1.057883e-01,1.116895e-01,1.1    
916     1.446344e-01,1.517909e-01,1.590700e-01,1.6    
917     1.961101e-01,2.034302e-01,2.106622e-01,2.1    
918     2.449747e-01,2.514023e-01,2.576744e-01,2.6    
919     2.866954e-01,2.920447e-01,2.972497e-01,3.0    
920     3.212540e-01,3.256827e-01,3.299983e-01,3.3    
921     3.500593e-01,3.538045e-01,3.574734e-01,3.6    
922     3.748540e-01,3.781939e-01,3.815172e-01,3.8    
923     3.986393e-01,4.024394e-01,4.065252e-01,4.1    
924     4.359366e-01,4.447861e-01,4.550284e-01,4.6    
925     5.326983e-01,5.547703e-01,5.795448e-01,6.0    
926     7.545360e-01,8.018700e-01,8.537886e-01,9.1    
927     1.173580e+00,1.245163e+00,1.317264e+00,1.3    
928     1.662267e+00,1.725924e+00,1.787763e+00,1.8    
929     2.078191e+00,2.134696e+00,2.191472e+00,2.2    
930     2.490835e+00,2.555683e+00,2.622564e+00,2.6    
931     2.984527e+00,3.060829e+00,3.137758e+00,3.2    
932     3.525611e+00,3.603556e+00,3.681827e+00,3.7    
933     4.084743e+00,4.169087e+00,4.255151e+00,4.3    
934     4.716965e+00,4.816561e+00,4.918816e+00,5.0    
935     5.472618e+00,5.592327e+00,5.715152e+00,5.8    
936     6.377929e+00,6.520625e+00,6.666852e+00,6.8    
937     7.453583e+00,7.622638e+00,7.795812e+00,7.9    
938     8.727236e+00,8.927450e+00,9.132601e+00,9.3    
939     1.023740e+01,1.047524e+01,1.071908e+01,1.0    
940     1.203469e+01,1.231843e+01,1.260954e+01,1.2    
941     1.418323e+01,1.452328e+01,1.487236e+01,1.5    
942     1.676306e+01,1.717234e+01,1.759273e+01,1.8    
943     1.987358e+01,2.036810e+01,2.087633e+01,2.1    
944     2.363808e+01,2.423775e+01,2.485434e+01,2.5    
945     2.820970e+01,2.893925e+01,2.968970e+01,3.0    
946     3.377889e+01,3.466908e+01,3.558515e+01,3.6    
947     4.058284e+01,4.167204e+01,4.279334e+01,4.3    
948     4.891749e+01,5.025360e+01,5.162956e+01,5.3    
949     5.915250e+01,6.079540e+01,6.248786e+01,6.4    
950     7.175028e+01,7.377494e+01,7.586129e+01,7.8    
951 static const G4double P25[nE]={                   
952     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
953     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
954     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
955     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
956     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
957     0.000000e+00,0.000000e+00,0.000000e+00,0.0    
958     2.695366e-05,1.456646e-04,2.940957e-04,4.7    
959     1.997124e-03,2.758357e-03,3.853132e-03,5.3    
960     2.139761e-02,3.019542e-02,4.204851e-02,5.7    
961     1.409425e-01,1.687811e-01,2.011671e-01,2.3    
962     4.578416e-01,5.296390e-01,6.068635e-01,6.8    
963     1.048148e+00,1.145963e+00,1.247793e+00,1.3    
964     1.850381e+00,1.996693e+00,2.153810e+00,2.3    
965     3.118309e+00,3.347312e+00,3.587475e+00,3.8    
966     4.928176e+00,5.218268e+00,5.513750e+00,5.8    
967     7.049470e+00,7.364768e+00,7.681941e+00,8.0    
968     9.288638e+00,9.613299e+00,9.938942e+00,1.0    
969     1.158282e+01,1.191526e+01,1.224919e+01,1.2    
970     1.394841e+01,1.429602e+01,1.464702e+01,1.5    
971     1.646782e+01,1.685103e+01,1.724405e+01,1.7    
972     1.946627e+01,2.000522e+01,2.060253e+01,2.1    
973     2.535341e+01,2.690469e+01,2.875518e+01,3.0    
974     4.424985e+01,4.902993e+01,5.455981e+01,6.0    
975     9.764114e+01,1.103048e+02,1.246205e+02,1.4    
976     2.215491e+02,2.452021e+02,2.697554e+02,2.9    
977     3.983095e+02,4.242881e+02,4.502985e+02,4.7    
978     5.840293e+02,6.125149e+02,6.420137e+02,6.7    
979     8.126858e+02,8.530658e+02,8.959874e+02,9.4    
980     1.150848e+03,1.209525e+03,1.270497e+03,1.3    
981     1.607330e+03,1.681358e+03,1.757972e+03,1.8    
982     2.190427e+03,2.289361e+03,2.393404e+03,2.5    
983     3.006011e+03,3.150290e+03,3.302954e+03,3.4    
984     4.210296e+03,4.424458e+03,4.650917e+03,4.8    
985     5.991817e+03,6.307078e+03,6.640026e+03,6.9    
986     8.605548e+03,9.066788e+03,9.553727e+03,1.0    
987     1.242741e+04,1.310199e+04,1.381436e+04,1.4    
988     1.802382e+04,1.901340e+04,2.005902e+04,2.1    
989     2.624943e+04,2.770738e+04,2.924888e+04,3.0    
990     3.839344e+04,4.055112e+04,4.283390e+04,4.5    
991     5.640180e+04,5.960874e+04,6.300357e+04,6.6    
992     8.321661e+04,8.800174e+04,9.306993e+04,9.8    
993     1.232946e+05,1.304601e+05,1.380532e+05,1.4    
994     1.834010e+05,1.941660e+05,2.055782e+05,2.1    
995     2.738267e+05,2.900474e+05,3.072504e+05,3.2    
996     4.102564e+05,4.347651e+05,4.607680e+05,4.8    
997     6.166422e+05,6.537685e+05,6.931719e+05,7.3    
998     9.296281e+05,9.860018e+05,1.045853e+06,1.1    
999     1.405371e+06,1.491162e+06,1.582272e+06,1.6    
1000 static const G4double P06[nE]={                  
1001     1.177941e-08,2.588247e-08,4.232501e-08,6.    
1002     1.622136e-07,1.942754e-07,2.292262e-07,2.    
1003     4.568609e-07,5.221504e-07,5.966563e-07,6.    
1004     1.117610e-06,1.253648e-06,1.407254e-06,1.    
1005     2.555265e-06,2.866400e-06,3.210941e-06,3.    
1006     5.771395e-06,6.471521e-06,7.274318e-06,8.    
1007     1.307801e-05,1.472581e-05,1.657663e-05,1.    
1008     2.993756e-05,3.371561e-05,3.801233e-05,4.    
1009     6.927898e-05,7.822982e-05,8.836866e-05,9.    
1010     1.638882e-04,1.858924e-04,2.113879e-04,2.    
1011     4.199591e-04,4.903319e-04,5.776505e-04,6.    
1012     1.545358e-03,1.876495e-03,2.234266e-03,2.    
1013     4.049256e-03,4.410120e-03,4.770247e-03,5.    
1014     6.493445e-03,6.811580e-03,7.118638e-03,7.    
1015     8.475244e-03,8.710011e-03,8.932790e-03,9.    
1016     9.878136e-03,1.003638e-02,1.018543e-02,1.    
1017     1.080940e-02,1.091309e-02,1.101074e-02,1.    
1018     1.142087e-02,1.148956e-02,1.155449e-02,1.    
1019     1.183119e-02,1.187846e-02,1.192350e-02,1.    
1020     1.211950e-02,1.215412e-02,1.218777e-02,1.    
1021     1.235014e-02,1.238405e-02,1.241978e-02,1.    
1022     1.265887e-02,1.272581e-02,1.280106e-02,1.    
1023     1.331759e-02,1.345191e-02,1.359785e-02,1.    
1024     1.452607e-02,1.475493e-02,1.499797e-02,1.    
1025     1.636002e-02,1.664033e-02,1.691510e-02,1.    
1026     1.812995e-02,1.833578e-02,1.853000e-02,1.    
1027     1.936343e-02,1.951030e-02,1.965297e-02,1.    
1028     2.033054e-02,2.046278e-02,2.059460e-02,2.    
1029     2.124276e-02,2.136797e-02,2.149083e-02,2.    
1030     2.206318e-02,2.216917e-02,2.227266e-02,2.    
1031     2.276069e-02,2.285400e-02,2.294633e-02,2.    
1032     2.339744e-02,2.348610e-02,2.357434e-02,2.    
1033     2.400939e-02,2.409515e-02,2.418046e-02,2.    
1034     2.459996e-02,2.468237e-02,2.476428e-02,2.    
1035     2.516612e-02,2.524495e-02,2.532328e-02,2.    
1036     2.570757e-02,2.578302e-02,2.585800e-02,2.    
1037     2.622645e-02,2.629891e-02,2.637098e-02,2.    
1038     2.672577e-02,2.679569e-02,2.686527e-02,2.    
1039     2.720856e-02,2.727634e-02,2.734385e-02,2.    
1040     2.767754e-02,2.774355e-02,2.780934e-02,2.    
1041     2.813508e-02,2.819963e-02,2.826399e-02,2.    
1042     2.858319e-02,2.864654e-02,2.870974e-02,2.    
1043     2.902362e-02,2.908601e-02,2.914827e-02,2.    
1044     2.945791e-02,2.951952e-02,2.958105e-02,2.    
1045     2.988738e-02,2.994842e-02,3.000938e-02,3.    
1046     3.031327e-02,3.037388e-02,3.043444e-02,3.    
1047     3.073665e-02,3.079699e-02,3.085730e-02,3.    
1048     3.115854e-02,3.121874e-02,3.127894e-02,3.    
1049 static const G4double P16[nE]={                  
1050     1.804258e-08,4.051712e-08,6.777154e-08,1.    
1051     2.930569e-07,3.599020e-07,4.355335e-07,5.    
1052     9.900430e-07,1.165832e-06,1.373705e-06,1.    
1053     3.002083e-06,3.470293e-06,4.017609e-06,4.    
1054     8.588820e-06,9.947777e-06,1.150413e-05,1.    
1055     2.439490e-05,2.825371e-05,3.282620e-05,3.    
1056     6.958951e-05,8.099809e-05,9.423242e-05,1.    
1057     2.003571e-04,2.331035e-04,2.715493e-04,3.    
1058     5.819333e-04,6.787875e-04,7.920050e-04,9.    
1059     1.727161e-03,2.023776e-03,2.378366e-03,2.    
1060     5.599801e-03,6.779798e-03,8.289933e-03,1.    
1061     2.693730e-02,3.383013e-02,4.151002e-02,4.    
1062     8.426714e-02,9.358336e-02,1.031704e-01,1.    
1063     1.534460e-01,1.636247e-01,1.737542e-01,1.    
1064     2.227163e-01,2.320190e-01,2.411202e-01,2.    
1065     2.833311e-01,2.910927e-01,2.986297e-01,3.    
1066     3.331003e-01,3.393933e-01,3.455017e-01,3.    
1067     3.735326e-01,3.786887e-01,3.837123e-01,3.    
1068     4.071056e-01,4.114928e-01,4.158014e-01,4.    
1069     4.362959e-01,4.402669e-01,4.442464e-01,4.    
1070     4.652924e-01,4.701007e-01,4.753240e-01,4.    
1071     5.138871e-01,5.256210e-01,5.392151e-01,5.    
1072     6.420599e-01,6.711473e-01,7.037187e-01,7.    
1073     9.317722e-01,9.929958e-01,1.060005e+00,1.    
1074     1.472082e+00,1.564697e+00,1.658261e+00,1.    
1075     2.110166e+00,2.194160e+00,2.275844e+00,2.    
1076     2.658844e+00,2.732878e+00,2.806997e+00,2.    
1077     3.192444e+00,3.274781e+00,3.359367e+00,3.    
1078     3.814868e+00,3.911152e+00,4.008524e+00,4.    
1079     4.504750e+00,4.605408e+00,4.706697e+00,4.    
1080     5.229531e+00,5.338975e+00,5.450581e+00,5.    
1081     6.047825e+00,6.176250e+00,6.307973e+00,6.    
1082     7.019380e+00,7.172786e+00,7.330069e+00,7.    
1083     8.177134e+00,8.359207e+00,8.545692e+00,8.    
1084     9.547688e+00,9.762753e+00,9.982984e+00,1.    
1085     1.116638e+01,1.142055e+01,1.168091e+01,1.    
1086     1.308209e+01,1.338354e+01,1.369254e+01,1.    
1087     1.535881e+01,1.571802e+01,1.608647e+01,1.    
1088     1.807750e+01,1.850756e+01,1.894900e+01,1.    
1089     2.133906e+01,2.185627e+01,2.238747e+01,2.    
1090     2.526871e+01,2.589326e+01,2.653505e+01,2.    
1091     3.002178e+01,3.077871e+01,3.155693e+01,3.    
1092     3.579106e+01,3.671150e+01,3.765827e+01,3.    
1093     4.281631e+01,4.393901e+01,4.509429e+01,4.    
1094     5.139612e+01,5.276937e+01,5.418302e+01,5.    
1095     6.190297e+01,6.358705e+01,6.532129e+01,6.    
1096     7.480211e+01,7.687239e+01,7.900504e+01,8.    
1097     9.067549e+01,9.322629e+01,9.585477e+01,9.    
1098 static const G4double P26[nE]={                  
1099     2.763953e-08,6.345949e-08,1.086410e-07,1.    
1100     5.332258e-07,6.726060e-07,8.362873e-07,1.    
1101     2.190769e-06,2.664122e-06,3.244155e-06,3.    
1102     8.346712e-06,9.958338e-06,1.190868e-05,1.    
1103     3.015277e-05,3.608891e-05,4.311988e-05,5.    
1104     1.081628e-04,1.294330e-04,1.554788e-04,1.    
1105     3.888571e-04,4.678518e-04,5.624921e-04,6.    
1106     1.407177e-03,1.691031e-03,2.035063e-03,2.    
1107     5.122493e-03,6.170607e-03,7.434979e-03,8.    
1108     1.903802e-02,2.303667e-02,2.796869e-02,3.    
1109     7.782102e-02,9.760857e-02,1.237277e-01,1.    
1110     4.837016e-01,6.271917e-01,7.920609e-01,9.    
1111     1.801311e+00,2.041842e+00,2.297083e+00,2.    
1112     3.766771e+00,4.092464e+00,4.426651e+00,4.    
1113     6.197200e+00,6.565854e+00,6.937695e+00,7.    
1114     8.826109e+00,9.206837e+00,9.587977e+00,9.    
1115     1.149595e+01,1.187788e+01,1.226005e+01,1.    
1116     1.417951e+01,1.456656e+01,1.495527e+01,1.    
1117     1.693687e+01,1.734407e+01,1.775627e+01,1.    
1118     1.990328e+01,2.035890e+01,2.082955e+01,2.    
1119     2.356264e+01,2.424460e+01,2.500821e+01,2.    
1120     3.123968e+01,3.329659e+01,3.575267e+01,3.    
1121     5.626792e+01,6.256728e+01,6.983744e+01,7.    
1122     1.259725e+02,1.423521e+02,1.608289e+02,1.    
1123     2.857326e+02,3.163352e+02,3.481975e+02,3.    
1124     5.166161e+02,5.508948e+02,5.852521e+02,6.    
1125     7.615940e+02,7.989157e+02,8.374251e+02,8.    
1126     1.057108e+03,1.108378e+03,1.162661e+03,1.    
1127     1.483376e+03,1.557421e+03,1.634595e+03,1.    
1128     2.065641e+03,2.161241e+03,2.260386e+03,2.    
1129     2.821563e+03,2.949940e+03,3.084860e+03,3.    
1130     3.877070e+03,4.063111e+03,4.259770e+03,4.    
1131     5.425268e+03,5.699711e+03,5.989703e+03,6.    
1132     7.703383e+03,8.105638e+03,8.530251e+03,8.    
1133     1.103351e+04,1.162029e+04,1.223954e+04,1.    
1134     1.589056e+04,1.674691e+04,1.765102e+04,1.    
1135     2.298963e+04,2.424390e+04,2.556894e+04,2.    
1136     3.340927e+04,3.525492e+04,3.720606e+04,3.    
1137     4.877559e+04,5.150442e+04,5.439111e+04,5.    
1138     7.154238e+04,7.559504e+04,7.988472e+04,8.    
1139     1.054184e+05,1.114616e+05,1.178617e+05,1.    
1140     1.560204e+05,1.650649e+05,1.746485e+05,1.    
1141     2.318726e+05,2.454544e+05,2.598521e+05,2.    
1142     3.459394e+05,3.663968e+05,3.880920e+05,4.    
1143     5.179768e+05,5.488770e+05,5.816596e+05,6.    
1144     7.781503e+05,8.249455e+05,8.746091e+05,9.    
1145     1.172603e+06,1.243642e+06,1.319059e+06,1.    
1146     1.772044e+06,1.880129e+06,1.994914e+06,2.    
1147 static const G4double P07[nE]={                  
1148     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1149     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1150     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1151     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1152     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1153     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1154     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1155     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1156     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1157     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1158     1.569824e-04,2.639676e-04,4.016252e-04,5.    
1159     2.302400e-03,3.161734e-03,4.237816e-03,5.    
1160     1.009134e-02,1.098015e-02,1.176593e-02,1.    
1161     1.464189e-02,1.507718e-02,1.548101e-02,1.    
1162     1.713770e-02,1.741048e-02,1.766689e-02,1.    
1163     1.873604e-02,1.891300e-02,1.907937e-02,1.    
1164     1.977430e-02,1.988986e-02,1.999879e-02,2.    
1165     2.045879e-02,2.053649e-02,2.061020e-02,2.    
1166     2.092857e-02,2.098394e-02,2.103706e-02,2.    
1167     2.127393e-02,2.131701e-02,2.135933e-02,2.    
1168     2.157041e-02,2.161596e-02,2.166444e-02,2.    
1169     2.199488e-02,2.208782e-02,2.219215e-02,2.    
1170     2.290275e-02,2.308599e-02,2.328443e-02,2.    
1171     2.453375e-02,2.483929e-02,2.516314e-02,2.    
1172     2.697807e-02,2.735392e-02,2.772372e-02,2.    
1173     2.937688e-02,2.965945e-02,2.992642e-02,3.    
1174     3.107052e-02,3.127081e-02,3.146466e-02,3.    
1175     3.237296e-02,3.254769e-02,3.272115e-02,3.    
1176     3.356878e-02,3.373266e-02,3.389385e-02,3.    
1177     3.465181e-02,3.479348e-02,3.493210e-02,3.    
1178     3.558802e-02,3.571351e-02,3.583761e-02,3.    
1179     3.644246e-02,3.656098e-02,3.667881e-02,3.    
1180     3.725808e-02,3.737196e-02,3.748515e-02,3.    
1181     3.804056e-02,3.814948e-02,3.825765e-02,3.    
1182     3.878758e-02,3.889140e-02,3.899452e-02,3.    
1183     3.949989e-02,3.959900e-02,3.969749e-02,3.    
1184     4.018101e-02,4.027603e-02,4.037051e-02,4.    
1185     4.083538e-02,4.092693e-02,4.101803e-02,4.    
1186     4.146726e-02,4.155592e-02,4.164422e-02,4.    
1187     4.208045e-02,4.216671e-02,4.225267e-02,4.    
1188     4.267816e-02,4.276246e-02,4.284650e-02,4.    
1189     4.326316e-02,4.334584e-02,4.342831e-02,4.    
1190     4.383781e-02,4.391918e-02,4.400038e-02,4.    
1191     4.440415e-02,4.448448e-02,4.456469e-02,4.    
1192     4.496398e-02,4.504353e-02,4.512297e-02,4.    
1193     4.551894e-02,4.559790e-02,4.567681e-02,4.    
1194     4.607046e-02,4.614905e-02,4.622760e-02,4.    
1195     4.661990e-02,4.669829e-02,4.677667e-02,4.    
1196 static const G4double P17[nE]={                  
1197     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1198     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1199     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1200     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1201     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1202     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1203     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1204     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1205     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1206     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1207     2.513681e-03,4.307796e-03,6.688706e-03,9.    
1208     4.349128e-02,6.138684e-02,8.449316e-02,1.    
1209     2.219850e-01,2.449239e-01,2.658362e-01,2.    
1210     3.493314e-01,3.632567e-01,3.765775e-01,3.    
1211     4.362941e-01,4.471031e-01,4.575777e-01,4.    
1212     5.053029e-01,5.139821e-01,5.223946e-01,5.    
1213     5.607851e-01,5.677982e-01,5.746128e-01,5.    
1214     6.060565e-01,6.118886e-01,6.175915e-01,6.    
1215     6.445170e-01,6.496560e-01,6.547378e-01,6.    
1216     6.795190e-01,6.844617e-01,6.894666e-01,6.    
1217     7.168454e-01,7.233059e-01,7.303926e-01,7.    
1218     7.837012e-01,7.999913e-01,8.188394e-01,8.    
1219     9.603004e-01,9.999814e-01,1.044271e+00,1.    
1220     1.351156e+00,1.432891e+00,1.522179e+00,1.    
1221     2.071307e+00,2.195490e+00,2.321415e+00,2.    
1222     2.936495e+00,3.051807e+00,3.164089e+00,3.    
1223     3.689814e+00,3.790775e+00,3.891483e+00,3.    
1224     4.408024e+00,4.516815e+00,4.628125e+00,4.    
1225     5.223773e+00,5.349792e+00,5.477539e+00,5.    
1226     6.134829e+00,6.269367e+00,6.405046e+00,6.    
1227     7.107768e+00,7.254950e+00,7.404963e+00,7.    
1228     8.205696e+00,8.377371e+00,8.553273e+00,8.    
1229     9.500452e+00,9.704160e+00,9.912848e+00,1.    
1230     1.103432e+01,1.127492e+01,1.152121e+01,1.    
1231     1.284257e+01,1.312581e+01,1.341573e+01,1.    
1232     1.497192e+01,1.530584e+01,1.564780e+01,1.    
1233     1.748656e+01,1.788187e+01,1.828699e+01,1.    
1234     2.047018e+01,2.094055e+01,2.142296e+01,2.    
1235     2.402840e+01,2.459093e+01,2.516825e+01,2.    
1236     2.829276e+01,2.896865e+01,2.966276e+01,3.    
1237     3.342632e+01,3.424188e+01,3.507988e+01,3.    
1238     3.963126e+01,4.061907e+01,4.163459e+01,4.    
1239     4.715847e+01,4.835904e+01,4.959386e+01,5.    
1240     5.631988e+01,5.778360e+01,5.928972e+01,6.    
1241     6.750389e+01,6.929358e+01,7.113583e+01,7.    
1242     8.119494e+01,8.338901e+01,8.564834e+01,8.    
1243     9.799814e+01,1.006946e+02,1.034722e+02,1.    
1244     1.186703e+02,1.219918e+02,1.254143e+02,1.    
1245 static const G4double P27[nE]={                  
1246     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1247     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1248     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1249     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1250     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1251     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1252     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1253     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1254     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1255     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1256     4.026430e-02,7.035354e-02,1.115368e-01,1.    
1257     8.253354e-01,1.198038e+00,1.694231e+00,2.    
1258     4.929838e+00,5.521904e+00,6.078501e+00,6.    
1259     8.507290e+00,8.952817e+00,9.392248e+00,9.    
1260     1.154897e+01,1.197730e+01,1.240524e+01,1.    
1261     1.453973e+01,1.496547e+01,1.539087e+01,1.    
1262     1.751579e+01,1.794142e+01,1.836777e+01,1.    
1263     2.052125e+01,2.095904e+01,2.140031e+01,2.    
1264     2.368188e+01,2.415886e+01,2.464505e+01,2.    
1265     2.724263e+01,2.780976e+01,2.840169e+01,2.    
1266     3.195970e+01,3.287602e+01,3.391207e+01,3.    
1267     4.252816e+01,4.538374e+01,4.878905e+01,5.    
1268     7.700256e+01,8.559606e+01,9.548163e+01,1.    
1269     1.710065e+02,1.928734e+02,2.174931e+02,2.    
1270     3.839498e+02,4.249834e+02,4.678667e+02,5.    
1271     6.971481e+02,7.442078e+02,7.914355e+02,8.    
1272     1.033470e+03,1.084366e+03,1.136689e+03,1.    
1273     1.430999e+03,1.498741e+03,1.570174e+03,1.    
1274     1.989542e+03,2.086454e+03,2.187703e+03,2.    
1275     2.758773e+03,2.886552e+03,3.019360e+03,3.    
1276     3.773646e+03,3.946289e+03,4.127639e+03,4.    
1277     5.189702e+03,5.438396e+03,5.701012e+03,5.    
1278     7.252683e+03,7.617117e+03,8.001887e+03,8.    
1279     1.027061e+04,1.080218e+04,1.136297e+04,1.    
1280     1.466398e+04,1.543675e+04,1.625197e+04,1.    
1281     2.105305e+04,2.217810e+04,2.336552e+04,2.    
1282     3.037128e+04,3.201606e+04,3.375326e+04,3.    
1283     4.402572e+04,4.644258e+04,4.899711e+04,5.    
1284     6.413682e+04,6.770616e+04,7.148143e+04,7.    
1285     9.390298e+04,9.919906e+04,1.048042e+05,1.    
1286     1.381568e+05,1.460482e+05,1.544049e+05,1.    
1287     2.042150e+05,2.160183e+05,2.285240e+05,2.    
1288     3.031790e+05,3.208942e+05,3.396723e+05,3.    
1289     4.519286e+05,4.786000e+05,5.068836e+05,5.    
1290     6.761823e+05,7.164532e+05,7.591751e+05,8.    
1291     1.015202e+06,1.076168e+06,1.140869e+06,1.    
1292     1.529038e+06,1.621562e+06,1.719787e+06,1.    
1293     2.309695e+06,2.450437e+06,2.599897e+06,2.    
1294 static const G4double P08[nE]={                  
1295     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1296     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1297     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1298     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1299     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1300     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1301     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1302     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1303     0.000000e+00,0.000000e+00,0.000000e+00,5.    
1304     1.811834e-04,2.505496e-04,3.370214e-04,4.    
1305     1.193306e-03,1.512903e-03,1.914433e-03,2.    
1306     6.198077e-03,7.742119e-03,9.530727e-03,1.    
1307     1.873750e-02,2.022106e-02,2.156321e-02,2.    
1308     2.660257e-02,2.735786e-02,2.805033e-02,2.    
1309     3.077178e-02,3.119783e-02,3.159205e-02,3.    
1310     3.316932e-02,3.341990e-02,3.365279e-02,3.    
1311     3.459791e-02,3.475094e-02,3.489419e-02,3.    
1312     3.549011e-02,3.558968e-02,3.568400e-02,3.    
1313     3.609170e-02,3.616308e-02,3.623182e-02,3.    
1314     3.653815e-02,3.659408e-02,3.664961e-02,3.    
1315     3.693479e-02,3.699783e-02,3.706528e-02,3.    
1316     3.752719e-02,3.765644e-02,3.780102e-02,3.    
1317     3.877367e-02,3.902186e-02,3.928968e-02,3.    
1318     4.095916e-02,4.136456e-02,4.179369e-02,4.    
1319     4.420281e-02,4.470553e-02,4.520218e-02,4.    
1320     4.744879e-02,4.783646e-02,4.820329e-02,4.    
1321     4.977449e-02,5.004805e-02,5.031198e-02,5.    
1322     5.153419e-02,5.176632e-02,5.199590e-02,5.    
1323     5.311040e-02,5.332567e-02,5.353771e-02,5.    
1324     5.454172e-02,5.473078e-02,5.491615e-02,5.    
1325     5.579619e-02,5.596472e-02,5.613134e-02,5.    
1326     5.694179e-02,5.710018e-02,5.725750e-02,5.    
1327     5.802884e-02,5.818008e-02,5.833029e-02,5.    
1328     5.906581e-02,5.920977e-02,5.935267e-02,5.    
1329     6.005165e-02,6.018840e-02,6.032417e-02,6.    
1330     6.098884e-02,6.111906e-02,6.124841e-02,6.    
1331     6.188297e-02,6.200757e-02,6.213144e-02,6.    
1332     6.274049e-02,6.286037e-02,6.297964e-02,6.    
1333     6.356743e-02,6.368338e-02,6.379883e-02,6.    
1334     6.436902e-02,6.448173e-02,6.459402e-02,6.    
1335     6.514969e-02,6.525973e-02,6.536943e-02,6.    
1336     6.591318e-02,6.602104e-02,6.612863e-02,6.    
1337     6.666270e-02,6.676880e-02,6.687467e-02,6.    
1338     6.740099e-02,6.750569e-02,6.761022e-02,6.    
1339     6.813049e-02,6.823412e-02,6.833761e-02,6.    
1340     6.885334e-02,6.895618e-02,6.905893e-02,6.    
1341     6.957150e-02,6.967381e-02,6.977607e-02,6.    
1342     7.028673e-02,7.038876e-02,7.049077e-02,7.    
1343 static const G4double P18[nE]={                  
1344     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1345     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1346     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1347     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1348     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1349     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1350     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1351     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1352     0.000000e+00,0.000000e+00,0.000000e+00,6.    
1353     2.271897e-03,3.207104e-03,4.409961e-03,5.    
1354     1.767709e-02,2.303624e-02,2.998068e-02,3.    
1355     1.124957e-01,1.446424e-01,1.830426e-01,2.    
1356     3.993962e-01,4.376870e-01,4.734083e-01,5.    
1357     6.197329e-01,6.438952e-01,6.667360e-01,6.    
1358     7.647338e-01,7.816149e-01,7.977194e-01,8.    
1359     8.680643e-01,8.803545e-01,8.921305e-01,9.    
1360     9.443113e-01,9.535980e-01,9.625599e-01,9.    
1361     1.003283e+00,1.010757e+00,1.018054e+00,1.    
1362     1.052537e+00,1.059162e+00,1.065739e+00,1.    
1363     1.097776e+00,1.104194e+00,1.110761e+00,1.    
1364     1.147772e+00,1.156712e+00,1.166573e+00,1.    
1365     1.241088e+00,1.263744e+00,1.289862e+00,1.    
1366     1.483446e+00,1.537190e+00,1.596963e+00,1.    
1367     2.006982e+00,2.115434e+00,2.233748e+00,2.    
1368     2.962731e+00,3.128831e+00,3.297953e+00,3.    
1369     4.134030e+00,4.292233e+00,4.446513e+00,4.    
1370     5.168447e+00,5.306340e+00,5.443454e+00,5.    
1371     6.138329e+00,6.282856e+00,6.430177e+00,6.    
1372     7.213302e+00,7.378833e+00,7.546883e+00,7.    
1373     8.417680e+00,8.597232e+00,8.778664e+00,8.    
1374     9.721535e+00,9.919203e+00,1.012061e+01,1.    
1375     1.119347e+01,1.142290e+01,1.165776e+01,1.    
1376     1.291891e+01,1.318945e+01,1.346640e+01,1.    
1377     1.495148e+01,1.526949e+01,1.559485e+01,1.    
1378     1.733768e+01,1.771076e+01,1.809249e+01,1.    
1379     2.013915e+01,2.057787e+01,2.102702e+01,2.    
1380     2.344009e+01,2.395847e+01,2.448959e+01,2.    
1381     2.734990e+01,2.796580e+01,2.859733e+01,2.    
1382     3.200638e+01,3.274207e+01,3.349698e+01,3.    
1383     3.758087e+01,3.846397e+01,3.937074e+01,4.    
1384     4.428566e+01,4.535037e+01,4.644426e+01,4.    
1385     5.238376e+01,5.367249e+01,5.499726e+01,5.    
1386     6.220151e+01,6.376694e+01,6.537690e+01,6.    
1387     7.414450e+01,7.605215e+01,7.801493e+01,8.    
1388     8.871771e+01,9.104924e+01,9.344911e+01,9.    
1389     1.065509e+02,1.094082e+02,1.123504e+02,1.    
1390     1.284306e+02,1.319411e+02,1.355571e+02,1.    
1391     1.553403e+02,1.596635e+02,1.641179e+02,1.    
1392 static const G4double P28[nE]={                  
1393     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1394     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1395     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1396     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1397     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1398     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1399     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1400     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1401     0.000000e+00,0.000000e+00,0.000000e+00,7.    
1402     2.852749e-02,4.113714e-02,5.787072e-02,7.    
1403     2.638146e-01,3.536865e-01,4.737998e-01,6.    
1404     2.066174e+00,2.735516e+00,3.560003e+00,4.    
1405     8.653933e+00,9.642298e+00,1.059309e+01,1.    
1406     1.485012e+01,1.562315e+01,1.637661e+01,1.    
1407     1.991230e+01,2.058123e+01,2.123916e+01,2.    
1408     2.438256e+01,2.498540e+01,2.558088e+01,2.    
1409     2.846739e+01,2.903101e+01,2.959170e+01,3.    
1410     3.237989e+01,3.294094e+01,3.350557e+01,3.    
1411     3.642781e+01,3.704271e+01,3.767193e+01,3.    
1412     4.102906e+01,4.176548e+01,4.254216e+01,4.    
1413     4.735458e+01,4.862267e+01,5.006425e+01,5.    
1414     6.210756e+01,6.607900e+01,7.079781e+01,7.    
1415     1.093979e+02,1.210367e+02,1.343782e+02,1.    
1416     2.352665e+02,2.642811e+02,2.969039e+02,3.    
1417     5.179023e+02,5.727873e+02,6.303818e+02,6.    
1418     9.421156e+02,1.006680e+03,1.071573e+03,1.    
1419     1.403917e+03,1.473430e+03,1.544668e+03,1.    
1420     1.940482e+03,2.030475e+03,2.125017e+03,2.    
1421     2.676334e+03,2.803633e+03,2.936826e+03,3.    
1422     3.693517e+03,3.864050e+03,4.041643e+03,4.    
1423     5.053744e+03,5.285606e+03,5.529085e+03,5.    
1424     6.952019e+03,7.284382e+03,7.635021e+03,8.    
1425     9.700911e+03,1.018492e+04,1.069554e+04,1.    
1426     1.369971e+04,1.440231e+04,1.514313e+04,1.    
1427     1.949694e+04,2.051483e+04,2.158819e+04,2.    
1428     2.790227e+04,2.938043e+04,3.094006e+04,3.    
1429     4.013379e+04,4.229066e+04,4.456817e+04,4.    
1430     5.802649e+04,6.119106e+04,6.453527e+04,6.    
1431     8.434440e+04,8.901242e+04,9.394902e+04,9.    
1432     1.232549e+05,1.301745e+05,1.374970e+05,1.    
1433     1.810526e+05,1.913548e+05,2.022635e+05,2.    
1434     2.672645e+05,2.826635e+05,2.989775e+05,3.    
1435     3.963420e+05,4.194409e+05,4.439240e+05,4.    
1436     5.902534e+05,6.250137e+05,6.618730e+05,7.    
1437     8.824628e+05,9.349258e+05,9.905789e+05,1.    
1438     1.324048e+06,1.403444e+06,1.487699e+06,1.    
1439     1.993117e+06,2.113574e+06,2.241448e+06,2.    
1440     3.009326e+06,3.192510e+06,3.387034e+06,3.    
1441 static const G4double P09[nE]={                  
1442     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1443     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1444     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1445     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1446     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1447     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1448     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1449     2.874685e-05,4.278170e-05,6.015953e-05,8.    
1450     2.225746e-04,2.816669e-04,3.547278e-04,4.    
1451     1.095439e-03,1.369207e-03,1.718647e-03,2.    
1452     5.283815e-03,6.581892e-03,8.166275e-03,1.    
1453     2.061095e-02,2.348657e-02,2.621326e-02,2.    
1454     3.665237e-02,3.817043e-02,3.954315e-02,4.    
1455     4.478859e-02,4.559645e-02,4.634342e-02,4.    
1456     4.934096e-02,4.981961e-02,5.026497e-02,5.    
1457     5.207454e-02,5.236721e-02,5.264095e-02,5.    
1458     5.377608e-02,5.396514e-02,5.414407e-02,5.    
1459     5.491664e-02,5.505189e-02,5.518220e-02,5.    
1460     5.577674e-02,5.588719e-02,5.599565e-02,5.    
1461     5.650990e-02,5.660948e-02,5.670951e-02,5.    
1462     5.723655e-02,5.735463e-02,5.748102e-02,5.    
1463     5.833443e-02,5.856813e-02,5.882715e-02,5.    
1464     6.052514e-02,6.094997e-02,6.140573e-02,6.    
1465     6.420419e-02,6.487694e-02,6.558787e-02,6.    
1466     6.959533e-02,7.044270e-02,7.128573e-02,7.    
1467     7.517973e-02,7.586367e-02,7.651297e-02,7.    
1468     7.929835e-02,7.978018e-02,8.024301e-02,8.    
1469     8.235029e-02,8.274289e-02,8.312891e-02,8.    
1470     8.498082e-02,8.533691e-02,8.568805e-02,8.    
1471     8.736369e-02,8.768234e-02,8.799567e-02,8.    
1472     8.949159e-02,8.977874e-02,9.006257e-02,9.    
1473     9.143968e-02,9.170781e-02,9.197375e-02,9.    
1474     9.327188e-02,9.352532e-02,9.377671e-02,9.    
1475     9.500306e-02,9.524227e-02,9.547949e-02,9.    
1476     9.663655e-02,9.686235e-02,9.708636e-02,9.    
1477     9.818069e-02,9.839467e-02,9.860712e-02,9.    
1478     9.964757e-02,9.985156e-02,1.000543e-01,1.    
1479     1.010497e-01,1.012454e-01,1.014400e-01,1.    
1480     1.023982e-01,1.025870e-01,1.027750e-01,1.    
1481     1.037025e-01,1.038858e-01,1.040682e-01,1.    
1482     1.049706e-01,1.051492e-01,1.053272e-01,1.    
1483     1.062090e-01,1.063838e-01,1.065581e-01,1.    
1484     1.074231e-01,1.075949e-01,1.077663e-01,1.    
1485     1.086179e-01,1.087872e-01,1.089563e-01,1.    
1486     1.097973e-01,1.099648e-01,1.101320e-01,1.    
1487     1.109651e-01,1.111312e-01,1.112971e-01,1.    
1488     1.121246e-01,1.122897e-01,1.124547e-01,1.    
1489     1.132786e-01,1.134432e-01,1.136078e-01,1.    
1490 static const G4double P19[nE]={                  
1491     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1492     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1493     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1494     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1495     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1496     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1497     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1498     2.349414e-04,3.566209e-04,5.121431e-04,7.    
1499     2.132778e-03,2.772367e-03,3.588379e-03,4.    
1500     1.281519e-02,1.650640e-02,2.136751e-02,2.    
1501     7.661025e-02,9.837506e-02,1.257751e-01,1.    
1502     3.636123e-01,4.234494e-01,4.819630e-01,5.    
1503     7.263609e-01,7.655421e-01,8.020772e-01,8.    
1504     9.544631e-01,9.803078e-01,1.004947e+00,1.    
1505     1.112933e+00,1.131898e+00,1.150092e+00,1.    
1506     1.230825e+00,1.245180e+00,1.259022e+00,1.    
1507     1.321726e+00,1.333200e+00,1.344393e+00,1.    
1508     1.397235e+00,1.407387e+00,1.417470e+00,1.    
1509     1.467817e+00,1.478068e+00,1.488444e+00,1.    
1510     1.542307e+00,1.553733e+00,1.565564e+00,1.    
1511     1.633991e+00,1.650739e+00,1.669215e+00,1.    
1512     1.806833e+00,1.847794e+00,1.894587e+00,1.    
1513     2.232374e+00,2.324368e+00,2.426086e+00,2.    
1514     3.113195e+00,3.293165e+00,3.489174e+00,3.    
1515     4.702026e+00,4.982009e+00,5.269084e+00,5.    
1516     6.718845e+00,6.997954e+00,7.271040e+00,7.    
1517     8.550792e+00,8.793666e+00,9.034113e+00,9.    
1518     1.023171e+01,1.047615e+01,1.072385e+01,1.    
1519     1.202490e+01,1.229872e+01,1.257701e+01,1.    
1520     1.403061e+01,1.433323e+01,1.463991e+01,1.    
1521     1.624276e+01,1.657956e+01,1.692265e+01,1.    
1522     1.874553e+01,1.913392e+01,1.953093e+01,1.    
1523     2.165319e+01,2.210656e+01,2.257002e+01,2.    
1524     2.504595e+01,2.557440e+01,2.611448e+01,2.    
1525     2.899935e+01,2.961537e+01,3.024518e+01,3.    
1526     3.361474e+01,3.433567e+01,3.507331e+01,3.    
1527     3.902978e+01,3.987847e+01,4.074762e+01,4.    
1528     4.542228e+01,4.642768e+01,4.745822e+01,4.    
1529     5.301540e+01,5.421353e+01,5.544261e+01,5.    
1530     6.208597e+01,6.352141e+01,6.499498e+01,6.    
1531     7.297649e+01,7.470440e+01,7.647932e+01,7.    
1532     8.611087e+01,8.819958e+01,9.034631e+01,9.    
1533     1.020148e+02,1.045491e+02,1.071551e+02,1.    
1534     1.213411e+02,1.244265e+02,1.276007e+02,1.    
1535     1.449027e+02,1.486706e+02,1.525485e+02,1.    
1536     1.737129e+02,1.783272e+02,1.830781e+02,1.    
1537     2.090368e+02,2.147025e+02,2.205380e+02,2.    
1538     2.524567e+02,2.594302e+02,2.666150e+02,2.    
1539 static const G4double P29[nE]={                  
1540     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1541     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1542     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1543     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1544     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1545     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1546     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1547     1.921474e-03,2.976502e-03,4.368459e-03,6.    
1548     2.057119e-02,2.749440e-02,3.660912e-02,4.    
1549     1.517623e-01,2.015348e-01,2.691639e-01,3.    
1550     1.126750e+00,1.491709e+00,1.965598e+00,2.    
1551     6.519588e+00,7.764801e+00,9.020575e+00,1.    
1552     1.475356e+01,1.576491e+01,1.673737e+01,1.    
1553     2.117302e+01,2.199989e+01,2.281269e+01,2.    
1554     2.671041e+01,2.746196e+01,2.820524e+01,2.    
1555     3.181409e+01,3.251822e+01,3.321818e+01,3.    
1556     3.668862e+01,3.738499e+01,3.808532e+01,3.    
1557     4.170647e+01,4.246862e+01,4.324879e+01,4.    
1558     4.752048e+01,4.847204e+01,4.946482e+01,5.    
1559     5.511729e+01,5.642837e+01,5.782761e+01,5.    
1560     6.672884e+01,6.910427e+01,7.180553e+01,7.    
1561     9.403866e+01,1.012187e+02,1.096725e+02,1.    
1562     1.769950e+02,1.969173e+02,2.196207e+02,2.    
1563     3.886448e+02,4.367924e+02,4.908381e+02,5.    
1564     8.585945e+02,9.511118e+02,1.048877e+03,1.    
1565     1.589651e+03,1.703560e+03,1.818425e+03,1.    
1566     2.407530e+03,2.529963e+03,2.654887e+03,2.    
1567     3.336794e+03,3.488995e+03,3.647955e+03,3.    
1568     4.563730e+03,4.774309e+03,4.994874e+03,5.    
1569     6.258253e+03,6.545669e+03,6.845863e+03,7.    
1570     8.566547e+03,8.961603e+03,9.376365e+03,9.    
1571     1.179388e+04,1.235652e+04,1.294923e+04,1.    
1572     1.642541e+04,1.723648e+04,1.809099e+04,1.    
1573     2.309921e+04,2.426670e+04,2.549643e+04,2.    
1574     3.270279e+04,3.438349e+04,3.615442e+04,3.    
1575     4.654932e+04,4.897832e+04,5.153970e+04,5.    
1576     6.661323e+04,7.014443e+04,7.387144e+04,7.    
1577     9.586609e+04,1.010320e+05,1.064891e+05,1.    
1578     1.387798e+05,1.463820e+05,1.544194e+05,1.    
1579     2.020912e+05,2.133387e+05,2.252383e+05,2.    
1580     2.959687e+05,3.126882e+05,3.303881e+05,3.    
1581     4.357933e+05,4.607512e+05,4.871874e+05,5.    
1582     6.448833e+05,6.822787e+05,7.219092e+05,7.    
1583     9.586691e+05,1.014890e+06,1.074498e+06,1.    
1584     1.431099e+06,1.515882e+06,1.605811e+06,1.    
1585     2.144486e+06,2.272704e+06,2.408756e+06,2.    
1586     3.224662e+06,3.419071e+06,3.625434e+06,3.    
1587     4.864335e+06,5.159824e+06,5.473583e+06,5.    
1588 static const G4double P010[nE]={                 
1589     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1590     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1591     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1592     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1593     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1594     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1595     1.420000e-05,4.701887e-05,8.694089e-05,1.    
1596     4.837658e-04,6.323004e-04,8.240919e-04,1.    
1597     2.792245e-03,3.550245e-03,4.501045e-03,5.    
1598     1.410796e-02,1.749772e-02,2.157807e-02,2.    
1599     5.345748e-02,6.167108e-02,6.986856e-02,7.    
1600     1.016834e-01,1.059752e-01,1.097284e-01,1.    
1601     1.228498e-01,1.246924e-01,1.263587e-01,1.    
1602     1.327952e-01,1.338059e-01,1.347478e-01,1.    
1603     1.386291e-01,1.392709e-01,1.398763e-01,1.    
1604     1.424548e-01,1.428983e-01,1.433229e-01,1.    
1605     1.452240e-01,1.455704e-01,1.459086e-01,1.    
1606     1.475078e-01,1.478144e-01,1.481178e-01,1.    
1607     1.495976e-01,1.498872e-01,1.501749e-01,1.    
1608     1.515759e-01,1.518503e-01,1.521247e-01,1.    
1609     1.535529e-01,1.538711e-01,1.542117e-01,1.    
1610     1.564842e-01,1.570896e-01,1.577520e-01,1.    
1611     1.619445e-01,1.629686e-01,1.640604e-01,1.    
1612     1.706592e-01,1.722273e-01,1.738801e-01,1.    
1613     1.831911e-01,1.851767e-01,1.871639e-01,1.    
1614     1.965384e-01,1.982195e-01,1.998234e-01,2.    
1615     2.067461e-01,2.079405e-01,2.090837e-01,2.    
1616     2.142048e-01,2.151398e-01,2.160530e-01,2.    
1617     2.203669e-01,2.211890e-01,2.219995e-01,2.    
1618     2.258871e-01,2.266330e-01,2.273689e-01,2.    
1619     2.309088e-01,2.315914e-01,2.322663e-01,2.    
1620     2.355358e-01,2.361701e-01,2.367983e-01,2.    
1621     2.398478e-01,2.404398e-01,2.410260e-01,2.    
1622     2.438703e-01,2.444223e-01,2.449688e-01,2.    
1623     2.476231e-01,2.481389e-01,2.486500e-01,2.    
1624     2.511383e-01,2.516232e-01,2.521042e-01,2.    
1625     2.544533e-01,2.549127e-01,2.553688e-01,2.    
1626     2.576034e-01,2.580418e-01,2.584775e-01,2.    
1627     2.606187e-01,2.610399e-01,2.614590e-01,2.    
1628     2.635238e-01,2.639311e-01,2.643365e-01,2.    
1629     2.663389e-01,2.667347e-01,2.671291e-01,2.    
1630     2.690805e-01,2.694670e-01,2.698523e-01,2.    
1631     2.717623e-01,2.721413e-01,2.725193e-01,2.    
1632     2.743960e-01,2.747690e-01,2.751411e-01,2.    
1633     2.769916e-01,2.773598e-01,2.777275e-01,2.    
1634     2.795578e-01,2.799225e-01,2.802868e-01,2.    
1635     2.821025e-01,2.824646e-01,2.828265e-01,2.    
1636     2.846324e-01,2.849930e-01,2.853535e-01,2.    
1637 static const G4double P110[nE]={                 
1638     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1639     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1640     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1641     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1642     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1643     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1644     9.591763e-05,3.231969e-04,6.087125e-04,9.    
1645     3.781121e-03,5.069066e-03,6.785649e-03,8.    
1646     2.645527e-02,3.465992e-02,4.527966e-02,5.    
1647     1.648770e-01,2.105763e-01,2.673302e-01,3.    
1648     7.582455e-01,8.958848e-01,1.037582e+00,1.    
1649     1.637162e+00,1.726448e+00,1.806975e+00,1.    
1650     2.113806e+00,2.161365e+00,2.205714e+00,2.    
1651     2.392761e+00,2.425097e+00,2.456165e+00,2.    
1652     2.596079e+00,2.621512e+00,2.646242e+00,2.    
1653     2.761415e+00,2.783170e+00,2.804641e+00,2.    
1654     2.909838e+00,2.930865e+00,2.952023e+00,2.    
1655     3.061616e+00,3.084631e+00,3.108110e+00,3.    
1656     3.233578e+00,3.260461e+00,3.287985e+00,3.    
1657     3.434812e+00,3.466295e+00,3.498755e+00,3.    
1658     3.684134e+00,3.729263e+00,3.779065e+00,3.    
1659     4.145355e+00,4.251473e+00,4.371124e+00,4.    
1660     5.204640e+00,5.426399e+00,5.670065e+00,5.    
1661     7.289825e+00,7.709327e+00,8.164993e+00,8.    
1662     1.098318e+01,1.163929e+01,1.231599e+01,1.    
1663     1.580777e+01,1.649383e+01,1.716840e+01,1.    
1664     2.034918e+01,2.095125e+01,2.154514e+01,2.    
1665     2.445438e+01,2.503650e+01,2.562246e+01,2.    
1666     2.865235e+01,2.928457e+01,2.992690e+01,3.    
1667     3.329988e+01,3.400829e+01,3.472857e+01,3.    
1668     3.852199e+01,3.932257e+01,4.013841e+01,4.    
1669     4.446606e+01,4.538488e+01,4.632254e+01,4.    
1670     5.130758e+01,5.236662e+01,5.344725e+01,5.    
1671     5.918909e+01,6.040851e+01,6.165285e+01,6.    
1672     6.826996e+01,6.967727e+01,7.111431e+01,7.    
1673     7.877533e+01,8.040918e+01,8.207925e+01,8.    
1674     9.101153e+01,9.292265e+01,9.487827e+01,9.    
1675     1.053723e+02,1.076245e+02,1.099316e+02,1.    
1676     1.223492e+02,1.250219e+02,1.277621e+02,1.    
1677     1.425502e+02,1.457410e+02,1.490151e+02,1.    
1678     1.667257e+02,1.705553e+02,1.744876e+02,1.    
1679     1.958027e+02,2.004204e+02,2.051650e+02,2.    
1680     2.309295e+02,2.365206e+02,2.422684e+02,2.    
1681     2.735316e+02,2.803262e+02,2.873146e+02,2.    
1682     3.253817e+02,3.336664e+02,3.421912e+02,3.    
1683     3.886893e+02,3.988214e+02,4.092516e+02,4.    
1684     4.662121e+02,4.786383e+02,4.914349e+02,5.    
1685     5.613974e+02,5.766763e+02,5.924161e+02,6.    
1686 static const G4double P210[nE]={                 
1687     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1688     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1689     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1690     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1691     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1692     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1693     6.479160e-04,2.222020e-03,4.264157e-03,6.    
1694     2.967511e-02,4.084384e-02,5.620899e-02,7.    
1695     2.531616e-01,3.419767e-01,4.606012e-01,6.    
1696     1.952294e+00,2.568440e+00,3.357900e+00,4.    
1697     1.093215e+01,1.323883e+01,1.568832e+01,1.    
1698     2.701003e+01,2.886770e+01,3.059561e+01,3.    
1699     3.778462e+01,3.901221e+01,4.019267e+01,4.    
1700     4.563906e+01,4.667362e+01,4.769852e+01,4.    
1701     5.275195e+01,5.375980e+01,5.477015e+01,5.    
1702     5.992444e+01,6.099162e+01,6.207742e+01,6.    
1703     6.791004e+01,6.918632e+01,7.051016e+01,7.    
1704     7.803498e+01,7.976279e+01,8.157970e+01,8.    
1705     9.223838e+01,9.473382e+01,9.736730e+01,1.    
1706     1.127846e+02,1.163971e+02,1.202362e+02,1.    
1707     1.443459e+02,1.507468e+02,1.580279e+02,1.    
1708     2.171791e+02,2.357799e+02,2.573961e+02,2.    
1709     4.234203e+02,4.714435e+02,5.258289e+02,5.    
1710     9.241665e+02,1.036396e+03,1.162036e+03,1.    
1711     2.016629e+03,2.233434e+03,2.463894e+03,2.    
1712     3.766953e+03,4.046947e+03,4.330688e+03,4.    
1713     5.794940e+03,6.098442e+03,6.406997e+03,6.    
1714     8.062854e+03,8.425307e+03,8.801331e+03,9.    
1715     1.093345e+04,1.141964e+04,1.192874e+04,1.    
1716     1.486078e+04,1.553362e+04,1.623867e+04,1.    
1717     2.031147e+04,2.125054e+04,2.223681e+04,2.    
1718     2.797583e+04,2.930685e+04,3.070674e+04,3.    
1719     3.887109e+04,4.076570e+04,4.275810e+04,4.    
1720     5.437120e+04,5.706526e+04,5.989850e+04,6.    
1721     7.642649e+04,8.026610e+04,8.430679e+04,8.    
1722     1.079389e+05,1.134438e+05,1.192429e+05,1.    
1723     1.532716e+05,1.612234e+05,1.696092e+05,1.    
1724     2.189820e+05,2.305546e+05,2.427716e+05,2.    
1725     3.149225e+05,3.318808e+05,3.497998e+05,3.    
1726     4.559143e+05,4.809159e+05,5.073551e+05,5.    
1727     6.642985e+05,7.013543e+05,7.405681e+05,7.    
1728     9.738297e+05,1.029007e+06,1.087434e+06,1.    
1729     1.435629e+06,1.518129e+06,1.605537e+06,1.    
1730     2.127300e+06,2.251108e+06,2.382345e+06,2.    
1731     3.166910e+06,3.353328e+06,3.551019e+06,3.    
1732     4.734475e+06,5.016015e+06,5.314705e+06,5.    
1733     7.105007e+06,7.531392e+06,7.983922e+06,8.    
1734     1.069945e+07,1.134686e+07,1.203421e+07,1.    
1735 static const G4double P011[nE]={                 
1736     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1737     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1738     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1739     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1740     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1741     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1742     4.819877e-04,6.963055e-04,9.598242e-04,1.    
1743     3.659508e-03,4.684845e-03,6.012414e-03,7.    
1744     1.946512e-02,2.442244e-02,3.048669e-02,3.    
1745     8.219187e-02,9.763501e-02,1.147626e-01,1.    
1746     2.058575e-01,2.199478e-01,2.320907e-01,2.    
1747     2.715525e-01,2.765412e-01,2.808940e-01,2.    
1748     2.962019e-01,2.983884e-01,3.003811e-01,3.    
1749     3.082963e-01,3.095866e-01,3.108058e-01,3.    
1750     3.160703e-01,3.169924e-01,3.178803e-01,3.    
1751     3.219164e-01,3.226615e-01,3.233911e-01,3.    
1752     3.268613e-01,3.275281e-01,3.281876e-01,3.    
1753     3.313894e-01,3.320112e-01,3.326266e-01,3.    
1754     3.355989e-01,3.361705e-01,3.367340e-01,3.    
1755     3.394248e-01,3.399435e-01,3.404607e-01,3.    
1756     3.431912e-01,3.438227e-01,3.445089e-01,3.    
1757     3.491345e-01,3.503420e-01,3.516483e-01,3.    
1758     3.597024e-01,3.616397e-01,3.636979e-01,3.    
1759     3.760363e-01,3.789488e-01,3.820118e-01,3.    
1760     3.991939e-01,4.028554e-01,4.065239e-01,4.    
1761     4.239442e-01,4.270940e-01,4.301065e-01,4.    
1762     4.431646e-01,4.454194e-01,4.475751e-01,4.    
1763     4.571727e-01,4.589098e-01,4.606012e-01,4.    
1764     4.685277e-01,4.700296e-01,4.715090e-01,4.    
1765     4.786094e-01,4.799758e-01,4.813257e-01,4.    
1766     4.878438e-01,4.891042e-01,4.903511e-01,4.    
1767     4.963917e-01,4.975623e-01,4.987206e-01,4.    
1768     5.043297e-01,5.054154e-01,5.064892e-01,5.    
1769     5.116834e-01,5.126883e-01,5.136823e-01,5.    
1770     5.184958e-01,5.194287e-01,5.203523e-01,5.    
1771     5.248375e-01,5.257097e-01,5.265740e-01,5.    
1772     5.307869e-01,5.316092e-01,5.324250e-01,5.    
1773     5.364158e-01,5.371974e-01,5.379738e-01,5.    
1774     5.417840e-01,5.425326e-01,5.432770e-01,5.    
1775     5.469403e-01,5.476620e-01,5.483803e-01,5.    
1776     5.519240e-01,5.526238e-01,5.533208e-01,5.    
1777     5.567669e-01,5.574489e-01,5.581286e-01,5.    
1778     5.614955e-01,5.621630e-01,5.628287e-01,5.    
1779     5.661318e-01,5.667877e-01,5.674422e-01,5.    
1780     5.706946e-01,5.713415e-01,5.719872e-01,5.    
1781     5.752004e-01,5.758403e-01,5.764794e-01,5.    
1782     5.796636e-01,5.802985e-01,5.809328e-01,5.    
1783     5.840970e-01,5.847286e-01,5.853599e-01,5.    
1784 static const G4double P111[nE]={                 
1785     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1786     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1787     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1788     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1789     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1790     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1791     3.107743e-03,4.592002e-03,6.476688e-03,8.    
1792     2.806967e-02,3.696045e-02,4.884261e-02,6.    
1793     1.832285e-01,2.368826e-01,3.046122e-01,3.    
1794     9.458741e-01,1.154034e+00,1.392213e+00,1.    
1795     2.780710e+00,3.016745e+00,3.226555e+00,3.    
1796     3.968634e+00,4.072420e+00,4.165812e+00,4.    
1797     4.523860e+00,4.580294e+00,4.633332e+00,4.    
1798     4.863513e+00,4.904792e+00,4.945010e+00,4.    
1799     5.134994e+00,5.171534e+00,5.207809e+00,5.    
1800     5.388348e+00,5.424899e+00,5.461794e+00,5.    
1801     5.654057e+00,5.694526e+00,5.735793e+00,5.    
1802     5.955296e+00,6.001974e+00,6.049594e+00,6.    
1803     6.301538e+00,6.354593e+00,6.408503e+00,6.    
1804     6.690414e+00,6.749930e+00,6.811096e+00,6.    
1805     7.165738e+00,7.255305e+00,7.355625e+00,7.    
1806     8.101065e+00,8.312698e+00,8.548675e+00,8.    
1807     1.014926e+01,1.056874e+01,1.102811e+01,1.    
1808     1.405626e+01,1.483537e+01,1.567986e+01,1.    
1809     2.088013e+01,2.208995e+01,2.333920e+01,2.    
1810     2.982889e+01,3.111433e+01,3.238135e+01,3.    
1811     3.838153e+01,3.951805e+01,4.063795e+01,4.    
1812     4.608939e+01,4.717095e+01,4.825627e+01,4.    
1813     5.382263e+01,5.497760e+01,5.615006e+01,5.    
1814     6.231082e+01,6.360849e+01,6.492971e+01,6.    
1815     7.191516e+01,7.339351e+01,7.490076e+01,7.    
1816     8.289618e+01,8.459176e+01,8.632091e+01,8.    
1817     9.548939e+01,9.743150e+01,9.941112e+01,1.    
1818     1.098962e+02,1.121162e+02,1.143792e+02,1.    
1819     1.263787e+02,1.289239e+02,1.315205e+02,1.    
1820     1.453293e+02,1.482676e+02,1.512688e+02,1.    
1821     1.672875e+02,1.707082e+02,1.742066e+02,1.    
1822     1.929462e+02,1.969619e+02,2.010732e+02,2.    
1823     2.231694e+02,2.279189e+02,2.327864e+02,2.    
1824     2.590228e+02,2.646774e+02,2.704774e+02,2.    
1825     3.018199e+02,3.085907e+02,3.155410e+02,3.    
1826     3.531820e+02,3.613300e+02,3.696996e+02,3.    
1827     4.151151e+02,4.249638e+02,4.350862e+02,4.    
1828     4.901079e+02,5.020589e+02,5.143485e+02,5.    
1829     5.812541e+02,5.958075e+02,6.107802e+02,6.    
1830     6.924079e+02,7.101869e+02,7.284863e+02,7.    
1831     8.283782e+02,8.501616e+02,8.725913e+02,8.    
1832     9.951752e+02,1.021937e+03,1.049502e+03,1.    
1833 static const G4double P211[nE]={                 
1834     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1835     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1836     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1837     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1838     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1839     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1840     2.006524e-02,3.034536e-02,4.382582e-02,6.    
1841     2.168665e-01,2.939658e-01,4.003237e-01,5.    
1842     1.745343e+00,2.326101e+00,3.082613e+00,4.    
1843     1.105116e+01,1.385719e+01,1.716964e+01,2.    
1844     3.837462e+01,4.232892e+01,4.595438e+01,4.    
1845     5.993661e+01,6.209595e+01,6.409989e+01,6.    
1846     7.249103e+01,7.394773e+01,7.535948e+01,7.    
1847     8.206649e+01,8.338723e+01,8.471403e+01,8.    
1848     9.158348e+01,9.303155e+01,9.451368e+01,9.    
1849     1.026050e+02,1.043980e+02,1.062640e+02,1.    
1850     1.169368e+02,1.193933e+02,1.219754e+02,1.    
1851     1.370529e+02,1.405572e+02,1.442423e+02,1.    
1852     1.656389e+02,1.705636e+02,1.757217e+02,1.    
1853     2.053134e+02,2.121426e+02,2.193769e+02,2.    
1854     2.655272e+02,2.782318e+02,2.928987e+02,3.    
1855     4.132549e+02,4.503497e+02,4.929799e+02,5.    
1856     8.116615e+02,9.025013e+02,1.005029e+03,1.    
1857     1.749612e+03,1.958046e+03,2.190894e+03,2.    
1858     3.767744e+03,4.167523e+03,4.592970e+03,5.    
1859     7.015188e+03,7.539806e+03,8.072749e+03,8.    
1860     1.083504e+04,1.140796e+04,1.198979e+04,1.    
1861     1.509210e+04,1.576553e+04,1.646199e+04,1.    
1862     2.037842e+04,2.126663e+04,2.219590e+04,2.    
1863     2.755151e+04,2.878401e+04,3.007731e+04,3.    
1864     3.757777e+04,3.931185e+04,4.113396e+04,4.    
1865     5.173679e+04,5.419304e+04,5.677458e+04,5.    
1866     7.178955e+04,7.526394e+04,7.891383e+04,8.    
1867     1.001189e+05,1.050235e+05,1.101763e+05,1.    
1868     1.401464e+05,1.470904e+05,1.543917e+05,1.    
1869     1.969859e+05,2.068858e+05,2.173072e+05,2.    
1870     2.783297e+05,2.925628e+05,3.075639e+05,3.    
1871     3.957283e+05,4.163613e+05,4.381322e+05,4.    
1872     5.665162e+05,5.966524e+05,6.284822e+05,6.    
1873     8.167405e+05,8.610475e+05,9.078852e+05,9.    
1874     1.185622e+06,1.251138e+06,1.320447e+06,1.    
1875     1.732364e+06,1.829725e+06,1.932793e+06,2.    
1876     2.546550e+06,2.691874e+06,2.845807e+06,3.    
1877     3.764073e+06,3.981840e+06,4.212628e+06,4.    
1878     5.591541e+06,5.919014e+06,6.266233e+06,6.    
1879     8.343775e+06,8.837797e+06,9.361836e+06,9.    
1880     1.250148e+07,1.324894e+07,1.404213e+07,1.    
1881     1.880006e+07,1.993402e+07,2.113780e+07,2.    
1882 static const G4double P012[nE]={                 
1883     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1884     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1885     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1886     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1887     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1888     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1889     2.064930e-03,3.284750e-03,4.790517e-03,6.    
1890     2.029405e-02,2.612967e-02,3.360287e-02,4.    
1891     1.048143e-01,1.292053e-01,1.582020e-01,1.    
1892     3.681689e-01,4.129651e-01,4.533231e-01,4.    
1893     5.846749e-01,6.003010e-01,6.136180e-01,6.    
1894     6.569272e-01,6.625246e-01,6.674590e-01,6.    
1895     6.854375e-01,6.881350e-01,6.906373e-01,6.    
1896     7.011385e-01,7.029615e-01,7.047201e-01,7.    
1897     7.127865e-01,7.142885e-01,7.157619e-01,7.    
1898     7.227802e-01,7.241260e-01,7.254552e-01,7.    
1899     7.318649e-01,7.330998e-01,7.343184e-01,7.    
1900     7.401512e-01,7.412621e-01,7.423534e-01,7.    
1901     7.475050e-01,7.484740e-01,7.494234e-01,7.    
1902     7.539251e-01,7.548030e-01,7.556892e-01,7.    
1903     7.607266e-01,7.619819e-01,7.633625e-01,7.    
1904     7.724363e-01,7.746990e-01,7.771113e-01,7.    
1905     7.915675e-01,7.949937e-01,7.986237e-01,8.    
1906     8.202389e-01,8.253131e-01,8.306399e-01,8.    
1907     8.603917e-01,8.667190e-01,8.730602e-01,8.    
1908     9.032762e-01,9.087677e-01,9.140286e-01,9.    
1909     9.369096e-01,9.408651e-01,9.446450e-01,9.    
1910     9.614123e-01,9.644302e-01,9.673624e-01,9.    
1911     9.810219e-01,9.835976e-01,9.861323e-01,9.    
1912     9.982862e-01,1.000627e+00,1.002941e+00,1.    
1913     1.014139e+00,1.016308e+00,1.018455e+00,1.    
1914     1.028859e+00,1.030874e+00,1.032867e+00,1.    
1915     1.042502e+00,1.044363e+00,1.046202e+00,1.    
1916     1.055077e+00,1.056790e+00,1.058483e+00,1.    
1917     1.066661e+00,1.068243e+00,1.069807e+00,1.    
1918     1.077388e+00,1.078859e+00,1.080316e+00,1.    
1919     1.087404e+00,1.088785e+00,1.090155e+00,1.    
1920     1.096843e+00,1.098150e+00,1.099449e+00,1.    
1921     1.105813e+00,1.107062e+00,1.108303e+00,1.    
1922     1.114404e+00,1.115605e+00,1.116799e+00,1.    
1923     1.122687e+00,1.123848e+00,1.125005e+00,1.    
1924     1.130718e+00,1.131848e+00,1.132973e+00,1.    
1925     1.138545e+00,1.139649e+00,1.140750e+00,1.    
1926     1.146208e+00,1.147291e+00,1.148371e+00,1.    
1927     1.153738e+00,1.154805e+00,1.155870e+00,1.    
1928     1.161166e+00,1.162220e+00,1.163272e+00,1.    
1929     1.168515e+00,1.169560e+00,1.170604e+00,1.    
1930     1.175809e+00,1.176847e+00,1.177885e+00,1.    
1931 static const G4double P112[nE]={                 
1932     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1933     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1934     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1935     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1936     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1937     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1938     1.353042e-02,2.197839e-02,3.274774e-02,4.    
1939     1.567549e-01,2.073542e-01,2.742402e-01,3.    
1940     9.845531e-01,1.248518e+00,1.572353e+00,1.    
1941     4.163437e+00,4.766978e+00,5.327878e+00,5.    
1942     7.314181e+00,7.575938e+00,7.806026e+00,8.    
1943     8.620696e+00,8.737148e+00,8.843020e+00,8.    
1944     9.263911e+00,9.333538e+00,9.400141e+00,9.    
1945     9.705909e+00,9.764235e+00,9.822253e+00,9.    
1946     1.011373e+01,1.017325e+01,1.023345e+01,1.    
1947     1.054767e+01,1.061369e+01,1.068090e+01,1.    
1948     1.103609e+01,1.111104e+01,1.118729e+01,1.    
1949     1.158702e+01,1.167041e+01,1.175485e+01,1.    
1950     1.219132e+01,1.228126e+01,1.237208e+01,1.    
1951     1.284374e+01,1.294447e+01,1.304928e+01,1.    
1952     1.370470e+01,1.388274e+01,1.408458e+01,1.    
1953     1.554557e+01,1.594212e+01,1.637787e+01,1.    
1954     1.924953e+01,1.999143e+01,2.080156e+01,2.    
1955     2.610582e+01,2.746321e+01,2.893183e+01,3.    
1956     3.793576e+01,4.002641e+01,4.218581e+01,4.    
1957     5.344347e+01,5.568452e+01,5.789724e+01,6.    
1958     6.841150e+01,7.040530e+01,7.236897e+01,7.    
1959     8.189187e+01,8.377082e+01,8.565227e+01,8.    
1960     9.524359e+01,9.722429e+01,9.923302e+01,1.    
1961     1.097786e+02,1.120016e+02,1.142666e+02,1.    
1962     1.262679e+02,1.288124e+02,1.314077e+02,1.    
1963     1.451783e+02,1.480970e+02,1.510722e+02,1.    
1964     1.668202e+02,1.701491e+02,1.735398e+02,1.    
1965     1.914545e+02,1.952384e+02,1.990926e+02,2.    
1966     2.194793e+02,2.237936e+02,2.281918e+02,2.    
1967     2.515305e+02,2.564865e+02,2.615454e+02,2.    
1968     2.884953e+02,2.942405e+02,3.001125e+02,3.    
1969     3.315166e+02,3.382361e+02,3.451123e+02,3.    
1970     3.820173e+02,3.899398e+02,3.980558e+02,4.    
1971     4.417506e+02,4.511576e+02,4.608034e+02,4.    
1972     5.128744e+02,5.241128e+02,5.356456e+02,5.    
1973     5.980497e+02,6.115475e+02,6.254086e+02,6.    
1974     7.005664e+02,7.168537e+02,7.335900e+02,7.    
1975     8.245030e+02,8.442381e+02,8.645283e+02,8.    
1976     9.749276e+02,9.989293e+02,1.023618e+03,1.    
1977     1.158151e+03,1.187440e+03,1.217582e+03,1.    
1978     1.382046e+03,1.417896e+03,1.454806e+03,1.    
1979     1.656447e+03,1.700453e+03,1.745775e+03,1.    
1980 static const G4double P212[nE]={                 
1981     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1982     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1983     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1984     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1985     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1986     0.000000e+00,0.000000e+00,0.000000e+00,0.    
1987     8.871441e-02,1.472267e-01,2.242564e-01,3.    
1988     1.218055e+00,1.656827e+00,2.255513e+00,3.    
1989     9.354437e+00,1.221137e+01,1.582824e+01,2.    
1990     4.786507e+01,5.599726e+01,6.379332e+01,7.    
1991     9.388937e+01,9.827447e+01,1.022502e+02,1.    
1992     1.176048e+02,1.200277e+02,1.222995e+02,1.    
1993     1.321723e+02,1.339697e+02,1.357426e+02,1.    
1994     1.446633e+02,1.465296e+02,1.484437e+02,1.    
1995     1.589966e+02,1.613556e+02,1.638152e+02,1.    
1996     1.779103e+02,1.811492e+02,1.845488e+02,1.    
1997     2.042690e+02,2.088184e+02,2.135892e+02,2.    
1998     2.410367e+02,2.472972e+02,2.538316e+02,2.    
1999     2.908809e+02,2.992293e+02,3.079187e+02,3.    
2000     3.574312e+02,3.689892e+02,3.813869e+02,3.    
2001     4.668260e+02,4.920818e+02,5.215910e+02,5.    
2002     7.572662e+02,8.267702e+02,9.054874e+02,9.    
2003     1.477006e+03,1.637665e+03,1.818480e+03,2.    
2004     3.122568e+03,3.485707e+03,3.890646e+03,4.    
2005     6.620672e+03,7.311514e+03,8.046928e+03,8.    
2006     1.224914e+04,1.316377e+04,1.409450e+04,1.    
2007     1.893518e+04,1.994025e+04,2.096045e+04,2.    
2008     2.637920e+04,2.754912e+04,2.875647e+04,3.    
2009     3.550403e+04,3.702724e+04,3.861931e+04,4.    
2010     4.778667e+04,4.989812e+04,5.211519e+04,5.    
2011     6.500194e+04,6.798661e+04,7.112399e+04,7.    
2012     8.938526e+04,9.361335e+04,9.805527e+04,1.    
2013     1.238442e+05,1.297996e+05,1.360510e+05,1.    
2014     1.722801e+05,1.806398e+05,1.894154e+05,1.    
2015     2.403312e+05,2.521018e+05,2.644689e+05,2.    
2016     3.364554e+05,3.531535e+05,3.707199e+05,3.    
2017     4.733811e+05,4.972852e+05,5.224649e+05,5.    
2018     6.702071e+05,7.047325e+05,7.411446e+05,7.    
2019     9.555652e+05,1.005835e+06,1.058907e+06,1.    
2020     1.372430e+06,1.446140e+06,1.524033e+06,1.    
2021     1.985444e+06,2.094187e+06,2.209195e+06,2.    
2022     2.892091e+06,3.053375e+06,3.224069e+06,3.    
2023     4.239756e+06,4.480087e+06,4.734597e+06,5.    
2024     6.251840e+06,6.611444e+06,6.992476e+06,7.    
2025     9.267762e+06,9.807836e+06,1.038038e+07,1.    
2026     1.380440e+07,1.461824e+07,1.548140e+07,1.    
2027     2.065053e+07,2.188068e+07,2.318592e+07,2.    
2028     3.101232e+07,3.287695e+07,3.485616e+07,3.    
2029 static const G4double P013[nE]={                 
2030     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2031     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2032     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2033     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2034     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2035     0.000000e+00,0.000000e+00,3.486588e-04,9.    
2036     5.154692e-03,7.030543e-03,9.343593e-03,1.    
2037     3.270264e-02,4.120096e-02,5.179545e-02,6.    
2038     1.402537e-01,1.647614e-01,1.907608e-01,2.    
2039     3.291280e-01,3.571385e-01,3.844637e-01,4.    
2040     4.931270e-01,5.074779e-01,5.198558e-01,5.    
2041     5.612855e-01,5.668161e-01,5.717392e-01,5.    
2042     5.901599e-01,5.930068e-01,5.956710e-01,5.    
2043     6.071155e-01,6.091467e-01,6.111183e-01,6.    
2044     6.203001e-01,6.220315e-01,6.237350e-01,6.    
2045     6.318973e-01,6.334668e-01,6.350169e-01,6.    
2046     6.424746e-01,6.439056e-01,6.453153e-01,6.    
2047     6.520247e-01,6.532950e-01,6.545406e-01,6.    
2048     6.603918e-01,6.614881e-01,6.625614e-01,6.    
2049     6.676680e-01,6.686761e-01,6.697018e-01,6.    
2050     6.757071e-01,6.772293e-01,6.789009e-01,6.    
2051     6.896938e-01,6.923421e-01,6.951548e-01,6.    
2052     7.118958e-01,7.158503e-01,7.200374e-01,7.    
2053     7.449331e-01,7.507700e-01,7.568949e-01,7.    
2054     7.910682e-01,7.983318e-01,8.056115e-01,8.    
2055     8.403238e-01,8.466393e-01,8.526920e-01,8.    
2056     8.790371e-01,8.835929e-01,8.879462e-01,8.    
2057     9.072423e-01,9.107112e-01,9.140799e-01,9.    
2058     9.297517e-01,9.327034e-01,9.356072e-01,9.    
2059     9.495258e-01,9.522065e-01,9.548567e-01,9.    
2060     9.676864e-01,9.701728e-01,9.726336e-01,9.    
2061     9.845594e-01,9.868688e-01,9.891530e-01,9.    
2062     1.000191e+00,1.002322e+00,1.004428e+00,1.    
2063     1.014583e+00,1.016542e+00,1.018478e+00,1.    
2064     1.027824e+00,1.029630e+00,1.031416e+00,1.    
2065     1.040068e+00,1.041746e+00,1.043408e+00,1.    
2066     1.051488e+00,1.053062e+00,1.054622e+00,1.    
2067     1.062238e+00,1.063727e+00,1.065205e+00,1.    
2068     1.072446e+00,1.073866e+00,1.075278e+00,1.    
2069     1.082215e+00,1.083579e+00,1.084937e+00,1.    
2070     1.091626e+00,1.092946e+00,1.094259e+00,1.    
2071     1.100747e+00,1.102030e+00,1.103308e+00,1.    
2072     1.109632e+00,1.110885e+00,1.112134e+00,1.    
2073     1.118326e+00,1.119555e+00,1.120780e+00,1.    
2074     1.126867e+00,1.128077e+00,1.129284e+00,1.    
2075     1.135289e+00,1.136484e+00,1.137677e+00,1.    
2076     1.143619e+00,1.144803e+00,1.145986e+00,1.    
2077     1.151884e+00,1.153061e+00,1.154237e+00,1.    
2078 static const G4double P113[nE]={                 
2079     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2080     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2081     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2082     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2083     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2084     0.000000e+00,0.000000e+00,1.995454e-03,5.    
2085     3.263520e-02,4.562653e-02,6.216957e-02,8.    
2086     2.489316e-01,3.226120e-01,4.174284e-01,5.    
2087     1.296801e+00,1.561942e+00,1.852224e+00,2.    
2088     3.551797e+00,3.929232e+00,4.309095e+00,4.    
2089     5.957193e+00,6.197595e+00,6.411465e+00,6.    
2090     7.191277e+00,7.306340e+00,7.411974e+00,7.    
2091     7.843477e+00,7.916960e+00,7.987876e+00,8.    
2092     8.321275e+00,8.386263e+00,8.451305e+00,8.    
2093     8.783185e+00,8.851798e+00,8.921398e+00,8.    
2094     9.286866e+00,9.363861e+00,9.442250e+00,9.    
2095     9.855482e+00,9.942335e+00,1.003054e+01,1.    
2096     1.049029e+01,1.058565e+01,1.068203e+01,1.    
2097     1.117772e+01,1.127947e+01,1.138215e+01,1.    
2098     1.191727e+01,1.203294e+01,1.215425e+01,1.    
2099     1.293604e+01,1.315194e+01,1.339632e+01,1.    
2100     1.513342e+01,1.559757e+01,1.610563e+01,1.    
2101     1.943082e+01,2.028712e+01,2.122158e+01,2.    
2102     2.733069e+01,2.889211e+01,3.058077e+01,3.    
2103     4.092260e+01,4.332263e+01,4.580166e+01,4.    
2104     5.873476e+01,6.131210e+01,6.385783e+01,6.    
2105     7.596408e+01,7.826048e+01,8.052200e+01,8.    
2106     9.148086e+01,9.364060e+01,9.580217e+01,9.    
2107     1.068062e+02,1.090760e+02,1.113772e+02,1.    
2108     1.234539e+02,1.259998e+02,1.285938e+02,1.    
2109     1.423440e+02,1.452603e+02,1.482350e+02,1.    
2110     1.640199e+02,1.673652e+02,1.707750e+02,1.    
2111     1.888164e+02,1.926284e+02,1.965104e+02,2.    
2112     2.170098e+02,2.213372e+02,2.257441e+02,2.    
2113     2.490408e+02,2.539681e+02,2.589904e+02,2.    
2114     2.856261e+02,2.912794e+02,2.970490e+02,3.    
2115     3.277707e+02,3.343169e+02,3.410068e+02,3.    
2116     3.767697e+02,3.844189e+02,3.922454e+02,4.    
2117     4.342361e+02,4.432473e+02,4.524776e+02,4.    
2118     5.021563e+02,5.128485e+02,5.238111e+02,5.    
2119     5.829747e+02,5.957407e+02,6.088401e+02,6.    
2120     6.797048e+02,6.950292e+02,7.107652e+02,7.    
2121     7.960713e+02,8.145545e+02,8.335460e+02,8.    
2122     9.366915e+02,9.590784e+02,9.820938e+02,1.    
2123     1.107302e+03,1.134519e+03,1.162515e+03,1.    
2124     1.315045e+03,1.348248e+03,1.382416e+03,1.    
2125     1.568829e+03,1.609460e+03,1.651289e+03,1.    
2126     1.879785e+03,1.929646e+03,1.980999e+03,2.    
2127 static const G4double P213[nE]={                 
2128     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2129     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2130     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2131     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2132     0.000000e+00,0.000000e+00,0.000000e+00,0.    
2133     0.000000e+00,0.000000e+00,1.142094e-02,3.    
2134     2.071675e-01,2.971476e-01,4.154746e-01,5.    
2135     1.911579e+00,2.550442e+00,3.399081e+00,4.    
2136     1.215838e+01,1.502709e+01,1.826833e+01,2.    
2137     3.918699e+01,4.427325e+01,4.955437e+01,5.    
2138     7.460021e+01,7.862768e+01,8.232331e+01,8.    
2139     9.703032e+01,9.942442e+01,1.016911e+02,1.    
2140     1.118190e+02,1.137160e+02,1.156037e+02,1.    
2141     1.253353e+02,1.274148e+02,1.295607e+02,1.    
2142     1.415800e+02,1.442993e+02,1.471432e+02,1.    
2143     1.635387e+02,1.673159e+02,1.712805e+02,1.    
2144     1.942219e+02,1.994937e+02,2.050125e+02,2.    
2145     2.365773e+02,2.437362e+02,2.511944e+02,2.    
2146     2.932671e+02,3.027119e+02,3.125358e+02,3.    
2147     3.687193e+02,3.819927e+02,3.963408e+02,4.    
2148     4.983109e+02,5.289362e+02,5.646649e+02,6.    
2149     8.447731e+02,9.261237e+02,1.017903e+03,1.    
2150     1.679620e+03,1.865051e+03,2.073617e+03,2.    
2151     3.575541e+03,3.993264e+03,4.458871e+03,4.    
2152     7.594502e+03,8.387576e+03,9.231842e+03,1.    
2153     1.405955e+04,1.511144e+04,1.618224e+04,1.    
2154     2.175592e+04,2.291353e+04,2.408848e+04,2.    
2155     3.032420e+04,3.166894e+04,3.305604e+04,3.    
2156     4.079722e+04,4.254275e+04,4.436670e+04,4.    
2157     5.486506e+04,5.728307e+04,5.982226e+04,6.    
2158     7.458700e+04,7.800776e+04,8.160383e+04,8.    
2159     1.025363e+05,1.073824e+05,1.124731e+05,1.    
2160     1.420175e+05,1.488372e+05,1.559945e+05,1.    
2161     1.974500e+05,2.070104e+05,2.170444e+05,2.    
2162     2.752272e+05,2.886704e+05,3.027922e+05,3.    
2163     3.849474e+05,4.039947e+05,4.240291e+05,4.    
2164     5.410568e+05,5.682940e+05,5.969807e+05,6.    
2165     7.652278e+05,8.045300e+05,8.459747e+05,8.    
2166     1.089941e+06,1.147119e+06,1.207478e+06,1.    
2167     1.563934e+06,1.647715e+06,1.736241e+06,1.    
2168     2.260499e+06,2.384024e+06,2.514653e+06,2.    
2169     3.290130e+06,3.473241e+06,3.667023e+06,3.    
2170     4.819852e+06,5.092583e+06,5.381389e+06,5.    
2171     7.102774e+06,7.510697e+06,7.942907e+06,8.    
2172     1.052339e+07,1.113582e+07,1.178504e+07,1.    
2173     1.566710e+07,1.658970e+07,1.756817e+07,1.    
2174     2.342715e+07,2.482132e+07,2.630054e+07,2.    
2175     3.516924e+07,3.728200e+07,3.952452e+07,4.    
2176 static const G4double* P0[nN]={P00,P01,P02,P0    
2177 static const G4double* P1[nN]={P10,P11,P12,P1    
2178 static const G4double* P2[nN]={P20,P21,P22,P2    
2179 //                                               
2180                                                  
2181 G4ElectroNuclearCrossSection::G4ElectroNuclea    
2182 currentN(0), currentZ(0), lastZ(0),lastE(0),     
2183 mNeut(neutron_mass_c2), mProt(proton_mass_c2)    
2184 {                                                
2185     SetForAllAtomsAndEnergies(true);             
2186     //Initialize caches                          
2187     lastUsedCacheEl = new cacheEl_t();           
2188     nistmngr = G4NistManager::Instance();        
2189                                                  
2190     for (G4int i=0;i<120;i++)                    
2191     {                                            
2192         cache.push_back(nullptr);                
2193     }                                            
2194                                                  
2195 }                                                
2196                                                  
2197 G4ElectroNuclearCrossSection::~G4ElectroNucle    
2198 {                                                
2199      std::vector<cacheEl_t*>::iterator it = c    
2200      while ( it != cache.end() )  /* Loop che    
2201      {                                           
2202          if ( *it ) {                            
2203              delete[] (*it)->J1; (*it)->J1 =     
2204              delete[] (*it)->J2; (*it)->J2 =     
2205              delete[] (*it)->J3; (*it)->J3 =     
2206              delete *it;                         
2207          }                                       
2208      ++it;                                       
2209      }                                           
2210      cache.clear();                              
2211      delete lastUsedCacheEl;                     
2212 }                                                
2213                                                  
2214 G4double G4ElectroNuclearCrossSection::HighEn    
2215 {                                                
2216     return ha*(lE*lE-lEMa2)-ab*(lE-lEMa)-cd*(    
2217 }                                                
2218                                                  
2219 G4double G4ElectroNuclearCrossSection::HighEn    
2220 {                                                
2221     return poc*((lE-1.)*E-le1)-ab*(E-EMa)+cd1    
2222 }                                                
2223                                                  
2224 G4double G4ElectroNuclearCrossSection::HighEn    
2225 {                                                
2226     return ha*((lE-.5)*E2-leh)-hab*(E2-EMa2)+    
2227 }                                                
2228                                                  
2229 G4double G4ElectroNuclearCrossSection::DFun(G    
2230 {                                                
2231     G4double y=G4Exp(x-lastG-lmel);              
2232     G4double flux=lastG*(2.-y*(2.-y))-1.;        
2233     return (poc*(x-pos)+shd*G4Exp(-reg*x))*fl    
2234 }                                                
2235                                                  
2236 G4double G4ElectroNuclearCrossSection::Fun(G4    
2237 {                                                
2238     // Integrated PhoNuc cross section           
2239     G4double dlg1=lastG+lastG-1.;                
2240     G4double lgoe=lastG/lastE;                   
2241     G4double HE2=HighEnergyJ2(x, G4Exp(x));      
2242     return dlg1*HighEnergyJ1(x)-lgoe*(HE2+HE2    
2243 }                                                
2244                                                  
2245 void                                             
2246 G4ElectroNuclearCrossSection::CrossSectionDes    
2247 {                                                
2248     outFile << "G4ElectroNuclearCrossSection     
2249     << "cross section for e- and e+ interacti    
2250     << "cross sections are retrieved from a t    
2251     << "generated using the equivalent photon    
2252     << "this approximation real gammas are pr    
2253     << "ones generated at the electromagnetic    
2254     << "section set is valid for incident ele    
2255     << "all energies.\n";                        
2256 }                                                
2257                                                  
2258 G4bool G4ElectroNuclearCrossSection::IsElemen    
2259 {                                                
2260   return true;                                   
2261 }                                                
2262                                                  
2263 G4double G4ElectroNuclearCrossSection::GetEle    
2264 {                                                
2265     const G4double Energy = aPart->GetKinetic    
2266                                                  
2267     if (Energy <= sLowEnergyLimit || ZZ >= 12    
2268                                                  
2269     if(ZZ!=lastZ)      // This nucleus was no    
2270     {                                            
2271         lastE    = 0.;                           
2272         lastG    = 0.;                           
2273         lastZ = ZZ;                              
2274                                                  
2275         //key to search in cache                 
2276         if(!cache[ZZ]){                          
2277             lastUsedCacheEl->J1 = new G4doubl    
2278             lastUsedCacheEl->J2 = new G4doubl    
2279             lastUsedCacheEl->J3 = new G4doubl    
2280             G4double Aa = nistmngr->GetAtomic    
2281             G4int N = (G4int)Aa - ZZ;            
2282             lastUsedCacheEl->F  = GetFunction    
2283             lastUsedCacheEl->H  = alop*Aa*(1.    
2284             lastUsedCacheEl->TH = ThresholdEn    
2285             cacheEl_t* new_el = new cacheEl_t    
2286             cache[ZZ] = new_el;                  
2287         }                                        
2288         else                                     
2289         { //found in cache                       
2290             const cacheEl_t& el = *(cache[ZZ]    
2291             lastUsedCacheEl->F  = el.F;          
2292             lastUsedCacheEl->TH = el.TH;         
2293             lastUsedCacheEl->H  = el.H;          
2294             lastUsedCacheEl->J1 = el.J1;         
2295             lastUsedCacheEl->J2 = el.J2;         
2296             lastUsedCacheEl->J3 = el.J3;         
2297         }                                        
2298     }                                            
2299     else                                         
2300     { //current isotope is the same as previo    
2301         if ( lastE == Energy ) return lastSig    
2302     }                                            
2303     //End of optimization: now lastUsedCacheE    
2304                                                  
2305     // ============================== NOW Cal    
2306     lastE=Energy;                          //    
2307                                                  
2308     if ( Energy <= lastUsedCacheEl->TH ) // c    
2309     {                                            
2310         lastSig=0.;                              
2311         return 0.;                               
2312     }                                            
2313                                                  
2314     G4double lE=G4Log(Energy);                   
2315                                                  
2316     lastG=lE-lmel;                         //    
2317     G4double dlg1=lastG+lastG-1.;                
2318     G4double lgoe=lastG/lastE;                   
2319     if(lE<lEMa) // Linear fit is made explici    
2320   {                                              
2321         G4double shift=(lE-lEMi)/dlnE;           
2322         G4int    blast=static_cast<G4int>(shi    
2323         if(blast<0)   blast=0;                   
2324         if(blast>=mLL) blast=mLL-1;              
2325         shift-=blast;                            
2326         lastL=blast+1;                           
2327         G4double YNi=dlg1*lastUsedCacheEl->J1    
2328         G4double YNj=dlg1*lastUsedCacheEl->J1    
2329         lastSig= YNi+shift*(YNj-YNi);            
2330         if(lastSig>YNj)lastSig=YNj;              
2331     }                                            
2332     else                                         
2333   {                                              
2334         lastL=mLL;                               
2335                                                  
2336         G4double term1=lastUsedCacheEl->J1[mL    
2337                                                  
2338         G4double term2=lastUsedCacheEl->J2[mL    
2339                                                  
2340         G4double En2 = Energy*Energy;            
2341         G4double term3=lastUsedCacheEl->J3[mL    
2342                                                  
2343         lastSig=dlg1*term1-lgoe*(term2+term2-    
2344     }                                            
2345                                                  
2346     if(lastSig<0.) lastSig = 0.;                 
2347                                                  
2348     return lastSig*millibarn;                    
2349 }                                                
2350                                                  
2351 // Gives the threshold energy for different n    
2352 G4double G4ElectroNuclearCrossSection::Thresh    
2353 {                                                
2354     // ---------                                 
2355                                                  
2356     G4int Aa=Z+N;                                
2357     if(Aa<1) return infEn;                       
2358     else if(Aa==1) return 134.9766; // Pi0 th    
2359                                                  
2360     G4double mT= 0.;                             
2361     if(G4NucleiProperties::IsInStableTable(Aa    
2362     // If it is not in the Table of Stable Nu    
2363     else return infEn;                           
2364     // ---------                                 
2365     G4double mP= infEn;                          
2366     //if(Z) mP= G4QPDGCode(111).GetNuclMass(Z    
2367     if(Z&&G4NucleiProperties::IsInStableTable    
2368     else return infEn;                           
2369     G4double mN= infEn;                          
2370     //if(N) mN= G4QPDGCode(111).GetNuclMass(Z    
2371     if(N&&G4NucleiProperties::IsInStableTable    
2372     else return infEn;                           
2373     G4double dP= mP+mProt-mT;                    
2374     G4double dN= mN+mNeut-mT;                    
2375     if(dP<dN)dN=dP;                              
2376     return dN;                                   
2377 }                                                
2378                                                  
2379                                                  
2380 // Calculate the functions for the G4Log(A)      
2381 G4int G4ElectroNuclearCrossSection::GetFuncti    
2382 {                                                
2383     // --------------------------------          
2384     G4int r=-1;                             /    
2385                                                  
2386     //AR-24Apr2018 Switch to allow transurani    
2387     const G4bool isHeavyElementAllowed = true    
2388     if(a<=.9999 || ( !isHeavyElementAllowed &    
2389     {                                            
2390         G4cout<<"***G4ElectroNuclearCrossSect    
2391         return r;                                
2392     }                                            
2393     G4int iA=static_cast<G4int>(a+.499);    /    
2394     G4double ai=iA;                              
2395     if(a!=ai) a=ai;                              
2396     for(G4int i=0; i<nN; i++)                    
2397     {                                            
2398         if(std::abs(a-A[i])<.0005)               
2399         {                                        
2400             for(G4int k=0; k<nE; k++)            
2401             {                                    
2402                 xx[k]=P0[i][k];                  
2403                 yy[k]=P1[i][k];                  
2404                 zz[k]=P2[i][k];                  
2405             }                                    
2406             r=LL[i];                             
2407         }                                        
2408         if(r<0)                                  
2409         {                                        
2410             G4int k=0;                           
2411             for(k=1; k<nN; k++)if(a<A[k]) bre    
2412             if(k<1) k=1;                         
2413             if(k>=nN) k=nN-1;                    
2414             G4int     k1=k-1;                    
2415             G4double  xi=A[k1];                  
2416             G4double   b=(a-xi)/(A[k]-xi);       
2417             for(G4int q=0; q<nE; q++)            
2418             {                                    
2419                 xi=P0[k1][q];                    
2420                 xx[q]=xi+(P0[k][q]-xi)*b;        
2421                 G4double yi=P1[k1][q];           
2422                 yy[q]=yi+(P1[k][q]-yi)*b;        
2423                 G4double zi=P2[k1][q];           
2424                 zz[q]=zi+(P2[k][q]-zi)*b;        
2425             }                                    
2426             r=LL[k];                             
2427             if(LL[k1]<r) r=LL[k1];               
2428         }                                        
2429     }                                            
2430     return r;                                    
2431 }                                                
2432                                                  
2433 G4double G4ElectroNuclearCrossSection::GetEqu    
2434 {                                                
2435     if(lastSig <= 0.0) { return 0.0; }  // VI    
2436     G4double phLE = 0.;                          
2437     G4double Y[nE] = {0.0};                      
2438                                                  
2439     G4double lastLE=lastG+lmel;   // recover     
2440     G4double dlg1=lastG+lastG-1.;                
2441     G4double lgoe=lastG/lastE;                   
2442     for (G4int i=lastUsedCacheEl->F;i<=lastL;    
2443         Y[i] = dlg1*lastUsedCacheEl->J1[i]-lg    
2444         if(Y[i] < 0.0) { Y[i] = 0.0; }           
2445     }                                            
2446     // Tempory IF of H.P.: delete it if the *    
2447     // show up M.K.                              
2448     if(lastSig>0.99*Y[lastL] && lastL<mLL &&     
2449     {                                            
2450         G4cerr << "*HP*G4ElNucCS::GetEqPhotE:    
2451         << ",l=" << lastL << ">" << mLL << G4    
2452         if(lastSig <= 0.0) { return 0.0; }  /    
2453     }                                            
2454     G4double ris = lastSig*G4UniformRand(); /    
2455     // is in the funct. region                   
2456                                                  
2457     if (ris < Y[lastL]) {               // Se    
2458         G4int j = lastUsedCacheEl->F;            
2459         G4double Yj = Y[j];               //     
2460         while (ris > Yj && j < lastL) {   //     
2461             j++;                                 
2462             Yj = Y[j];                      /    
2463         }                                        
2464         G4int j1 = j-1;                          
2465         G4double Yi = Y[j1];              //     
2466         phLE = lEMi + (j1 + (ris-Yi)/(Yj-Yi)     
2467     } else {                            // Se    
2468         if (lastL < mLL) G4cerr << "**G4EleNu    
2469             << lastSig << ",Y=" << Y[lastL] <    
2470         G4double f = (ris-Y[lastL])/lastUsedC    
2471         phLE=SolveTheEquation(f);           /    
2472     }                                            
2473                                                  
2474     if (phLE>lastLE) {                           
2475         G4cerr << "***G4ElectroNuclearCS::Get    
2476         << currentZ << ", lpE" << phLE << ">l    
2477         << lastSig << ",rndSig=" << ris << ",    
2478         << lastL << ",Y=" << Y[lastL] << G4en    
2479         if(lastLE<7.2) phLE=G4Log(G4Exp(lastL    
2480         else phLE=7.;                            
2481     }                                            
2482     return G4Exp(phLE);                          
2483 }                                                
2484                                                  
2485                                                  
2486 G4double G4ElectroNuclearCrossSection::SolveT    
2487 {                                                
2488     G4double lastLE=lastG+lmel;                  
2489     G4double topLim=lastLE-.001;                 
2490     G4double rE=EMa/G4Exp(lastLE);               
2491     G4double x=lEMa+f/phte/(lastG*(2.-rE*(2.-    
2492     if(x>topLim) x=topLim;                       
2493     for(G4int i=0; i<imax; i++)                  
2494     {                                            
2495         G4double fx=Fun(x);                      
2496         G4double df=DFun(x);                     
2497         G4double d=(f-fx)/df;                    
2498         x=x+d;                                   
2499         if(x>=lastLE)                            
2500         {                                        
2501             G4cerr<<"*G4ElNCS::SolveTheEq:*Co    
2502             <<",fx="<<fx<<",df="<<df<<",A(Z="    
2503             x=topLim;                            
2504         }                                        
2505         if(std::abs(d)<eps) break;               
2506         if(i+1>=imax) G4cerr<<"*G4ElNucCS::So    
2507             <<lastLE<<",Z="<<currentZ<<", N="    
2508     }                                            
2509     return x;                                    
2510 }                                                
2511                                                  
2512 G4double G4ElectroNuclearCrossSection::GetEqu    
2513 {                                                
2514     if(lastG <= 0.0 || lastE <= 0.0) { return    
2515     if(lastSig <= 0.0) { return 0.0; }  // VI    
2516     G4double y=nu/lastE;                    /    
2517     if(y>=1.-1./(lastG+lastG)) return 0.;   /    
2518     G4double y2=y*y;                        /    
2519     G4double ye=1.-y;                       /    
2520     G4double Qi2=mel2*y2/ye;                /    
2521     G4double Qa2=4*lastE*lastE*ye;          /    
2522     G4double iar=Qi2/Qa2;                   /    
2523     G4double Dy=ye+.5*y2;                   /    
2524     G4double Py=ye/Dy;                      /    
2525     G4double ePy=1.-G4Exp(Py);                   
2526     G4double Uy=Py*(1.-iar);                /    
2527     G4double Fy=(ye+ye)*(1.+ye)*iar/y2;     /    
2528     G4double fr=iar/(1.-ePy*iar);           /    
2529     if(Fy<=-fr)                                  
2530     {                                            
2531         return 0.;                               
2532     }                                            
2533     G4double LyQa2=G4Log(Fy+fr);                 
2534     G4bool cond=true;                            
2535     G4int maxTry=3;                              
2536     G4int cntTry=0;                              
2537     G4double Q2=Qi2;                             
2538     while(cond&&cntTry<maxTry)             //    
2539     {                                            
2540         G4double R=G4UniformRand();              
2541         Q2=Qi2*(ePy+1./(G4Exp(R*LyQa2-(1.-R)*    
2542         cntTry++;                                
2543         cond = Q2>1878.*nu;                      
2544     }                                            
2545     if(Q2<Qi2)                                   
2546     {                                            
2547         return Qi2;                              
2548     }                                            
2549     if(Q2>Qa2)                                   
2550     {                                            
2551         return Qa2;                              
2552     }                                            
2553     return Q2;                                   
2554 }                                                
2555                                                  
2556 G4double G4ElectroNuclearCrossSection::GetVir    
2557 {                                                
2558     if(nu <= 0.0 || Q2 <= 0.0) { return 0.0;     
2559     //G4double x=Q2/dM/nu;                  /    
2560     G4double K=nu-Q2/dM;                    /    
2561     if(K <= 0.) // VI                            
2562     {                                            
2563         return 0.;                               
2564     }                                            
2565     G4double lK=G4Log(K);                        
2566     G4double x=1.-K/nu;                     /    
2567     G4double GD=1.+Q2/Q02;                  /    
2568     G4double b=G4Exp(bp*(lK-blK0));              
2569     G4double c=G4Exp(cp*(lK-clK0));              
2570     G4double r=.5*G4Log(Q2+nu*nu)-lK;            
2571     G4double ef=G4Exp(r*(b-c*r*r));              
2572     return (1.-x)*ef/GD/GD;                      
2573 }                                                
2574