| Geant4 Cross Reference |
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 // The lust update: M.V. Kossov, CERN/ITEP(Mos
27 //
28 // G4 Physics class: G4PhotoNuclearCrossSectio
29 // Created: M.V. Kossov, CERN/ITEP(Moscow), 10
30 // The last update: M.V. Kossov, CERN/ITEP (Mo
31
32
33 #include <iostream>
34 #include <fstream>
35
36 #include "G4PhotoNuclearCrossSection.hh"
37 #include "G4SystemOfUnits.hh"
38
39 static const G4double THmin=2.; // mi
40 static const G4double dE=1.; // st
41 static const G4int nL=105; // A#
42 // (each MeV from 2 to 106)
43 static const G4double Emin=THmin+(nL-1)*dE; //
44 static const G4double Emax=50000.; // ma
45 static const G4int nH=224; // A#
46 static const G4double milE=G4Log(Emin); // Lo
47 // the HighE part
48 static const G4double malE=G4Log(Emax); // Hi
49 // (each 2.75 percent)
50 static const G4double dlE=(malE-milE)/(nH-1);
51 // in the HighE part
52 //
53 //static const G4double shd=1.075-.0023*G4Log(
54 static const G4double shd=1.0734;
55 static const G4double shc=0.072;
56 static const G4double poc=0.0375;
57 static const G4double pos=16.5;
58 static const G4double reg=.11;
59
60 static const G4double infEn = 9.e27;
61
62 static const G4int nLA=49; // A#
63 static const G4double LA[nLA]={
64 2., 4., 6., 7., 9., 12., 1
65 24., 27., 28., 32., 34., 40.,
66 63.5, 65.4, 76., 82., 107.9, 112.4, 1
67 165., 168., 174., 178., 180., 181.,
68 209., 232., 235., 238., 239.};
69
70 static const G4int nHA=14; // A#of HResonanc
71 static const G4double HA[nHA]={
72 1.,2.,3.,4.,6.,7.,9.,12.,16.,27.,63.5,118.
73
74 // If the cross section approximation formula
75
76 static const G4double SL0[nL]={
77 7.094260e-01,1.532987e+00,2.449381e+00,2.7
78 2.128172e+00,1.780549e+00,1.506934e+00,1.2
79 9.953850e-01,8.879274e-01,7.995356e-01,7.2
80 6.104038e-01,5.645786e-01,5.247229e-01,4.8
81 4.315429e-01,4.070560e-01,3.850576e-01,3.6
82 3.307971e-01,3.158133e-01,3.020711e-01,2.8
83 2.669563e-01,2.569336e-01,2.476099e-01,2.3
84 2.231848e-01,2.160475e-01,2.093390e-01,2.0
85 1.914383e-01,1.861152e-01,1.810725e-01,1.7
86 1.674254e-01,1.633120e-01,1.593914e-01,1.5
87 1.486616e-01,1.453926e-01,1.422615e-01,1.3
88 1.336147e-01,1.309573e-01,1.284017e-01,1.2
89 1.212914e-01,1.190904e-01,1.169666e-01,1.1
90 1.110206e-01,1.091688e-01,1.073770e-01,1.0
91 1.023336e-01,1.007548e-01,9.922335e-02,9.7
92 9.489316e-02,9.353161e-02,9.220814e-02,9.0
93 8.845106e-02,8.726514e-02,8.611027e-02,8.4
94 8.282039e-02,8.177841e-02,8.076208e-02,7.9
95 7.785794e-02,7.693536e-02,7.603421e-02,7.5
96 7.345216e-02,7.262971e-02,7.182534e-02,7.1
97 6.951498e-02,6.877732e-02,6.805505e-02,6.7
98
99 static const G4double SL1[nL]={
100 2.017310e-04,9.866847e-04,3.081371e-03,7.4
101 2.873865e-02,4.915763e-02,7.909925e-02,1.2
102 2.563778e-01,3.574504e-01,4.874600e-01,6.5
103 1.109763e+00,1.413389e+00,1.768398e+00,2.1
104 2.960166e+00,3.267863e+00,3.467401e+00,3.5
105 3.480379e+00,3.369351e+00,3.235278e+00,3.0
106 2.793315e+00,2.649362e+00,2.511013e+00,2.3
107 2.136980e+00,2.026587e+00,1.923127e+00,1.8
108 1.651181e+00,1.572129e+00,1.498250e+00,1.4
109 1.304111e+00,1.247474e+00,1.194383e+00,1.1
110 1.053842e+00,1.012488e+00,9.735481e-01,9.3
111 8.695413e-01,8.386509e-01,8.094332e-01,7.8
112 7.307199e-01,7.071392e-01,6.847445e-01,6.6
113 6.239481e-01,6.055964e-01,5.881058e-01,5.7
114 5.403069e-01,5.257844e-01,5.119013e-01,4.9
115 4.737459e-01,4.620893e-01,4.509173e-01,4.4
116 4.200659e-01,4.105971e-01,4.015025e-01,3.9
117 3.762885e-01,3.685206e-01,3.610467e-01,3.5
118 3.402609e-01,3.338384e-01,3.276510e-01,3.2
119 3.104042e-01,3.050648e-01,2.999168e-01,2.9
120 2.855516e-01,2.811014e-01,2.768107e-01,2.7
121
122 static const G4double SL2[nL]={
123 4.776434e-03,2.412116e-02,7.595870e-02,1.8
124 6.466818e-01,9.877908e-01,1.324697e+00,1.5
125 1.715038e+00,1.669943e+00,1.614318e+00,1.5
126 1.751419e+00,1.984310e+00,2.321229e+00,2.7
127 3.366107e+00,3.454042e+00,3.371971e+00,3.1
128 2.626333e+00,2.354348e+00,2.103834e+00,1.8
129 1.507070e+00,1.354935e+00,1.221899e+00,1.1
130 9.132844e-01,8.339468e-01,7.637380e-01,7.0
131 5.964196e-01,5.520797e-01,5.122769e-01,4.7
132 4.149083e-01,3.884171e-01,3.643543e-01,3.4
133 3.042291e-01,2.875255e-01,2.722111e-01,2.5
134 2.333017e-01,2.223151e-01,2.121693e-01,2.0
135 1.860586e-01,1.785946e-01,1.716644e-01,1.6
136 1.536595e-01,1.484670e-01,1.436278e-01,1.3
137 1.309760e-01,1.273073e-01,1.238809e-01,1.2
138 1.148962e-01,1.122861e-01,1.098477e-01,1.0
139 1.034596e-01,1.016087e-01,9.988340e-02,9.8
140 9.539245e-02,9.410337e-02,9.290919e-02,9.1
141 8.984997e-02,8.899108e-02,8.820664e-02,8.7
142 8.627010e-02,8.575488e-02,8.530112e-02,8.4
143 8.429132e-02,8.406729e-02,8.389779e-02,8.3
144
145 static const G4double SL3[nL]={
146 1.375991e-04,6.420490e-04,2.009594e-03,5.0
147 2.408187e-02,5.091978e-02,1.151175e-01,2.9
148 1.635125e+00,1.931572e+00,2.185333e+00,2.7
149 3.632210e+00,3.708366e+00,3.594398e+00,3.4
150 3.149899e+00,3.091216e+00,3.075568e+00,3.0
151 3.162196e+00,3.198373e+00,3.224873e+00,3.2
152 3.212852e+00,3.176382e+00,3.125768e+00,3.0
153 2.911700e+00,2.827186e+00,2.739409e+00,2.6
154 2.470840e+00,2.383021e+00,2.297216e+00,2.2
155 2.055557e+00,1.980911e+00,1.909346e+00,1.8
156 1.712877e+00,1.653235e+00,1.596357e+00,1.5
157 1.441186e+00,1.394230e+00,1.349469e+00,1.3
158 1.227254e+00,1.190196e+00,1.154820e+00,1.1
159 1.057908e+00,1.028411e+00,1.000196e+00,9.7
160 9.225887e-01,8.988681e-01,8.761312e-01,8.5
161 8.133424e-01,7.940736e-01,7.755685e-01,7.5
162 7.242788e-01,7.084823e-01,6.932862e-01,6.7
163 6.510355e-01,6.379834e-01,6.254100e-01,6.1
164 5.903633e-01,5.795116e-01,5.690472e-01,5.5
165 5.398295e-01,5.307700e-01,5.220296e-01,5.1
166 4.976128e-01,4.900419e-01,4.827400e-01,4.7
167
168 static const G4double SL4[nL]={
169 1.531367e-04,6.750684e-04,2.023434e-03,4.8
170 1.816857e-02,3.094217e-02,4.965477e-02,7.6
171 1.614108e-01,2.270208e-01,3.153403e-01,4.3
172 8.886525e-01,1.345605e+00,2.121366e+00,3.2
173 5.172459e+00,5.243522e+00,5.175754e+00,5.1
174 5.151144e+00,5.108382e+00,5.025027e+00,4.9
175 4.621981e+00,4.464473e+00,4.303590e+00,4.1
176 3.822656e+00,3.667551e+00,3.516631e+00,3.3
177 3.094556e+00,2.965180e+00,2.841688e+00,2.7
178 2.506035e+00,2.405305e+00,2.309830e+00,2.2
179 2.052608e+00,1.975802e+00,1.903066e+00,1.8
180 1.707024e+00,1.648361e+00,1.592714e+00,1.5
181 1.442123e+00,1.396846e+00,1.353788e+00,1.3
182 1.236668e+00,1.201266e+00,1.167510e+00,1.1
183 1.075223e+00,1.047183e+00,1.020382e+00,9.9
184 9.467696e-01,9.243013e-01,9.027797e-01,8.8
185 8.434307e-01,8.252457e-01,8.077958e-01,7.9
186 7.595173e-01,7.446784e-01,7.304196e-01,7.1
187 6.908733e-01,6.786917e-01,6.669756e-01,6.5
188 6.344356e-01,6.244022e-01,6.147497e-01,6.0
189 5.879438e-01,5.796850e-01,5.717461e-01,5.6
190
191 static const G4double SL5[nL]={
192 1.905569e-04,7.771730e-04,2.250919e-03,5.2
193 1.969996e-02,3.365091e-02,5.440813e-02,8.4
194 1.866020e-01,2.707115e-01,3.912405e-01,5.7
195 1.317340e+00,2.143911e+00,3.657987e+00,6.3
196 1.571664e+01,1.840405e+01,1.776700e+01,1.5
197 1.138076e+01,9.874227e+00,8.700723e+00,7.7
198 6.458855e+00,5.969695e+00,5.556515e+00,5.2
199 4.609287e+00,4.358030e+00,4.129172e+00,3.9
200 3.545861e+00,3.378977e+00,3.223486e+00,3.0
201 2.815610e+00,2.696573e+00,2.584914e+00,2.4
202 2.288930e+00,2.201736e+00,2.119606e+00,2.0
203 1.900204e+00,1.835064e+00,1.773480e+00,1.7
204 1.607794e+00,1.558252e+00,1.511256e+00,1.4
205 1.384013e+00,1.345726e+00,1.309299e+00,1.2
206 1.210133e+00,1.180134e+00,1.151523e+00,1.1
207 1.073277e+00,1.049504e+00,1.026785e+00,1.0
208 9.644275e-01,9.454166e-01,9.272213e-01,9.0
209 8.771474e-01,8.618447e-01,8.471837e-01,8.3
210 8.067717e-01,7.944065e-01,7.825554e-01,7.7
211 7.498817e-01,7.398883e-01,7.303156e-01,7.2
212 7.039741e-01,6.959417e-01,6.882635e-01,6.8
213
214 static const G4double SL6[nL]={
215 2.222448e-04,8.620556e-04,2.444896e-03,5.7
216 2.205349e-02,3.918281e-02,6.696997e-02,1.1
217 2.959155e-01,4.753435e-01,7.596938e-01,1.2
218 3.111254e+00,5.043478e+00,8.232698e+00,1.3
219 2.638709e+01,2.859878e+01,2.728600e+01,2.4
220 1.931943e+01,1.710886e+01,1.522482e+01,1.3
221 1.113950e+01,1.015995e+01,9.314220e+00,8.5
222 7.356396e+00,6.844921e+00,6.386042e+00,5.9
223 5.256885e+00,4.946267e+00,4.662228e+00,4.4
224 3.942438e+00,3.739345e+00,3.551711e+00,3.3
225 3.067782e+00,2.928947e+00,2.799729e+00,2.6
226 2.461965e+00,2.363783e+00,2.271852e+00,2.1
227 2.028885e+00,1.957488e+00,1.890304e+00,1.8
228 1.711104e+00,1.657972e+00,1.607769e+00,1.5
229 1.472853e+00,1.432558e+00,1.394357e+00,1.3
230 1.291070e+00,1.260042e+00,1.230549e+00,1.2
231 1.150425e+00,1.126246e+00,1.103218e+00,1.0
232 1.040438e+00,1.021433e+00,1.003308e+00,9.8
233 9.537847e-01,9.387662e-01,9.244342e-01,9.1
234 8.852523e-01,8.733710e-01,8.620378e-01,8.5
235 8.311102e-01,8.217603e-01,8.128613e-01,8.0
236
237 static const G4double SL7[nL]={
238 2.400132e-04,9.082999e-04,2.545511e-03,5.9
239 2.297237e-02,4.117454e-02,7.124517e-02,1.2
240 3.306145e-01,5.401144e-01,8.769596e-01,1.4
241 3.727247e+00,6.087430e+00,9.967584e+00,1.6
242 2.968929e+01,3.091057e+01,2.878613e+01,2.5
243 1.981778e+01,1.747914e+01,1.551277e+01,1.3
244 1.131384e+01,1.031654e+01,9.457362e+00,8.7
245 7.472543e+00,6.954497e+00,6.489768e+00,6.0
246 5.346149e+00,5.031499e+00,4.743750e+00,4.4
247 4.014464e+00,3.808676e+00,3.618553e+00,3.4
248 3.128239e+00,2.987593e+00,2.856703e+00,2.7
249 2.514665e+00,2.415276e+00,2.322236e+00,2.2
250 2.076456e+00,2.004275e+00,1.936374e+00,1.8
251 1.755391e+00,1.701773e+00,1.651134e+00,1.6
252 1.515169e+00,1.474604e+00,1.436168e+00,1.3
253 1.332368e+00,1.301226e+00,1.271645e+00,1.2
254 1.191399e+00,1.167222e+00,1.144214e+00,1.1
255 1.081599e+00,1.062680e+00,1.044655e+00,1.0
256 9.955079e-01,9.806377e-01,9.664643e-01,9.5
257 9.278146e-01,9.161277e-01,9.049967e-01,8.9
258 8.747192e-01,8.655994e-01,8.569368e-01,8.4
259
260 static const G4double SL8[nL]={
261 2.590923e-04,9.573672e-04,2.651275e-03,6.1
262 2.396311e-02,4.335926e-02,7.599430e-02,1.3
263 3.685911e-01,6.084138e-01,9.922345e-01,1.5
264 4.001018e+00,6.212589e+00,9.507804e+00,1.4
265 2.598385e+01,2.841920e+01,2.739643e+01,2.4
266 1.941582e+01,1.718791e+01,1.530626e+01,1.3
267 1.126539e+01,1.029864e+01,9.462554e+00,8.7
268 7.516956e+00,7.005488e+00,6.545306e+00,6.1
269 4.805315e+00,4.541037e+00,4.297946e+00,4.0
270 3.334802e+00,3.182394e+00,3.040697e+00,2.9
271 2.463791e+00,2.369986e+00,2.282085e+00,2.1
272 1.916681e+00,1.856017e+00,1.798827e+00,1.7
273 1.557271e+00,1.516539e+00,1.477966e+00,1.4
274 1.313142e+00,1.285040e+00,1.258344e+00,1.2
275 1.143384e+00,1.123646e+00,1.104863e+00,1.0
276 1.023662e+00,1.009685e+00,9.963805e-01,9.8
277 9.389503e-01,9.291078e-01,9.197629e-01,9.1
278
279 static const G4double SL9[nL]={
280 3.243985e-04,1.122034e-03,3.000932e-03,6.8
281 9.887958e-02,1.966468e-01,4.282973e-01,1.0
282 1.162472e+01,1.124054e+01,1.202416e+01,1.4
283 2.059083e+01,1.993672e+01,1.873926e+01,1.7
284 1.230272e+01,1.135820e+01,1.052046e+01,9.7
285 7.451996e+00,6.995366e+00,6.576679e+00,6.1
286 4.931401e+00,4.673850e+00,4.435420e+00,4.2
287 3.478323e+00,3.325356e+00,3.182848e+00,3.0
288 2.600432e+00,2.505518e+00,2.416558e+00,2.3
289 2.046943e+00,1.985682e+00,1.927976e+00,1.8
290 1.685032e+00,1.644242e+00,1.605677e+00,1.5
291 1.441781e+00,1.414020e+00,1.387711e+00,1.3
292 1.275266e+00,1.256130e+00,1.237976e+00,1.2
293 1.160265e+00,1.147042e+00,1.134507e+00,1.1
294 1.081106e+00,1.072098e+00,1.063597e+00,1.0
295
296 static const G4double SL10[nL]={
297 4.311217e-04,1.384716e-03,3.549518e-03,7.9
298 1.266167e-01,2.409191e-01,4.501490e-01,8.2
299 7.790746e+00,1.287033e+01,1.909053e+01,2.3
300 2.536786e+01,2.330746e+01,2.113183e+01,1.9
301 1.304262e+01,1.200459e+01,1.109667e+01,1.0
302 7.837784e+00,7.358628e+00,6.920153e+00,6.5
303 5.203837e+00,4.936070e+00,4.688412e+00,4.4
304 3.696964e+00,3.539013e+00,3.392026e+00,3.2
305 2.793474e+00,2.696364e+00,2.605493e+00,2.5
306 2.229875e+00,2.168007e+00,2.109859e+00,2.0
307 1.866773e+00,1.826299e+00,1.788145e+00,1.7
308 1.627479e+00,1.600555e+00,1.575133e+00,1.5
309 1.467726e+00,1.449689e+00,1.432656e+00,1.4
310 1.360776e+00,1.348747e+00,1.337409e+00,1.3
311 1.289997e+00,1.282179e+00,1.274863e+00,1.2
312
313 static const G4double SL11[nL]={
314 4.614524e-04,1.458509e-03,3.702639e-03,8.3
315 1.367187e-01,2.630019e-01,4.950763e-01,9.0
316 7.867272e+00,1.250247e+01,1.878669e+01,2.5
317 2.689359e+01,2.438858e+01,2.190927e+01,1.9
318 1.327514e+01,1.220222e+01,1.126790e+01,1.0
319 7.940129e+00,7.453592e+00,7.008848e+00,6.6
320 5.271842e+00,5.001344e+00,4.751274e+00,4.5
321 3.751336e+00,3.592227e+00,3.444221e+00,3.3
322 2.842215e+00,2.744679e+00,2.653450e+00,2.5
323 2.276915e+00,2.215004e+00,2.156852e+00,2.1
324 1.914217e+00,1.873910e+00,1.835942e+00,1.8
325 1.676451e+00,1.649799e+00,1.624659e+00,1.6
326 1.518758e+00,1.501034e+00,1.484317e+00,1.4
327 1.414022e+00,1.402307e+00,1.391282e+00,1.3
328 1.345386e+00,1.337862e+00,1.330834e+00,1.3
329
330 static const G4double SL12[nL]={
331 5.615148e-04,1.700309e-03,4.203181e-03,9.3
332 1.726852e-01,3.430025e-01,6.623201e-01,1.2
333 1.070034e+01,1.656745e+01,2.430795e+01,3.3
334 3.426056e+01,2.971854e+01,2.578645e+01,2.2
335 1.434242e+01,1.307807e+01,1.199938e+01,1.1
336 8.310522e+00,7.788614e+00,7.314751e+00,6.8
337 5.489325e+00,5.208212e+00,4.949006e+00,4.7
338 3.918886e+00,3.755935e+00,3.604610e+00,3.4
339 2.991993e+00,2.893250e+00,2.801049e+00,2.7
340 2.422448e+00,2.360566e+00,2.302559e+00,2.2
341 2.062033e+00,2.022365e+00,1.985093e+00,1.9
342 1.829709e+00,1.803970e+00,1.779763e+00,1.7
343 1.678712e+00,1.661974e+00,1.646241e+00,1.6
344 1.580780e+00,1.570002e+00,1.559901e+00,1.5
345 1.518397e+00,1.511701e+00,1.505485e+00,1.4
346
347 static const G4double SL13[nL]={
348 5.979521e-04,1.787895e-03,4.384312e-03,9.7
349 1.867226e-01,3.746609e-01,7.299098e-01,1.3
350 1.267423e+01,2.032656e+01,3.102537e+01,4.2
351 3.635559e+01,3.114336e+01,2.680306e+01,2.3
352 1.463175e+01,1.331842e+01,1.220266e+01,1.1
353 8.419927e+00,7.888645e+00,7.407021e+00,6.9
354 5.557733e+00,5.273725e+00,5.012027e+00,4.7
355 3.973720e+00,3.809747e+00,3.657550e+00,3.5
356 3.042294e+00,2.943289e+00,2.850895e+00,2.7
357 2.472142e+00,2.410358e+00,2.352483e+00,2.2
358 2.113006e+00,2.073607e+00,2.036619e+00,2.0
359 1.882814e+00,1.857411e+00,1.833544e+00,1.8
360 1.734206e+00,1.717807e+00,1.702410e+00,1.6
361 1.638561e+00,1.628089e+00,1.618286e+00,1.6
362 1.578167e+00,1.571726e+00,1.565757e+00,1.5
363
364 static const G4double SL14[nL]={
365 7.595609e-04,2.174487e-03,5.184472e-03,1.1
366 2.534357e-01,5.242273e-01,1.043266e+00,1.9
367 1.696937e+01,2.517085e+01,3.466133e+01,4.3
368 3.547368e+01,3.049881e+01,2.637019e+01,2.3
369 1.475365e+01,1.348118e+01,1.239409e+01,1.1
370 8.657732e+00,8.126345e+00,7.643355e+00,7.2
371 5.779927e+00,5.492917e+00,5.228379e+00,4.9
372 4.179304e+00,4.013924e+00,3.860561e+00,3.7
373 3.242733e+00,3.143782e+00,3.051603e+00,2.9
374 2.675929e+00,2.615087e+00,2.558239e+00,2.5
375 2.324889e+00,2.286859e+00,2.251272e+00,2.2
376 2.104744e+00,2.080815e+00,2.058418e+00,2.0
377 1.966234e+00,1.951205e+00,1.937150e+00,1.9
378 1.879545e+00,1.870216e+00,1.861520e+00,1.8
379 1.826362e+00,1.820798e+00,1.815670e+00,1.8
380
381 static const G4double SL15[nL]={
382 8.500963e-04,2.390172e-03,5.632030e-03,1.2
383 2.941224e-01,6.161245e-01,1.237476e+00,2.3
384 2.037233e+01,3.017649e+01,4.195065e+01,5.4
385 4.707792e+01,3.899572e+01,3.258608e+01,2.7
386 1.631941e+01,1.471906e+01,1.338580e+01,1.2
387 9.054659e+00,8.469380e+00,7.943327e+00,7.4
388 5.960212e+00,5.660481e+00,5.385398e+00,5.1
389 4.304884e+00,4.135994e+00,3.979725e+00,3.8
390 3.353748e+00,3.254064e+00,3.161362e+00,3.0
391 2.785372e+00,2.724800e+00,2.668301e+00,2.6
392 2.437529e+00,2.400129e+00,2.365194e+00,2.3
393 2.222116e+00,2.198888e+00,2.177188e+00,2.1
394 2.088348e+00,2.073945e+00,2.060501e+00,2.0
395 2.005648e+00,1.996806e+00,1.988574e+00,1.9
396 1.955397e+00,1.950164e+00,1.945345e+00,1.9
397
398 static const G4double SL16[nL]={
399 1.161977e-03,3.130797e-03,7.178175e-03,1.5
400 4.493705e-01,9.672255e-01,1.976461e+00,3.8
401 2.929890e+01,3.906811e+01,4.690664e+01,5.2
402 4.880825e+01,4.096475e+01,3.429691e+01,2.9
403 1.707654e+01,1.539807e+01,1.400363e+01,1.2
404 9.501359e+00,8.896578e+00,8.353847e+00,7.8
405 6.316575e+00,6.010193e+00,5.729483e+00,5.4
406 4.632774e+00,4.462505e+00,4.305337e+00,4.1
407 3.680544e+00,3.581977e+00,3.490615e+00,3.4
408 3.123820e+00,3.065443e+00,3.011214e+00,2.9
409 2.792457e+00,2.757506e+00,2.725010e+00,2.6
410 2.593694e+00,2.572682e+00,2.553138e+00,2.5
411 2.474040e+00,2.461354e+00,2.449544e+00,2.4
412 2.401579e+00,2.393854e+00,2.386654e+00,2.3
413 2.357400e+00,2.352718e+00,2.348380e+00,2.3
414
415 static const G4double SL17[nL]={
416 2.137065e-03,5.442007e-03,1.210645e-02,2.7
417 1.057471e+00,2.354951e+00,4.918482e+00,9.6
418 6.860085e+01,9.419933e+01,1.267075e+02,1.5
419 7.155618e+01,5.654375e+01,4.555533e+01,3.7
420 2.039992e+01,1.816677e+01,1.635312e+01,1.4
421 1.081664e+01,1.011536e+01,9.495083e+00,8.9
422 7.243487e+00,6.915590e+00,6.617778e+00,6.3
423 5.480788e+00,5.308639e+00,5.150985e+00,5.0
424 4.538429e+00,4.444301e+00,4.357799e+00,4.2
425 4.019090e+00,3.966677e+00,3.918419e+00,3.8
426 3.728406e+00,3.698788e+00,3.671436e+00,3.6
427 3.562470e+00,3.545177e+00,3.529084e+00,3.5
428 3.463103e+00,3.452209e+00,3.441916e+00,3.4
429 3.397582e+00,3.389836e+00,3.382382e+00,3.3
430 3.348730e+00,3.342620e+00,3.336699e+00,3.3
431
432 static const G4double SL18[nL]={
433 2.220534e-03,5.640053e-03,1.253572e-02,2.8
434 1.115760e+00,2.485562e+00,5.183559e+00,1.0
435 6.345164e+01,7.777111e+01,8.950804e+01,9.3
436 4.641759e+01,3.856198e+01,3.257293e+01,2.7
437 1.742802e+01,1.589540e+01,1.460538e+01,1.3
438 1.031123e+01,9.718498e+00,9.183826e+00,8.6
439 7.161876e+00,6.857372e+00,6.578785e+00,6.3
440 5.499048e+00,5.333567e+00,5.181614e+00,5.0
441 4.587960e+00,4.496341e+00,4.412068e+00,4.3
442 4.081496e+00,4.030275e+00,3.983101e+00,3.9
443 3.797233e+00,3.768237e+00,3.741451e+00,3.7
444 3.634571e+00,3.617565e+00,3.601721e+00,3.5
445 3.536470e+00,3.525629e+00,3.515361e+00,3.5
446 3.470787e+00,3.462928e+00,3.455342e+00,3.4
447 3.420792e+00,3.414463e+00,3.408314e+00,3.4
448
449 static const G4double SL19[nL]={
450 2.305897e-03,5.842654e-03,1.297593e-02,2.9
451 1.179632e+00,2.638978e+00,5.539887e+00,1.0
452 8.715375e+01,1.188643e+02,1.303680e+02,1.1
453 4.768315e+01,3.944920e+01,3.322398e+01,2.8
454 1.762968e+01,1.606532e+01,1.475137e+01,1.3
455 1.039872e+01,9.800670e+00,9.261781e+00,8.7
456 7.228815e+00,6.923346e+00,6.644047e+00,6.3
457 5.563316e+00,5.397976e+00,5.246237e+00,5.1
458 4.654358e+00,4.563177e+00,4.479354e+00,4.4
459 4.151059e+00,4.100274e+00,4.053523e+00,4.0
460 3.869510e+00,3.840823e+00,3.814323e+00,3.7
461 3.708509e+00,3.691642e+00,3.675912e+00,3.6
462 3.610872e+00,3.600004e+00,3.589688e+00,3.5
463 3.544574e+00,3.536552e+00,3.528789e+00,3.5
464 3.493151e+00,3.486572e+00,3.480165e+00,3.4
465
466 static const G4double SL20[nL]={
467 2.545914e-03,6.412659e-03,1.422001e-02,3.3
468 1.354715e+00,3.029540e+00,6.323258e+00,1.2
469 6.796031e+01,7.669870e+01,8.176394e+01,8.7
470 5.574957e+01,4.551980e+01,3.772808e+01,3.1
471 1.882555e+01,1.702523e+01,1.553424e+01,1.4
472 1.075719e+01,1.012157e+01,9.553019e+00,9.0
473 7.441691e+00,7.128756e+00,6.843578e+00,6.5
474 5.748752e+00,5.582537e+00,5.430325e+00,5.2
475 4.839973e+00,4.749573e+00,4.666616e+00,4.5
476 4.343241e+00,4.293452e+00,4.247676e+00,4.2
477 4.067961e+00,4.039982e+00,4.014130e+00,3.9
478 3.910607e+00,3.894002e+00,3.878468e+00,3.8
479 3.813448e+00,3.802402e+00,3.791850e+00,3.7
480 3.744778e+00,3.736225e+00,3.727891e+00,3.7
481 3.688887e+00,3.681555e+00,3.674374e+00,3.6
482
483 static const G4double SL21[nL]={
484 2.564250e-03,6.456227e-03,1.431544e-02,3.3
485 1.368585e+00,3.061075e+00,6.389076e+00,1.2
486 6.804530e+01,7.622389e+01,8.041323e+01,8.5
487 6.319441e+01,5.115788e+01,4.188960e+01,3.4
488 1.981637e+01,1.779045e+01,1.613068e+01,1.4
489 1.095233e+01,1.028156e+01,9.685270e+00,9.1
490 7.499706e+00,7.179292e+00,6.888008e+00,6.6
491 5.775573e+00,5.607438e+00,5.453629e+00,5.3
492 4.858517e+00,4.767585e+00,4.684186e+00,4.6
493 4.359506e+00,4.309575e+00,4.263683e+00,4.2
494 4.083632e+00,4.055616e+00,4.029732e+00,4.0
495 3.926087e+00,3.909459e+00,3.893901e+00,3.8
496 3.828735e+00,3.817652e+00,3.807062e+00,3.7
497 3.759759e+00,3.751154e+00,3.742764e+00,3.7
498 3.703458e+00,3.696062e+00,3.688816e+00,3.6
499
500 static const G4double SL22[nL]={
501 3.007427e-03,7.510236e-03,1.663782e-02,3.9
502 1.719543e+00,3.868385e+00,8.112125e+00,1.5
503 9.555037e+01,1.153230e+02,1.159518e+02,9.9
504 4.371769e+01,3.676918e+01,3.143378e+01,2.7
505 1.759555e+01,1.614984e+01,1.492332e+01,1.3
506 1.076997e+01,1.018862e+01,9.662970e+00,9.1
507 7.669857e+00,7.369937e+00,7.095886e+00,6.8
508 6.039585e+00,5.879011e+00,5.732019e+00,5.5
509 5.163281e+00,5.076520e+00,4.997014e+00,4.9
510 4.688303e+00,4.640952e+00,4.597451e+00,4.5
511 4.426712e+00,4.400054e+00,4.375369e+00,4.3
512 4.275425e+00,4.259104e+00,4.243719e+00,4.2
513 4.177536e+00,4.165907e+00,4.154667e+00,4.1
514 4.102746e+00,4.092980e+00,4.083363e+00,4.0
515 4.037123e+00,4.028220e+00,4.019441e+00,4.0
516
517 static const G4double SL23[nL]={
518 3.202591e-03,7.975022e-03,1.767003e-02,4.1
519 1.881188e+00,4.239060e+00,8.896849e+00,1.7
520 1.037149e+02,1.232794e+02,1.201643e+02,1.0
521 4.398670e+01,3.699713e+01,3.163451e+01,2.7
522 1.773851e+01,1.628791e+01,1.505750e+01,1.3
523 1.089372e+01,1.031145e+01,9.785155e+00,9.3
524 7.791975e+00,7.492525e+00,7.219066e+00,6.9
525 6.167007e+00,6.007447e+00,5.861497e+00,5.7
526 5.298080e+00,5.212351e+00,5.133851e+00,5.0
527 4.829614e+00,4.783020e+00,4.740226e+00,4.7
528 4.572181e+00,4.545889e+00,4.521511e+00,4.4
529 4.422252e+00,4.405903e+00,4.390437e+00,4.3
530 4.323132e+00,4.311147e+00,4.299511e+00,4.2
531 4.245096e+00,4.234743e+00,4.224516e+00,4.2
532 4.174947e+00,4.165339e+00,4.155849e+00,4.1
533
534 static const G4double SL24[nL]={
535 4.424391e-03,1.089365e-02,2.425454e-02,5.9
536 3.000059e+00,6.849073e+00,1.459387e+01,2.9
537 1.455690e+02,1.285470e+02,1.059779e+02,8.6
538 4.051569e+01,3.473537e+01,3.020519e+01,2.6
539 1.794276e+01,1.660309e+01,1.545498e+01,1.4
540 1.148642e+01,1.092239e+01,1.041132e+01,9.9
541 8.471504e+00,8.180659e+00,7.915525e+00,7.6
542 6.902666e+00,6.750545e+00,6.611879e+00,6.4
543 6.081868e+00,6.002071e+00,5.929208e+00,5.8
544 5.648227e+00,5.605225e+00,5.565653e+00,5.5
545 5.408150e+00,5.382864e+00,5.359136e+00,5.3
546 5.258062e+00,5.240405e+00,5.223341e+00,5.2
547 5.144141e+00,5.129087e+00,5.114183e+00,5.0
548 5.040923e+00,5.026397e+00,5.011893e+00,4.9
549 4.939879e+00,4.925668e+00,4.911569e+00,4.8
550
551 static const G4double SL25[nL]={
552 5.218262e-03,1.279812e-02,2.863691e-02,7.1
553 3.785207e+00,8.620686e+00,1.819011e+01,3.5
554 1.795448e+02,1.653851e+02,1.332773e+02,1.0
555 4.486079e+01,3.794689e+01,3.263665e+01,2.8
556 1.881372e+01,1.736415e+01,1.613394e+01,1.5
557 1.197431e+01,1.139498e+01,1.087253e+01,1.0
558 8.911366e+00,8.620462e+00,8.356070e+00,8.1
559 7.353992e+00,7.204745e+00,7.069027e+00,6.9
560 6.553303e+00,6.476044e+00,6.405559e+00,6.3
561 6.133589e+00,6.091756e+00,6.053132e+00,6.0
562 5.897073e+00,5.871438e+00,5.847158e+00,5.8
563 5.740528e+00,5.721253e+00,5.702421e+00,5.6
564 5.612426e+00,5.594882e+00,5.577398e+00,5.5
565 5.490177e+00,5.472703e+00,5.455220e+00,5.4
566 5.368096e+00,5.350876e+00,5.333791e+00,5.3
567
568 static const G4double SL26[nL]={
569 9.533418e-03,2.324917e-02,5.364098e-02,1.4
570 8.890991e+00,2.014541e+01,4.138069e+01,7.5
571 2.305942e+02,2.457159e+02,2.095925e+02,1.6
572 5.903474e+01,4.855591e+01,4.078885e+01,3.4
573 2.206960e+01,2.028240e+01,1.880303e+01,1.7
574 1.409475e+01,1.347801e+01,1.293034e+01,1.2
575 1.094945e+01,1.066638e+01,1.041174e+01,1.0
576 9.470694e+00,9.333803e+00,9.209956e+00,9.0
577 8.741745e+00,8.671018e+00,8.605956e+00,8.5
578 8.344033e+00,8.300834e+00,8.259766e+00,8.2
579 8.076498e+00,8.042705e+00,8.009465e+00,7.9
580 7.847632e+00,7.815563e+00,7.783460e+00,7.7
581 7.621469e+00,7.588688e+00,7.555781e+00,7.5
582 7.389908e+00,7.356640e+00,7.323419e+00,7.2
583 7.159594e+00,7.127657e+00,7.096141e+00,7.0
584
585 static const G4double SL27[nL]={
586 1.043535e-02,2.545247e-02,5.908485e-02,1.6
587 1.012392e+01,2.297013e+01,4.720890e+01,8.6
588 2.968858e+02,2.729073e+02,2.110699e+02,1.5
589 5.807038e+01,4.795957e+01,4.044610e+01,3.4
590 2.223556e+01,2.048595e+01,1.903566e+01,1.7
591 1.440731e+01,1.379991e+01,1.326049e+01,1.2
592 1.131048e+01,1.103217e+01,1.078194e+01,1.0
593 9.858006e+00,9.723600e+00,9.601935e+00,9.4
594 9.140195e+00,9.069912e+00,9.005027e+00,8.9
595 8.740293e+00,8.695865e+00,8.653370e+00,8.6
596 8.460418e+00,8.424256e+00,8.388530e+00,8.3
597 8.212851e+00,8.177798e+00,8.142660e+00,8.1
598 7.965007e+00,7.929047e+00,7.892960e+00,7.8
599 7.711389e+00,7.675061e+00,7.638818e+00,7.6
600 7.460641e+00,7.426023e+00,7.391901e+00,7.3
601
602 static const G4double SL28[nL]={
603 1.177612e-02,2.873855e-02,6.729897e-02,1.8
604 1.203191e+01,2.730811e+01,5.598968e+01,1.0
605 2.688536e+02,2.132293e+02,1.608280e+02,1.2
606 5.063436e+01,4.279795e+01,3.684661e+01,3.2
607 2.174808e+01,2.021937e+01,1.893580e+01,1.7
608 1.471984e+01,1.415200e+01,1.364501e+01,1.3
609 1.179407e+01,1.152799e+01,1.128840e+01,1.1
610 1.040057e+01,1.027077e+01,1.015302e+01,1.0
611 9.701922e+00,9.632220e+00,9.567466e+00,9.5
612 9.297506e+00,9.251017e+00,9.206173e+00,9.1
613 8.997851e+00,8.958023e+00,8.918474e+00,8.8
614 8.721877e+00,8.682378e+00,8.642738e+00,8.6
615 8.442103e+00,8.401522e+00,8.360828e+00,8.3
616 8.156696e+00,8.116003e+00,8.075461e+00,8.0
617 7.877018e+00,7.838644e+00,7.800883e+00,7.7
618
619 static const G4double SL29[nL]={
620 1.365967e-02,3.337537e-02,7.906748e-02,2.2
621 1.480860e+01,3.341772e+01,6.741315e+01,1.1
622 3.008511e+02,2.541744e+02,1.921484e+02,1.4
623 5.584793e+01,4.667796e+01,3.980970e+01,3.4
624 2.290273e+01,2.125263e+01,1.988051e+01,1.8
625 1.547709e+01,1.489750e+01,1.438287e+01,1.3
626 1.252665e+01,1.226263e+01,1.202543e+01,1.1
627 1.114933e+01,1.102117e+01,1.090467e+01,1.0
628 1.045353e+01,1.038251e+01,1.031602e+01,1.0
629 1.003150e+01,9.981059e+00,9.931961e+00,9.8
630 9.698622e+00,9.653189e+00,9.607882e+00,9.5
631 9.380825e+00,9.335018e+00,9.289031e+00,9.2
632 9.056561e+00,9.009673e+00,8.962717e+00,8.9
633 8.728269e+00,8.681773e+00,8.635535e+00,8.5
634 8.410514e+00,8.367269e+00,8.324807e+00,8.2
635
636 static const G4double SL30[nL]={
637 2.103117e-02,5.172350e-02,1.273601e-01,3.8
638 2.709930e+01,6.005843e+01,1.151063e+02,1.8
639 2.282780e+02,2.025860e+02,1.607669e+02,1.2
640 5.328718e+01,4.543843e+01,3.947900e+01,3.4
641 2.445368e+01,2.295694e+01,2.170826e+01,2.0
642 1.768163e+01,1.715075e+01,1.667954e+01,1.6
643 1.498116e+01,1.473901e+01,1.452086e+01,1.4
644 1.370056e+01,1.357624e+01,1.346135e+01,1.3
645 1.298788e+01,1.290703e+01,1.282916e+01,1.2
646 1.246828e+01,1.239939e+01,1.233100e+01,1.2
647 1.199175e+01,1.192379e+01,1.185567e+01,1.1
648 1.151237e+01,1.144324e+01,1.137400e+01,1.1
649 1.102773e+01,1.095880e+01,1.089011e+01,1.0
650 1.055244e+01,1.048657e+01,1.042144e+01,1.0
651 1.011014e+01,1.005149e+01,9.994313e+00,9.9
652
653 static const G4double SL31[nL]={
654 2.164664e-02,5.326850e-02,1.315360e-01,3.9
655 2.822812e+01,6.251535e+01,1.195506e+02,1.8
656 2.417968e+02,1.964670e+02,1.512314e+02,1.1
657 5.170376e+01,4.435567e+01,3.874082e+01,3.4
658 2.440774e+01,2.296126e+01,2.175087e+01,2.0
659 1.782157e+01,1.730033e+01,1.683705e+01,1.6
660 1.516228e+01,1.492277e+01,1.470678e+01,1.4
661 1.389184e+01,1.376774e+01,1.365284e+01,1.3
662 1.317653e+01,1.309467e+01,1.301566e+01,1.2
663 1.264768e+01,1.257715e+01,1.250708e+01,1.2
664 1.215890e+01,1.208909e+01,1.201911e+01,1.1
665 1.166655e+01,1.159559e+01,1.152454e+01,1.1
666 1.116955e+01,1.109896e+01,1.102863e+01,1.0
667 1.068331e+01,1.061602e+01,1.054952e+01,1.0
668 1.023205e+01,1.017232e+01,1.011413e+01,1.0
669
670 static const G4double SL32[nL]={
671 2.258863e-02,5.563670e-02,1.379665e-01,4.1
672 2.998579e+01,6.634936e+01,1.265293e+02,1.9
673 2.792187e+02,2.183968e+02,1.646189e+02,1.2
674 5.375330e+01,4.588048e+01,3.990576e+01,3.5
675 2.486877e+01,2.337590e+01,2.213188e+01,2.1
676 1.813222e+01,1.760647e+01,1.714012e+01,1.6
677 1.546043e+01,1.522073e+01,1.500457e+01,1.4
678 1.418745e+01,1.406248e+01,1.394655e+01,1.3
679 1.346262e+01,1.337878e+01,1.329768e+01,1.3
680 1.291767e+01,1.284450e+01,1.277173e+01,1.2
681 1.240949e+01,1.233681e+01,1.226397e+01,1.2
682 1.189720e+01,1.182346e+01,1.174967e+01,1.1
683 1.138148e+01,1.130838e+01,1.123558e+01,1.1
684 1.087876e+01,1.080935e+01,1.074079e+01,1.0
685 1.041413e+01,1.035280e+01,1.029310e+01,1.0
686
687 static const G4double SL33[nL]={
688 2.454062e-02,6.055745e-02,1.514382e-01,4.6
689 3.362441e+01,7.394409e+01,1.387856e+02,2.0
690 2.124190e+02,1.687167e+02,1.317569e+02,1.0
691 4.945432e+01,4.295680e+01,3.792803e+01,3.3
692 2.479430e+01,2.344031e+01,2.230232e+01,2.1
693 1.857355e+01,1.807474e+01,1.763054e+01,1.7
694 1.601607e+01,1.578341e+01,1.557289e+01,1.5
695 1.476805e+01,1.464303e+01,1.452637e+01,1.4
696 1.403079e+01,1.394334e+01,1.385828e+01,1.3
697 1.345466e+01,1.337620e+01,1.329802e+01,1.3
698 1.290750e+01,1.282905e+01,1.275042e+01,1.2
699 1.235509e+01,1.227578e+01,1.219646e+01,1.2
700 1.180186e+01,1.172374e+01,1.164603e+01,1.1
701 1.126633e+01,1.119272e+01,1.112009e+01,1.1
702 1.077532e+01,1.071087e+01,1.064824e+01,1.0
703
704 static const G4double SL34[nL]={
705 2.555084e-02,6.311099e-02,1.584856e-01,4.8
706 3.561279e+01,7.827264e+01,1.466546e+02,2.2
707 2.658424e+02,2.035783e+02,1.541383e+02,1.1
708 5.281838e+01,4.541822e+01,3.977013e+01,3.5
709 2.541809e+01,2.398190e+01,2.278352e+01,2.1
710 1.892070e+01,1.841172e+01,1.795992e+01,1.7
711 1.632757e+01,1.609322e+01,1.588122e+01,1.5
712 1.506950e+01,1.494290e+01,1.482456e+01,1.4
713 1.431870e+01,1.422884e+01,1.414128e+01,1.4
714 1.372397e+01,1.364260e+01,1.356147e+01,1.3
715 1.315592e+01,1.307446e+01,1.299284e+01,1.2
716 1.258282e+01,1.250066e+01,1.241854e+01,1.2
717 1.201058e+01,1.192993e+01,1.184976e+01,1.1
718 1.145863e+01,1.138292e+01,1.130827e+01,1.1
719 1.095457e+01,1.088860e+01,1.082454e+01,1.0
720
721 static const G4double SL35[nL]={
722 2.764032e-02,6.840683e-02,1.732180e-01,5.4
723 3.965730e+01,8.650314e+01,1.589997e+02,2.3
724 1.892896e+02,1.529510e+02,1.218148e+02,9.8
725 4.866270e+01,4.259935e+01,3.787106e+01,3.4
726 2.537244e+01,2.407192e+01,2.297720e+01,2.2
727 1.937927e+01,1.889649e+01,1.846614e+01,1.8
728 1.689546e+01,1.666738e+01,1.646026e+01,1.6
729 1.565683e+01,1.552939e+01,1.540955e+01,1.5
730 1.488843e+01,1.479431e+01,1.470218e+01,1.4
731 1.425863e+01,1.417154e+01,1.408459e+01,1.3
732 1.364917e+01,1.356169e+01,1.347407e+01,1.3
733 1.303472e+01,1.294688e+01,1.285916e+01,1.2
734 1.242457e+01,1.233891e+01,1.225383e+01,1.2
735 1.184003e+01,1.176019e+01,1.168155e+01,1.1
736 1.131034e+01,1.124140e+01,1.117458e+01,1.1
737
738 static const G4double SL36[nL]={
739 2.908478e-02,7.207879e-02,1.835222e-01,5.7
740 4.268370e+01,9.316524e+01,1.715413e+02,2.5
741 3.599490e+02,2.689394e+02,1.969985e+02,1.4
742 5.959877e+01,5.047252e+01,4.364148e+01,3.8
743 2.698235e+01,2.539504e+01,2.408733e+01,2.2
744 1.999559e+01,1.947149e+01,1.900901e+01,1.8
745 1.735437e+01,1.711766e+01,1.690324e+01,1.6
746 1.607388e+01,1.594222e+01,1.581828e+01,1.5
747 1.527688e+01,1.517863e+01,1.508233e+01,1.4
748 1.461738e+01,1.452595e+01,1.443465e+01,1.4
749 1.397765e+01,1.388593e+01,1.379410e+01,1.3
750 1.333447e+01,1.324275e+01,1.315122e+01,1.3
751 1.269862e+01,1.260959e+01,1.252123e+01,1.2
752 1.209230e+01,1.200973e+01,1.192845e+01,1.1
753 1.154574e+01,1.147488e+01,1.140627e+01,1.1
754
755 static const G4double SL37[nL]={
756 2.982256e-02,7.395762e-02,1.888215e-01,5.9
757 4.418698e+01,9.625560e+01,1.763759e+02,2.5
758 3.605764e+02,2.693717e+02,1.973299e+02,1.4
759 5.978476e+01,5.065157e+01,4.381637e+01,3.8
760 2.715812e+01,2.557296e+01,2.426760e+01,2.3
761 2.018745e+01,1.966529e+01,1.920456e+01,1.8
762 1.755549e+01,1.731924e+01,1.710506e+01,1.6
763 1.627392e+01,1.614137e+01,1.601639e+01,1.5
764 1.546791e+01,1.536794e+01,1.526984e+01,1.5
765 1.479512e+01,1.470162e+01,1.460824e+01,1.4
766 1.414075e+01,1.404695e+01,1.395306e+01,1.3
767 1.348342e+01,1.338977e+01,1.329635e+01,1.3
768 1.283483e+01,1.274413e+01,1.265413e+01,1.2
769 1.221773e+01,1.213380e+01,1.205122e+01,1.1
770 1.166286e+01,1.159105e+01,1.152158e+01,1.1
771
772 static const G4double SL38[nL]={
773 3.019534e-02,7.490782e-02,1.915089e-01,6.0
774 4.516020e+01,9.894848e+01,1.839628e+02,2.8
775 3.029151e+02,2.212396e+02,1.647273e+02,1.2
776 5.515881e+01,4.738384e+01,4.148198e+01,3.6
777 2.668305e+01,2.522895e+01,2.402196e+01,2.3
778 2.017852e+01,1.967755e+01,1.923363e+01,1.8
779 1.763002e+01,1.739830e+01,1.718776e+01,1.6
780 1.636599e+01,1.623415e+01,1.610960e+01,1.5
781 1.556054e+01,1.546006e+01,1.536137e+01,1.5
782 1.488273e+01,1.478833e+01,1.469401e+01,1.4
783 1.422164e+01,1.412685e+01,1.403197e+01,1.3
784 1.355748e+01,1.346290e+01,1.336855e+01,1.3
785 1.290264e+01,1.281111e+01,1.272032e+01,1.2
786 1.228022e+01,1.219562e+01,1.211239e+01,1.2
787 1.172125e+01,1.164898e+01,1.157908e+01,1.1
788
789 static const G4double SL39[nL]={
790 3.132934e-02,7.780169e-02,1.997199e-01,6.3
791 4.750616e+01,1.037394e+02,1.913355e+02,2.8
792 2.761193e+02,2.038392e+02,1.536530e+02,1.1
793 5.380369e+01,4.649042e+01,4.090667e+01,3.6
794 2.675392e+01,2.534816e+01,2.417843e+01,2.3
795 2.043239e+01,1.994127e+01,1.950544e+01,1.9
796 1.792446e+01,1.769479e+01,1.748570e+01,1.7
797 1.666414e+01,1.653123e+01,1.640532e+01,1.6
798 1.584607e+01,1.574304e+01,1.564166e+01,1.5
799 1.514829e+01,1.505077e+01,1.495332e+01,1.4
800 1.446512e+01,1.436719e+01,1.426920e+01,1.4
801 1.377966e+01,1.368219e+01,1.358501e+01,1.3
802 1.310574e+01,1.301172e+01,1.291849e+01,1.2
803 1.246724e+01,1.238063e+01,1.229548e+01,1.2
804 1.189599e+01,1.182235e+01,1.175118e+01,1.1
805
806 static const G4double SL40[nL]={
807 3.209841e-02,7.976722e-02,2.053207e-01,6.5
808 4.932232e+01,1.081304e+02,2.014268e+02,3.1
809 2.780145e+02,2.044508e+02,1.540079e+02,1.1
810 5.399013e+01,4.667044e+01,4.108296e+01,3.6
811 2.693231e+01,2.552882e+01,2.436152e+01,2.3
812 2.062710e+01,2.013787e+01,1.970371e+01,1.9
813 1.812768e+01,1.789832e+01,1.768932e+01,1.7
814 1.686515e+01,1.673120e+01,1.660410e+01,1.6
815 1.603707e+01,1.593219e+01,1.582891e+01,1.5
816 1.532528e+01,1.522563e+01,1.512602e+01,1.5
817 1.462708e+01,1.452704e+01,1.442695e+01,1.4
818 1.392729e+01,1.382789e+01,1.372881e+01,1.3
819 1.324062e+01,1.314493e+01,1.305008e+01,1.2
820 1.259142e+01,1.250348e+01,1.241705e+01,1.2
821 1.201207e+01,1.193753e+01,1.186553e+01,1.1
822
823 static const G4double SL41[nL]={
824 3.651675e-02,9.110302e-02,2.379739e-01,7.7
825 5.938433e+01,1.302605e+02,2.436285e+02,3.9
826 3.113977e+02,2.250866e+02,1.673725e+02,1.2
827 5.665354e+01,4.882987e+01,4.290028e+01,3.8
828 2.812414e+01,2.668871e+01,2.550153e+01,2.4
829 2.175080e+01,2.126446e+01,2.083343e+01,2.0
830 1.926597e+01,1.903585e+01,1.882512e+01,1.8
831 1.797832e+01,1.783729e+01,1.770239e+01,1.7
832 1.708753e+01,1.697174e+01,1.685722e+01,1.6
833 1.629429e+01,1.618244e+01,1.607062e+01,1.5
834 1.551112e+01,1.539924e+01,1.528743e+01,1.5
835 1.473164e+01,1.462157e+01,1.451203e+01,1.4
836 1.397483e+01,1.387005e+01,1.376636e+01,1.3
837 1.326746e+01,1.317233e+01,1.307902e+01,1.2
838 1.264478e+01,1.256554e+01,1.248927e+01,1.2
839
840 static const G4double SL42[nL]={
841 3.967023e-02,9.923729e-02,2.617528e-01,8.5
842 6.608688e+01,1.423929e+02,2.556972e+02,3.8
843 4.319374e+02,3.026257e+02,2.171501e+02,1.6
844 6.421712e+01,5.439270e+01,4.709187e+01,4.1
845 2.962908e+01,2.801379e+01,2.669475e+01,2.5
846 2.265018e+01,2.214060e+01,2.169159e+01,2.1
847 2.007380e+01,1.983699e+01,1.961986e+01,1.9
848 1.874067e+01,1.859262e+01,1.845046e+01,1.8
849 1.779617e+01,1.767200e+01,1.754898e+01,1.7
850 1.694269e+01,1.682214e+01,1.670165e+01,1.6
851 1.609985e+01,1.597979e+01,1.585993e+01,1.5
852 1.526585e+01,1.514856e+01,1.503195e+01,1.4
853 1.446186e+01,1.435101e+01,1.424143e+01,1.4
854 1.371597e+01,1.361614e+01,1.351835e+01,1.3
855 1.306540e+01,1.298324e+01,1.290436e+01,1.2
856
857 static const G4double SL43[nL]={
858 4.090181e-02,1.024236e-01,2.711450e-01,8.9
859 7.002434e+01,1.540039e+02,2.909962e+02,4.8
860 3.262939e+02,2.342891e+02,1.734786e+02,1.3
861 5.825243e+01,5.021882e+01,4.414930e+01,3.9
862 2.914791e+01,2.770777e+01,2.652069e+01,2.5
863 2.279664e+01,2.231596e+01,2.188986e+01,2.1
864 2.033214e+01,2.010070e+01,1.988760e+01,1.9
865 1.901461e+01,1.886590e+01,1.872267e+01,1.8
866 1.805846e+01,1.793169e+01,1.780594e+01,1.7
867 1.718477e+01,1.706112e+01,1.693751e+01,1.6
868 1.632027e+01,1.619719e+01,1.607435e+01,1.5
869 1.546605e+01,1.534607e+01,1.522683e+01,1.5
870 1.464448e+01,1.453138e+01,1.441961e+01,1.4
871 1.388429e+01,1.378272e+01,1.368328e+01,1.3
872 1.322351e+01,1.314030e+01,1.306050e+01,1.2
873
874 static const G4double SL44[nL]={
875 4.170472e-02,1.045035e-01,2.772928e-01,9.1
876 7.022927e+01,1.487823e+02,2.562178e+02,3.4
877 2.830064e+02,2.087710e+02,1.575446e+02,1.2
878 5.615029e+01,4.876829e+01,4.314714e+01,3.8
879 2.904621e+01,2.767098e+01,2.653307e+01,2.5
880 2.293024e+01,2.246069e+01,2.204342e+01,2.1
881 2.050838e+01,2.027875e+01,2.006688e+01,1.9
882 1.919382e+01,1.904422e+01,1.889988e+01,1.8
883 1.822798e+01,1.809938e+01,1.797173e+01,1.7
884 1.734050e+01,1.721479e+01,1.708912e+01,1.6
885 1.646175e+01,1.633671e+01,1.621193e+01,1.6
886 1.559442e+01,1.547270e+01,1.535176e+01,1.5
887 1.476153e+01,1.464698e+01,1.453380e+01,1.4
888 1.399216e+01,1.388949e+01,1.378900e+01,1.3
889 1.332489e+01,1.324103e+01,1.316065e+01,1.3
890
891 static const G4double SL45[nL]={
892 5.274378e-02,1.333110e-01,3.641511e-01,1.2
893 9.639532e+01,1.976332e+02,3.177821e+02,3.9
894 2.539062e+02,1.908284e+02,1.468253e+02,1.1
895 5.636581e+01,4.951281e+01,4.426627e+01,4.0
896 3.103504e+01,2.974522e+01,2.867873e+01,2.7
897 2.529440e+01,2.484848e+01,2.444965e+01,2.4
898 2.294188e+01,2.270688e+01,2.248667e+01,2.2
899 2.153584e+01,2.136504e+01,2.119811e+01,2.1
900 2.039825e+01,2.024210e+01,2.008657e+01,1.9
901 1.931412e+01,1.916036e+01,1.900686e+01,1.8
902 1.824507e+01,1.809433e+01,1.794431e+01,1.7
903 1.720795e+01,1.706403e+01,1.692142e+01,1.6
904 1.623135e+01,1.609859e+01,1.596783e+01,1.5
905 1.534810e+01,1.523193e+01,1.511872e+01,1.5
906 1.460399e+01,1.451294e+01,1.442648e+01,1.4
907
908 static const G4double SL46[nL]={
909 5.429151e-02,1.373796e-01,3.766506e-01,1.2
910 1.001858e+02,2.044221e+02,3.254972e+02,3.9
911 2.490584e+02,1.878782e+02,1.450486e+02,1.1
912 5.634628e+01,4.957827e+01,4.439170e+01,4.0
913 3.129630e+01,3.001910e+01,2.896297e+01,2.8
914 2.560848e+01,2.516547e+01,2.476879e+01,2.4
915 2.326261e+01,2.302646e+01,2.280470e+01,2.2
916 2.184156e+01,2.166762e+01,2.149737e+01,2.1
917 2.067932e+01,2.051935e+01,2.035997e+01,2.0
918 1.956833e+01,1.941080e+01,1.925358e+01,1.9
919 1.847398e+01,1.831986e+01,1.816653e+01,1.8
920 1.741470e+01,1.726790e+01,1.712250e+01,1.6
921 1.641963e+01,1.628456e+01,1.615157e+01,1.6
922 1.552206e+01,1.540424e+01,1.528948e+01,1.5
923 1.476885e+01,1.467703e+01,1.458997e+01,1.4
924
925 static const G4double SL47[nL]={
926 5.586443e-02,1.415214e-01,3.894297e-01,1.3
927 1.523387e+01,4.273995e+01,1.039652e+02,2.1
928 3.922061e+02,3.966237e+02,4.015659e+02,3.7
929 2.186308e+02,1.651709e+02,1.280396e+02,1.0
930 6.954857e+01,5.943717e+01,5.185713e+01,4.6
931 3.832984e+01,3.568426e+01,3.359747e+01,3.1
932 2.947046e+01,2.854456e+01,2.776093e+01,2.7
933 2.599335e+01,2.553984e+01,2.513469e+01,2.4
934 2.413690e+01,2.385890e+01,2.360137e+01,2.3
935 2.292175e+01,2.271830e+01,2.252324e+01,2.2
936 2.197463e+01,2.180035e+01,2.162898e+01,2.1
937 2.112668e+01,2.096180e+01,2.079771e+01,2.0
938 2.030859e+01,2.014627e+01,1.998426e+01,1.9
939 1.950011e+01,1.933948e+01,1.917932e+01,1.9
940 1.870239e+01,1.854486e+01,1.838819e+01,1.8
941 1.792421e+01,1.777185e+01,1.762079e+01,1.7
942 1.717627e+01,1.703129e+01,1.688806e+01,1.6
943 1.646987e+01,1.633466e+01,1.620173e+01,1.6
944 1.581792e+01,1.569546e+01,1.557599e+01,1.5
945 1.523747e+01,1.513195e+01,1.503048e+01,1.4
946 1.475318e+01,1.467081e+01,1.459405e+01,1.4
947
948 static const G4double SL48[nL]={
949 5.639434e-02,1.429184e-01,3.937524e-01,1.3
950 1.543710e+01,4.331539e+01,1.053565e+02,2.1
951 4.008312e+02,4.175491e+02,4.294249e+02,3.7
952 2.137731e+02,1.618452e+02,1.258669e+02,1.0
953 6.893285e+01,5.903288e+01,5.159623e+01,4.5
954 3.828743e+01,3.567806e+01,3.361799e+01,3.1
955 2.953782e+01,2.862094e+01,2.784434e+01,2.7
956 2.608997e+01,2.563908e+01,2.523594e+01,2.4
957 2.424142e+01,2.396383e+01,2.370644e+01,2.3
958 2.302601e+01,2.282197e+01,2.262619e+01,2.2
959 2.207484e+01,2.189950e+01,2.172701e+01,2.1
960 2.122110e+01,2.105496e+01,2.088961e+01,2.0
961 2.039663e+01,2.023303e+01,2.006974e+01,1.9
962 1.958183e+01,1.941998e+01,1.925861e+01,1.9
963 1.877819e+01,1.861953e+01,1.846176e+01,1.8
964 1.799463e+01,1.784127e+01,1.768923e+01,1.7
965 1.724194e+01,1.709608e+01,1.695201e+01,1.6
966 1.653144e+01,1.639550e+01,1.626187e+01,1.6
967 1.587614e+01,1.575310e+01,1.563310e+01,1.5
968 1.529319e+01,1.518730e+01,1.508550e+01,1.4
969 1.480751e+01,1.472502e+01,1.464821e+01,1.4
970
971 static const G4double SH0[nH]={
972 1.718841e-05,1.912141e-05,2.128656e-05,2.3
973 3.873597e-05,4.577051e-05,5.661516e-05,7.5
974 5.948094e-04,1.184449e-03,2.411855e-03,4.9
975 5.419455e-02,7.777948e-02,1.011811e-01,1.2
976 2.057743e-01,2.307930e-01,2.589428e-01,2.9
977 4.481828e-01,4.887166e-01,5.240358e-01,5.5
978 5.206446e-01,4.856647e-01,4.465759e-01,4.0
979 2.745971e-01,2.511726e-01,2.314485e-01,2.1
980 1.790291e-01,1.764601e-01,1.763488e-01,1.7
981 2.141163e-01,2.289594e-01,2.451064e-01,2.6
982 2.851743e-01,2.756664e-01,2.618121e-01,2.4
983 2.017331e-01,2.014838e-01,2.059886e-01,2.1
984 2.023557e-01,1.911989e-01,1.808918e-01,1.7
985 1.541608e-01,1.524546e-01,1.513189e-01,1.5
986 1.480828e-01,1.474096e-01,1.467148e-01,1.4
987 1.433514e-01,1.427339e-01,1.421336e-01,1.4
988 1.393046e-01,1.387609e-01,1.382221e-01,1.3
989 1.355952e-01,1.350837e-01,1.345775e-01,1.3
990 1.321338e-01,1.316644e-01,1.312019e-01,1.3
991 1.289978e-01,1.285792e-01,1.281681e-01,1.2
992 1.262246e-01,1.258580e-01,1.254987e-01,1.2
993 1.238082e-01,1.234908e-01,1.231801e-01,1.2
994 1.217239e-01,1.214515e-01,1.211852e-01,1.2
995 1.199423e-01,1.197109e-01,1.194850e-01,1.1
996 1.184365e-01,1.182425e-01,1.180536e-01,1.1
997 1.171834e-01,1.170239e-01,1.168690e-01,1.1
998 1.161639e-01,1.160364e-01,1.159132e-01,1.1
999 1.153622e-01,1.152646e-01,1.151712e-01,1.1
1000 1.147652e-01,1.146960e-01,1.146307e-01,1.
1001 1.143621e-01,1.143198e-01,1.142812e-01,1.
1002 1.141437e-01,1.141271e-01,1.141140e-01,1.
1003 1.141019e-01,1.141099e-01,1.141214e-01,1.
1004
1005 static const G4double SH1[nH]={
1006 6.668702e-02,6.471599e-02,6.280838e-02,6.
1007 5.422715e-02,5.275790e-02,5.145368e-02,5.
1008 5.988460e-02,7.281273e-02,9.420106e-02,1.
1009 2.385444e-01,2.610207e-01,2.835546e-01,3.
1010 4.304707e-01,4.707493e-01,5.154414e-01,5.
1011 7.879273e-01,8.402160e-01,8.843244e-01,9.
1012 8.836029e-01,8.406235e-01,7.905189e-01,7.
1013 5.417618e-01,5.037798e-01,4.708677e-01,4.
1014 3.766827e-01,3.705478e-01,3.683091e-01,3.
1015 4.108509e-01,4.276609e-01,4.452665e-01,4.
1016 4.860977e-01,4.762589e-01,4.610644e-01,4.
1017 3.693549e-01,3.587356e-01,3.525376e-01,3.
1018 3.303977e-01,3.210857e-01,3.125479e-01,3.
1019 2.883070e-01,2.858170e-01,2.837185e-01,2.
1020 2.769203e-01,2.760078e-01,2.751677e-01,2.
1021 2.713259e-01,2.705376e-01,2.697270e-01,2.
1022 2.653104e-01,2.643629e-01,2.633994e-01,2.
1023 2.584370e-01,2.574337e-01,2.564319e-01,2.
1024 2.515094e-01,2.505518e-01,2.496056e-01,2.
1025 2.450708e-01,2.442063e-01,2.433567e-01,2.
1026 2.393357e-01,2.385771e-01,2.378336e-01,2.
1027 2.343390e-01,2.336837e-01,2.330425e-01,2.
1028 2.300425e-01,2.294824e-01,2.289352e-01,2.
1029 2.263859e-01,2.259122e-01,2.254502e-01,2.
1030 2.233091e-01,2.229137e-01,2.225289e-01,2.
1031 2.207590e-01,2.204350e-01,2.201206e-01,2.
1032 2.186908e-01,2.184324e-01,2.181831e-01,2.
1033 2.170676e-01,2.168702e-01,2.166812e-01,2.
1034 2.158589e-01,2.157185e-01,2.155860e-01,2.
1035 2.150393e-01,2.149526e-01,2.148734e-01,2.
1036 2.145873e-01,2.145516e-01,2.145230e-01,2.
1037 2.144847e-01,2.144976e-01,2.145173e-01,2.
1038
1039 static const G4double SH2[nH]={
1040 1.542383e-01,1.519749e-01,1.500571e-01,1.
1041 1.469146e-01,1.482202e-01,1.505946e-01,1.
1042 1.913930e-01,2.082703e-01,2.286959e-01,2.
1043 3.752810e-01,4.110019e-01,4.485617e-01,4.
1044 6.764587e-01,7.332866e-01,7.946123e-01,8.
1045 1.143552e+00,1.207044e+00,1.260057e+00,1.
1046 1.264525e+00,1.214592e+00,1.155254e+00,1.
1047 8.454562e-01,7.961371e-01,7.530613e-01,7.
1048 6.293909e-01,6.215310e-01,6.186705e-01,6.
1049 6.619737e-01,6.759222e-01,6.883771e-01,6.
1050 6.786831e-01,6.594410e-01,6.363628e-01,6.
1051 5.122400e-01,4.921277e-01,4.744756e-01,4.
1052 4.199996e-01,4.145245e-01,4.103489e-01,4.
1053 4.023771e-01,4.022854e-01,4.024166e-01,4.
1054 4.038884e-01,4.040156e-01,4.040242e-01,4.
1055 4.020532e-01,4.012654e-01,4.003614e-01,3.
1056 3.944278e-01,3.930263e-01,3.915747e-01,3.
1057 3.838335e-01,3.822342e-01,3.806303e-01,3.
1058 3.726773e-01,3.711210e-01,3.695814e-01,3.
1059 3.621847e-01,3.607727e-01,3.593848e-01,3.
1060 3.528157e-01,3.515768e-01,3.503629e-01,3.
1061 3.446616e-01,3.435935e-01,3.425488e-01,3.
1062 3.376660e-01,3.367553e-01,3.358660e-01,3.
1063 3.317265e-01,3.309579e-01,3.302084e-01,3.
1064 3.267373e-01,3.260964e-01,3.254728e-01,3.
1065 3.226036e-01,3.220781e-01,3.215681e-01,3.
1066 3.192451e-01,3.188245e-01,3.184182e-01,3.
1067 3.165949e-01,3.162708e-01,3.159598e-01,3.
1068 3.145978e-01,3.143631e-01,3.141406e-01,3.
1069 3.132080e-01,3.130567e-01,3.129170e-01,3.
1070 3.123871e-01,3.123143e-01,3.122524e-01,3.
1071 3.121025e-01,3.121040e-01,3.121157e-01,3.
1072
1073 static const G4double SH3[nH]={
1074 2.629521e-01,2.526203e-01,2.431577e-01,2.
1075 2.107267e-01,2.081600e-01,2.075143e-01,2.
1076 2.514498e-01,2.737640e-01,3.014354e-01,3.
1077 5.060360e-01,5.563139e-01,6.091512e-01,6.
1078 9.259290e-01,1.003412e+00,1.086287e+00,1.
1079 1.544923e+00,1.626692e+00,1.695047e+00,1.
1080 1.714557e+00,1.654820e+00,1.582408e+00,1.
1081 1.186603e+00,1.120659e+00,1.062094e+00,1.
1082 8.805089e-01,8.646577e-01,8.547494e-01,8.
1083 8.728198e-01,8.833474e-01,8.930526e-01,9.
1084 8.812964e-01,8.626545e-01,8.396586e-01,8.
1085 7.038504e-01,6.797291e-01,6.579596e-01,6.
1086 5.854904e-01,5.773026e-01,5.707205e-01,5.
1087 5.543723e-01,5.531226e-01,5.522016e-01,5.
1088 5.492791e-01,5.485986e-01,5.478012e-01,5.
1089 5.417006e-01,5.400585e-01,5.382913e-01,5.
1090 5.279674e-01,5.256850e-01,5.233556e-01,5.
1091 5.113224e-01,5.088925e-01,5.064694e-01,5.
1092 4.946129e-01,4.923172e-01,4.900525e-01,4.
1093 4.792477e-01,4.771973e-01,4.751850e-01,4.
1094 4.657009e-01,4.639189e-01,4.621746e-01,4.
1095 4.540052e-01,4.524784e-01,4.509863e-01,4.
1096 4.440254e-01,4.427295e-01,4.414644e-01,4.
1097 4.355835e-01,4.344928e-01,4.334295e-01,4.
1098 4.285079e-01,4.275996e-01,4.267157e-01,4.
1099 4.226493e-01,4.219041e-01,4.211809e-01,4.
1100 4.178831e-01,4.172851e-01,4.167071e-01,4.
1101 4.141068e-01,4.136430e-01,4.131974e-01,4.
1102 4.112360e-01,4.108956e-01,4.105721e-01,4.
1103 4.092013e-01,4.089753e-01,4.087651e-01,4.
1104 4.079445e-01,4.078256e-01,4.077213e-01,4.
1105 4.074170e-01,4.073987e-01,4.073943e-01,4.
1106
1107 static const G4double SH4[nH]={
1108 7.405778e-02,7.529642e-02,7.695159e-02,7.
1109 9.636447e-02,1.043620e-01,1.148402e-01,1.
1110 2.353712e-01,2.797917e-01,3.322864e-01,3.
1111 6.910366e-01,7.754430e-01,8.630582e-01,9.
1112 1.368302e+00,1.487447e+00,1.613306e+00,1.
1113 2.288176e+00,2.406638e+00,2.506475e+00,2.
1114 2.569112e+00,2.493232e+00,2.398030e+00,2.
1115 1.841504e+00,1.743297e+00,1.654408e+00,1.
1116 1.357542e+00,1.325684e+00,1.301743e+00,1.
1117 1.272221e+00,1.276770e+00,1.281590e+00,1.
1118 1.260999e+00,1.242178e+00,1.218585e+00,1.
1119 1.065486e+00,1.035376e+00,1.007278e+00,9.
1120 9.049702e-01,8.919552e-01,8.810082e-01,8.
1121 8.487605e-01,8.452151e-01,8.421882e-01,8.
1122 8.303629e-01,8.280477e-01,8.256249e-01,8.
1123 8.112975e-01,8.079728e-01,8.045118e-01,8.
1124 7.856415e-01,7.816556e-01,7.776315e-01,7.
1125 7.573474e-01,7.533262e-01,7.493353e-01,7.
1126 7.300255e-01,7.263208e-01,7.226752e-01,7.
1127 7.053878e-01,7.021244e-01,6.989265e-01,6.
1128 6.839097e-01,6.810975e-01,6.783474e-01,6.
1129 6.654984e-01,6.631024e-01,6.607622e-01,6.
1130 6.498624e-01,6.478361e-01,6.458588e-01,6.
1131 6.366764e-01,6.349747e-01,6.333162e-01,6.
1132 6.256428e-01,6.242270e-01,6.228493e-01,6.
1133 6.165095e-01,6.153472e-01,6.142189e-01,6.
1134 6.090686e-01,6.081335e-01,6.072289e-01,6.
1135 6.031506e-01,6.024211e-01,6.017194e-01,6.
1136 5.986171e-01,5.980756e-01,5.975596e-01,5.
1137 5.953544e-01,5.949863e-01,5.946420e-01,5.
1138 5.932682e-01,5.930615e-01,5.928769e-01,5.
1139 5.922798e-01,5.922242e-01,5.921894e-01,5.
1140
1141 static const G4double SH5[nH]={
1142 4.659776e-01,4.476902e-01,4.309775e-01,4.
1143 3.750402e-01,3.713091e-01,3.712933e-01,3.
1144 4.625322e-01,5.066943e-01,5.609019e-01,6.
1145 9.505815e-01,1.043997e+00,1.141265e+00,1.
1146 1.705950e+00,1.839605e+00,1.980828e+00,2.
1147 2.736990e+00,2.869544e+00,2.981348e+00,3.
1148 3.052863e+00,2.967066e+00,2.858644e+00,2.
1149 2.211463e+00,2.094818e+00,1.988432e+00,1.
1150 1.622712e+00,1.580840e+00,1.547984e+00,1.
1151 1.487081e+00,1.487447e+00,1.488669e+00,1.
1152 1.457487e+00,1.437429e+00,1.412571e+00,1.
1153 1.249221e+00,1.216200e+00,1.185043e+00,1.
1154 1.067888e+00,1.052343e+00,1.039061e+00,1.
1155 9.978077e-01,9.928816e-01,9.885710e-01,9.
1156 9.712366e-01,9.679192e-01,9.645088e-01,9.
1157 9.453386e-01,9.410538e-01,9.366353e-01,9.
1158 9.130459e-01,9.081377e-01,9.032012e-01,8.
1159 8.785340e-01,8.736773e-01,8.688658e-01,8.
1160 8.456872e-01,8.412565e-01,8.369010e-01,8.
1161 8.162992e-01,8.124188e-01,8.086186e-01,8.
1162 7.908030e-01,7.874715e-01,7.842150e-01,7.
1163 7.690169e-01,7.661858e-01,7.634213e-01,7.
1164 7.505555e-01,7.481653e-01,7.458334e-01,7.
1165 7.350099e-01,7.330050e-01,7.310511e-01,7.
1166 7.220140e-01,7.203469e-01,7.187247e-01,7.
1167 7.112604e-01,7.098919e-01,7.085633e-01,7.
1168 7.024976e-01,7.013959e-01,7.003300e-01,6.
1169 6.955214e-01,6.946606e-01,6.938324e-01,6.
1170 6.901658e-01,6.895248e-01,6.889136e-01,6.
1171 6.862947e-01,6.858561e-01,6.854452e-01,6.
1172 6.837958e-01,6.835452e-01,6.833203e-01,6.
1173 6.825752e-01,6.825004e-01,6.824498e-01,6.
1174
1175 static const G4double SH6[nH]={
1176 5.445765e-01,5.259720e-01,5.092147e-01,4.
1177 4.589754e-01,4.583113e-01,4.624465e-01,4.
1178 6.024632e-01,6.642824e-01,7.390624e-01,8.
1179 1.261606e+00,1.384524e+00,1.511770e+00,1.
1180 2.237925e+00,2.407239e+00,2.585320e+00,2.
1181 3.530524e+00,3.696566e+00,3.837865e+00,3.
1182 3.957835e+00,3.858084e+00,3.728662e+00,3.
1183 2.918748e+00,2.767352e+00,2.627710e+00,2.
1184 2.129397e+00,2.068133e+00,2.018030e+00,1.
1185 1.897116e+00,1.889649e+00,1.884158e+00,1.
1186 1.832895e+00,1.810468e+00,1.783286e+00,1.
1187 1.602978e+00,1.565190e+00,1.528974e+00,1.
1188 1.386391e+00,1.366358e+00,1.348893e+00,1.
1189 1.291059e+00,1.283522e+00,1.276765e+00,1.
1190 1.248754e+00,1.243470e+00,1.238117e+00,1.
1191 1.209386e+00,1.203202e+00,1.196890e+00,1.
1192 1.163967e+00,1.157235e+00,1.150494e+00,1.
1193 1.117164e+00,1.110656e+00,1.104223e+00,1.
1194 1.073400e+00,1.067535e+00,1.061776e+00,1.
1195 1.034624e+00,1.029523e+00,1.024532e+00,1.
1196 1.001182e+00,9.968228e-01,9.925643e-01,9.
1197 9.727164e-01,9.690235e-01,9.654188e-01,9.
1198 9.486577e-01,9.455463e-01,9.425115e-01,9.
1199 9.284330e-01,9.258264e-01,9.232864e-01,9.
1200 9.115411e-01,9.093748e-01,9.072668e-01,9.
1201 8.975665e-01,8.957878e-01,8.940607e-01,8.
1202 8.861720e-01,8.847382e-01,8.833507e-01,8.
1203 8.770849e-01,8.759618e-01,8.748805e-01,8.
1204 8.700849e-01,8.692443e-01,8.684421e-01,8.
1205 8.649917e-01,8.644108e-01,8.638654e-01,8.
1206 8.616572e-01,8.613169e-01,8.610097e-01,8.
1207 8.599577e-01,8.598422e-01,8.597575e-01,8.
1208
1209 static const G4double SH7[nH]={
1210 6.601789e-01,6.429790e-01,6.279582e-01,6.
1211 5.941728e-01,5.994893e-01,6.111752e-01,6.
1212 8.300428e-01,9.191103e-01,1.025252e+00,1.
1213 1.744040e+00,1.909484e+00,2.079667e+00,2.
1214 3.034365e+00,3.253973e+00,3.484139e+00,3.
1215 4.701980e+00,4.918458e+00,5.105403e+00,5.
1216 5.320522e+00,5.205479e+00,5.049727e+00,4.
1217 4.008169e+00,3.804492e+00,3.614203e+00,3.
1218 2.908433e+00,2.816062e+00,2.738112e+00,2.
1219 2.518816e+00,2.498248e+00,2.481393e+00,2.
1220 2.393491e+00,2.366348e+00,2.334725e+00,2.
1221 2.128833e+00,2.084617e+00,2.041630e+00,2.
1222 1.864390e+00,1.838007e+00,1.814532e+00,1.
1223 1.731819e+00,1.720184e+00,1.709544e+00,1.
1224 1.664220e+00,1.655719e+00,1.647190e+00,1.
1225 1.602933e+00,1.593696e+00,1.584347e+00,1.
1226 1.536603e+00,1.526999e+00,1.517424e+00,1.
1227 1.470568e+00,1.461496e+00,1.452549e+00,1.
1228 1.409919e+00,1.401847e+00,1.393931e+00,1.
1229 1.356733e+00,1.349767e+00,1.342956e+00,1.
1230 1.311158e+00,1.305234e+00,1.299449e+00,1.
1231 1.272529e+00,1.267527e+00,1.262647e+00,1.
1232 1.239975e+00,1.235770e+00,1.231670e+00,1.
1233 1.212660e+00,1.209142e+00,1.205714e+00,1.
1234 1.189869e+00,1.186947e+00,1.184104e+00,1.
1235 1.171019e+00,1.168619e+00,1.166288e+00,1.
1236 1.155638e+00,1.153701e+00,1.151825e+00,1.
1237 1.143348e+00,1.141826e+00,1.140360e+00,1.
1238 1.133844e+00,1.132699e+00,1.131605e+00,1.
1239 1.126879e+00,1.126079e+00,1.125326e+00,1.
1240 1.122249e+00,1.121768e+00,1.121331e+00,1.
1241 1.119786e+00,1.119603e+00,1.119460e+00,1.
1242
1243 static const G4double SH8[nH]={
1244 8.326620e-01,8.187860e-01,8.074156e-01,7.
1245 8.004391e-01,8.143676e-01,8.366683e-01,8.
1246 1.161963e+00,1.287041e+00,1.434275e+00,1.
1247 2.403252e+00,2.621760e+00,2.845287e+00,3.
1248 4.083048e+00,4.365291e+00,4.660621e+00,4.
1249 6.228396e+00,6.512543e+00,6.762334e+00,6.
1250 7.136586e+00,7.009905e+00,6.827389e+00,6.
1251 5.500809e+00,5.228183e+00,4.970100e+00,4.
1252 3.977736e+00,3.841183e+00,3.723389e+00,3.
1253 3.360729e+00,3.320633e+00,3.286663e+00,3.
1254 3.140973e+00,3.106044e+00,3.067232e+00,3.
1255 2.825830e+00,2.773655e+00,2.722401e+00,2.
1256 2.502745e+00,2.468369e+00,2.437227e+00,2.
1257 2.321429e+00,2.304110e+00,2.288023e+00,2.
1258 2.217898e+00,2.204739e+00,2.191611e+00,2.
1259 2.124976e+00,2.111368e+00,2.097687e+00,2.
1260 2.028956e+00,2.015316e+00,2.001766e+00,1.
1261 1.936048e+00,1.923419e+00,1.910988e+00,1.
1262 1.852064e+00,1.840954e+00,1.830074e+00,1.
1263 1.779107e+00,1.769589e+00,1.760290e+00,1.
1264 1.716968e+00,1.708913e+00,1.701051e+00,1.
1265 1.664517e+00,1.657738e+00,1.651125e+00,1.
1266 1.620437e+00,1.614751e+00,1.609207e+00,1.
1267 1.583519e+00,1.578768e+00,1.574140e+00,1.
1268 1.552749e+00,1.548805e+00,1.544967e+00,1.
1269 1.527303e+00,1.524063e+00,1.520915e+00,1.
1270 1.506525e+00,1.503906e+00,1.501370e+00,1.
1271 1.489891e+00,1.487827e+00,1.485838e+00,1.
1272 1.476979e+00,1.475418e+00,1.473924e+00,1.
1273 1.467449e+00,1.466346e+00,1.465306e+00,1.
1274 1.461021e+00,1.460342e+00,1.459721e+00,1.
1275 1.457464e+00,1.457179e+00,1.456948e+00,1.
1276
1277 static const G4double SH9[nH]={
1278 1.425410e+00,1.421381e+00,1.420064e+00,1.
1279 1.477877e+00,1.512528e+00,1.560992e+00,1.
1280 2.165508e+00,2.386174e+00,2.641752e+00,2.
1281 4.250887e+00,4.602541e+00,4.959837e+00,5.
1282 6.915215e+00,7.358674e+00,7.822640e+00,8.
1283 1.031441e+01,1.078267e+01,1.120652e+01,1.
1284 1.208752e+01,1.195370e+01,1.172359e+01,1.
1285 9.710864e+00,9.256615e+00,8.816313e+00,8.
1286 7.018988e+00,6.753498e+00,6.518232e+00,6.
1287 5.721783e+00,5.620790e+00,5.532858e+00,5.
1288 5.197938e+00,5.136693e+00,5.073880e+00,5.
1289 4.727591e+00,4.654788e+00,4.582647e+00,4.
1290 4.257712e+00,4.203068e+00,4.152199e+00,4.
1291 3.946978e+00,3.913492e+00,3.881699e+00,3.
1292 3.737984e+00,3.710781e+00,3.683763e+00,3.
1293 3.549640e+00,3.522927e+00,3.496284e+00,3.
1294 3.365257e+00,3.339727e+00,3.314496e+00,3.
1295 3.193692e+00,3.170733e+00,3.148206e+00,3.
1296 3.042254e+00,3.022416e+00,3.003026e+00,2.
1297 2.912654e+00,2.895854e+00,2.879463e+00,2.
1298 2.803363e+00,2.789258e+00,2.775505e+00,2.
1299 2.711744e+00,2.699939e+00,2.688431e+00,2.
1300 2.635115e+00,2.625250e+00,2.615637e+00,2.
1301 2.571141e+00,2.562918e+00,2.554909e+00,2.
1302 2.517912e+00,2.511092e+00,2.504456e+00,2.
1303 2.473909e+00,2.468302e+00,2.462856e+00,2.
1304 2.437930e+00,2.433387e+00,2.428985e+00,2.
1305 2.409023e+00,2.405425e+00,2.401953e+00,2.
1306 2.386432e+00,2.383683e+00,2.381049e+00,2.
1307 2.369546e+00,2.367569e+00,2.365697e+00,2.
1308 2.357868e+00,2.356600e+00,2.355429e+00,2.
1309 2.350989e+00,2.350378e+00,2.349856e+00,2.
1310
1311 static const G4double SH10[nH]={
1312 3.918292, 3.904931, 3.893792, 3.886847, 3
1313 3.973155, 4.052892, 4.173448, 4.346251, 4
1314 5.767472, 6.325587, 6.952077, 7.631192, 8
1315 10.58224, 11.34205, 12.11228, 12.89914, 1
1316 16.35126, 17.31704, 18.32759, 19.37949, 2
1317 2.377554e+01,2.481942e+01,2.578236e+01,2.
1318 2.821866e+01,2.807489e+01,2.771701e+01,2.
1319 2.369877e+01,2.269021e+01,2.168502e+01,2.
1320 1.728387e+01,1.658334e+01,1.594670e+01,1.
1321 1.361791e+01,1.329360e+01,1.300539e+01,1.
1322 1.194464e+01,1.178007e+01,1.162372e+01,1.
1323 1.090421e+01,1.076726e+01,1.063254e+01,1.
1324 1.000757e+01,9.894853e+00,9.786587e+00,9.
1325 9.304761e+00,9.217824e+00,9.133076e+00,9.
1326 8.730205e+00,8.652092e+00,8.574465e+00,8.
1327 8.192917e+00,8.118146e+00,8.044035e+00,7.
1328 7.686057e+00,7.617475e+00,7.550037e+00,7.
1329 7.231395e+00,7.171554e+00,7.113040e+00,7.
1330 6.840281e+00,6.789624e+00,6.740229e+00,6.
1331 6.511433e+00,6.469143e+00,6.427951e+00,6.
1332 6.237554e+00,6.202409e+00,6.168183e+00,6.
1333 6.010017e+00,5.980819e+00,5.952382e+00,5.
1334 5.820930e+00,5.796658e+00,5.773017e+00,5.
1335 5.663751e+00,5.643583e+00,5.623944e+00,5.
1336 5.533276e+00,5.516567e+00,5.500309e+00,5.
1337 5.425436e+00,5.411684e+00,5.398319e+00,5.
1338 5.337056e+00,5.325867e+00,5.315016e+00,5.
1339 5.265647e+00,5.256711e+00,5.248076e+00,5.
1340 5.209243e+00,5.202315e+00,5.195656e+00,5.
1341 5.166273e+00,5.161154e+00,5.156279e+00,5.
1342 5.135467e+00,5.131996e+00,5.128750e+00,5.
1343 5.115788e+00,5.113833e+00,5.112084e+00,5.
1344
1345 static const G4double SH11[nH]={
1346 7.590321e+00,7.509120e+00,7.439927e+00,7.
1347 7.378718e+00,7.443759e+00,7.577042e+00,7.
1348 8.587142e+00,9.191438e+00,9.948696e+00,1.
1349 1.303902e+01,1.425810e+01,1.552108e+01,1.
1350 1.938727e+01,2.068629e+01,2.200102e+01,2.
1351 2.616717e+01,2.767219e+01,2.925247e+01,3.
1352 3.448491e+01,3.638304e+01,3.833817e+01,4.
1353 4.426974e+01,4.613359e+01,4.785354e+01,4.
1354 5.234188e+01,5.215066e+01,5.158403e+01,5.
1355 4.466144e+01,4.285783e+01,4.103934e+01,3.
1356 3.281973e+01,3.146624e+01,3.022230e+01,2.
1357 2.552116e+01,2.484197e+01,2.423289e+01,2.
1358 2.198543e+01,2.164944e+01,2.133824e+01,2.
1359 2.002347e+01,1.979144e+01,1.956643e+01,1.
1360 1.853145e+01,1.834084e+01,1.815514e+01,1.
1361 1.728639e+01,1.712128e+01,1.695797e+01,1.
1362 1.615706e+01,1.599876e+01,1.584098e+01,1.
1363 1.506367e+01,1.491163e+01,1.476115e+01,1.
1364 1.403806e+01,1.390031e+01,1.376511e+01,1.
1365 1.312955e+01,1.301078e+01,1.289481e+01,1.
1366 1.235641e+01,1.225680e+01,1.215977e+01,1.
1367 1.171172e+01,1.162914e+01,1.154878e+01,1.
1368 1.117817e+01,1.110991e+01,1.104348e+01,1.
1369 1.073701e+01,1.068053e+01,1.062555e+01,1.
1370 1.037170e+01,1.032488e+01,1.027930e+01,1.
1371 1.006876e+01,1.002992e+01,9.992112e+00,9.
1372 9.817586e+00,9.785425e+00,9.754129e+00,9.
1373 9.609955e+00,9.583456e+00,9.557697e+00,9.
1374 9.439481e+00,9.417856e+00,9.396874e+00,9.
1375 9.301186e+00,9.283817e+00,9.267013e+00,9.
1376 9.191149e+00,9.177546e+00,9.164447e+00,9.
1377 9.106252e+00,9.096024e+00,9.086249e+00,9.
1378 9.043985e+00,9.036813e+00,9.030055e+00,9.
1379 9.002302e+00,8.997926e+00,8.993930e+00,8.
1380
1381 static const G4double SH12[nH]={
1382 1.274173e+01,1.261154e+01,1.253680e+01,1.
1383 1.286155e+01,1.324712e+01,1.382513e+01,1.
1384 1.694815e+01,1.845745e+01,2.015574e+01,2.
1385 2.587848e+01,2.784823e+01,2.980938e+01,3.
1386 3.568546e+01,3.769757e+01,3.977278e+01,4.
1387 4.658588e+01,4.910951e+01,5.177874e+01,5.
1388 6.067137e+01,6.389725e+01,6.721344e+01,7.
1389 7.721799e+01,8.034325e+01,8.322076e+01,8.
1390 8.942639e+01,9.041346e+01,9.077151e+01,9.
1391 8.811385e+01,8.614235e+01,8.376122e+01,8.
1392 7.511689e+01,7.203142e+01,6.896584e+01,6.
1393 6.036126e+01,5.778900e+01,5.538996e+01,5.
1394 4.925181e+01,4.754501e+01,4.599451e+01,4.
1395 4.217515e+01,4.114200e+01,4.021009e+01,3.
1396 3.791594e+01,3.728600e+01,3.670900e+01,3.
1397 3.522478e+01,3.479282e+01,3.438448e+01,3.
1398 3.326834e+01,3.292429e+01,3.259109e+01,3.
1399 3.164275e+01,3.133978e+01,3.104162e+01,3.
1400 3.016670e+01,2.987885e+01,2.959177e+01,2.
1401 2.873069e+01,2.844289e+01,2.815459e+01,2.
1402 2.728821e+01,2.699984e+01,2.671229e+01,2.
1403 2.585875e+01,2.557869e+01,2.530156e+01,2.
1404 2.449130e+01,2.422931e+01,2.397182e+01,2.
1405 2.322821e+01,2.299044e+01,2.275787e+01,2.
1406 2.209182e+01,2.188038e+01,2.167419e+01,2.
1407 2.108662e+01,2.090087e+01,2.072003e+01,2.
1408 2.020604e+01,2.004387e+01,1.988611e+01,1.
1409 1.943817e+01,1.929694e+01,1.915958e+01,1.
1410 1.876959e+01,1.864662e+01,1.852699e+01,1.
1411 1.818725e+01,1.808008e+01,1.797579e+01,1.
1412 1.767946e+01,1.758593e+01,1.749489e+01,1.
1413 1.723608e+01,1.715436e+01,1.707480e+01,1.
1414 1.684856e+01,1.677711e+01,1.670754e+01,1.
1415 1.650973e+01,1.644726e+01,1.638646e+01,1.
1416 1.621363e+01,1.615908e+01,1.610601e+01,1.
1417 1.595530e+01,1.590779e+01,1.586160e+01,1.
1418 1.573061e+01,1.568940e+01,1.564937e+01,1.
1419 1.553612e+01,1.550058e+01,1.546611e+01,1.
1420 1.536892e+01,1.533854e+01,1.530912e+01,1.
1421 1.522657e+01,1.520089e+01,1.517611e+01,1.
1422 1.510699e+01,1.508564e+01,1.506512e+01,1.
1423 1.500840e+01,1.499106e+01,1.497449e+01,1.
1424 1.492928e+01,1.491568e+01,1.490279e+01,1.
1425 1.486832e+01,1.485821e+01,1.484876e+01,1.
1426 1.482437e+01,1.481753e+01,1.481133e+01,1.
1427
1428 static const G4double SH13[nH]={
1429 1.444282e+01,1.433200e+01,1.430197e+01,1.
1430 1.497946e+01,1.557548e+01,1.641373e+01,1.
1431 2.051932e+01,2.236282e+01,2.436683e+01,2.
1432 3.078116e+01,3.292713e+01,3.505686e+01,3.
1433 4.148940e+01,4.372398e+01,4.604543e+01,4.
1434 5.375059e+01,5.662368e+01,5.966783e+01,6.
1435 6.982292e+01,7.350621e+01,7.729059e+01,8.
1436 8.869083e+01,9.224647e+01,9.551803e+01,9.
1437 1.025696e+02,1.036932e+02,1.041041e+02,1.
1438 1.011066e+02,9.887372e+01,9.617299e+01,9.
1439 8.633954e+01,8.281934e+01,7.931645e+01,7.
1440 6.945404e+01,6.649625e+01,6.373321e+01,6.
1441 5.664119e+01,5.466241e+01,5.286187e+01,5.
1442 4.841316e+01,4.720643e+01,4.611690e+01,4.
1443 4.343327e+01,4.269734e+01,4.202440e+01,4.
1444 4.030251e+01,3.980491e+01,3.933629e+01,3.
1445 3.806404e+01,3.767398e+01,3.729683e+01,3.
1446 3.622434e+01,3.588127e+01,3.554314e+01,3.
1447 3.454699e+01,3.421782e+01,3.388880e+01,3.
1448 3.289811e+01,3.256591e+01,3.223270e+01,3.
1449 3.122945e+01,3.089506e+01,3.056146e+01,3.
1450 2.957065e+01,2.924544e+01,2.892361e+01,2.
1451 2.798266e+01,2.767845e+01,2.737947e+01,2.
1452 2.651624e+01,2.624029e+01,2.597040e+01,2.
1453 2.519768e+01,2.495244e+01,2.471333e+01,2.
1454 2.403210e+01,2.381680e+01,2.360721e+01,2.
1455 2.301165e+01,2.282379e+01,2.264106e+01,2.
1456 2.212233e+01,2.195881e+01,2.179978e+01,2.
1457 2.134838e+01,2.120607e+01,2.106765e+01,2.
1458 2.067458e+01,2.055060e+01,2.042998e+01,2.
1459 2.008725e+01,1.997909e+01,1.987382e+01,1.
1460 1.957458e+01,1.948010e+01,1.938812e+01,1.
1461 1.912659e+01,1.904399e+01,1.896358e+01,1.
1462 1.873493e+01,1.866273e+01,1.859245e+01,1.
1463 1.839266e+01,1.832961e+01,1.826826e+01,1.
1464 1.809402e+01,1.803908e+01,1.798567e+01,1.
1465 1.783418e+01,1.778651e+01,1.774019e+01,1.
1466 1.760914e+01,1.756801e+01,1.752810e+01,1.
1467 1.741553e+01,1.738032e+01,1.734624e+01,1.
1468 1.725050e+01,1.722071e+01,1.719194e+01,1.
1469 1.711164e+01,1.708682e+01,1.706295e+01,1.
1470 1.699688e+01,1.697667e+01,1.695733e+01,1.
1471 1.690446e+01,1.688852e+01,1.687340e+01,1.
1472 1.683285e+01,1.682090e+01,1.680973e+01,1.
1473 1.678072e+01,1.677252e+01,1.676505e+01,1.
1474
1475 static const G4double* SL[nLA]={
1476 SL0, SL1, SL2, SL3, SL4, SL5, SL6, SL7, S
1477 SL14,SL15,SL16,SL17,SL18,SL19,SL20,SL21,S
1478 SL28,SL29,SL30,SL31,SL32,SL33,SL34,SL35,S
1479 SL42,SL43,SL44,SL45,SL46,SL47,SL48};
1480
1481 static const G4double* SH[nHA]={
1482 SH0,SH1,SH2,SH3,SH4,SH5,SH6,SH7,SH8,SH9,S
1483
1484
1485 // Bug 2092 fix (ejc3) Init deuteron_GDR ptr
1486 G4PhotoNuclearCrossSection::G4PhotoNuclearCro
1487 : G4VCrossSectionDataSet(Default_Name()),
1488 mNeut(G4NucleiProperties::GetNuclearMass(1
1489 mProt(G4NucleiProperties::GetNuclearMass(1
1490 {
1491 SetForAllAtomsAndEnergies(true);
1492 nistmngr = G4NistManager::Instance();
1493
1494 GDR.resize(120, nullptr);
1495 HEN.resize(120, nullptr);
1496 spA.resize(120, 0.0);
1497 eTH.resize(120, 0.0);
1498 }
1499
1500
1501 G4PhotoNuclearCrossSection::~G4PhotoNuclearCr
1502 {
1503 for (auto const & ptr : GDR) { delete [] pt
1504 for (auto const & ptr : HEN) { delete [] pt
1505
1506 delete [] deuteron_GDR;
1507 delete [] deuteron_HR;
1508 delete [] triton_GDR;
1509 delete [] triton_HR;
1510 delete [] he3_GDR;
1511 delete [] he3_HR;
1512 }
1513
1514
1515 void
1516 G4PhotoNuclearCrossSection::CrossSectionDescr
1517 {
1518 outFile << "G4PhotoNuclearCrossSection prov
1519 << "cross section for photon intera
1520 << "cross section is a parameteriza
1521 << "all incident gamma energies.\n"
1522 }
1523
1524
1525 // Allow D, T, 3He targets
1526 G4bool
1527 G4PhotoNuclearCrossSection::IsIsoApplicable(c
1528 G
1529 c
1530 c
1531 {
1532 // explicitly allow deuterium, tritium, and
1533 return ((Z == 1 && A == 2) || (Z == 1 && A
1534 }
1535
1536
1537 G4bool
1538 G4PhotoNuclearCrossSection::IsElementApplicab
1539
1540 {
1541 return true;
1542 }
1543
1544
1545 // Get cross sections for deuterium, tritium
1546 G4double
1547 G4PhotoNuclearCrossSection::GetIsoCrossSectio
1548
1549
1550
1551
1552 {
1553 return ComputeIsoXSection(aPart->GetKinetic
1554 }
1555
1556 G4double
1557 G4PhotoNuclearCrossSection::ComputeIsoXSectio
1558 {
1559 if (Energy <= THmin) { return 0.0; }
1560 G4double sigma;
1561 G4double lE;
1562 if (Z == 1 && A == 2) {
1563 // init the XS table if need be
1564 if (deuteron_GDR == nullptr) {
1565 deuteron_TH = ThresholdEnergy(1,1); //
1566 deuteron_GDR = new G4double[nL]; //
1567 for (G4int i = 0; i < nL; ++i) deuteron
1568 deuteron_HR = new G4double[nH];
1569 for (G4int i = 0; i < nH; ++i) deuteron
1570 deuteron_SP = 1; //
1571 }
1572 if (Energy < deuteron_TH) {
1573 sigma = 0.;
1574 } else if (Energy < Emin) { // GDR r
1575 sigma = EquLinearFit(Energy,nL,THmin,dE
1576 } else if (Energy < Emax) { // High
1577 lE = G4Log(Energy);
1578 sigma = EquLinearFit(lE,nH,milE,dlE,deu
1579 } else { // Very
1580 lE = G4Log(Energy);
1581 sigma = deuteron_SP*(poc*(lE-pos)+shd*s
1582 }
1583
1584 } else if (Z == 1 && A == 3) {
1585 if (triton_GDR == nullptr) {
1586 triton_TH = ThresholdEnergy(1,2);
1587 triton_GDR = new G4double[nL]; //
1588 for (G4int i = 0; i < nL; ++i) triton_G
1589 triton_HR = new G4double[nH];
1590 for (G4int i = 0; i < nH; ++i) triton_H
1591 triton_SP = 1;
1592 }
1593 if (Energy < triton_TH) {
1594 sigma = 0.;
1595 } else if (Energy < Emin) { // GDR r
1596 sigma = EquLinearFit(Energy,nL,THmin,dE
1597 } else if (Energy < Emax) { // High
1598 lE = G4Log(Energy);
1599 sigma = EquLinearFit(lE,nH,milE,dlE,tri
1600 } else {
1601 lE = G4Log(Energy);
1602 sigma = triton_SP*(poc*(lE-pos)+shd*std
1603 }
1604
1605 } else if (Z == 2 && A == 3) {
1606 if (he3_GDR == nullptr) {
1607 he3_TH = ThresholdEnergy(2,1);
1608 he3_GDR = new G4double[nL]; // same as
1609 for (G4int i = 0 ; i < nL ; ++i) he3_GD
1610 he3_HR = new G4double[nH];
1611 for (G4int i = 0 ; i < nH ; ++i) he3_HR
1612 he3_SP = 2;
1613 }
1614 if (Energy < he3_TH) {
1615 sigma = 0.;
1616 } else if (Energy < Emin) { // GDR r
1617 sigma = EquLinearFit(Energy,nL,THmin,dE
1618 } else if (Energy < Emax) { // High
1619 lE = G4Log(Energy);
1620 sigma = EquLinearFit(lE,nH,milE,dlE,he3
1621 } else {
1622 lE = G4Log(Energy);
1623 sigma = he3_SP*(poc*(lE-pos)+shd*std::e
1624 }
1625
1626 } else {
1627 return ComputeElementXSection(Energy, Z);
1628 }
1629
1630 if(sigma < 0.) sigma = 0.;
1631 return sigma*millibarn;
1632 }
1633
1634
1635 // The main member function giving the gamma-
1636 // (E in MeV, CS in mb)
1637 G4double
1638 G4PhotoNuclearCrossSection::GetElementCrossSe
1639
1640 {
1641 return ComputeElementXSection(aPart->GetKin
1642 }
1643
1644 G4double
1645 G4PhotoNuclearCrossSection::ComputeElementXSe
1646 {
1647 if (Energy <= THmin) { return 0.0; }
1648 G4int ZZ = std::min(Z, nL - 1);
1649
1650 G4double sigma=0.;
1651
1652 if(ZZ!=lastZ) // Otherwise the set of p
1653 {
1654 lastZ = ZZ; // The last Z
1655 if(GDR[ZZ]) // Calculated nuclei in D
1656 {
1657 lastGDR=GDR[ZZ]; // Pointer to
1658 lastHEN=HEN[ZZ]; // Pointer to
1659 lastTH =eTH[ZZ]; // Energy Thre
1660 lastSP =spA[ZZ]; // Shadowing c
1661 }
1662 else
1663 {
1664 G4double Aa = nistmngr->GetAtomic
1665 G4int N = static_cast<G4int>(Aa)
1666
1667 G4double lnA=G4Log(Aa); // The nu
1668 if(Aa==1.) lastSP=1.; // The
1669 else lastSP=Aa*(1.-shc*lnA);
1670 lastTH=ThresholdEnergy(ZZ, N); //
1671 lastGDR = new G4double[nL]; // A
1672 // GDR cross sections
1673 lastHEN = new G4double[nH]; // A
1674 // HEN cross sections
1675 G4int er=GetFunctions(Aa,lastGDR,
1676 // filling of the functions
1677 if(er<1) G4cerr << "***G4PhotoNuc
1678 << Aa << " failed" << G4endl;
1679
1680 GDR[ZZ]=lastGDR; // added GDR
1681 HEN[ZZ]=lastHEN; // added HEN
1682 eTH[ZZ]=lastTH; // Threshold
1683 spA[ZZ]=lastSP; // Pomeron S
1684 } // End of creation of the new set o
1685 }// End of parameters udate
1686
1687 //
1688 // =================== NOW the Magic Form
1689 //
1690 if (Energy<lastTH)
1691 {
1692 lastE=Energy;
1693 lastSig=0.;
1694 return 0.;
1695 }
1696 else if (Energy<Emin) // GDR region (ap
1697 {
1698 sigma=EquLinearFit(Energy,nL,THmin,dE
1699 }
1700 else if (Energy<Emax)
1701 {
1702 G4double lE=G4Log(Energy);
1703 sigma=EquLinearFit(lE,nH,milE,dlE,las
1704 }
1705 else // UHE region (calcula
1706 {
1707 G4double lE=G4Log(Energy);
1708 sigma=lastSP*(poc*(lE-pos)+shd*G4Exp(
1709 }
1710 // End of "sigma" calculation
1711
1712 sigma = std::max(sigma, 0.);
1713 return sigma*millibarn;
1714 }
1715
1716 // Threshold energy for nuclei: explicitly ca
1717 // min of p- and n-threshold for heavier targ
1718 G4double G4PhotoNuclearCrossSection::Threshol
1719 {
1720 G4int A = Z + N;
1721 if (A < 1) return infEn;
1722
1723 // Thresholds for p, d, t, 3He in lab frame
1724 else if (A == 1) return 144.6821;
1725 else if (Z == 1 && N == 1) return 2.2263;
1726 else if (Z == 1 && N == 2) return 6.2650;
1727 else if (Z == 2 && N == 1) return 5.4994;
1728
1729 G4double mT = 0.;
1730 if (G4NucleiProperties::IsInStableTable(A,Z
1731 mT = G4NucleiProperties::GetNuclearMass(A
1732 } else {
1733 return infEn;
1734 }
1735
1736 G4double mP = infEn;
1737 if (Z && G4NucleiProperties::IsInStableTabl
1738 mP = G4NucleiProperties::GetNuclearMass(A
1739 }
1740 G4double mN = infEn;
1741 if (N && G4NucleiProperties::IsInStableTabl
1742 mN = G4NucleiProperties::GetNuclearMass(A
1743 }
1744
1745 G4double dP = mP + mProt - mT;
1746 G4double dN = mN + mNeut - mT;
1747 if (dP < dN) dN = dP;
1748 return dN;
1749 }
1750
1751 //
1752 // Linear fit for YN[N] tabulated (from X0 w
1753 //
1754 G4double
1755 G4PhotoNuclearCrossSection::EquLinearFit(G4do
1756 cons
1757 cons
1758 {
1759 if(DX<=0. || N<2)
1760 {
1761 G4cout<<"***G4PhotoNuclearCrossSection:
1762 <<DX<<", N="<<N<<", X0="<<X0<<", X="<<X
1763 return Y[0];
1764 }
1765 G4int N2=N-2;
1766 G4double d=(X-X0)/DX;
1767 G4int j=static_cast<int>(d);
1768 if (j<0) j=0;
1769 else if(j>N2) j=N2;
1770 d-=j; // excess
1771 G4double yi=Y[j];
1772 G4double sigma=yi+(Y[j+1]-yi)*d;
1773 return sigma;
1774 }
1775
1776 // Calculate the functions for the G4Log(A)
1777
1778 G4int
1779 G4PhotoNuclearCrossSection::GetFunctions(G4do
1780 {
1781 if(a<=.9)
1782 {
1783 G4cout << "***G4PhotoNuclearCS::GetFu
1784 << "(?). No CS returned!" << G
1785 return -1;
1786 }
1787 G4int r=0; // Low channe
1788 for(G4int i=0; i<nLA; ++i) if(std::abs(a-LA
1789 {
1790 for(G4int k=0; k<nL; ++k) y[k]=SL[i][
1791 r=1; // Flag
1792 }
1793 G4int h=0;
1794 for(G4int j=0; j<nHA; ++j) if(std::abs(a-HA
1795 {
1796 for(G4int k=0; k<nH; ++k) z[k]=SH[j][
1797 h=1; // Flag
1798 }
1799 if(0 == r) // GDR part is not fil
1800 {
1801 G4int k=0; // !! To be good for dif
1802 for(k=1; k<nLA; ++k) if(a<LA[k]) brea
1803 if(k<1) k=1; // Extr
1804 if(k>=nLA) k=nLA-1; // Extr
1805 G4int k1=k-1;
1806 G4double xi=LA[k1];
1807 G4double b=(a-xi)/(LA[k]-xi);
1808 for(G4int q=0; q<nL; ++q)
1809 {
1810 if(a>1.5)
1811 {
1812 G4double yi=SL[k1][q];
1813 y[q]=yi+(SL[k][q]-yi)*b;
1814 }
1815 else y[q]=0.;
1816 }
1817 r=1;
1818 }
1819 if(0 == h) // High Energy part is not fi
1820 {
1821 G4int k=0;
1822 for(k=1; k<nHA; ++k) if(a<HA[k]) brea
1823 if(k<1) k=1; // Extrapolation
1824 if(k>=nHA) k=nHA-1; // Extrapolation
1825 G4int k1=k-1;
1826 G4double xi=HA[k1];
1827 G4double b=(a-xi)/(HA[k]-xi);
1828 for(G4int q=0; q<nH; ++q)
1829 {
1830 G4double zi=SH[k1][q];
1831 z[q]=zi+(SH[k][q]-zi)*b;
1832 }
1833 h=1;
1834 }
1835 return r*h;
1836 }
1837