00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054 #include "G4HadronCrossSections.hh"
00055 #include "G4ios.hh"
00056 #include "G4SystemOfUnits.hh"
00057 #include "G4HadronicException.hh"
00058 #include "G4Pow.hh"
00059
00060
00061 G4HadronCrossSections* G4HadronCrossSections::theInstance = 0;
00062
00063
00064
00065
00066
00067
00068
00069 G4float G4HadronCrossSections::plab[TSIZE] = {
00070 0.00000E+00 , 0.10000 , 0.15000 , 0.20000 , 0.25000 ,
00071 0.30000 , 0.35000 , 0.40000 , 0.45000 , 0.50000 ,
00072 0.55000 , 0.60000 , 0.65000 , 0.70000 , 0.75000 ,
00073 0.80000 , 0.85000 , 0.90000 , 0.95000 , 1.0000 ,
00074 1.1000 , 1.2000 , 1.3000 , 1.4000 , 1.5000 ,
00075 1.6000 , 1.8000 , 2.0000 , 2.2000 , 2.4000 ,
00076 2.6000 , 2.8000 , 3.0000 , 4.0000 , 5.0000 ,
00077 6.0000 , 8.0000 , 10.000 , 20.000 , 100.00 ,
00078 1000.0
00079 };
00080
00081
00082
00083
00084
00085 G4float G4HadronCrossSections::csel[NPARTS][TSIZE] = {
00086 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00087 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00088 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00089 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00090 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00091 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00092 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00093 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00094 0.00000E+00},
00095
00096 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00097 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00098 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00099 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00100 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00101 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00102 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00103 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00104 0.00000E+00},
00105
00106 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00107 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00108 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00109 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00110 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00111 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00112 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00113 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00114 0.00000E+00},
00115
00116 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00117 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00118 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00119 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00120 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00121 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00122 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00123 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00124 0.00000E+00},
00125
00126 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00127 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00128 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00129 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00130 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00131 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00132 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00133 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00134 0.00000E+00},
00135
00136 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00137 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00138 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00139 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00140 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00141 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00142 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00143 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00144 0.00000E+00},
00145
00146
00147 {0.00000E+00, 6.0000 , 20.000 , 71.000 , 155.00 ,
00148 195.00 , 130.00 , 78.000 , 60.000 , 32.000 ,
00149 23.500 , 18.500 , 15.000 , 12.500 , 10.000 ,
00150 9.1000 , 8.6000 , 8.8000 , 9.5000 , 10.600 ,
00151 13.000 , 15.500 , 17.100 , 17.200 , 16.200 ,
00152 15.000 , 12.300 , 10.200 , 9.0000 , 8.0000 ,
00153 7.3000 , 6.8000 , 6.5000 , 5.8000 , 5.4000 ,
00154 5.2000 , 5.0000 , 4.9000 , 3.8000 , 3.2000 ,
00155 3.5000},
00156
00157 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00158 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00159 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00160 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00161 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00162 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00163 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00164 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00165 0.00000E+00},
00166
00167
00168 {0.00000E+00, 1.0000 , 3.0000 , 8.0000 , 18.000 ,
00169 25.000 , 27.500 , 12.300 , 10.600 , 11.000 ,
00170 12.500 , 14.500 , 17.000 , 19.400 , 19.800 ,
00171 16.800 , 14.000 , 14.800 , 20.000 , 26.100 ,
00172 19.500 , 15.000 , 12.800 , 11.500 , 10.500 ,
00173 9.8000 , 8.8000 , 8.2000 , 7.8000 , 7.5000 ,
00174 7.2000 , 7.0000 , 6.8000 , 6.1000 , 5.7000 ,
00175 5.4000 , 4.9000 , 4.6000 , 4.0000 , 3.3000 ,
00176 3.5000},
00177
00178 {10.000 , 11.200 , 11.300 , 11.400 , 11.500 ,
00179 11.600 , 11.800 , 12.000 , 12.100 , 12.200 ,
00180 12.300 , 12.400 , 12.500 , 12.500 , 12.500 ,
00181 12.400 , 12.300 , 12.200 , 12.000 , 11.800 ,
00182 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 ,
00183 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 ,
00184 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 ,
00185 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 ,
00186 3.0000},
00187
00188 {10.000 , 11.200 , 11.300 , 11.400 , 11.500 ,
00189 11.600 , 11.800 , 12.000 , 12.100 , 12.200 ,
00190 12.300 , 12.400 , 12.500 , 12.500 , 12.500 ,
00191 12.400 , 12.300 , 12.200 , 12.000 , 11.800 ,
00192 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 ,
00193 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 ,
00194 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 ,
00195 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 ,
00196 3.0000},
00197
00198 {160.83 , 82.800 , 58.575 , 43.683 , 34.792 ,
00199 28.650 , 24.367 , 20.917 , 18.192 , 16.300 ,
00200 14.608 , 13.017 , 12.250 , 11.700 , 12.017 ,
00201 14.075 , 15.842 , 16.433 , 16.042 , 15.008 ,
00202 12.575 , 10.708 , 9.2000 , 8.0167 , 7.2833 ,
00203 7.0750 , 6.6333 , 6.1250 , 5.6583 , 5.2750 ,
00204 4.9333 , 4.6250 , 4.4583 , 3.7333 , 3.3833 ,
00205 3.1833 , 2.9833 , 2.7500 , 2.3667 , 2.2000 ,
00206 2.6000},
00207
00208 {300.00 , 140.00 , 97.000 , 70.000 , 55.000 ,
00209 45.000 , 37.000 , 31.000 , 26.000 , 23.000 ,
00210 20.000 , 17.000 , 15.500 , 14.500 , 14.700 ,
00211 18.500 , 22.000 , 23.000 , 22.500 , 20.700 ,
00212 16.500 , 14.000 , 11.500 , 9.6000 , 8.6000 ,
00213 8.5000 , 8.3000 , 7.6000 , 7.0000 , 6.4000 ,
00214 5.9000 , 5.5000 , 5.3000 , 4.4000 , 4.1000 ,
00215 3.9000 , 3.7000 , 3.3000 , 2.6000 , 2.5000 ,
00216 3.0000},
00217
00218
00219 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00220 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00221 20.000 , 20.000 , 20.500 , 21.000 , 22.000 ,
00222 23.000 , 24.000 , 24.000 , 24.400 , 24.500 ,
00223 25.000 , 25.500 , 26.000 , 26.500 , 27.000 ,
00224 27.000 , 26.000 , 23.000 , 21.500 , 20.000 ,
00225 19.000 , 18.000 , 17.000 , 13.000 , 11.500 ,
00226 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 ,
00227 7.5000},
00228
00229 {200.00 , 163.00 , 141.00 , 120.00 , 111.00 ,
00230 99.500 , 92.500 , 86.500 , 82.000 , 78.000 ,
00231 74.000 , 71.000 , 67.500 , 65.000 , 62.500 ,
00232 59.700 , 58.100 , 56.300 , 54.700 , 52.700 ,
00233 50.000 , 48.400 , 47.000 , 46.000 , 45.200 ,
00234 42.800 , 39.200 , 36.300 , 32.800 , 30.400 ,
00235 28.100 , 26.300 , 24.500 , 19.250 , 16.840 ,
00236 14.600 , 12.340 , 11.210 , 8.8500 , 7.5000 ,
00237 7.5000} ,
00238
00239 {4200.0 , 440.00 , 420.00 , 400.00 , 230.00 ,
00240 160.00 , 105.00 , 80.000 , 62.000 , 50.000 ,
00241 45.000 , 41.000 , 38.000 , 36.000 , 35.000 ,
00242 34.000 , 33.000 , 32.000 , 31.500 , 31.000 ,
00243 30.500 , 30.000 , 29.500 , 29.000 , 28.500 ,
00244 28.000 , 26.000 , 23.000 , 21.500 , 20.000 ,
00245 19.000 , 18.000 , 17.000 , 13.000 , 11.500 ,
00246 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 ,
00247 7.5000} ,
00248
00249 {185.88 , 133.23 , 119.37 , 102.86 , 93.102 ,
00250 82.752 , 76.205 , 71.008 , 67.366 , 64.096 ,
00251 60.891 , 58.501 , 55.735 , 53.773 , 51.839 ,
00252 49.671 , 48.485 , 47.045 , 45.803 , 44.306 ,
00253 42.623 , 41.786 , 41.115 , 40.630 , 40.129 ,
00254 38.242 , 35.233 , 32.662 , 29.639 , 27.573 ,
00255 25.536 , 23.948 , 22.356 , 17.723 , 15.614 ,
00256 13.653 , 11.675 , 10.653 , 8.6198 , 7.4464 ,
00257 7.4821} ,
00258
00259 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00260 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00261 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
00262 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
00263 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
00264 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
00265 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
00266 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
00267 6.1145} ,
00268
00269 {157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
00270 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
00271 34.674 , 33.105 , 31.712 , 30.685 , 29.613 ,
00272 28.602 , 28.336 , 28.075 , 27.786 , 27.215 ,
00273 26.380 , 26.146 , 25.108 , 24.783 , 24.360 ,
00274 23.219 , 21.431 , 20.095 , 18.382 , 17.267 ,
00275 16.100 , 15.175 , 14.271 , 11.573 , 10.305 ,
00276 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 ,
00277 6.0774} ,
00278
00279 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00280 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00281 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
00282 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
00283 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
00284 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
00285 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
00286 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
00287 6.1145} ,
00288
00289 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00290 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00291 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00292 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00293 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00294 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00295 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00296 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00297 0.00000E+00},
00298
00299 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00300 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00301 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
00302 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
00303 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
00304 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
00305 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
00306 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
00307 6.1145} ,
00308
00309 {185.88 , 133.23 , 119.37 , 102.86 , 93.102 ,
00310 82.752 , 76.205 , 71.008 , 67.366 , 64.096 ,
00311 60.891 , 58.104 , 55.241 , 53.140 , 50.934 ,
00312 48.660 , 47.566 , 46.585 , 45.581 , 44.003 ,
00313 41.134 , 39.374 , 36.878 , 35.523 , 34.503 ,
00314 32.334 , 29.365 , 27.370 , 24.705 , 22.921 ,
00315 21.229 , 19.879 , 18.559 , 14.625 , 12.758 ,
00316 11.041 , 9.3440 , 8.5484 , 6.7104 , 6.0000 ,
00317 6.1131} ,
00318
00319 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00320 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00321 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00322 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00323 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00324 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00325 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00326 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00327 0.00000E+00},
00328
00329 {157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
00330 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
00331 34.674 , 33.105 , 31.712 , 30.685 , 29.613 ,
00332 28.602 , 28.336 , 28.075 , 27.786 , 27.215 ,
00333 26.380 , 26.146 , 25.108 , 24.783 , 24.360 ,
00334 23.219 , 21.431 , 20.095 , 18.382 , 17.267 ,
00335 16.100 , 15.175 , 14.271 , 11.573 , 10.305 ,
00336 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 ,
00337 6.0774} ,
00338
00339 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00340 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00341 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
00342 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
00343 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
00344 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
00345 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
00346 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
00347 4.7289} ,
00348
00349 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00350 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00351 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
00352 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
00353 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
00354 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
00355 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
00356 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
00357 4.7289} ,
00358
00359 {157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
00360 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
00361 34.674 , 32.708 , 31.218 , 30.052 , 28.707 ,
00362 27.591 , 27.417 , 27.615 , 27.564 , 26.913 ,
00363 24.891 , 23.734 , 20.871 , 19.677 , 18.734 ,
00364 17.311 , 15.563 , 14.803 , 13.448 , 12.615 ,
00365 11.794 , 11.106 , 10.474 , 8.4745 , 7.4498 ,
00366 6.5350 , 5.6835 , 5.3300 , 4.3406 , 4.4464 ,
00367 4.7083} ,
00368
00369 {143.53 , 43.935 , 54.462 , 51.429 , 39.407 ,
00370 32.510 , 27.321 , 24.532 , 23.465 , 22.383 ,
00371 21.566 , 20.209 , 19.453 , 18.825 , 18.046 ,
00372 17.562 , 17.802 , 18.360 , 18.667 , 18.519 ,
00373 17.514 , 17.120 , 14.985 , 14.306 , 13.663 ,
00374 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 ,
00375 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 ,
00376 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 ,
00377 4.6905} ,
00378
00379 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00380 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00381 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00382 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00383 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00384 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00385 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00386 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00387 0.00000E+00},
00388
00389 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00390 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00391 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00392 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00393 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00394 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00395 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00396 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00397 0.00000E+00},
00398
00399 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00400 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00401 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00402 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00403 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00404 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00405 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00406 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00407 0.00000E+00},
00408
00409 {1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
00410 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
00411 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
00412 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
00413 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
00414 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
00415 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
00416 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
00417 4.7289} ,
00418
00419 {143.53 , 43.935 , 54.462 , 51.429 , 39.407 ,
00420 32.510 , 27.321 , 24.532 , 23.465 , 22.383 ,
00421 21.566 , 20.209 , 19.453 , 18.825 , 18.046 ,
00422 17.562 , 17.802 , 18.360 , 18.667 , 18.519 ,
00423 17.514 , 17.120 , 14.985 , 14.306 , 13.663 ,
00424 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 ,
00425 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 ,
00426 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 ,
00427 4.6905} ,
00428
00429 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00430 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00431 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00432 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00433 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00434 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00435 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00436 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00437 0.00000E+00}
00438 };
00439
00440
00441
00442
00443
00444 G4float G4HadronCrossSections::csin[NPARTS][TSIZE] = {
00445
00446 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00447 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00448 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00449 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00450 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00451 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00452 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00453 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00454 0.00000E+00},
00455
00456 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00457 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00458 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00459 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00460 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00461 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00462 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00463 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00464 0.00000E+00},
00465
00466 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00467 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00468 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00469 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00470 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00471 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00472 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00473 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00474 0.00000E+00},
00475
00476 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00477 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00478 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00479 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00480 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00481 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00482 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00483 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00484 0.00000E+00},
00485
00486 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00487 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00488 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00489 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00490 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00491 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00492 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00493 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00494 0.00000E+00},
00495
00496 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00497 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00498 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00499 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00500 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00501 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00502 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00503 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00504 0.00000E+00},
00505
00506
00507 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00508 0.00000E+00, 0.00000E+00, 0.50000 , 1.2000 , 1.7000 ,
00509 2.2500 , 3.0000 , 3.6000 , 4.5000 , 5.4000 ,
00510 6.3000 , 8.6000 , 9.0000 , 10.000 , 11.500 ,
00511 14.000 , 17.000 , 19.500 , 22.000 , 24.000 ,
00512 21.500 , 18.500 , 19.000 , 20.500 , 22.200 ,
00513 23.000 , 23.300 , 23.000 , 21.000 , 20.500 ,
00514 20.200 , 20.100 , 20.000 , 20.000 , 20.000 ,
00515 21.000},
00516
00517
00518 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00519 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00520 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00521 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00522 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00523 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00524 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00525 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00526 0.00000E+00},
00527
00528
00529 {0.00000E+00, 3.0000 , 9.2000 , 20.500 , 36.500 ,
00530 45.000 , 28.000 , 19.500 , 15.500 , 14.200 ,
00531 15.500 , 17.500 , 20.000 , 23.000 , 26.000 ,
00532 20.000 , 23.000 , 26.500 , 32.000 , 35.000 ,
00533 28.500 , 22.000 , 22.500 , 23.500 , 24.000 ,
00534 24.500 , 26.000 , 27.500 , 27.500 , 27.000 ,
00535 26.500 , 25.500 , 25.000 , 23.000 , 22.500 ,
00536 22.200 , 22.000 , 22.000 , 21.200 , 20.700 ,
00537 21.000} ,
00538
00539 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00540 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00541 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00542 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 ,
00543 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 ,
00544 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 ,
00545 12.200 , 12.400 , 12.600 , 13.200 , 13.500 ,
00546 13.700 , 14.000 , 14.200 , 14.500 , 16.400 ,
00547 17.000} ,
00548
00549 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00550 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00551 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00552 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 ,
00553 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 ,
00554 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 ,
00555 12.200 , 12.400 , 12.600 , 13.200 , 13.500 ,
00556 13.700 , 14.000 , 14.200 , 14.500 , 16.400 ,
00557 17.000} ,
00558
00559 {266.67 , 133.33 , 83.333 , 57.083 , 44.500 ,
00560 33.250 , 24.583 , 20.833 , 18.333 , 16.083 ,
00561 15.625 , 15.083 , 14.833 , 15.083 , 15.833 ,
00562 17.042 , 18.958 , 20.758 , 22.533 , 22.825 ,
00563 21.250 , 18.567 , 17.767 , 18.100 , 19.933 ,
00564 20.783 , 21.225 , 21.000 , 20.558 , 20.258 ,
00565 20.017 , 19.767 , 19.600 , 19.183 , 18.850 ,
00566 18.575 , 18.350 , 18.175 , 17.808 , 17.558 ,
00567 19.250} ,
00568
00569 {400.00 , 200.00 , 120.00 , 81.000 , 62.000 ,
00570 47.000 , 35.000 , 28.000 , 24.000 , 21.000 ,
00571 19.500 , 19.000 , 18.800 , 19.000 , 20.000 ,
00572 21.000 , 23.000 , 25.000 , 27.000 , 27.500 ,
00573 25.500 , 22.000 , 20.800 , 21.000 , 23.000 ,
00574 24.000 , 24.000 , 23.800 , 23.000 , 22.500 ,
00575 22.000 , 21.600 , 21.400 , 21.000 , 20.500 ,
00576 20.200 , 19.800 , 19.500 , 18.600 , 17.500 ,
00577 20.000} ,
00578
00579
00580 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00581 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00582 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00583 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 ,
00584 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 ,
00585 22.000 , 23.500 , 24.800 , 25.800 , 26.500 ,
00586 27.000 , 27.500 , 28.000 , 30.000 , 31.000 ,
00587 32.000 , 32.500 , 32.500 , 33.000 , 33.500 ,
00588 33.500} ,
00589
00590 {1500.0 , 1160.0 , 310.00 , 230.00 , 178.00 ,
00591 153.00 , 134.00 , 124.00 , 113.00 , 106.00 ,
00592 101.00 , 96.000 , 92.000 , 89.000 , 87.000 ,
00593 84.000 , 81.000 , 78.500 , 76.500 , 75.000 ,
00594 72.000 , 70.000 , 68.000 , 64.500 , 63.000 ,
00595 62.000 , 61.000 , 59.500 , 58.500 , 56.500 ,
00596 56.500 , 56.000 , 55.500 , 52.000 , 50.000 ,
00597 48.000 , 45.000 , 44.000 , 39.200 , 34.500 ,
00598 34.500} ,
00599
00600 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00601 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00602 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00603 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 ,
00604 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 ,
00605 22.000 , 23.500 , 24.800 , 25.800 , 26.500 ,
00606 27.000 , 27.500 , 28.000 , 30.000 , 31.000 ,
00607 32.000 , 32.500 , 32.500 , 33.000 , 33.500 ,
00608 34.000} ,
00609
00610 {1394.1 , 948.17 , 262.43 , 197.14 , 149.30 ,
00611 127.25 , 110.39 , 101.79 , 92.834 , 87.104 ,
00612 83.109 , 79.099 , 75.965 , 73.627 , 72.161 ,
00613 69.889 , 67.595 , 65.595 , 64.057 , 63.054 ,
00614 61.377 , 60.434 , 59.485 , 56.970 , 55.931 ,
00615 55.398 , 54.827 , 53.538 , 52.861 , 51.247 ,
00616 51.344 , 50.992 , 50.644 , 47.876 , 46.358 ,
00617 44.887 , 42.577 , 41.815 , 38.180 , 34.254 ,
00618 34.418} ,
00619
00620 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00621 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00622 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00623 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
00624 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
00625 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
00626 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
00627 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
00628 27.719} ,
00629
00630 {1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
00631 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
00632 47.326 , 44.762 , 43.222 , 42.015 , 41.221 ,
00633 40.244 , 39.504 , 39.145 , 38.860 , 38.731 ,
00634 37.987 , 37.814 , 36.326 , 34.750 , 33.953 ,
00635 33.635 , 33.349 , 32.938 , 32.785 , 32.092 ,
00636 32.373 , 32.312 , 32.329 , 31.261 , 30.597 ,
00637 30.073 , 29.228 , 29.182 , 27.683 , 27.107 ,
00638 27.956} ,
00639
00640 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00641 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00642 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00643 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
00644 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
00645 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
00646 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
00647 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
00648 27.719} ,
00649
00650 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00651 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00652 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00653 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00654 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00655 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00656 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00657 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00658 0.00000E+00},
00659
00660 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00661 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00662 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00663 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
00664 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
00665 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
00666 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
00667 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
00668 27.719} ,
00669
00670 {1394.1 , 948.17 , 262.43 , 197.14 , 149.30 ,
00671 127.25 , 110.39 , 101.79 , 92.834 , 87.104 ,
00672 83.109 , 78.563 , 75.292 , 72.760 , 70.900 ,
00673 68.467 , 66.314 , 64.955 , 63.746 , 62.623 ,
00674 59.233 , 56.946 , 53.355 , 49.810 , 48.090 ,
00675 46.839 , 45.695 , 44.863 , 44.062 , 42.599 ,
00676 42.684 , 42.328 , 42.041 , 39.508 , 37.880 ,
00677 36.299 , 34.075 , 33.553 , 29.723 , 27.600 ,
00678 28.120} ,
00679
00680 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00681 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00682 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00683 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00684 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00685 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00686 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00687 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00688 0.00000E+00},
00689
00690 {1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
00691 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
00692 47.326 , 44.762 , 43.222 , 42.015 , 41.221 ,
00693 40.244 , 39.504 , 39.145 , 38.860 , 38.731 ,
00694 37.987 , 37.814 , 36.326 , 34.750 , 33.953 ,
00695 33.635 , 33.349 , 32.938 , 32.785 , 32.092 ,
00696 32.373 , 32.312 , 32.329 , 31.261 , 30.597 ,
00697 30.073 , 29.228 , 29.182 , 27.683 , 27.107 ,
00698 27.956} ,
00699
00700 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00701 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00702 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00703 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
00704 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
00705 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
00706 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
00707 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
00708 21.438} ,
00709
00710 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00711 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00712 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00713 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
00714 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
00715 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
00716 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
00717 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
00718 21.438} ,
00719
00720 {1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
00721 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
00722 47.326 , 44.225 , 42.549 , 41.148 , 39.960 ,
00723 38.822 , 38.223 , 38.505 , 38.549 , 38.301 ,
00724 35.843 , 34.326 , 30.196 , 27.590 , 26.112 ,
00725 25.076 , 24.217 , 24.264 , 23.985 , 23.445 ,
00726 23.713 , 23.647 , 23.726 , 22.892 , 22.119 ,
00727 21.485 , 20.726 , 20.921 , 19.226 , 20.454 ,
00728 21.658} ,
00729
00730 {1076.5 , 312.66 , 119.74 , 98.571 , 63.193 ,
00731 49.990 , 39.579 , 35.168 , 32.335 , 30.417 ,
00732 29.434 , 27.325 , 26.514 , 25.775 , 25.120 ,
00733 24.711 , 24.818 , 25.600 , 26.106 , 26.355 ,
00734 25.220 , 24.760 , 21.681 , 20.060 , 19.044 ,
00735 18.474 , 18.044 , 18.301 , 18.347 , 18.192 ,
00736 18.557 , 18.639 , 18.870 , 18.769 , 18.478 ,
00737 18.372 , 18.302 , 18.735 , 18.206 , 20.207 ,
00738 21.576} ,
00739
00740 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00741 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00742 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00743 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00744 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00745 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00746 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00747 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00748 0.00000E+00},
00749
00750 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00751 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00752 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00753 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00754 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00755 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00756 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00757 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00758 0.00000E+00},
00759
00760 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00761 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00762 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00763 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00764 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00765 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00766 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00767 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00768 0.00000E+00},
00769
00770 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00771 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00772 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00773 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
00774 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
00775 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
00776 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
00777 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
00778 21.438} ,
00779
00780 {1076.5 , 312.66 , 119.74 , 98.571 , 63.193 ,
00781 49.990 , 39.579 , 35.168 , 32.335 , 30.417 ,
00782 29.434 , 27.325 , 26.514 , 25.775 , 25.120 ,
00783 24.711 , 24.818 , 25.600 , 26.106 , 26.355 ,
00784 25.220 , 24.760 , 21.681 , 20.060 , 19.044 ,
00785 18.474 , 18.044 , 18.301 , 18.347 , 18.192 ,
00786 18.557 , 18.639 , 18.870 , 18.769 , 18.478 ,
00787 18.372 , 18.302 , 18.735 , 18.206 , 20.207 ,
00788 21.576} ,
00789
00790 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00791 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00792 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00793 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00794 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00795 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00796 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00797 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
00798 0.00000E+00}
00799
00800 };
00801
00802
00803
00804
00805
00806 G4float G4HadronCrossSections::cspiel[3][TSIZE] = {
00807
00808
00809 {0.00000E+00, 350.00 , 580.00 , 600.00 , 550.00 ,
00810 450.00 , 410.00 , 370.00 , 340.00 , 230.00 ,
00811 220.00 , 205.00 , 180.00 , 155.00 , 145.00 ,
00812 140.00 , 160.00 , 195.00 , 235.00 , 250.00 ,
00813 270.00 , 280.00 , 300.00 , 300.00 , 290.00 ,
00814 285.00 , 265.00 , 240.00 , 230.00 , 222.00 ,
00815 204.00 , 196.00 , 190.00 , 170.00 , 170.00 ,
00816 160.00 , 150.00 , 140.00 , 120.00 , 80.000 ,
00817 80.000},
00818
00819
00820 {0.00000E+00, 700.00 , 1000.0 , 1200.0 , 1300.0 ,
00821 1300.0 , 1250.0 , 1250.0 , 1100.0 , 1000.0 ,
00822 940.00 , 740.00 , 700.00 , 670.00 , 660.00 ,
00823 670.00 , 680.00 , 700.00 , 735.00 , 800.00 ,
00824 810.00 , 820.00 , 820.00 , 810.00 , 800.00 ,
00825 800.00 , 700.00 , 600.00 , 500.00 , 470.00 ,
00826 440.00 , 410.00 , 380.00 , 330.00 , 330.00 ,
00827 330.00 , 330.00 , 330.00 , 285.00 , 240.00 ,
00828 240.00},
00829
00830
00831 {0.00000E+00, 1700.0 , 2200.0 , 2200.0 , 1800.0 ,
00832 1300.0 , 1200.0 , 900.00 , 900.00 , 1000.0 ,
00833 1100.0 , 1300.0 , 1400.0 , 1420.0 , 1490.0 ,
00834 1560.0 , 1580.0 , 1690.0 , 1795.0 , 2000.0 ,
00835 2070.0 , 2140.0 , 2050.0 , 2010.0 , 1970.0 ,
00836 1880.0 , 1690.0 , 1500.0 , 1420.0 , 1390.0 ,
00837 1350.0 , 1360.0 , 1370.0 , 1280.0 , 1290.0 ,
00838 1295.0 , 1250.0 , 1200.0 , 1050.0 , 900.00 ,
00839 900.00}
00840 };
00841
00842
00843
00844
00845
00846 G4float G4HadronCrossSections::cspiin[3][TSIZE] = {
00847
00848
00849 {0.00000E+00, 200.00 , 320.00 , 500.00 , 600.00 ,
00850 600.00 , 590.00 , 530.00 , 510.00 , 470.00 ,
00851 430.00 , 425.00 , 420.00 , 425.00 , 425.00 ,
00852 430.00 , 430.00 , 435.00 , 435.00 , 440.00 ,
00853 430.00 , 430.00 , 420.00 , 420.00 , 420.00 ,
00854 415.00 , 415.00 , 410.00 , 410.00 , 408.00 ,
00855 406.00 , 404.00 , 400.00 , 380.00 , 340.00 ,
00856 340.00 , 340.00 , 340.00 , 340.00 , 340.00 ,
00857 340.00} ,
00858
00859
00860 {0.00000E+00, 400.00 , 800.00 , 1000.0 , 1100.0 ,
00861 1200.0 , 1150.0 , 1050.0 , 1000.0 , 900.00 ,
00862 860.00 , 860.00 , 850.00 , 850.00 , 840.00 ,
00863 830.00 , 820.00 , 810.00 , 805.00 , 800.00 ,
00864 800.00 , 800.00 , 800.00 , 800.00 , 800.00 ,
00865 800.00 , 800.00 , 800.00 , 800.00 , 780.00 ,
00866 760.00 , 740.00 , 720.00 , 720.00 , 700.00 ,
00867 690.00 , 680.00 , 670.00 , 665.00 , 660.00 ,
00868 660.00} ,
00869
00870
00871 {0.00000E+00, 1000.0 , 1900.0 , 2600.0 , 2900.0 ,
00872 3000.0 , 2800.0 , 2600.0 , 2500.0 , 2300.0 ,
00873 2200.0 , 2000.0 , 1900.0 , 1880.0 , 1860.0 ,
00874 1840.0 , 1820.0 , 1810.0 , 1805.0 , 1800.0 ,
00875 1780.0 , 1760.0 , 1750.0 , 1740.0 , 1730.0 ,
00876 1720.0 , 1710.0 , 1700.0 , 1680.0 , 1660.0 ,
00877 1650.0 , 1640.0 , 1630.0 , 1620.0 , 1610.0 ,
00878 1605.0 , 1600.0 , 1600.0 , 1550.0 , 1500.0 ,
00879 1500.0}
00880 };
00881
00882
00883
00884
00885
00886 G4float G4HadronCrossSections::cspnel[3][TSIZE] = {
00887
00888
00889 {2100.0 , 1800.0 , 1500.0 , 1050.0 , 900.00 ,
00890 950.00 , 800.00 , 650.00 , 570.00 , 390.00 ,
00891 300.00 , 240.00 , 230.00 , 230.00 , 220.00 ,
00892 220.00 , 225.00 , 225.00 , 240.00 , 240.00 ,
00893 290.00 , 330.00 , 335.00 , 350.00 , 355.00 ,
00894 370.00 , 350.00 , 330.00 , 310.00 , 290.00 ,
00895 270.00 , 265.00 , 260.00 , 230.00 , 210.00 ,
00896 210.00 , 200.00 , 200.00 , 190.00 , 180.00 ,
00897 180.00},
00898
00899
00900 {3800.0 , 2900.0 , 1850.0 , 1550.0 , 1450.0 ,
00901 1520.0 , 1460.0 , 1300.0 , 1140.0 , 880.00 ,
00902 700.00 , 620.00 , 540.00 , 560.00 , 460.00 ,
00903 460.00 , 470.00 , 470.00 , 480.00 , 480.00 ,
00904 580.00 , 600.00 , 610.00 , 620.00 , 620.00 ,
00905 620.00 , 590.00 , 580.00 , 460.00 , 440.00 ,
00906 420.00 , 400.00 , 480.00 , 430.00 , 380.00 ,
00907 380.00 , 380.00 , 380.00 , 380.00 , 380.00 ,
00908 380.00},
00909
00910
00911 {7000.0 , 6000.0 , 4500.0 , 3350.0 , 2700.0 ,
00912 3000.0 , 3550.0 , 3970.0 , 3280.0 , 2490.0 ,
00913 2100.0 , 1510.0 , 1440.0 , 1370.0 , 1370.0 ,
00914 1370.0 , 1400.0 , 1400.0 , 1420.0 , 1420.0 ,
00915 1440.0 , 1460.0 , 1460.0 , 1450.0 , 1450.0 ,
00916 1470.0 , 1400.0 , 1400.0 , 1380.0 , 1370.0 ,
00917 1360.0 , 1350.0 , 1340.0 , 1330.0 , 1320.0 ,
00918 1310.0 , 1305.0 , 1300.0 , 1300.0 , 1300.0 ,
00919 1300.0}
00920 };
00921
00922
00923
00924
00925
00926 G4float G4HadronCrossSections::cspnin[3][TSIZE] = {
00927
00928
00929 {0.00000E+00, 200.00 , 400.00 , 800.00 , 800.00 ,
00930 550.00 , 500.00 , 450.00 , 430.00 , 410.00 ,
00931 400.00 , 390.00 , 380.00 , 370.00 , 370.00 ,
00932 370.00 , 365.00 , 365.00 , 360.00 , 360.00 ,
00933 360.00 , 360.00 , 365.00 , 370.00 , 375.00 ,
00934 380.00 , 400.00 , 410.00 , 420.00 , 430.00 ,
00935 440.00 , 440.00 , 440.00 , 440.00 , 440.00 ,
00936 440.00 , 440.00 , 440.00 , 440.00 , 440.00 ,
00937 440.00} ,
00938
00939
00940 {0.00000E+00, 400.00 , 950.00 , 1050.0 , 1050.0 ,
00941 980.00 , 940.00 , 900.00 , 860.00 , 820.00 ,
00942 800.00 , 780.00 , 760.00 , 740.00 , 740.00 ,
00943 740.00 , 730.00 , 730.00 , 720.00 , 720.00 ,
00944 720.00 , 720.00 , 730.00 , 740.00 , 750.00 ,
00945 760.00 , 800.00 , 820.00 , 820.00 , 820.00 ,
00946 820.00 , 820.00 , 820.00 , 820.00 , 820.00 ,
00947 820.00 , 820.00 , 820.00 , 820.00 , 820.00 ,
00948 820.00} ,
00949
00950
00951 {0.00000E+00, 0.00000E+00, 500.00 , 1450.0 , 1700.0 ,
00952 1800.0 , 1750.0 , 1730.0 , 1720.0 , 1710.0 ,
00953 1700.0 , 1690.0 , 1660.0 , 1630.0 , 1630.0 ,
00954 1630.0 , 1600.0 , 1600.0 , 1580.0 , 1580.0 ,
00955 1580.0 , 1580.0 , 1600.0 , 1630.0 , 1650.0 ,
00956 1670.0 , 1760.0 , 1800.0 , 1800.0 , 1800.0 ,
00957 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 ,
00958 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 ,
00959 1800.0}
00960 };
00961
00962
00963
00964
00965 G4float G4HadronCrossSections::elab[NELAB] = {
00966 0.10000E-03, 0.20000E-03, 0.30000E-03, 0.40000E-03, 0.50000E-03,
00967 0.70000E-03, 0.10000E-02, 0.20000E-02, 0.30000E-02, 0.40000E-02,
00968 0.50000E-02, 0.70000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01,
00969 0.25000E-01, 0.32700E-01
00970 };
00971
00972
00973
00974
00975
00976
00977 G4float G4HadronCrossSections::cnlwat[NCNLW] = {
00978 1.0000 , 16.000 , 27.000 , 56.000 , 59.000 ,
00979 64.000 , 91.000 , 112.00 , 119.00 , 127.00 ,
00980 137.00 , 181.00 , 207.00 , 209.00 , 238.00
00981 };
00982
00983 G4float G4HadronCrossSections::cnlwel[NCNLW][NELAB] = {
00984 {6000.0 , 5500.0 , 5200.0 , 4900.0 , 4800.0 ,
00985 4400.0 , 4000.0 , 2900.0 , 2200.0 , 1800.0 ,
00986 1400.0 , 1100.0 , 900.00 , 700.00 , 600.00 ,
00987 560.00 , 520.00} ,
00988 {5400.0 , 5050.0 , 4800.0 , 4600.0 , 4399.0 ,
00989 4090.0 , 3700.0 , 2600.0 , 1950.0 , 1600.0 ,
00990 1300.0 , 900.00 , 700.00 , 800.00 , 1050.0 ,
00991 1250.0 , 1320.0} ,
00992 {5500.0 , 5150.0 , 4900.0 , 4699.0 , 4490.0 ,
00993 4150.0 , 3750.0 , 2790.0 , 2100.0 , 1650.0 ,
00994 1300.0 , 950.00 , 800.00 , 860.00 , 1000.0 ,
00995 1090.0 , 1080.0} ,
00996 {5499.0 , 4970.0 , 4450.0 , 4080.0 , 3750.0 ,
00997 3380.0 , 2900.0 , 2400.0 , 2380.0 , 2350.0 ,
00998 2300.0 , 2100.0 , 1720.0 , 1370.0 , 1200.0 ,
00999 1060.0 , 870.00} ,
01000 {5399.0 , 4710.0 , 4180.0 , 3760.0 , 3460.0 ,
01001 3150.0 , 2730.0 , 2270.0 , 1850.0 , 1850.0 ,
01002 2130.0 , 2330.0 , 2120.0 , 1640.0 , 1310.0 ,
01003 1100.0 , 1050.0} ,
01004 {5099.0 , 4405.0 , 3825.0 , 3455.0 , 3125.0 ,
01005 2695.0 , 2350.0 , 1850.0 , 1580.0 , 1820.0 ,
01006 2050.0 , 2210.0 , 2000.0 , 1590.0 , 1310.0 ,
01007 1120.0 , 1040.0} ,
01008 {6290.0 , 5960.0 , 5640.0 , 5370.0 , 5150.0 ,
01009 4800.0 , 4250.0 , 3150.0 , 2470.0 , 2100.0 ,
01010 2230.0 , 2420.0 , 2450.0 , 2050.0 , 1760.0 ,
01011 1550.0 , 1330.0} ,
01012 {6885.0 , 6650.0 , 6350.0 , 6150.0 , 6000.0 ,
01013 5700.0 , 5360.0 , 4250.0 , 2800.0 , 1870.0 ,
01014 1810.0 , 1820.0 , 2170.0 , 2450.0 , 2150.0 ,
01015 1700.0 , 1390.0} ,
01016 {6600.0 , 6500.0 , 6400.0 , 6249.0 , 6190.0 ,
01017 5950.0 , 5520.0 , 4250.0 , 2750.0 , 1900.0 ,
01018 1850.0 , 1950.0 , 2340.0 , 2800.0 , 2540.0 ,
01019 2100.0 , 1760.0} ,
01020 {7400.0 , 7200.0 , 6999.0 , 6840.0 , 6655.0 ,
01021 6320.0 , 5820.0 , 4400.0 , 2850.0 , 2000.0 ,
01022 1800.0 , 1800.0 , 2150.0 , 2600.0 , 2350.0 ,
01023 1950.0 , 2100.0} ,
01024 {7900.0 , 7700.0 , 7499.0 , 7390.0 , 7202.0 ,
01025 6810.0 , 6360.0 , 4920.0 , 3450.0 , 2600.0 ,
01026 2200.0 , 1950.0 , 2300.0 , 2800.0 , 2650.0 ,
01027 2250.0 , 2050.0} ,
01028 {7900.0 , 7750.0 , 7699.0 , 7590.0 , 7450.0 ,
01029 7200.0 , 6850.0 , 5650.0 , 4400.0 , 3700.0 ,
01030 3400.0 , 2800.0 , 2700.0 , 3100.0 , 3250.0 ,
01031 3100.0 , 2750.0} ,
01032 {6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 ,
01033 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 ,
01034 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 ,
01035 3500.0 , 2980.0} ,
01036 {6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 ,
01037 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 ,
01038 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 ,
01039 3500.0 , 2980.0} ,
01040 {6600.0 , 6350.0 , 6100.0 , 5899.0 , 5690.0 ,
01041 5300.0 , 4850.0 , 4450.0 , 5650.0 , 5700.0 ,
01042 4950.0 , 3850.0 , 3050.0 , 3050.0 , 3460.0 ,
01043 3650.0 , 3340.0}
01044 };
01045
01046 G4float G4HadronCrossSections::cnlwin[NCNLW][NELAB] = {
01047 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01048 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01049 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01050 0.00000E+00, 0.00000E+00},
01051 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 ,
01052 10.000 , 50.000 , 100.00 , 200.00 , 300.00 ,
01053 400.00 , 600.00 , 700.00 , 750.00 , 700.00 ,
01054 700.00 , 680.00} ,
01055 {0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
01056 50.000 , 100.00 , 260.00 , 450.00 , 600.00 ,
01057 700.00 , 800.00 , 900.00 , 940.00 , 900.00 ,
01058 860.00 , 820.00} ,
01059 {1.0000 , 80.000 , 200.00 , 320.00 , 400.00 ,
01060 520.00 , 700.00 , 1000.0 , 1120.0 , 1200.0 ,
01061 1200.0 , 1200.0 , 1180.0 , 1130.0 , 1100.0 ,
01062 1090.0 , 1080.0} ,
01063 {1.0000 , 90.000 , 220.00 , 340.00 , 420.00 ,
01064 550.00 , 720.00 , 1080.0 , 1300.0 , 1400.0 ,
01065 1420.0 , 1420.0 , 1380.0 , 1260.0 , 1190.0 ,
01066 1150.0 , 1100.0} ,
01067 {1.0000 , 95.000 , 225.00 , 345.00 , 425.00 ,
01068 555.00 , 750.00 , 1150.0 , 1500.0 , 1680.0 ,
01069 1700.0 , 1690.0 , 1550.0 , 1360.0 , 1240.0 ,
01070 1180.0 , 1120.0} ,
01071 {10.000 , 140.00 , 260.00 , 380.00 , 450.00 ,
01072 600.00 , 750.00 , 1200.0 , 1580.0 , 1800.0 ,
01073 1820.0 , 1830.0 , 1800.0 , 1750.0 , 1690.0 ,
01074 1650.0 , 1620.0} ,
01075 {15.000 , 150.00 , 300.00 , 400.00 , 500.00 ,
01076 650.00 , 840.00 , 1500.0 , 2100.0 , 2130.0 ,
01077 2140.0 , 2130.0 , 2080.0 , 2000.0 , 1950.0 ,
01078 1900.0 , 1860.0} ,
01079 {0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
01080 150.00 , 380.00 , 1000.0 , 1650.0 , 2100.0 ,
01081 2100.0 , 2100.0 , 2060.0 , 1950.0 , 1860.0 ,
01082 1800.0 , 1740.0} ,
01083 {0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 45.000 ,
01084 180.00 , 380.00 , 1050.0 , 1900.0 , 2300.0 ,
01085 2300.0 , 2200.0 , 2150.0 , 2000.0 , 1900.0 ,
01086 1800.0 , 1750.0} ,
01087 {0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 48.000 ,
01088 190.00 , 390.00 , 1080.0 , 2000.0 , 2400.0 ,
01089 2400.0 , 2300.0 , 2200.0 , 2100.0 , 1950.0 ,
01090 1850.0 , 1800.0} ,
01091 {0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 50.000 ,
01092 200.00 , 400.00 , 1100.0 , 2100.0 , 2500.0 ,
01093 2500.0 , 2450.0 , 2300.0 , 2100.0 , 2000.0 ,
01094 1900.0 , 1850.0} ,
01095 {0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
01096 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 ,
01097 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 ,
01098 2200.0 , 2120.0} ,
01099 {0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
01100 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 ,
01101 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 ,
01102 2200.0 , 2120.0} ,
01103 {0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
01104 100.00 , 400.00 , 950.00 , 1600.0 , 2200.0 ,
01105 2550.0 , 2750.0 , 2700.0 , 2600.0 , 2540.0 ,
01106 2450.0 , 2360.0}
01107 };
01108
01109 G4float G4HadronCrossSections::cscap[100] = {
01110 6.0000 , 5.7000 , 5.5000 , 5.3000 , 5.2000 ,
01111 5.1000 , 5.0000 , 4.9000 , 4.8000 , 4.8000 ,
01112 4.8000 , 4.8000 , 4.8000 , 4.8000 , 4.8000 ,
01113 4.8000 , 4.9000 , 5.0000 , 5.2000 , 5.5000 ,
01114 6.0000 , 6.7000 , 7.5000 , 8.5000 , 10.000 ,
01115 12.000 , 14.500 , 19.000 , 26.500 , 40.000 ,
01116 75.000 , 120.00 , 180.00 , 260.00 , 360.00 ,
01117 330.00 , 60.000 , 7.0000 , 9.5000 , 20.000 ,
01118 75.000 , 140.00 , 250.00 , 360.00 , 480.00 ,
01119 580.00 , 590.00 , 500.00 , 300.00 , 100.00 ,
01120 200.00 , 300.00 , 400.00 , 470.00 , 500.00 ,
01121 430.00 , 100.00 , 20.000 , 22.000 , 40.000 ,
01122 560.00 , 950.00 , 1000.0 , 1000.0 , 1000.0 ,
01123 990.00 , 920.00 , 860.00 , 790.00 , 740.00 ,
01124 650.00 , 600.00 , 540.00 , 470.00 , 440.00 ,
01125 390.00 , 360.00 , 340.00 , 320.00 , 310.00 ,
01126 280.00 , 2.0000 , 2.5000 , 6.0000 , 13.000 ,
01127 38.000 , 65.000 , 140.00 , 280.00 , 300.00 ,
01128 430.00 , 580.00 , 650.00 , 800.00 , 920.00 ,
01129 1100.0 , 1250.0 , 1400.0 , 1550.0 , 1700.0
01130 };
01131
01132
01133
01134
01135
01136
01137 G4float G4HadronCrossSections::ekfiss[NFISS] = {
01138 0.10000E-03, 0.20000E-03, 0.30000E-03, 0.50000E-03, 0.70000E-03,
01139 0.10000E-02, 0.15000E-02, 0.20000E-02, 0.30000E-02, 0.50000E-02,
01140 0.70000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01, 0.50000E-01,
01141 0.10000 , 0.20000 , 0.30000 , 0.40000 , 0.50000 ,
01142 1000.0
01143 };
01144
01145 G4float G4HadronCrossSections::csfiss[4][NFISS] = {
01146 {2600.0 , 2300.0 , 2300.0 , 2100.0 , 2000.0 ,
01147 1950.0 , 1930.0 , 1900.0 , 1800.0 , 1600.0 ,
01148 2100.0 , 2300.0 , 0.00000E+00, 0.00000E+00, 0.00000E+00,
01149 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01150 0.00000E+00},
01151 {1850.0 , 1400.0 , 1300.0 , 1150.0 , 1100.0 ,
01152 1200.0 , 1250.0 , 1300.0 , 1250.0 , 1150.0 ,
01153 1600.0 , 1900.0 , 0.00000E+00, 0.00000E+00, 0.00000E+00,
01154 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01155 0.00000E+00},
01156 {1700.0 , 1650.0 , 1650.0 , 1700.0 , 1700.0 ,
01157 1800.0 , 1900.0 , 2000.0 , 1950.0 , 1800.0 ,
01158 2150.0 , 2450.0 , 0.00000E+00, 0.00000E+00, 0.00000E+00,
01159 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01160 0.00000E+00},
01161 {0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
01162 0.00000E+00, 250.00 , 550.00 , 550.00 , 550.00 ,
01163 550.00 , 550.00 , 1000.0 , 1400.0 , 1600.0 ,
01164 1500.0 , 1400.0 , 1300.0 , 1200.0 , 1100.0 ,
01165 1000.0}
01166 };
01167
01168 G4float G4HadronCrossSections::alpha[NPARTS] = {
01169 0.7,0.7,0.7,0.7,0.7,0.7,
01170 0.75,0.75,0.75,
01171 0.76,0.76,0.76,0.76,
01172 0.685,0.63,0.685,0.63,0.685,0.63,
01173 0.685,0.685,0.685,0.63,0.63,0.63,0.685,0.685,0.63,0.63,
01174 0.7,0.7,0.7,0.685,0.63,0.7
01175 };
01176
01177 G4float G4HadronCrossSections::alphac[TSIZE] = {
01178 1.2,1.2,1.2,1.15,0.90,0.91,0.98,1.06,1.10,1.11,
01179 1.10,1.08,1.05,1.01,0.985,0.962,0.945,0.932,
01180 0.925,0.920,0.920,0.921,0.922,0.923,0.928,0.931,
01181 0.940,0.945,0.950,0.955,0.958,0.962,0.965,0.976,
01182 0.982,0.988,0.992,1.010,1.020,1.030,1.040
01183 };
01184
01185 G4float G4HadronCrossSections::partel[NPARTS] = {
01186 0.,0.,0.,0.,0.,0.,
01187 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,
01188 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,
01189 1.,1.,1.,1.,1.,1.,1.,1.,1.
01190 };
01191 G4float G4HadronCrossSections::partin[NPARTS] = {
01192 0.,0.,0.,0.,0.,0.,
01193 1.00,0.00,1.05,1.20,1.35,1.30,1.20,1.00,1.30,
01194 1.00,1.30,1.00,1.30,1.00,1.00,1.00,1.30,1.30,1.30,
01195 1.00,1.00,1.30,1.30,1.00,1.,1.,1.,1.3,1.
01196 };
01197
01198
01199 G4int G4HadronCrossSections::icorr[NPARTS] = {
01200 1,1,1,1,1,1,1,1,1,1,1,1,1,1,
01201 0,1,0,1,0,1,1,1,0,0,0,1,1,0,0,1,1,1,1,0,0
01202 };
01203
01204
01205 G4int G4HadronCrossSections::intrc[NPARTS] = {
01206 0,0,0,0,0,0,
01207 1,0,1,1,1,1,1,1,1,1,1,1,1,1,0,2,2,0,
01208 1,1,1,1,1,1,1,1,1,1,0
01209 };
01210
01211 G4float G4HadronCrossSections::csa[4] = {1., 27.00, 63.54, 207.19};
01212
01213 G4int G4HadronCrossSections::ipart2[7] = {9, 8, 7, 11, 10, 13, 12};
01214
01215 G4bool G4HadronCrossSections::correctInelasticNearZero = 0;
01216
01217 G4HadronCrossSections::G4HadronCrossSections()
01218 : prevParticleDefinition(0), prevElement(0), prevZZ(0), prevAA(0),
01219 prevKineticEnergy(DBL_MAX), lastEkx(0.), lastEkxPower(0.), verboseLevel(0)
01220 {
01221 siginelastic = sigelastic = 0.0;
01222 g4pow = G4Pow::GetInstance();
01223 }
01224
01225 G4HadronCrossSections::~G4HadronCrossSections()
01226 {}
01227
01228 G4HadronCrossSections* G4HadronCrossSections::Instance()
01229 {
01230 if (!theInstance) {
01231 static G4HadronCrossSections xsection;
01232 theInstance = &xsection;
01233 }
01234 return theInstance;
01235 }
01236
01237 G4bool G4HadronCrossSections::IsApplicable(const G4DynamicParticle* aParticle)
01238 {
01239 return (GetParticleCode(aParticle) > 0);
01240 }
01241
01242 G4double
01243 G4HadronCrossSections::GetInelasticCrossSection(const G4DynamicParticle* particle,
01244 G4int ZZ, G4int AA)
01245 {
01246 if (particle->GetDefinition() != prevParticleDefinition ||
01247 particle->GetKineticEnergy() != prevKineticEnergy ||
01248 ZZ != prevZZ || AA != prevAA)
01249 { CalcScatteringCrossSections(particle, ZZ, AA); }
01250 if (verboseLevel > 1) {
01251 G4cout << " x-section inelastic(mb)= "
01252 << siginelastic/millibarn << G4endl;
01253 }
01254
01255 return siginelastic;
01256 }
01257
01258
01259 G4double
01260 G4HadronCrossSections::GetElasticCrossSection(const G4DynamicParticle* particle,
01261 G4int ZZ, G4int AA)
01262 {
01263 if (particle->GetDefinition() != prevParticleDefinition ||
01264 particle->GetKineticEnergy() != prevKineticEnergy ||
01265 ZZ != prevZZ || AA != prevAA)
01266 { CalcScatteringCrossSections(particle, ZZ, AA); }
01267
01268 if (verboseLevel > 1) {
01269 G4cout << " x-section elastic(mb)= "
01270 << sigelastic/millibarn << G4endl;
01271 }
01272 return sigelastic;
01273 }
01274
01275
01276
01277
01278
01279
01280
01281 void
01282 G4HadronCrossSections::CalcScatteringCrossSections(
01283 const G4DynamicParticle* aParticle,
01284 G4int ZZ, G4int AA)
01285 {
01286 G4double sigel, sigin, sigtot;
01287 G4double xsecel, xsecin=0;
01288 xsecel = 0;
01289 G4double dx, dy, rc, rca, b;
01290 G4double crel, crin;
01291 G4double xspiel, xspiin;
01292
01293 G4int ipart = GetParticleCode(aParticle);
01294 prevAA = AA;
01295 prevZZ = ZZ;
01296 prevParticleDefinition = aParticle->GetDefinition();
01297 prevKineticEnergy = aParticle->GetKineticEnergy();
01298
01299 if (verboseLevel > 1) {
01300 G4cout << "G4HadronCrossSections::ScatteringCrossSections: "
01301 << aParticle->GetDefinition()->GetParticleName()
01302 << " E(MeV)= " << aParticle->GetKineticEnergy()/MeV
01303 << " Z= " << ZZ << " A= " << AA << G4endl;
01304 }
01305
01306
01307
01308 if (ipart >= 30 && ipart <= 32) {
01309
01310 G4double apart=0;
01311 if (ipart == 30) apart = g4pow->Z13(2);
01312 else if (ipart == 31) apart = g4pow->Z13(3);
01313 else if (ipart == 32) apart = g4pow->Z13(4);
01314
01315 G4double term = apart + g4pow->Z13(AA);
01316 sigin = 49.*term*term;
01317
01318 siginelastic = sigin*millibarn;
01319 if(prevKineticEnergy < 6*MeV) siginelastic = 0;
01320 sigelastic = 0.;
01321 return;
01322 }
01323
01324 G4double ek = prevKineticEnergy/GeV;
01325
01326
01327
01328 if (ipart == 16 && ek <= 0.0327) {
01329
01330
01331
01332
01333
01334
01335
01336
01337
01338
01339
01340
01341 G4int je1 = 0;
01342 G4int je2 = NELAB - 1;
01343 do {
01344 G4int midBin = (je1 + je2)/2;
01345 if (ek < elab[midBin])
01346 je2 = midBin;
01347 else
01348 je1 = midBin;
01349 } while (je2 - je1 > 1);
01350
01351 G4double delab = elab[je2] - elab[je1];
01352
01353
01354
01355
01356
01357
01358
01359
01360
01361
01362
01363
01364 G4int ja1 = 0;
01365 G4int ja2 = NCNLW - 1;
01366 do {
01367 G4int midBin = (ja1 + ja2)/2;
01368 if (AA < cnlwat[midBin])
01369 ja2 = midBin;
01370 else
01371 ja1 = midBin;
01372 } while (ja2 - ja1 > 1);
01373
01374 G4double dnlwat = cnlwat[ja2] - cnlwat[ja1];
01375
01376
01377
01378 dy = cnlwel[ja1][je2] - cnlwel[ja1][je1];
01379 G4double rce = dy/delab;
01380
01381 dy = cnlwel[ja2][je1] - cnlwel[ja1][je1];
01382 rca = dy/dnlwat;
01383 b = cnlwel[ja1][je1] - rce*elab[je1] - rca*cnlwat[ja1];
01384 sigelastic = rce*ek + rca*AA + b;
01385
01386
01387 dy = cnlwin[ja1][je2] - cnlwin[ja1][je1];
01388 rce = dy/delab;
01389
01390 dy = cnlwin[ja2][je1] - cnlwin[ja1][je1];
01391 rca = dy/dnlwat;
01392 b = cnlwin[ja1][je1] - rce*elab[je1] - rca*cnlwat[ja1];
01393 siginelastic = rce*ek + rca*AA + b;
01394
01395 sigelastic = sigelastic*millibarn;
01396 siginelastic = siginelastic*millibarn;
01397 return;
01398 }
01399
01400
01401
01402
01403 G4double p = aParticle->GetTotalMomentum()/GeV;
01404
01405
01406
01407
01408
01409
01410
01411
01412
01413
01414 G4int je1 = 0;
01415 G4int je2 = TSIZE - 1;
01416 do {
01417 G4int midBin = (je1 + je2)/2;
01418 if (p < plab[midBin])
01419 je2 = midBin;
01420 else
01421 je1 = midBin;
01422 } while (je2 - je1 > 1);
01423
01424 G4int ipart1 = ipart - 1;
01425
01426
01427 dx = plab[je2] - plab[je1];
01428
01429 dy = csel[ipart1][je2] - csel[ipart1][je1];
01430 rc = dy/dx;
01431 b = csel[ipart1][je1] - rc*plab[je1];
01432 sigel = rc*p + b;
01433
01434 dy = csin[ipart1][je2] - csin[ipart1][je1];
01435 rc = dy/dx;
01436 b = csin[ipart1][je1] - rc*plab[je1];
01437 sigin = rc*p + b;
01438 if (verboseLevel > 1) {
01439 G4cout << "sigel " << sigel << G4endl;
01440 G4cout << "sigin " << sigin << G4endl;
01441 }
01442 G4double alph = alpha[ipart1];
01443 if (ipart < 14) {
01444 dy = alphac[je2] - alphac[je1];
01445 rc = dy/dx;
01446 b = alphac[je1] - rc*plab[je1];
01447 G4double corfac = rc*p + b;
01448 alph = alph*corfac;
01449 G4int ipart3 = ipart2[ipart - 7];
01450 ipart3 = ipart3 - 1;
01451
01452 dy = csel[ipart3][je2] - csel[ipart3][je1];
01453 rc = dy/dx;
01454 b = csel[ipart3][je1] - rc*plab[je1];
01455 xsecel = rc*p + b;
01456
01457 dy = csin[ipart3][je2] - csin[ipart3][je1];
01458 rc = dy/dx;
01459 b = csin[ipart3][je1] - rc*plab[je1];
01460 xsecin = rc*p + b;
01461 }
01462
01463
01464
01465 if (AA >= 2) {
01466
01467 crel = 1.;
01468 crin = 1.;
01469
01470 G4int i = 2;
01471 if (AA < 50) i = 1;
01472 else if (AA > 100) i = 3;
01473
01474
01475 if (ipart == 14 || ipart == 16) {
01476 dy = cspnel[i - 1][je2] - cspnel[i - 1][je1];
01477 rc = dy/dx;
01478 b = cspnel[i - 1][je1] - rc*plab[je1];
01479 xsecel = rc*p + b;
01480 dy = cspnin[i - 1][je2] - cspnin[i - 1][je1];
01481 rc = dy/dx;
01482 b = cspnin[i - 1][je1] - rc*plab[je1];
01483 xsecin = rc*p + b;
01484
01485
01486
01487
01488
01489 if (correctInelasticNearZero && je1 == 0 && i <= 3) {
01490 G4double m0 = aParticle->GetMass()/GeV;
01491 G4double T = std::sqrt(m0*m0 + p*p) - m0;
01492 dx = std::sqrt(m0*m0 + plab[1]*plab[1]) - m0;
01493 rc = dy/dx;
01494 xsecin = rc*T + b;
01495 }
01496
01497 if (sigel >= 0.001)
01498 crel = xsecel/(0.36*sigel*std::pow(G4double(csa[i]), 1.17));
01499 sigtot = sigel + sigin;
01500 if (sigtot >= 0.001)
01501 crin = xsecin/(sigtot*std::pow(G4double(csa[i]), alph));
01502 }
01503
01504 else if (ipart < 15) {
01505
01506
01507 G4double wgch = 0.5;
01508 if (AA < 20) { wgch = 0.5 + 0.5*std::exp(-G4double(AA - 1)); }
01509 sigel = wgch*sigel + (1. - wgch)*xsecel;
01510 sigin = wgch*sigin + (1. - wgch)*xsecin;
01511
01512
01513 if (ipart < 10) {
01514 dy = cspiel[i - 1][je2] - cspiel[i - 1][je1];
01515 rc = dy/dx;
01516 b = cspiel[i - 1][je1] - rc*plab[je1];
01517 xspiel = rc*p + b;
01518 dy = cspiin[i - 1][je2] - cspiin[i - 1][je1];
01519 rc = dy/dx;
01520 b = cspiin[i - 1][je1] - rc*plab[je1];
01521 xspiin = rc*p + b;
01522 if (verboseLevel > 1) {
01523 G4cout << "xspiel " << xspiel << G4endl;
01524 G4cout << "xspiin " << xspiin << G4endl;
01525 }
01526 if (sigel >= 0.001)
01527 crel = xspiel/(0.36*sigel*std::pow(G4double(csa[i]),1.17));
01528 sigtot = sigel + sigin;
01529 if (sigtot >= 0.001)
01530 crin = xspiin/(sigtot*std::pow(G4double(csa[i]), alph));
01531 }
01532 }
01533
01534
01535 sigin = crin*(sigin + sigel)*g4pow->powZ(AA, alph);
01536 sigel = crel*0.36*sigel*g4pow->powZ(AA, 1.17);
01537 sigel = sigel*partel[ipart1];
01538 sigin = sigin*partin[ipart1];
01539 }
01540
01541
01542 G4double corh = 1.;
01543 if (p > 100.) corh = 0.1085736156*std::log(p) + 0.5;
01544
01545 sigel = corh*sigel;
01546 sigin = corh*sigin;
01547
01548 sigelastic = sigel*millibarn;
01549 siginelastic = sigin*millibarn;
01550
01551 return;
01552 }
01553
01554
01555 G4double
01556 G4HadronCrossSections::GetCaptureCrossSection(const G4DynamicParticle* aParticle,
01557 G4int ZZ)
01558 {
01559 if (GetParticleCode(aParticle) != 16) return 0.;
01560 G4double ek = aParticle->GetKineticEnergy()/GeV;
01561 if (ek > 0.0327) { return 0.; }
01562
01563 G4double ekx = std::max(ek, 1.e-9);
01564 if (ekx != lastEkx) {
01565 lastEkx = ekx;
01566 lastEkxPower = std::pow(ekx*1.e6, 0.577);
01567 }
01568
01569 G4int izno = ZZ;
01570 if (izno > 100) izno = 100;
01571 izno = izno - 1;
01572 G4double sigcap = 11.12*cscap[izno]/lastEkxPower;
01573
01574
01575 sigcap = sigcap*millibarn;
01576 return sigcap;
01577 }
01578
01579
01580 G4double
01581 G4HadronCrossSections::GetFissionCrossSection(const G4DynamicParticle* aParticle,
01582 G4int ZZ, G4int AA)
01583 {
01584 if (AA < 230) return 0;
01585
01586 G4double ek = aParticle->GetKineticEnergy()/GeV;
01587
01588 G4int ie1 = 0;
01589 G4int ie2 = NFISS - 1;
01590 do {
01591 G4int midBin = (ie1 + ie2)/2;
01592 if (ek < ekfiss[midBin])
01593 ie2 = midBin;
01594 else
01595 ie1 = midBin;
01596 } while (ie2 - ie1 > 1);
01597 G4int i = ie2;
01598 if (ek < ekfiss[0]) i = 0;
01599
01600 G4int j = 4;
01601 if (ek <= 0.01) {
01602 if (ZZ == 92 && AA == 233) j = 1;
01603 else if (ZZ == 92 && AA == 235) j = 2;
01604 else if (ZZ == 94 && AA == 239) j = 3;
01605 }
01606
01607 G4double z43ba;
01608 if (j == 4) {
01609 z43ba = ZZ*g4pow->Z13(ZZ)/G4double(AA);
01610 z43ba = std::max(-67. + 38.7*z43ba, 0.);
01611 } else {
01612 z43ba = 1.;
01613 }
01614 j = j - 1;
01615
01616 G4double sigfiss = csfiss[j][i]*z43ba;
01617 sigfiss = sigfiss*millibarn;
01618 return sigfiss;
01619 }
01620
01621
01622 G4int
01623 G4HadronCrossSections::GetParticleCode(const G4DynamicParticle* aParticle)
01624 {
01625
01626
01627
01628 G4int ipart;
01629
01630 switch( aParticle->GetPDGcode()) {
01631 case 111:
01632 ipart = 8;
01633 break;
01634 case 211:
01635 ipart = 7;
01636 break;
01637 case -211:
01638 ipart = 9;
01639 break;
01640 case 2112:
01641 ipart = 16;
01642 break;
01643 case 2212:
01644 ipart = 14;
01645 break;
01646 case 321:
01647 ipart = 10;
01648 break;
01649 case -321:
01650 ipart = 13;
01651 break;
01652 case 130:
01653 ipart = 12;
01654 break;
01655 case 310:
01656 ipart = 11;
01657 break;
01658 case 1000010020:
01659 ipart = 30;
01660 break;
01661 case 1000010030:
01662 ipart = 31;
01663 break;
01664 case 1000020040:
01665 ipart = 32;
01666 break;
01667 case 3122:
01668 ipart = 18;
01669 break;
01670 case -2112:
01671 ipart = 17;
01672 break;
01673 case -2212:
01674 ipart = 15;
01675 break;
01676 case -3122:
01677 ipart = 19;
01678 break;
01679 case 3222:
01680 ipart = 20;
01681 break;
01682 case 3212:
01683 ipart = 21;
01684 break;
01685 case 3112:
01686 ipart = 22;
01687 break;
01688 case 3322:
01689 ipart = 26;
01690 break;
01691 case 3312:
01692 ipart = 27;
01693 break;
01694 case 3334:
01695 ipart = 33;
01696 break;
01697 case -3222:
01698 ipart = 23;
01699 break;
01700 case -3212:
01701 ipart = 24;
01702 break;
01703 case -3112:
01704 ipart = 25;
01705 break;
01706 case -3322:
01707 ipart = 28;
01708 break;
01709 case -3312:
01710 ipart = 29;
01711 break;
01712 case -3334:
01713 ipart = 34;
01714 break;
01715 default:
01716 throw G4HadronicException(__FILE__, __LINE__,
01717 "G4HadronCrossSections::GetParticleCode: unsupported particle");
01718 return 0;
01719 }
01720
01721 return ipart;
01722 }