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 #include "G4RPGPionInelastic.hh"
00030 #include "G4SystemOfUnits.hh"
00031 #include "Randomize.hh"
00032
00033 G4RPGPionInelastic::G4RPGPionInelastic(const G4String& modelName)
00034 :G4RPGInelastic(modelName)
00035 {
00036 SetMinEnergy( 0.0 );
00037 SetMaxEnergy( 30.*GeV );
00038
00039
00040
00041
00042
00043 G4int i, k, j;
00044 G4int start, stop;
00045
00046 for (j = 0; j < 8; j++) {
00047 start = pipPindex[j][0];
00048 stop = pipPindex[j][1] + 1;
00049 for (k = 0; k < 30; k++) {
00050 t32_dSigma_dMult[j][k] = 0.0;
00051 for (i = start; i < stop; i++) t32_dSigma_dMult[j][k] += pipPCrossSections[i][k];
00052 }
00053
00054 start = pimPindex[j][0];
00055 stop = pimPindex[j][1] + 1;
00056 for (k = 0; k < 30; k++) {
00057 t12_dSigma_dMult[j][k] = 0.0;
00058 for (i = start; i < stop; i++) t12_dSigma_dMult[j][k] += pimPCrossSections[i][k];
00059 }
00060 }
00061
00062
00063 for (k = 0; k < 30; k++) {
00064 pipPtot[k] = 0.0;
00065 pimPtot[k] = 0.0;
00066 for (j = 0; j < 8; j++) {
00067 pipPtot[k] += t32_dSigma_dMult[j][k];
00068 pimPtot[k] += t12_dSigma_dMult[j][k];
00069 }
00070 }
00071
00072
00073 }
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099 G4int G4RPGPionInelastic::GetMultiplicityT12(G4double KE) const
00100 {
00101 G4double multint(0.);
00102 std::vector<G4double> sigma;
00103
00104 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00105 G4int k = epair.first;
00106 G4double fraction = epair.second;
00107
00108 for(G4int j = 0; j < 8; j++) {
00109 multint = t12_dSigma_dMult[j][k]
00110 + fraction*(t12_dSigma_dMult[j][k+1] - t12_dSigma_dMult[j][k]);
00111 sigma.push_back(multint);
00112 }
00113
00114 return sampleFlat(sigma) + 2;
00115 }
00116
00117
00118 G4int G4RPGPionInelastic::GetMultiplicityT32(G4double KE) const
00119 {
00120 G4double multint(0.);
00121 std::vector<G4double> sigma;
00122
00123 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00124 G4int k = epair.first;
00125 G4double fraction = epair.second;
00126
00127 for (G4int j = 0; j < 8; j++) {
00128 multint = t32_dSigma_dMult[j][k]
00129 + fraction*(t32_dSigma_dMult[j][k+1] - t32_dSigma_dMult[j][k]);
00130 sigma.push_back(multint);
00131 }
00132
00133 return sampleFlat(sigma) + 2;
00134 }
00135
00136
00137 std::vector<G4int>
00138 G4RPGPionInelastic::GetFSPartTypesForT32(G4int mult, G4double KE, G4int tzindex) const
00139 {
00140 G4int i;
00141 G4double sigint(0.);
00142 std::vector<G4double> sigma;
00143
00144 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00145 G4int k = epair.first;
00146 G4double fraction = epair.second;
00147
00148 G4int start = pipPindex[mult-2][0];
00149 G4int stop = pipPindex[mult-2][1];
00150
00151 for(i = start; i < stop; i++) {
00152 sigint = pipPCrossSections[i][k]
00153 + fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[i][k]);
00154 sigma.push_back(sigint);
00155 }
00156
00157 G4int channel = sampleFlat(sigma);
00158
00159 std::vector<G4int> kinds;
00160
00161 if (mult == 2) {
00162 for(i = 0; i < mult; i++) kinds.push_back(T32_2bfs[tzindex][channel][i]);
00163 } else if (mult == 3) {
00164 for(i = 0; i < mult; i++) kinds.push_back(T32_3bfs[tzindex][channel][i]);
00165 } else if (mult == 4) {
00166 for(i = 0; i < mult; i++) kinds.push_back(T32_4bfs[tzindex][channel][i]);
00167 } else if (mult == 5) {
00168 for(i = 0; i < mult; i++) kinds.push_back(T32_5bfs[tzindex][channel][i]);
00169 } else if (mult == 6) {
00170 for(i = 0; i < mult; i++) kinds.push_back(T32_6bfs[tzindex][channel][i]);
00171 } else if (mult == 7) {
00172 for(i = 0; i < mult; i++) kinds.push_back(T32_7bfs[tzindex][channel][i]);
00173 } else if (mult == 8) {
00174 for(i = 0; i < mult; i++) kinds.push_back(T32_8bfs[tzindex][channel][i]);
00175 } else if (mult == 9) {
00176 for(i = 0; i < mult; i++) kinds.push_back(T32_9bfs[tzindex][channel][i]);
00177 } else {
00178 G4cout << " Illegal multiplicity " << G4endl;
00179 }
00180
00181 return kinds;
00182 }
00183
00184
00185 std::vector<G4int>
00186 G4RPGPionInelastic::GetFSPartTypesForT12(G4int mult, G4double KE, G4int tzindex) const
00187 {
00188 G4int i;
00189 G4double sigint(0.);
00190 std::vector<G4double> sigma;
00191
00192 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00193 G4int k = epair.first;
00194 G4double fraction = epair.second;
00195
00196 G4int start = pimPindex[mult-2][0];
00197 G4int stop = pimPindex[mult-2][1];
00198
00199 for(i = start; i < stop; i++) {
00200 sigint = pimPCrossSections[i][k]
00201 + fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
00202 sigma.push_back(sigint);
00203 }
00204
00205 G4int channel = sampleFlat(sigma);
00206
00207 std::vector<G4int> kinds;
00208
00209 if (mult == 2) {
00210 for(i = 0; i < mult; i++) kinds.push_back(T12_2bfs[tzindex][channel][i]);
00211 } else if (mult == 3) {
00212 for(i = 0; i < mult; i++) kinds.push_back(T12_3bfs[tzindex][channel][i]);
00213 } else if (mult == 4) {
00214 for(i = 0; i < mult; i++) kinds.push_back(T12_4bfs[tzindex][channel][i]);
00215 } else if (mult == 5) {
00216 for(i = 0; i < mult; i++) kinds.push_back(T12_5bfs[tzindex][channel][i]);
00217 } else if (mult == 6) {
00218 for(i = 0; i < mult; i++) kinds.push_back(T12_6bfs[tzindex][channel][i]);
00219 } else if (mult == 7) {
00220 for(i = 0; i < mult; i++) kinds.push_back(T12_7bfs[tzindex][channel][i]);
00221 } else if (mult == 8) {
00222 for(i = 0; i < mult; i++) kinds.push_back(T12_8bfs[tzindex][channel][i]);
00223 } else if (mult == 9) {
00224 for(i = 0; i < mult; i++) kinds.push_back(T12_9bfs[tzindex][channel][i]);
00225 } else {
00226 G4cout << " Illegal multiplicity " << G4endl;
00227 }
00228
00229 return kinds;
00230 }
00231
00232
00234
00235
00236
00237
00239
00240
00241 G4double G4RPGPionInelastic::pipPtot[30];
00242
00243
00244 G4double G4RPGPionInelastic::t32_dSigma_dMult[8][30];
00245
00246 const G4int G4RPGPionInelastic::pipPindex[8][2] =
00247 {{0, 1}, {2, 8}, {9,23}, {24,47}, {48,52}, {53,58}, {59,65}, {66,73}};
00248
00249
00250
00251
00252 const G4int G4RPGPionInelastic::T32_2bfs[2][2][2] =
00253 {{{pro,pip}, {sp,kp}},
00254
00255 {{neu,pim}, {sm,k0}}};
00256
00257 const G4int G4RPGPionInelastic::T32_3bfs[2][7][3] =
00258 {{{pro,pip,pi0}, {neu,pip,pip}, {sp,kp,pi0}, {sp,k0,pip},
00259 {s0,kp,pip}, {lam,kp,pip}, {pro,kp,k0b}},
00260
00261 {{neu,pim,pi0}, {pro,pim,pim}, {sm,k0,pi0}, {sm,kp,pim},
00262 {s0,k0,pim}, {lam,k0,pim}, {neu,k0,km}}};
00263
00264 const G4int G4RPGPionInelastic::T32_4bfs[2][15][4] =
00265 {{{pro,pip,pip,pim},{pro,pip,pi0,pi0},{neu,pip,pip,pi0},
00266 {sp,kp,pip,pim}, {sp,kp,pi0,pi0}, {sp,k0,pip,pi0},
00267 {s0,k0,pip,pip}, {s0,kp,pip,pi0}, {lam,kp,pip,pi0},
00268 {lam,k0,pip,pip}, {sm,kp,pip,pip}, {pro,pip,kp,km},
00269 {pro,pip,k0,k0b}, {pro,pi0,kp,k0b}, {neu,pip,kp,k0b}},
00270
00271 {{neu,pip,pim,pim},{neu,pim,pi0,pi0},{pro,pim,pim,pi0},
00272 {sm,k0,pip,pim}, {sm,k0,pi0,pi0}, {sm,kp,pim,pi0},
00273 {s0,kp,pim,pim}, {s0,k0,pim,pi0}, {lam,k0,pim,pi0},
00274 {lam,kp,pim,pim}, {sp,k0,pim,pim}, {neu,pim,k0,k0b},
00275 {neu,pim,kp,km}, {neu,pi0,k0,km}, {pro,pim,k0,km}}};
00276
00277 const G4int G4RPGPionInelastic::T32_5bfs[2][24][5] =
00278 {{{pro,pip,pip,pim,pi0}, {pro,pip,pi0,pi0,pi0}, {neu,pip,pip,pip,pim},
00279 {neu,pip,pip,pi0,pi0}, {sp,kp,pip,pim,pi0}, {sp,kp,pi0,pi0,pi0},
00280 {sp,k0,pip,pip,pim}, {sp,k0,pip,pi0,pi0}, {lam,k0,pip,pip,pi0},
00281 {lam,kp,pip,pip,pim}, {lam,kp,pip,pi0,pi0}, {s0,kp,pip,pip,pim},
00282 {s0,kp,pip,pi0,pi0}, {s0,k0,pip,pip,pi0}, {sm,kp,pip,pip,pi0},
00283 {sm,k0,pip,pip,pip}, {pro,pip,pim,kp,k0b}, {pro,pip,pip,k0,km},
00284 {pro,pip,pi0,kp,km}, {pro,pip,pi0,k0,k0b}, {pro,pi0,pi0,kp,k0b},
00285 {neu,pip,pip,kp,km}, {neu,pip,pip,k0,k0b}, {neu,pip,pi0,kp,k0b}},
00286
00287 {{neu,pip,pim,pim,pi0}, {neu,pim,pi0,pi0,pi0}, {pro,pip,pim,pim,pim},
00288 {pro,pim,pim,pi0,pi0}, {sm,k0,pip,pim,pi0}, {sm,k0,pi0,pi0,pi0},
00289 {sm,kp,pip,pim,pim}, {sm,kp,pim,pi0,pi0}, {lam,kp,pim,pim,pi0},
00290 {lam,k0,pip,pim,pim}, {lam,k0,pim,pi0,pi0}, {s0,k0,pip,pim,pim},
00291 {s0,k0,pim,pi0,pi0}, {s0,kp,pim,pim,pi0}, {sp,k0,pim,pim,pi0},
00292 {sp,kp,pim,pim,pim}, {neu,pip,pim,k0,km}, {neu,pim,pim,kp,k0b},
00293 {neu,pim,pi0,k0,k0b}, {neu,pim,pi0,kp,km}, {neu,pi0,pi0,k0,km},
00294 {pro,pim,pim,k0,k0b}, {pro,pim,pim,kp,km}, {pro,pim,pi0,k0,km}}};
00295
00296 const G4int G4RPGPionInelastic::T32_6bfs[2][5][6] =
00297 {{{pro,pip,pip,pip,pim,pim}, {pro,pip,pip,pim,pi0,pi0},
00298 {pro,pip,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0},
00299 {neu,pip,pip,pip,pim,pi0}},
00300
00301 {{neu,pip,pip,pim,pim,pim}, {neu,pip,pim,pim,pi0,pi0},
00302 {neu,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0},
00303 {pro,pip,pim,pim,pim,pi0}}};
00304
00305 const G4int G4RPGPionInelastic::T32_7bfs[2][6][7] =
00306 {{{pro,pip,pip,pip,pim,pim,pi0}, {pro,pip,pip,pim,pi0,pi0,pi0},
00307 {pro,pip,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim},
00308 {neu,pip,pip,pip,pim,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0}},
00309
00310 {{neu,pip,pip,pim,pim,pim,pi0}, {neu,pip,pim,pim,pi0,pi0,pi0},
00311 {neu,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pim},
00312 {pro,pip,pim,pim,pim,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0}}};
00313
00314 const G4int G4RPGPionInelastic::T32_8bfs[2][7][8] =
00315 {{{pro,pip,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pip,pim,pim,pi0,pi0},
00316 {pro,pip,pip,pim,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0},
00317 {neu,pip,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pip,pim,pi0,pi0,pi0},
00318 {neu,pip,pip,pi0,pi0,pi0,pi0,pi0}},
00319
00320 {{neu,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pip,pim,pim,pim,pi0,pi0},
00321 {neu,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00322 {pro,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pim,pim,pim,pi0,pi0,pi0},
00323 {pro,pim,pim,pi0,pi0,pi0,pi0,pi0}}};
00324
00325 const G4int G4RPGPionInelastic::T32_9bfs[2][8][9] =
00326 {{{pro,pip,pip,pip,pip,pim,pim,pim,pi0}, {pro,pip,pip,pip,pim,pim,pi0,pi0,pi0},
00327 {pro,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00328 {neu,pip,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pip,pim,pim,pi0,pi0},
00329 {neu,pip,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
00330
00331 {{neu,pip,pip,pip,pim,pim,pim,pim,pi0}, {neu,pip,pip,pim,pim,pim,pi0,pi0,pi0},
00332 {neu,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00333 {pro,pip,pip,pip,pim,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pim,pi0,pi0},
00334 {pro,pip,pim,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
00335
00336
00337
00338
00339
00340
00341
00342
00343
00344
00345
00346
00347
00348
00349
00350
00351 const G4float G4RPGPionInelastic::pipPCrossSections[74][30] = {
00352
00353
00354
00355
00356 { 0.00, 1.20, 2.50, 3.80, 5.00, 7.00, 9.00, 15.0, 30.0, 64.0,
00357 130.0,190.0,130.0, 55.7, 27.2, 14.0, 8.50, 13.0, 18.0, 11.0,
00358 8.50, 7.00, 6.20, 5.60, 5.00, 4.50, 4.20, 4.00, 3.80, 3.60 },
00359
00360
00361 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00362 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.16, 0.60, 0.32,
00363 0.19, 0.1, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00364
00365
00366
00367
00368 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00369 0.0, 0.0, 0.0, 0.2, 0.6, 2.4, 8.8, 10.0, 12.0, 6.2,
00370 4.00, 2.40, 1.69, 1.10, 0.73, 0.49, 0.41, 0.31, 0.24, 0.15 },
00371
00372
00373 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00374 0.0, 0.0, 0.0, 0.10, 0.20, 0.60, 1.50, 2.30, 3.60, 3.00,
00375 2.30, 1.70, 1.30, 0.95, 0.69, 0.46, 0.38, 0.27, 0.20, 0.15 },
00376
00377
00378 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00379 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
00380 0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00381
00382
00383 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00384 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
00385 0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00386
00387
00388 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00389 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
00390 0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00391
00392
00393 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00394 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
00395 0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00396
00397
00398 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00399 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
00400 0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00401
00402
00403
00404
00405 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00406 0.0, 0.0, 0.0, 0.0, 0.0, 0.06, 0.20, 0.78, 2.20, 3.20,
00407 3.50, 3.10, 2.70, 2.30, 2.00, 1.50, 1.40, 1.20, 1.00, 0.90 },
00408
00409
00410 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00411 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
00412 2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
00413
00414
00415 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00416 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
00417 2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
00418
00419
00420 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00421 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00422 0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
00423
00424
00425 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00426 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00427 0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
00428
00429
00430 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00431 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,
00432 0.015,0.06, 0.06, 0.05, 0.04,0.032,0.028, 0.02,0.017,0.014 },
00433
00434
00435 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00436 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00437 0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
00438
00439
00440 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00441 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00442 0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
00443
00444
00445 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00446 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00447 0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
00448
00449
00450 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00451 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00452 0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
00453
00454
00455 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00456 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,
00457 0.02, 0.16, 0.18, 0.13, 0.09, 0.06,0.045, 0.03,0.025, 0.02 },
00458
00459
00460 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00461 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00462 0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05,0.045, 0.04 },
00463
00464
00465 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00466 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00467 0.005, 0.05, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05,0.045, 0.04 },
00468
00469
00470 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00471 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00472 0.01, 0.05,0.075, 0.08,0.075, 0.07, 0.06, 0.05,0.045, 0.04 },
00473
00474
00475 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00476 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00477 0.01, 0.05, 0.07,0.065,0.055,0.048, 0.04, 0.03,0.027,0.023 },
00478
00479
00480
00481
00482 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00483 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.05, 0.30, 2.00,
00484 3.20, 3.70, 3.00, 2.50, 2.10, 1.60, 1.40, 1.10, 0.89, 0.70 },
00485
00486
00487 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00488 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.03, 0.21, 1.40,
00489 2.20, 2.60, 2.10, 1.80, 1.50, 1.10, 1.00, 0.80, 0.62, 0.50 },
00490
00491
00492 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00493 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007, 0.02, 0.05, 0.19,
00494 0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
00495
00496
00497 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00498 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007, 0.02, 0.05, 0.19,
00499 0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
00500
00501
00502 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00503 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00504 0.003,0.013, 0.05, 0.09, 0.08, 0.06, 0.05, 0.04,0.036,0.03 },
00505
00506
00507 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00508 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00509 0.001,0.006, 0.02, 0.04, 0.04, 0.03,0.025, 0.02, 0.02, 0.01 },
00510
00511
00512 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00513 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00514 0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
00515
00516
00517 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00518 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00519 0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
00520
00521
00522 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00523 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00524 0.002,0.013, 0.04,0.052,0.059,0.053, 0.05,0.043,0.037, 0.03 },
00525
00526
00527 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00528 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00529 0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
00530
00531
00532 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00533 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00534 0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
00535
00536
00537 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00538 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00539 0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00540
00541
00542 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00543 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00544 0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00545
00546
00547 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00548 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00549 0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00550
00551
00552 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00553 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00554 0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
00555
00556
00557 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00558 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00559 0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
00560
00561
00562 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00563 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00564 0.001,0.005, 0.02,0.065, 0.08, 0.07, 0.06,0.054,0.048, 0.04 },
00565
00566
00567 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00568 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00569 0.001,0.005, 0.02,0.045, 0.05,0.047, 0.04,0.033, 0.03,0.026 },
00570
00571
00572 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00573 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00574 0.001, 0.02, 0.06, 0.06,0.054,0.048,0.042,0.038,0.035, 0.03 },
00575
00576
00577 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00578 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00579 0.001, 0.01, 0.04, 0.05, 0.05, 0.04, 0.04,0.035,0.032, 0.03 },
00580
00581
00582 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00583 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00584 0.001,0.005, 0.01, 0.03, 0.04,0.035, 0.03,0.025,0.025, 0.02 },
00585
00586
00587 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00588 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00589 0.0, 0.02, 0.06, 0.05,0.042,0.038,0.035, 0.03,0.027,0.022 },
00590
00591
00592 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00593 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00594 0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
00595
00596
00597 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00598 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00599 0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
00600
00601
00602
00603
00604 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00605 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
00606 0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00607
00608
00609 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00610 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
00611 0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00612
00613
00614 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00615 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
00616 0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00617
00618
00619 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00620 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
00621 0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00622
00623
00624 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00625 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
00626 0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00627
00628
00629
00630
00631 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00632 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00633 0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00634
00635
00636 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00637 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00638 0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00639
00640
00641 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00642 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00643 0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00644
00645
00646 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00647 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00648 0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00649
00650
00651 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00652 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00653 0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00654
00655
00656 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00657 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00658 0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00659
00660
00661
00662
00663 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00664 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00665 0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00666
00667
00668 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00669 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00670 0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00671
00672
00673 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00674 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00675 0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00676
00677
00678 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00679 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00680 0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00681
00682
00683 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00684 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00685 0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00686
00687
00688 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00689 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00690 0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00691
00692
00693 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00694 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00695 0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00696
00697
00698
00699
00700 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00701 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00702 0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00703
00704
00705 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00706 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00707 0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00708
00709
00710 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00711 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00712 0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00713
00714
00715 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00716 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00717 0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00718
00719
00720 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00721 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00722 0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00723
00724
00725 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00726 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00727 0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00728
00729
00730 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00731 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00732 0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00733
00734
00735 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00736 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00737 0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 } };
00738
00740
00741
00742
00743
00745
00746
00747 G4double G4RPGPionInelastic::pimPtot[30];
00748
00749
00750 G4double G4RPGPionInelastic::t12_dSigma_dMult[8][30];
00751
00752 const G4int G4RPGPionInelastic::pimPindex[8][2] =
00753 {{0, 4}, {5,17}, {18,39}, {40,70}, {71,76}, {77,83}, {84,91}, {92,100}};
00754
00755
00756
00757
00758 const G4int G4RPGPionInelastic::T12_2bfs[2][5][2] =
00759 {{{pro,pim}, {neu,pi0}, {lam,k0}, {s0,k0}, {sm,kp}},
00760
00761 {{neu,pip}, {pro,pi0}, {lam,kp}, {s0,kp}, {sp,k0}}};
00762
00763 const G4int G4RPGPionInelastic::T12_3bfs[2][13][3] =
00764 {{{pro,pim,pi0}, {neu,pip,pim}, {neu,pi0,pi0}, {lam,k0,pi0},
00765 {lam,kp,pim}, {sm,k0,pip}, {sm,kp,pi0}, {sp,k0,pim},
00766 {s0,kp,pim}, {s0,k0,pi0}, {pro,k0,km}, {neu,kp,km},
00767 {neu,k0,k0b}},
00768
00769 {{neu,pip,pi0}, {pro,pip,pim}, {pro,pi0,pi0}, {lam,kp,pi0},
00770 {lam,k0,pip}, {sp,kp,pim}, {sp,k0,pi0}, {sm,kp,pip},
00771 {s0,k0,pip}, {s0,kp,pi0}, {neu,kp,k0b}, {pro,k0,k0b},
00772 {pro,kp,km}}};
00773
00774 const G4int G4RPGPionInelastic::T12_4bfs[2][22][4] =
00775 {{{pro,pip,pim,pim}, {pro,pim,pi0,pi0}, {neu,pip,pim,pi0},
00776 {neu,pi0,pi0,pi0}, {lam,k0,pip,pim}, {lam,k0,pi0,pi0},
00777 {lam,kp,pim,pi0}, {s0,k0,pip,pim}, {s0,k0,pi0,pi0},
00778 {s0,kp,pim,pi0}, {sp,kp,pim,pim}, {sp,k0,pim,pi0},
00779 {sm,kp,pip,pim}, {sm,kp,pi0,pi0}, {sm,k0,pip,pi0},
00780 {pro,pim,kp,km}, {pro,pim,k0,k0b}, {pro,pi0,k0,km},
00781 {neu,pip,k0,km}, {neu,pi0,k0,k0b}, {neu,pi0,kp,km},
00782 {neu,pim,kp,k0b}},
00783
00784 {{neu,pip,pip,pim}, {neu,pip,pi0,pi0}, {pro,pip,pim,pi0},
00785 {pro,pi0,pi0,pi0}, {lam,kp,pip,pim}, {lam,kp,pi0,pi0},
00786 {lam,k0,pip,pi0}, {s0,kp,pip,pim}, {s0,kp,pi0,pi0},
00787 {s0,k0,pip,pi0}, {sm,k0,pip,pip}, {sm,kp,pip,pi0},
00788 {sp,k0,pip,pim}, {sp,k0,pi0,pi0}, {sp,kp,pim,pi0},
00789 {neu,pip,k0,k0b}, {neu,pip,kp,km}, {neu,pi0,kp,k0b},
00790 {pro,pim,kp,k0b}, {pro,pi0,kp,km}, {pro,pi0,k0,k0b},
00791 {pro,pip,k0,km}}};
00792
00793 const G4int G4RPGPionInelastic::T12_5bfs[2][31][5] =
00794 {{{pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}, {neu,pip,pip,pim,pim},
00795 {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0}, {lam,k0,pip,pim,pi0},
00796 {lam,kp,pim,pi0,pi0}, {lam,kp,pip,pim,pim}, {lam,k0,pi0,pi0,pi0},
00797 {s0,kp,pip,pim,pim}, {s0,kp,pim,pi0,pi0}, {s0,k0,pip,pim,pi0},
00798 {s0,k0,pi0,pi0,pi0}, {sp,k0,pip,pim,pim}, {sp,k0,pim,pi0,pi0},
00799 {sp,kp,pim,pim,pi0}, {sm,k0,pip,pip,pim}, {sm,k0,pip,pi0,pi0},
00800 {sm,kp,pip,pim,pi0}, {sm,kp,pi0,pi0,pi0}, {pro,pim,pi0,kp,km},
00801 {pro,pim,pi0,k0,k0b}, {pro,pip,pim,k0,km}, {pro,pi0,pi0,k0,km},
00802 {pro,pim,pim,kp,k0b}, {neu,pip,pim,kp,km}, {neu,pip,pim,k0,k0b},
00803 {neu,pip,pi0,k0,km}, {neu,pim,pi0,kp,k0b}, {neu,pi0,pi0,k0,k0b},
00804 {neu,pi0,pi0,kp,km}},
00805
00806 {{neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}, {pro,pip,pip,pim,pim},
00807 {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0}, {lam,kp,pip,pim,pi0},
00808 {lam,k0,pip,pi0,pi0}, {lam,k0,pip,pip,pim}, {lam,kp,pi0,pi0,pi0},
00809 {s0,k0,pip,pip,pim}, {s0,k0,pip,pi0,pi0}, {s0,kp,pip,pim,pi0},
00810 {s0,kp,pi0,pi0,pi0}, {sm,kp,pip,pip,pim}, {sm,kp,pip,pi0,pi0},
00811 {sm,k0,pip,pip,pi0}, {sp,kp,pip,pim,pim}, {sp,kp,pim,pi0,pi0},
00812 {sp,k0,pip,pim,pi0}, {sp,k0,pi0,pi0,pi0}, {neu,pip,pi0,k0,k0b},
00813 {neu,pip,pi0,kp,km}, {neu,pip,pim,kp,k0b}, {neu,pi0,pi0,kp,k0b},
00814 {neu,pip,pip,k0,km}, {pro,pip,pim,k0,k0b}, {pro,pip,pim,kp,km},
00815 {pro,pim,pi0,kp,k0b}, {pro,pip,pi0,k0,km}, {pro,pi0,pi0,kp,km},
00816 {pro,pi0,pi0,k0,k0b}}};
00817
00818 const G4int G4RPGPionInelastic::T12_6bfs[2][6][6] =
00819 {{{pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0},
00820 {pro,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pim,pim,pi0},
00821 {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0}},
00822
00823 {{neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0},
00824 {neu,pip,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pi0},
00825 {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0}}};
00826
00827 const G4int G4RPGPionInelastic::T12_7bfs[2][7][7] =
00828 {{{pro,pip,pip,pim,pim,pim,pi0}, {pro,pip,pim,pim,pi0,pi0,pi0},
00829 {pro,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim},
00830 {neu,pip,pip,pim,pim,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0},
00831 {neu,pi0,pi0,pi0,pi0,pi0,pi0}},
00832
00833 {{neu,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pim,pi0,pi0,pi0},
00834 {neu,pip,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim},
00835 {pro,pip,pip,pim,pim,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0},
00836 {pro,pi0,pi0,pi0,pi0,pi0,pi0}}};
00837
00838 const G4int G4RPGPionInelastic::T12_8bfs[2][8][8] =
00839 {{{pro,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
00840 {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00841 {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
00842 {neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0}},
00843
00844 {{neu,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
00845 {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
00846 {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
00847 {pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0}}};
00848
00849 const G4int G4RPGPionInelastic::T12_9bfs[2][9][9] =
00850 {{{pro,pip,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0},
00851 {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00852 {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00853 {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0,pi0},
00854 {neu,pip,pip,pip,pip,pim,pim,pim,pim}},
00855
00856 {{neu,pip,pip,pip,pip,pim,pim,pim,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0},
00857 {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00858 {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00859 {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0,pi0},
00860 {pro,pip,pip,pip,pip,pim,pim,pim,pim}}};
00861
00862
00863
00864
00865
00866
00867
00868
00869
00870
00871
00872
00873
00874
00875
00876 const G4float G4RPGPionInelastic::pimPCrossSections[101][30] = {
00877
00878
00879
00880
00881
00882 { 0.0, 1.1, 1.2, 1.4, 1.5, 1.8, 2.0, 3.0, 3.4, 7.0,
00883 14.0, 24.0, 14.7, 10.5, 11.8, 20.0, 14.0, 25.0, 12.0, 9.5,
00884 8.0, 7.0, 6.0, 5.7, 5.0, 4.6, 4.3, 4.0, 3.8, 3.7 },
00885
00886
00887 { 0.0, 2.4, 2.8, 3.3, 4.5, 5.7, 6.3, 9.0, 11.0, 17.0,
00888 30.0, 43.0, 30.0, 16.5, 11.0, 7.0, 4.3, 5.0, 2.0, 0.9,
00889 0.5, 0.24, 0.15,0.094,0.061,0.048,0.035,0.023,0.018,0.014 },
00890
00891
00892 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00893 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.65, 0.28, 0.16,
00894 0.13,0.075, 0.05,0.032,0.022,0.015,0.011,0.008,0.006,0.004 },
00895
00896
00897 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00898 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.28, 0.18, 0.11,
00899 0.091,0.055,0.037,0.025,0.018,0.012,0.008,0.005,0.004,0.003 },
00900
00901
00902 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00903 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2, 0.24, 0.09,
00904 0.04,0.012,0.004,0.002, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 },
00905
00906
00907
00908
00909 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00910 0.0, 0.0, 0.05, 0.18, 0.86, 4.4, 5.2, 6.6, 5.4, 4.4,
00911 3.5, 2.5, 2.0, 1.4, 0.97, 0.68, 0.55, 0.36, 0.3, 0.22 },
00912
00913
00914 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00915 0.0, 0.0, 0.09, 0.9, 3.3, 6.6, 8.9, 8.2, 8.8, 6.6,
00916 5.2, 3.8, 2.9, 1.9, 1.3, 0.9, 0.75, 0.38, 0.31, 0.22 },
00917
00918
00919 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00920 0.0, 0.0, 0.22, 0.74, 1.8, 2.7, 3.0, 2.5, 1.3, 0.64,
00921 0.41, 0.24, 0.15, 0.10,0.065,0.042, 0.0, 0.0, 0.0, 0.0 },
00922
00923
00924 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00925 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.049, 0.18,
00926 0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01,0.007,0.004 },
00927
00928
00929 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00930 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.035, 0.12,
00931 0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01,0.007,0.004 },
00932
00933
00934 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00935 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10,
00936 0.13, 0.06, 0.03, 0.02, 0.01,0.003,0.001, 0.0, 0.0, 0.0 },
00937
00938
00939 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00940 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
00941 0.05, 0.03, 0.02, 0.01,0.005,0.002,0.001, 0.0, 0.0, 0.0 },
00942
00943
00944 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00945 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
00946 0.06, 0.04, 0.02, 0.01, 0.01,0.005,0.003,0.001, 0.0, 0.0 },
00947
00948
00949 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00950 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
00951 0.07, 0.04, 0.03, 0.02, 0.02, 0.01,0.004,0.002,0.001, 0.0 },
00952
00953
00954 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00955 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.09,
00956 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01,0.005,0.003,0.002 },
00957
00958
00959 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00960 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03,
00961 0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00962
00963
00964 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00965 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
00966 0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01, 0.0, 0.0, 0.0 },
00967
00968
00969 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00970 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06,
00971 0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0 },
00972
00973
00974
00975
00976 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00977 0.0, 0.0, 0.0, 0.0, 0.01, 0.05, 0.2, 0.48, 0.88, 1.8,
00978 2.0, 2.2, 2.0, 1.9, 1.6, 1.4, 1.3, 1.1, 1.05, 1.0 },
00979
00980
00981 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00982 0.0, 0.0, 0.0, 0.0, 0.01, 0.06, 0.18, 0.42, 1.0, 1.9,
00983 2.4, 2.4, 2.0, 1.8, 1.6, 1.4, 1.3, 1.2, 1.1, 1.0 },
00984
00985
00986 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00987 0.0, 0.0, 0.0, 0.0, 0.1, 0.3, 0.5, 1.0, 2.0, 2.0,
00988 1.0, 0.8, 0.6, 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1 },
00989
00990
00991 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00992 0.0, 0.0, 0.0, 0.0, 0.1, 0.55, 1.0, 0.87, 0.7, 0.6,
00993 0.5, 0.4, 0.35, 0.29, 0.21, 0.14, 0.1, 0.06, 0.03, 0.0 },
00994
00995
00996 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
00997 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,
00998 0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
00999
01000
01001 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01002 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,
01003 0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
01004
01005
01006 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01007 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.009,
01008 0.06, 0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02 },
01009
01010
01011 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01012 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,
01013 0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
01014
01015
01016 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01017 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,
01018 0.04,0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
01019
01020
01021 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01022 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01023 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01024
01025
01026 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01027 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01028 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01029
01030
01031 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01032 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01033 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01034
01035
01036 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01037 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01038 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01039
01040
01041 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01042 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01043 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01044
01045
01046 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01047 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
01048 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01049
01050
01051 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01052 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01053 0.003,0.06, 0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05 },
01054
01055
01056 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01057 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01058 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01059
01060
01061 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01062 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01063 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01064
01065
01066 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01067 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01068 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01069
01070
01071 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01072 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01073 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01074
01075
01076 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01077 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01078 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01079
01080
01081 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01082 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01083 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01084
01085
01086
01087
01088 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01089 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,0.022, 0.10, 0.70,
01090 1.3, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62 },
01091
01092
01093 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01094 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,0.022, 0.10, 0.70,
01095 1.3, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62 },
01096
01097
01098 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01099 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,0.029, 0.10, 0.30,
01100 0.56, 0.93, 1.2, 1.2, 1.2, 0.94, 0.74, 0.53, 0.40, 0.30 },
01101
01102
01103 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01104 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.006, 0.02, 0.09, 0.63,
01105 1.2, 1.7, 2.0, 1.8, 1.5, 1.3, 1.1, 0.80, 0.70, 0.60 },
01106
01107
01108 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01109 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.006, 0.02, 0.08, 0.60,
01110 1.0, 1.5, 2.0, 1.6, 1.4, 1.1, 1.0, 0.70, 0.60, 0.50 },
01111
01112
01113 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01114 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01115 0.007, 0.05, 0.07, 0.08, 0.08, 0.07,0.063, 0.06,0.053,0.048 },
01116
01117
01118 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01119 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01120 0.003,0.025,0.035, 0.04, 0.04,0.035, 0.03, 0.03,0.025,0.024 },
01121
01122
01123 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01124 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01125 0.002, 0.01,0.036, 0.04, 0.04,0.033,0.031,0.029,0.025,0.022 },
01126
01127
01128 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01129 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01130 0.001,0.005,0.018, 0.02, 0.02,0.017,0.016,0.015,0.012,0.011 },
01131
01132
01133 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01134 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01135 0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
01136
01137
01138 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01139 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01140 0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
01141
01142
01143 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01144 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01145 0.003, 0.02, 0.03, 0.04, 0.04,0.035, 0.03, 0.03,0.025, 0.02 },
01146
01147
01148 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01149 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01150 0.001, 0.01,0.015, 0.02, 0.02,0.017,0.015,0.015,0.012, 0.01 },
01151
01152
01153 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01154 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01155 0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
01156
01157
01158 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01159 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01160 0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
01161
01162
01163 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01164 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01165 0.001,0.002,0.006,0.009,0.008,0.007,0.006,0.005,0.005,0.004 },
01166
01167
01168 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01169 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01170 0.002,0.005,0.015,0.014,0.012,0.009,0.008,0.006,0.005,0.004 },
01171
01172
01173 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01174 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01175 0.001,0.002,0.008,0.007,0.006,0.005,0.004,0.003,0.003,0.002 },
01176
01177
01178 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01179 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01180 0.002,0.007,0.023,0.025,0.021,0.018,0.015,0.013,0.010,0.009 },
01181
01182
01183 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01184 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01185 0.001,0.004,0.011,0.013,0.011,0.009,0.008,0.007,0.006,0.005 },
01186
01187
01188 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01189 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01190 0.0, 0.01, 0.07, 0.07, 0.06,0.055, 0.05,0.047,0.042, 0.04 },
01191
01192
01193 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01194 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01195 0.0,0.003, 0.02, 0.02,0.018,0.017,0.015,0.012,0.011, 0.01 },
01196
01197
01198 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01199 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01200 0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
01201
01202
01203 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01204 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01205 0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
01206
01207
01208 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01209 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01210 0.0,0.007, 0.04, 0.05, 0.04,0.035, 0.03, 0.03,0.027,0.025 },
01211
01212
01213 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01214 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01215 0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01216
01217
01218 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01219 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01220 0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01221
01222
01223 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01224 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01225 0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01226
01227
01228 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01229 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01230 0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01231
01232
01233 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01234 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01235 0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01236
01237
01238 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01239 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01240 0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01241
01242
01243
01244
01245 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01246 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.04,
01247 0.06, 0.11, 0.16, 0.22, 0.31, 0.34, 0.3, 0.24, 0.19, 0.16 },
01248
01249
01250 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01251 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
01252 0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
01253
01254
01255 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01256 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
01257 0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
01258
01259
01260 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01261 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
01262 0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
01263
01264
01265 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01266 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
01267 0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
01268
01269
01270 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01271 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
01272 0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
01273
01274
01275
01276
01277 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01278 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014, 0.03,
01279 0.05, 0.11, 0.19, 0.37, 0.67, 0.75, 0.7, 0.59, 0.52, 0.47 },
01280
01281
01282 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01283 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
01284 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
01285
01286
01287 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01288 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
01289 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
01290
01291
01292 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01293 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,0.006,
01294 0.012,0.029,0.062, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13 },
01295
01296
01297 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01298 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
01299 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
01300
01301
01302 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01303 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
01304 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
01305
01306
01307 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01308 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
01309 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
01310
01311
01312
01313
01314 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01315 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01316 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01317
01318
01319 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01320 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01321 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01322
01323
01324 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01325 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01326 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01327
01328
01329 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01330 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01331 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01332
01333
01334 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01335 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01336 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01337
01338
01339 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01340 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01341 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01342
01343
01344 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01345 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01346 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01347
01348
01349 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01350 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01351 0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01352
01353
01354
01355
01356 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01357 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01358 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01359
01360
01361 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01362 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01363 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01364
01365
01366 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01367 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01368 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01369
01370
01371 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01372 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01373 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01374
01375
01376 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01377 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01378 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01379
01380
01381 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01382 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01383 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01384
01385
01386 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01387 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01388 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01389
01390
01391 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01392 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01393 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01394
01395
01396 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01397 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
01398 0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 } };
01399
01400
01401