Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4RPGNucleonInelastic.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27
28#include "G4RPGNucleonInelastic.hh"
29#include "G4SystemOfUnits.hh"
30#include "Randomize.hh"
31
32G4RPGNucleonInelastic::G4RPGNucleonInelastic(const G4String& modelName)
33 :G4RPGInelastic(modelName)
34{
35 SetMinEnergy( 0.0 );
36 SetMaxEnergy( 30.*GeV );
37
38 // Initialize t1_dSigma_dMult, t0_dSigma_dMult,
39 // nucleon-nucleon inelastic cross sections for a given multiplicity
40 // for |T_z| = 1 and 0, respectively
41
42 G4int i, j, k;
43 G4int start, stop;
44
45 for (j = 0; j < 8; j++) {
46 start = pPindex[j][0];
47 stop = pPindex[j][1] + 1;
48 for (k = 0; k < 30; k++) {
49 t1_dSigma_dMult[j][k] = 0.0;
50 for (i = start; i < stop; i++) t1_dSigma_dMult[j][k] += pPCrossSections[i][k];
51 }
52
53 start = pNindex[j][0];
54 stop = pNindex[j][1] + 1;
55 for (k = 0; k < 30; k++) {
56 t0_dSigma_dMult[j][k] = 0.0;
57 for (i = start; i < stop; i++) t0_dSigma_dMult[j][k] += pNCrossSections[i][k];
58 }
59 }
60
61 // Initialize total cross section array
62
63 for (k = 0; k < 30; k++) {
64 pPtot[k] = 0.0;
65 pNtot[k] = 0.0;
66 for (j = 0; j < 8; j++) {
67 pPtot[k] += t1_dSigma_dMult[j][k];
68 pNtot[k] += t0_dSigma_dMult[j][k];
69 }
70 }
71
72 // printCrossSections();
73}
74
75/*
76void G4RPGNucleonInelastic::printCrossSections() const
77{
78 G4cout << " pp total cross sections (mb) " << G4endl;
79 for (G4int i = 0; i < 5; i++) {
80 G4int istart = i*6;
81 G4int istop = istart + 6;
82 for (G4int t = istart; t < istop; t++) G4cout << pPtot[t] << " " ;
83 G4cout << G4endl;
84 }
85 G4cout << G4endl;
86
87 G4cout << " pn total cross sections (mb) " << G4endl;
88 for (G4int i = 0; i < 5; i++) {
89 G4int istart = i*6;
90 G4int istop = istart + 6;
91 for (G4int t = istart; t < istop; t++) G4cout << pNtot[t] << " " ;
92 G4cout << G4endl;
93 }
94}
95*/
96
97
98G4int G4RPGNucleonInelastic::GetMultiplicityT0(G4double KE) const
99{
100 G4double multint(0.);
101 std::vector<G4double> sigma;
102
103 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
104 G4int k = epair.first;
105 G4double fraction = epair.second;
106
107 for(G4int j = 0; j < 8; j++) {
108 multint = t0_dSigma_dMult[j][k]
109 + fraction*(t0_dSigma_dMult[j][k+1] - t0_dSigma_dMult[j][k]);
110 sigma.push_back(multint);
111 }
112
113 return sampleFlat(sigma) + 2;
114}
115
116
117G4int G4RPGNucleonInelastic::GetMultiplicityT1(G4double KE) const
118{
119 G4double multint(0.);
120 std::vector<G4double> sigma;
121
122 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
123 G4int k = epair.first;
124 G4double fraction = epair.second;
125
126 for(G4int j = 0; j < 8; j++) {
127 multint = t1_dSigma_dMult[j][k]
128 + fraction*(t1_dSigma_dMult[j][k+1] - t1_dSigma_dMult[j][k]);
129 sigma.push_back(multint);
130 }
131
132 return sampleFlat(sigma) + 2;
133}
134
135
136std::vector<G4int>
137G4RPGNucleonInelastic::GetFSPartTypesForT0(G4int mult, G4double KE) const
138{
139 G4int i;
140 G4double sigint(0.);
141 std::vector<G4double> sigma;
142
143 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
144 G4int k = epair.first;
145 G4double fraction = epair.second;
146
147 G4int start = pNindex[mult-2][0];
148 G4int stop = pNindex[mult-2][1];
149
150 for(i = start; i < stop; i++) {
151 sigint = pNCrossSections[i][k]
152 + fraction*(pNCrossSections[i][k+1] - pNCrossSections[i][k]);
153 sigma.push_back(sigint);
154 }
155
156 G4int channel = sampleFlat(sigma);
157
158 std::vector<G4int> kinds;
159
160 if (mult == 2) {
161 for(i = 0; i < mult; i++) kinds.push_back(T0_2bfs[channel][i]);
162 } else if (mult == 3) {
163 for(i = 0; i < mult; i++) kinds.push_back(T0_3bfs[channel][i]);
164 } else if (mult == 4) {
165 for(i = 0; i < mult; i++) kinds.push_back(T0_4bfs[channel][i]);
166 } else if (mult == 5) {
167 for(i = 0; i < mult; i++) kinds.push_back(T0_5bfs[channel][i]);
168 } else if (mult == 6) {
169 for(i = 0; i < mult; i++) kinds.push_back(T0_6bfs[channel][i]);
170 } else if (mult == 7) {
171 for(i = 0; i < mult; i++) kinds.push_back(T0_7bfs[channel][i]);
172 } else if (mult == 8) {
173 for(i = 0; i < mult; i++) kinds.push_back(T0_8bfs[channel][i]);
174 } else if (mult == 9) {
175 for(i = 0; i < mult; i++) kinds.push_back(T0_9bfs[channel][i]);
176 } else {
177 G4cout << " Illegal multiplicity " << G4endl;
178 }
179
180 return kinds;
181}
182
183
184std::vector<G4int>
185G4RPGNucleonInelastic::GetFSPartTypesForT1(G4int mult, G4double KE,
186 G4int tzindex) const
187{
188 G4int i;
189 G4double sigint(0.);
190 std::vector<G4double> sigma;
191
192 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
193 G4int k = epair.first;
194 G4double fraction = epair.second;
195
196 G4int start = pPindex[mult-2][0];
197 G4int stop = pPindex[mult-2][1];
198
199 for(i = start; i < stop; i++) {
200 sigint = pPCrossSections[i][k]
201 + fraction*(pPCrossSections[i][k+1] - pPCrossSections[i][k]);
202 sigma.push_back(sigint);
203 }
204
205 G4int channel = sampleFlat(sigma);
206
207 std::vector<G4int> kinds;
208
209 if (mult == 2) {
210 for(i = 0; i < mult; i++) kinds.push_back(T1_2bfs[tzindex][channel][i]);
211 } else if (mult == 3) {
212 for(i = 0; i < mult; i++) kinds.push_back(T1_3bfs[tzindex][channel][i]);
213 } else if (mult == 4) {
214 for(i = 0; i < mult; i++) kinds.push_back(T1_4bfs[tzindex][channel][i]);
215 } else if (mult == 5) {
216 for(i = 0; i < mult; i++) kinds.push_back(T1_5bfs[tzindex][channel][i]);
217 } else if (mult == 6) {
218 for(i = 0; i < mult; i++) kinds.push_back(T1_6bfs[tzindex][channel][i]);
219 } else if (mult == 7) {
220 for(i = 0; i < mult; i++) kinds.push_back(T1_7bfs[tzindex][channel][i]);
221 } else if (mult == 8) {
222 for(i = 0; i < mult; i++) kinds.push_back(T1_8bfs[tzindex][channel][i]);
223 } else if (mult == 9) {
224 for(i = 0; i < mult; i++) kinds.push_back(T1_9bfs[tzindex][channel][i]);
225 } else {
226 G4cout << " Illegal multiplicity " << G4endl;
227 }
228
229 return kinds;
230}
231
232///////////////////////////////////////////////////////////////////////////////
233// //
234// p p and n n (|Tz| = 1) cross sections //
235// and final state particle types //
236// //
237///////////////////////////////////////////////////////////////////////////////
238
239// Total p p cross section as a function of kinetic energy
240G4ThreadLocal G4double G4RPGNucleonInelastic::pPtot[30];
241
242// p p multiplicities as a function of kinetic energy
243G4ThreadLocal G4double G4RPGNucleonInelastic::t1_dSigma_dMult[8][30];
244
245const G4int G4RPGNucleonInelastic::pPindex[8][2] =
246 {{0, 0}, {1, 6}, {7,24}, {25,56}, {57,63}, {64,71}, {72,81}, {82,92}};
247
248// Outgoing particle types of a given multiplicity
249// T1_nbfs = final state types for p p and n n
250
251const G4int G4RPGNucleonInelastic::T1_2bfs[2][1][2] =
252 {{{pro,pro}},
253
254 {{neu,neu}}};
255
256const G4int G4RPGNucleonInelastic::T1_3bfs[2][6][3] =
257 {{{pro,pro,pi0}, {pro,neu,pip}, {pro,lam,kp},
258 {pro,s0,kp}, {pro,sp,k0}, {neu,sp,kp}},
259
260 {{neu,neu,pi0}, {pro,neu,pim}, {neu,lam,k0},
261 {neu,s0,k0}, {neu,sm,kp}, {pro,sm,k0}}};
262
263const G4int G4RPGNucleonInelastic::T1_4bfs[2][18][4] =
264 {{{pro,pro,pip,pim},{pro,neu,pip,pi0},{pro,pro,pi0,pi0},
265 {neu,neu,pip,pip},{pro,lam,kp,pi0}, {pro,lam,k0,pip},
266 {neu,lam,kp,pip}, {neu,s0,kp,pip}, {pro,s0,kp,pi0},
267 {pro,s0,k0,pip}, {pro,sm,kp,pip}, {pro,sp,k0,pi0},
268 {neu,sp,k0,pip}, {pro,sp,kp,pim}, {neu,sp,kp,pi0},
269 {pro,pro,k0,k0b}, {pro,pro,kp,km}, {pro,neu,kp,k0b}},
270
271 {{neu,neu,pip,pim},{pro,neu,pim,pi0},{neu,neu,pi0,pi0},
272 {pro,pro,pim,pim},{neu,lam,k0,pi0}, {neu,lam,kp,pim},
273 {pro,lam,k0,pim}, {pro,s0,k0,pim}, {neu,s0,k0,pi0},
274 {neu,s0,kp,pim}, {neu,sp,k0,pim}, {neu,sm,kp,pi0},
275 {pro,sm,kp,pim}, {neu,sm,k0,pip}, {pro,sm,k0,pi0},
276 {neu,neu,kp,km}, {neu,neu,k0,k0b}, {pro,neu,k0,km}}};
277
278const G4int G4RPGNucleonInelastic::T1_5bfs[2][32][5] =
279 {{{pro,pro,pip,pim,pi0},{pro,pro,pi0,pi0,pi0},{pro,neu,pip,pip,pim},
280 {pro,neu,pip,pi0,pi0},{neu,neu,pip,pip,pi0},{pro,lam,kp,pip,pim},
281 {pro,lam,kp,pi0,pi0}, {pro,lam,k0,pip,pi0}, {pro,s0,kp,pip,pim},
282 {pro,s0,kp,pi0,pi0}, {pro,s0,k0,pip,pi0}, {pro,sp,k0,pip,pim},
283 {pro,sp,k0,pi0,pi0}, {pro,sp,kp,pim,pi0}, {pro,sm,kp,pip,pi0},
284 {pro,sm,k0,pip,pip}, {neu,lam,kp,pip,pi0}, {neu,lam,k0,pip,pip},
285 {neu,s0,kp,pip,pi0}, {neu,s0,k0,pip,pip}, {neu,sp,k0,pip,pi0},
286 {neu,sp,kp,pip,pim}, {neu,sp,kp,pi0,pi0}, {neu,sm,kp,pip,pip},
287 {pro,pro,pip,k0,km}, {pro,pro,pim,kp,k0b}, {pro,pro,pi0,k0,k0b},
288 {pro,pro,pi0,kp,km}, {pro,neu,pip,k0,k0b}, {pro,neu,pip,kp,km},
289 {pro,neu,pi0,kp,k0b}, {neu,neu,pip,kp,k0b}},
290
291 {{neu,neu,pip,pim,pi0},{neu,neu,pi0,pi0,pi0},{pro,neu,pip,pim,pim},
292 {pro,neu,pim,pi0,pi0},{pro,pro,pim,pim,pi0},{neu,lam,k0,pip,pim},
293 {neu,lam,k0,pi0,pi0}, {neu,lam,kp,pim,pi0}, {neu,s0,k0,pip,pim},
294 {neu,s0,k0,pi0,pi0}, {neu,s0,kp,pim,pi0}, {neu,sm,kp,pip,pim},
295 {neu,sm,kp,pi0,pi0}, {neu,sm,k0,pip,pi0}, {neu,sp,k0,pim,pi0},
296 {neu,sp,kp,pim,pim}, {pro,lam,k0,pim,pi0}, {pro,lam,kp,pim,pim},
297 {pro,s0,k0,pim,pi0}, {pro,s0,kp,pim,pim}, {pro,sm,kp,pim,pi0},
298 {pro,sm,k0,pip,pim}, {pro,sm,k0,pi0,pi0}, {pro,sp,k0,pim,pim},
299 {neu,neu,pim,kp,k0b}, {neu,neu,pip,k0,km}, {neu,neu,pi0,kp,km},
300 {neu,neu,pi0,k0,k0b}, {pro,neu,pim,kp,km}, {pro,neu,pim,k0,k0b},
301 {pro,neu,pi0,k0,km}, {pro,pro,pim,k0,km}}};
302
303const G4int G4RPGNucleonInelastic::T1_6bfs[2][7][6] =
304 {{{pro,pro,pip,pip,pim,pim},{pro,pro,pip,pim,pi0,pi0},
305 {pro,pro,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pi0},
306 {pro,neu,pip,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim},
307 {neu,neu,pip,pip,pi0,pi0}},
308
309 {{neu,neu,pip,pip,pim,pim},{neu,neu,pip,pim,pi0,pi0},
310 {neu,neu,pi0,pi0,pi0,pi0},{pro,neu,pip,pim,pim,pi0},
311 {pro,neu,pim,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pim},
312 {pro,pro,pim,pim,pi0,pi0}}};
313
314const G4int G4RPGNucleonInelastic::T1_7bfs[2][8][7] =
315 {{{pro,pro,pip,pip,pim,pim,pi0},{pro,pro,pip,pim,pi0,pi0,pi0},
316 {pro,pro,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pip,pim,pim},
317 {pro,neu,pip,pip,pim,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0},
318 {neu,neu,pip,pip,pip,pim,pi0},{neu,neu,pip,pip,pi0,pi0,pi0}},
319
320 {{neu,neu,pip,pip,pim,pim,pi0},{neu,neu,pip,pim,pi0,pi0,pi0},
321 {neu,neu,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pim,pim},
322 {pro,neu,pip,pim,pim,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0},
323 {pro,pro,pip,pim,pim,pim,pi0},{pro,pro,pim,pim,pi0,pi0,pi0}}};
324
325const G4int G4RPGNucleonInelastic::T1_8bfs[2][10][8] =
326 {{{pro,pro,pip,pip,pip,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pi0,pi0},
327 {pro,pro,pip,pim,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0},
328 {pro,neu,pip,pip,pip,pim,pim,pi0},{pro,neu,pip,pip,pim,pi0,pi0,pi0},
329 {pro,neu,pip,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pip,pim,pim},
330 {neu,neu,pip,pip,pip,pim,pi0,pi0},{neu,neu,pip,pip,pi0,pi0,pi0,pi0}},
331
332 {{neu,neu,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pim,pim,pi0,pi0},
333 {neu,neu,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0},
334 {pro,neu,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pim,pim,pi0,pi0,pi0},
335 {pro,neu,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim,pim},
336 {pro,pro,pip,pim,pim,pim,pi0,pi0},{pro,pro,pim,pim,pi0,pi0,pi0,pi0}}};
337
338const G4int G4RPGNucleonInelastic::T1_9bfs[2][11][9] =
339{{{pro,pro,pip,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pip,pim,pim,pi0,pi0,pi0},
340 {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
341 {pro,neu,pip,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pip,pim,pim,pi0,pi0},
342 {pro,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
343 {neu,neu,pip,pip,pip,pip,pim,pim,pi0},{neu,neu,pip,pip,pip,pim,pi0,pi0,pi0},
344 {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}},
345
346 {{neu,neu,pip,pip,pip,pim,pim,pim,pi0},{neu,neu,pip,pip,pim,pim,pi0,pi0,pi0},
347 {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
348 {pro,neu,pip,pip,pip,pim,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pim,pi0,pi0},
349 {pro,neu,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
350 {pro,pro,pip,pip,pim,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pim,pi0,pi0,pi0},
351 {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}}};
352
353//
354// Cross sections (in mb) for p p -> 2-9 body final states
355//
356// first index: 0: channels for mult = 2
357// 1-6: channels for mult = 3
358// 7-24: channels for mult = 4
359// 25-56: channels for mult = 5
360// 57-63: channels for mult = 6
361// 64-71: channels for mult = 7
362// 72-81: channels for mult = 8
363// 82-92: channels for mult = 9
364//
365// second index: kinetic energy
366//
367
368const G4float G4RPGNucleonInelastic::pPCrossSections[93][30] = {
369//
370// multiplicity 2 (1 channel)
371//
372// p p (n n)
373 { 0.0f,330.0f,240.0f,160.0f,110.0f, 85.0f, 63.0f, 44.0f, 33.0f, 28.0f,
374 25.0f, 24.0f, 23.0f, 23.0f, 26.3f, 26.1f, 25.0f, 23.5f, 21.0f, 18.0f,
375 16.0f, 14.3f, 12.5f, 11.2f, 10.3f, 9.6f, 9.0f, 8.5f, 8.0f, 7.7f },
376//
377// multiplicity 3 (6 channels)
378//
379// p p pi0 (n n pi0)
380 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
381 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.4f, 4.0f, 4.3f, 4.0f, 4.0f,
382 3.6f, 3.0f, 2.8f, 2.5f, 1.7f, 1.3f, 1.1f, 1.0f, 0.9f, 0.85f },
383
384// p n pi+ (p n pi-)
385 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
386 0.0f, 0.0f, 0.0f, 0.0f, 0.7f, 4.5f, 15.0f, 19.1f, 18.0f, 16.0f,
387 13.0f, 10.0f, 8.2f, 6.0f, 4.3f, 3.3f, 2.6f, 2.0f, 1.65f, 1.4f },
388
389// p L K+ (n L K0)
390 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
391 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.012f,
392 0.03f, 0.06f, 0.06f, 0.055f,0.05f, 0.047f,0.043f,0.04f, 0.037f,0.033f },
393
394// p S0 K+ (n S0 K0)
395 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
396 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
397 0.006f,0.02f, 0.027f,0.026f,0.021f,0.018f,0.015f,0.011f,0.009f,0.007f },
398
399// p S+ K0 (n S- K+)
400 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
401 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
402 0.013f,0.025f,0.03f, 0.029f,0.027f,0.026f,0.024f,0.022f,0.021f,0.019f },
403
404// n S+ K+ (p S- K0)
405 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
406 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
407 0.015f,0.06f, 0.07f, 0.065f,0.05f, 0.04f, 0.033f,0.026f,0.02f, 0.015f },
408//
409// multiplicity 4 (18 channels)
410//
411// p p pi+ pi- (n n pi+ pi-)
412 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
413 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.05f, 0.6f, 1.9f,
414 2.8f, 3.0f, 3.0f, 2.8f, 2.5f, 2.1f, 1.9f, 1.6f, 1.4f, 1.2f },
415
416// p n pi+ pi0 (p n pi- pi0)
417 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
418 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.05f, 0.6f, 3.5f,
419 4.0f, 3.9f, 3.5f, 3.1f, 2.8f, 2.4f, 2.2f, 1.9f, 1.7f, 1.5f },
420
421// p p pi0 pi0 (n n pi0 pi0)
422 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
423 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.02f, 0.24f, 0.76f,
424 1.1f, 1.2f, 1.2f, 1.1f, 1.0f, 0.84f, 0.76f, 0.64f, 0.56f, 0.48f },
425
426// n n pi+ pi+ (p p pi- pi-)
427 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
428 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.02f, 0.24f, 1.4f,
429 1.6f, 1.6f, 1.4f, 1.2f, 1.1f, 1.0f, 0.88f, 0.76f, 0.68f, 0.6f },
430
431// L K+ p pi0 (L K0 n pi0)
432 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
433 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
434 0.002f,0.015f,0.05f, 0.06f, 0.052f,0.042f,0.037f,0.029f,0.025f,0.020f },
435
436// L K0 p pi+ (L K+ n pi-)
437 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
438 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
439 0.002f,0.015f,0.06f, 0.086f,0.09f, 0.082f,0.072f,0.06f, 0.051f,0.043f },
440
441// L K+ n pi+ (L K0 p pi-)
442 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
443 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
444 0.002f,0.01f, 0.04f, 0.08f, 0.066f,0.058f,0.05f, 0.04f, 0.035f,0.03f },
445
446// S0 K+ n pi+ (S0 K0 p pi-)
447 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
448 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
449 0.0f, 0.003f,0.02f, 0.03f, 0.03f, 0.025f,0.02f, 0.015f,0.011f,0.01f },
450
451// S0 K+ p pi0 (S0 K0 n pi0)
452 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
453 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
454 0.0f, 0.005f,0.02f, 0.025f,0.022f,0.02f, 0.015f,0.01f, 0.008f,0.007f },
455
456// S0 K0 p pi+ (S0 K+ n pi-)
457 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
458 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
459 0.0f, 0.003f,0.012f,0.04f, 0.037f,0.03f, 0.027f,0.022f,0.019f,0.016f },
460
461// S- K+ p pi+ (S+ K0 n pi-)
462 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
463 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
464 0.0f, 0.004f,0.016f,0.037f,0.031f,0.028f,0.023f,0.02f, 0.017f,0.014f },
465
466// S+ K0 p pi0 (S- K+ n pi0)
467 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
468 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
469 0.0f, 0.003f,0.007f,0.021f,0.021f,0.019f,0.017f,0.014f,0.012f,0.01f },
470
471// S+ K0 n pi+ (S- K+ p pi-)
472 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
473 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
474 0.0f, 0.003f,0.009f,0.036f,0.052f,0.043f,0.038f,0.03f, 0.026f,0.02f },
475
476// S+ K+ p pi- (S- K0 n pi+)
477 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
478 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
479 0.0f, 0.004f,0.012f,0.038f,0.037f,0.03f, 0.026f,0.02f, 0.017f,0.014f },
480
481// S+ K+ n pi0 (S- K0 p pi0)
482 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
483 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
484 0.0f, 0.004f,0.012f,0.038f,0.037f,0.03f, 0.026f,0.02f, 0.017f,0.014f },
485
486// p p K0 K0bar (n n K+ K-)
487 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
488 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
489 0.0f, 0.001f,0.003f,0.009f,0.02f, 0.02f, 0.017f,0.014f,0.012f,0.009f },
490
491// p p K+ K- (n n K0 K0bar)
492 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
493 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
494 0.0f, 0.001f,0.003f,0.009f,0.02f, 0.02f, 0.017f,0.014f,0.012f,0.009f },
495
496// p n K+ K0bar (p n K0 K-)
497 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
498 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
499 0.0f, 0.002f,0.007f,0.029f,0.024f,0.02f, 0.017f,0.014f,0.012f,0.009f },
500//
501// multiplicity 5 (32 channels)
502//
503// p p pi+ pi- pi0 (n n pi+ pi- pi0)
504 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
505 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.06f,
506 0.4f, 1.1f, 1.8f, 2.4f, 2.4f, 2.2f, 2.0f, 1.7f, 1.5f, 1.3f },
507
508// p p pi0 pi0 pi0 (n n pi0 pi0 pi0)
509 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
510 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.003f,0.018f,
511 0.12f, 0.33f, 0.54f, 0.72f, 0.72f, 0.66f, 0.6f, 0.51f, 0.45f, 0.39f },
512
513// p n pi+ pi+ pi- (p n pi+ pi- pi-)
514 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
515 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.12f, 0.26f,
516 0.7f, 1.6f, 2.4f, 2.6f, 2.3f, 2.0f, 1.8f, 1.6f, 1.4f, 1.2f },
517
518// p n pi+ pi0 pi0 (p n pi- pi0 pi0)
519 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
520 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.036f,0.078f,
521 0.21f, 0.48f, 0.72f, 0.78f, 0.69f, 0.6f, 0.54f, 0.48f, 0.42f, 0.36f },
522
523// n n pi+ pi+ pi0 (p p pi- pi- pi0)
524 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
525 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.036f,
526 0.24f, 0.66f, 1.08f, 1.44f, 1.44f, 1.32f, 1.2f, 1.0f, 0.9f, 0.78f },
527
528// p L K+ pi+ pi- (n L K0 pi+ pi-)
529 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
530 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
531 0.0f, 0.0f, 0.002f,0.010f,0.040f,0.045f,0.040f,0.035f,0.030f,0.020f },
532
533// p L K+ pi0 pi0 (n L K0 pi0 pi0)
534 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
535 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
536 0.0f, 0.0f, 0.001f,0.005f,0.020f,0.022f,0.020f,0.017f,0.015f,0.010f },
537
538// p L K0 pi+ pi0 (n L K+ pi- pi0)
539 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
540 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
541 0.0f, 0.0f, 0.002f,0.010f,0.040f,0.045f,0.040f,0.035f,0.030f,0.020f },
542
543// p S0 K+ pi+ pi- (n S0 K0 pi+ pi-)
544 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
545 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
546 0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
547
548// p S0 K+ pi0 pi0 (n S0 K0 pi0 pi0)
549 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
550 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
551 0.0f, 0.0f, 0.001f,0.004f,0.015f,0.017f,0.015f,0.014f,0.011f,0.009f },
552
553// p S0 K0 pi+ pi0 (n S0 K+ pi- pi0)
554 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
555 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
556 0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
557
558// p S+ K0 pi+ pi- (n S- K+ pi+ pi-)
559 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
560 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
561 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
562
563// p S+ K0 pi0 pi0 (n S- K+ pi0 pi0)
564 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
565 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
566 0.0f, 0.0f, 0.001f,0.002f,0.009f,0.020f,0.013f,0.014f,0.011f,0.009f },
567
568// p S+ K+ pi- pi0 (n S- K0 pi+ pi0)
569 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
570 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
571 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
572
573// p S- K+ pi+ pi0 (n S+ K0 pi- pi0)
574 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
575 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
576 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
577
578// p S- K0 pi+ pi+ (n S+ K+ pi- pi-)
579 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
580 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
581 0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
582
583// n L K+ pi+ pi0 (p L K0 pi- pi0)
584 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
585 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
586 0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
587
588// n L K0 pi+ pi+ (p L K+ pi- pi-)
589 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
590 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
591 0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
592
593// n S0 K+ pi+ pi0 (p S0 K0 pi- pi0)
594 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
595 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
596 0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
597
598// n S0 K0 pi+ pi+ (p S0 K+ pi- pi-)
599 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
600 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
601 0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
602
603// n S+ K0 pi+ pi0 (p S- K+ pi- pi0)
604 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
605 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
606 0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
607
608// n S+ K+ pi+ pi- (p S- K0 pi+ pi-)
609 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
610 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
611 0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
612
613// n S+ K+ pi0 pi0 (p S- K0 pi0 pi0)
614 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
615 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
616 0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
617
618// n S- K+ pi+ pi+ (p S+ K0 pi- pi-)
619 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
620 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
621 0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
622
623// p p pi+ K0 K- (n n pi- K+ K0bar)
624 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
625 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
626 0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
627
628// p p pi- K+ K0bar (n n pi+ K0 K-)
629 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
630 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
631 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
632
633// p p pi0 K0 K0bar (n n pi0 K+ K-)
634 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
635 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
636 0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
637
638// p p pi0 K+ K- (n n pi0 K0 K0bar)
639 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
640 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
641 0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
642
643// p n pi+ K0 K0bar (p n pi- K+ K-)
644 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
645 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
646 0.0f, 0.0f, 0.002f,0.012f,0.035f,0.058f,0.050f,0.033f,0.023f,0.016f },
647
648// p n pi+ K+ K- (p n pi- K0 K0bar)
649 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
650 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
651 0.0f, 0.0f, 0.002f,0.012f,0.035f,0.058f,0.050f,0.033f,0.023f,0.016f },
652
653// p n pi0 K+ K0bar (p n pi0 K0 K-)
654 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
655 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
656 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
657
658// n n pi+ K+ K0bar (p p pi- K0 K-)
659 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
660 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
661 0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
662//
663// multiplicity 6 (7 channels)
664//
665// p p pi+ pi+ pi- pi- (n n pi+ pi+ pi- pi-)
666 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
667 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
668 0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
669
670// p p pi+ pi- pi0 pi0 (n n pi+ pi- pi0 pi0)
671 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
672 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
673 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
674
675// p p pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0)
676 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
677 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
678 0.01f, 0.02f, 0.05f, 0.1f, 0.13f, 0.12f, 0.11f, 0.1f, 0.1f, 0.09f },
679
680// p n pi+ pi+ pi- pi0 (p n pi+ pi- pi- pi0)
681 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
682 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
683 0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
684
685// p n pi+ pi0 pi0 pi0 (p n pi- pi0 pi0 pi0)
686 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
687 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
688 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
689
690// n n pi+ pi+ pi+ pi- (p p pi+ pi- pi- pi-)
691 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
692 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
693 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
694
695// n n pi+ pi+ pi0 pi0 (p p pi- pi- pi0 pi0)
696 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
697 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
698 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
699//
700// multiplicity 7 (8 channels)
701//
702// p p pi+ pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi- pi0)
703 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
704 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
705 0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
706
707// p p pi+ pi- pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0)
708 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
709 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
710 0.0f, 0.0f, 0.036f,0.096f,0.30f, 0.42f, 0.42f, 0.42f, 0.40f, 0.37f },
711
712// p p pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0)
713 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
714 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
715 0.0f, 0.0f, 0.018f,0.048f,0.14f, 0.20f, 0.22f, 0.20f, 0.19f, 0.18f },
716
717// p n pi+ pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi- pi-)
718 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
719 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
720 0.0f, 0.0f, 0.06f, 0.19f, 0.31f, 0.41f, 0.44f, 0.47f, 0.45f, 0.45f },
721
722// p n pi+ pi+ pi- pi0 pi0 (p n pi+ pi- pi- pi0 pi0)
723 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
724 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
725 0.0f, 0.0f, 0.036f,0.12f, 0.18f, 0.24f, 0.26f, 0.23f, 0.28f, 0.26f },
726
727// p n pi+ pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0)
728 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
729 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
730 0.0f, 0.0f, 0.024f,0.06f, 0.08f, 0.12f, 0.13f, 0.14f, 0.13f, 0.13f },
731
732// n n pi+ pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi- pi0)
733 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
734 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
735 0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
736
737// n n pi+ pi+ pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0)
738 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
739 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
740 0.0f, 0.0f, 0.036f,0.096f,0.30f, 0.42f, 0.42f, 0.41f, 0.40f, 0.37f },
741//
742// multiplicity 8 (10 channels)
743//
744// p p pi+ pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi- pi-)
745 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
746 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
747 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f,0.27f, 0.30f, 0.27f, 0.24f },
748
749// p p pi+ pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0)
750 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
751 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
752 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f,0.27f, 0.30f, 0.27f, 0.24f },
753
754// p p pi+ pi- pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0)
755 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
756 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
757 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
758
759// p p pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0)
760 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
761 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
762 0.0f, 0.0f, 0.003f,0.01f, 0.024f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
763
764// p n pi+ pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi- pi0)
765 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
766 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
767 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
768
769// p n pi+ pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0)
770 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
771 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
772 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
773
774// p n pi+ pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0)
775 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
776 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
777 0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
778
779// n n pi+ pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi- pi-)
780 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
781 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
782 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
783
784// n n pi+ pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0)
785 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
786 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
787 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
788
789// n n pi+ pi+ pi0 pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0 pi0)
790 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
791 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
792 0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
793//
794// multiplicity 9 (11 channels)
795//
796// p p pi+ pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi- pi0)
797 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
798 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
799 0.0f, 0.002f,0.008f,0.025f,0.074f,0.11f, 0.14f, 0.15f, 0.15f, 0.15f },
800
801// p p pi+ pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0 pi0)
802 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
803 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
804 0.0f, 0.001f,0.006f,0.020f,0.058f,0.091f,0.11f, 0.12f, 0.12f, 0.12f },
805
806// p p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0 pi0)
807 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
808 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
809 0.0f, 0.001f,0.004f,0.012f,0.035f,0.055f,0.065f,0.07f, 0.07f, 0.07f },
810
811// p p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
812 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
813 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
814 0.0f, 0.0f, 0.0f, 0.006f,0.027f,0.032f,0.04f, 0.042f,0.042f,0.042f },
815
816// p n pi+ pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi- pi-)
817 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
818 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
819 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
820
821// p n pi+ pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi- pi0 pi0)
822 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
823 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
824 0.0f, 0.003f,0.008f,0.026f,0.078f,0.20f, 0.25f, 0.29f, 0.29f, 0.29f },
825
826// p n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0 pi0)
827 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
828 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
829 0.0f, 0.001f,0.005f,0.016f,0.047f,0.12f, 0.15f, 0.17f, 0.17f, 0.17f },
830
831// p n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0 pi0)
832 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
833 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
834 0.0f, 0.001f,0.003f,0.009f,0.029f,0.07f, 0.094f,0.10f, 0.10f, 0.10f },
835
836// n n pi+ pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi- pi0)
837 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
838 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
839 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
840
841// n n pi+ pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0 pi0)
842 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
843 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
844 0.0f, 0.001f,0.005f,0.015f,0.047f,0.12f, 0.15f, 0.17f, 0.17f, 0.17f },
845
846// n n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p p pi+ pi- pi0 pi0 pi0 pi0 pi0)
847 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
848 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
849 0.0f, 0.001f,0.003f,0.009f,0.029f,0.07f, 0.094f,0.10f, 0.10f, 0.10f }};
850
851///////////////////////////////////////////////////////////////////////////////
852// //
853// p n and n p (|Tz| = 0) cross sections //
854// and final state particle types //
855// //
856///////////////////////////////////////////////////////////////////////////////
857
858// Total p n cross section as a function of kinetic energy
859G4ThreadLocal G4double G4RPGNucleonInelastic::pNtot[30];
860
861// p n multiplicities as a function of kinetic energy
862G4ThreadLocal G4double G4RPGNucleonInelastic::t0_dSigma_dMult[8][30];
863
864const G4int G4RPGNucleonInelastic::pNindex[8][2] =
865 {{0, 0}, {1,9}, {10,31}, {32,69}, {70,76}, {77,85}, {86,95}, {96,107}};
866// first index: 0: channels for mult = 2
867// 1-9: channels for mult = 3
868// 10-31: channels for mult = 4
869// 32-69: channels for mult = 5
870// 70-76: channels for mult = 6
871// 77-85: channels for mult = 7
872// 86-95: channels for mult = 8
873// 96-107: channels for mult = 9
874
875// Outgoing particle types of a given multiplicity
876// T0_nbfs = final state types for p n and n p
877
878const G4int G4RPGNucleonInelastic::T0_2bfs[1][2] =
879 {{pro,neu}};
880
881const G4int G4RPGNucleonInelastic::T0_3bfs[9][3] =
882 {{pro,pro,pim},{pro,neu,pi0},{neu,neu,pip},{pro,lam,k0},
883 {pro,s0,k0}, {pro,sm,kp}, {neu,lam,kp}, {neu,s0,kp},
884 {neu,sp,k0}};
885
886const G4int G4RPGNucleonInelastic::T0_4bfs[22][4] =
887 {{pro,neu,pip,pim},{pro,pro,pim,pi0},{pro,neu,pi0,pi0},
888 {neu,neu,pip,pi0},{pro,lam,kp,pim}, {pro,s0,kp,pim},
889 {pro,lam,k0,pi0}, {pro,s0,k0,pi0}, {pro,sp,k0,pim},
890 {pro,sm,kp,pi0}, {pro,sm,k0,pip}, {neu,lam,kp,pi0},
891 {neu,lam,k0,pip}, {neu,sp,kp,pim}, {neu,sp,k0,pi0},
892 {neu,s0,kp,pi0}, {neu,s0,k0,pip}, {neu,sm,kp,pip},
893 {pro,neu,kp,km}, {pro,neu,k0,k0b}, {pro,pro,k0,km},
894 {neu,neu,kp,k0b}};
895
896const G4int G4RPGNucleonInelastic::T0_5bfs[38][5] =
897 {{pro,neu,pip,pim,pi0},{pro,neu,pi0,pi0,pi0},{pro,pro,pip,pim,pim},
898 {pro,pro,pim,pi0,pi0},{neu,neu,pip,pip,pim},{neu,neu,pip,pi0,pi0},
899 {pro,lam,kp,pim,pi0}, {pro,lam,k0,pip,pim}, {pro,lam,k0,pi0,pi0},
900 {pro,s0,k0,pip,pim}, {pro,s0,k0,pi0,pi0}, {pro,s0,kp,pim,pi0},
901 {pro,sp,kp,pim,pim}, {pro,sp,k0,pim,pi0}, {pro,sm,k0,pip,pi0},
902 {pro,sm,kp,pip,pim}, {pro,sm,kp,pi0,pi0}, {neu,lam,kp,pip,pim},
903 {neu,lam,kp,pi0,pi0}, {neu,lam,k0,pip,pi0}, {neu,s0,kp,pip,pim},
904 {neu,s0,kp,pi0,pi0}, {neu,s0,k0,pip,pi0}, {neu,sp,k0,pip,pim},
905 {neu,sp,k0,pi0,pi0}, {neu,sp,kp,pim,pi0}, {neu,sm,kp,pip,pi0},
906 {neu,sm,k0,pip,pip}, {pro,neu,kp,km,pi0}, {pro,neu,k0,k0b,pi0},
907 {pro,neu,k0,km,pip}, {pro,neu,kp,k0b,pim}, {pro,pro,k0,k0b,pim},
908 {pro,pro,kp,km,pim}, {pro,pro,k0,km,pi0}, {neu,neu,kp,km,pip},
909 {neu,neu,k0,k0b,pip}, {neu,neu,kp,k0b,pi0}};
910
911const G4int G4RPGNucleonInelastic::T0_6bfs[7][6] =
912 {{pro,neu,pip,pip,pim,pim},{pro,neu,pip,pim,pi0,pi0},
913 {pro,neu,pi0,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pi0},
914 {pro,pro,pim,pi0,pi0,pi0},{neu,neu,pip,pip,pim,pi0},
915 {neu,neu,pip,pi0,pi0,pi0}};
916
917const G4int G4RPGNucleonInelastic::T0_7bfs[9][7] =
918 {{pro,neu,pip,pip,pim,pim,pi0},{pro,neu,pip,pim,pi0,pi0,pi0},
919 {pro,neu,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim},
920 {pro,pro,pip,pim,pim,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0},
921 {neu,neu,pip,pip,pip,pim,pim},{neu,neu,pip,pip,pim,pi0,pi0},
922 {neu,neu,pip,pi0,pi0,pi0,pi0}};
923
924const G4int G4RPGNucleonInelastic::T0_8bfs[10][8] =
925{{pro,neu,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pi0,pi0},
926 {pro,neu,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0},
927 {pro,pro,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pi0,pi0,pi0},
928 {pro,pro,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim,pim,pi0},
929 {neu,neu,pip,pip,pim,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0}};
930
931const G4int G4RPGNucleonInelastic::T0_9bfs[12][9] =
932{{pro,neu,pip,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pip,pim,pim,pi0,pi0,pi0},
933 {pro,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
934 {pro,pro,pip,pip,pip,pim,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pim,pi0,pi0},
935 {pro,pro,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
936 {neu,neu,pip,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pip,pim,pim,pi0,pi0},
937 {neu,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}};
938
939//
940// Cross sections (in mb) for p n -> 2-9 body final states
941//
942// first index: 0: channels for mult = 2
943// 1-9: channels for mult = 3
944// 10-31: channels for mult = 4
945// 32-69: channels for mult = 5
946// 70-76: channels for mult = 6
947// 77-85: channels for mult = 7
948// 86-95: channels for mult = 8
949// 96-107: channels for mult = 9
950//
951// second index: kinetic energy
952//
953const G4float G4RPGNucleonInelastic::pNCrossSections[108][30] = {
954//
955// multiplicity 2 (1 channel)
956//
957// p n (p n)
958 { 0.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f,
959 44.0f, 42.0f, 40.0f, 35.0f, 31.0f, 27.0f, 23.0f, 19.0f, 17.0f, 15.5f,
960 14.0f, 13.0f, 12.0f, 11.0f, 10.0f, 9.5f, 9.0f, 8.5f, 8.0f, 7.7f },
961//
962// multiplicity 3 (9 channels)
963//
964// p p pi- (n n pi+)
965 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
966 0.0f, 0.0f, 0.0f, 0.0f, 0.25f, 0.9f, 1.75f, 2.3f, 2.8f, 2.8f,
967 2.2f, 1.9f, 1.6f, 1.35f, 1.1f, 0.95f, 0.8f, 0.7f, 0.6f, 0.53f },
968
969// p n pi0 (p n pi0)
970 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
971 0.0f, 0.0f, 0.0f, 0.0f, 1.8f, 4.7f, 8.3f, 11.3f, 12.0f, 10.2f,
972 8.2f, 6.0f, 4.9f, 3.6f, 2.5f, 2.0f, 1.6f, 1.2f, 1.0f, 0.08f },
973
974// n n pi+ (p p pi-)
975 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
976 0.0f, 0.0f, 0.0f, 0.0f, 0.95f, 2.4f, 4.2f, 5.6f, 6.1f, 5.1f,
977 4.1f, 3.0f, 2.5f, 1.8f, 1.2f, 1.0f, 0.8f, 0.6f, 0.5f, 0.41f},
978
979// p L K0 (n L K+)
980 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
981 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
982 0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
983
984// p S0 K0 (n S0 K+)
985 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
986 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
987 0.003f,0.010f,0.016f,0.020f,0.016f,0.015f,0.014f,0.013f,0.011f,0.010f},
988
989// p S- K+ (n S+ K0)
990 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
991 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
992 0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
993
994// n L K+ (p L K0)
995 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
996 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
997 0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
998
999// n S0 K+ (p S0 K0)
1000 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1001 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1002 0.003f,0.010f,0.016f,0.020f,0.016f,0.015f,0.014f,0.013f,0.011f,0.010f},
1003
1004// n S+ K0 (p S- K+)
1005 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1006 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1007 0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
1008//
1009// multiplicity 4 (22 channels)
1010//
1011// p n pi+ pi- (p n pi+ pi-)
1012 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1013 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.12f, 0.38f, 1.1f, 3.5f,
1014 5.9f, 5.9f, 5.1f, 4.2f, 3.7f, 3.0f, 2.6f, 2.1f, 1.8f, 1.4f },
1015
1016// p p pi- pi0 (n n pi+ pi0)
1017 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1018 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.03f, 0.1f, 0.24f, 0.55f,
1019 1.2f, 1.5f, 1.45f, 1.25f, 1.0f, 0.9f, 0.8f, 0.7f, 0.6f, 0.53f },
1020
1021// p n pi0 pi0 (p n pi0 pi0)
1022 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1023 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.07f, 0.24f, 0.66f, 2.1f,
1024 3.6f, 3.6f, 3.1f, 2.5f, 2.2f, 1.8f, 1.5f, 1.2f, 1.1f, 0.84f },
1025
1026// n n pi+ pi0 (p p pi- pi0)
1027 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1028 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.03f, 0.1f, 0.24f, 0.55f,
1029 1.2f, 1.5f, 1.45f, 1.25f, 1.0f, 0.9f, 0.8f, 0.7f, 0.6f, 0.53f },
1030
1031// p L K+ pi- (n L K0 pi+)
1032 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1033 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1034 0.0f, 0.002f,0.006f,0.018f,0.037f,0.036f,0.033f,0.030f,0.028f,0.023f },
1035
1036// p S0 K+ pi- (n S0 K0 pi+)
1037 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1038 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1039 0.0f, 0.0f, 0.002f,0.007f,0.021f,0.025f,0.022f,0.020f,0.018f,0.017f },
1040
1041// p L K0 pi0 (n L K+ pi0)
1042 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1043 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1044 0.0f, 0.001f,0.002f,0.005f,0.009f,0.009f,0.008f,0.007f,0.007f,0.006f },
1045
1046// p S0 K0 pi0 (n S0 K+ pi0)
1047 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1048 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1049 0.0f, 0.0f, 0.001f,0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.004f },
1050
1051// p S+ K0 pi- (n S- K+ pi+)
1052 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1053 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1054 0.0f, 0.0015f,0.004f,0.013f,0.02f,0.016f,0.013f,0.01f, 0.009f,0.007f },
1055
1056// p S- K+ pi0 (n S+ K0 pi0)
1057 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1058 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1059 0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1060
1061// p S- K0 pi+ (n S+ K+ pi-)
1062 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1063 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1064 0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1065
1066// n L K+ pi0 (p L K0 pi0)
1067 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1068 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1069 0.0f, 0.002f,0.006f,0.018f,0.037f,0.036f,0.033f,0.030f,0.028f,0.023f },
1070
1071// n L K0 pi+ (p L K+ pi-)
1072 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1073 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1074 0.0f, 0.001f,0.002f,0.005f,0.009f,0.009f,0.008f,0.007f,0.007f,0.006f },
1075
1076// n S+ K+ pi- (p S- K0 pi+)
1077 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1078 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1079 0.0f, 0.0f, 0.002f,0.007f,0.021f,0.025f,0.022f,0.020f,0.018f,0.017f },
1080
1081// n S+ K0 pi0 (p S- K+ pi0)
1082 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1083 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1084 0.0f, 0.0f, 0.001f,0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.004f },
1085
1086// n S0 K+ pi0 (p S0 K0 pi0)
1087 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1088 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1089 0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1090
1091// n S0 K0 pi+ (p S0 K+ pi-)
1092 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1093 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1094 0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1095
1096// n S- K+ pi+ (p S+ K0 pi-)
1097 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1098 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1099 0.0f, 0.0015f,0.004f,0.013f,0.02f, 0.016f,0.013f,0.01f, 0.009f,0.007f },
1100
1101// p n K+ K- (p n K0 K0bar)
1102 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1103 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1104 0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1105
1106// p n K0 K0bar (p n K+ K-)
1107 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1108 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1109 0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1110
1111// p p K0 K- (n n K+ K0bar)
1112 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1113 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1114 0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1115
1116// n n K+ K0bar (p p K0 K-)
1117 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1118 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1119 0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1120//
1121// multiplicity 5 (38 channels)
1122//
1123// p n pi+ pi- pi0 (p n pi+ pi- pi0)
1124 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1125 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.04f,
1126 0.3f, 0.82f, 1.35f, 1.8f, 1.8f, 1.65f, 1.5f, 1.28f, 1.12f, 0.98f },
1127
1128// p n pi0 pi0 pi0 (p n pi0 pi0 pi0)
1129 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1130 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.004f,0.022f,
1131 0.15f, 0.41f, 0.68f, 0.9f, 0.9f, 0.82f, 0.75f, 0.64f, 0.55f, 0.49f },
1132
1133// p p pi+ pi- pi- (n n pi+ pi+ pi-)
1134 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1135 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.09f, 0.2f,
1136 0.52f, 1.2f, 1.8f, 2.0f, 1.7f, 1.5f, 1.35f, 1.2f, 1.05f, 0.9f },
1137
1138// p p pi- pi0 pi0 (n n pi+ pi0 pi0)
1139 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1140 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.04f, 0.1f,
1141 0.26f, 0.6f, 0.9f, 0.98f, 0.86f, 0.75f, 0.68f, 0.6f, 0.52f, 0.45f },
1142
1143// n n pi+ pi+ pi- (p p pi+ pi- pi-)
1144 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1145 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.04f,
1146 0.3f, 0.82f, 1.35f, 1.8f, 1.8f, 1.65f, 1.5f, 1.28f, 1.12f, 0.98f },
1147
1148// n n pi+ pi0 pi0 (p p pi- pi0 pi0)
1149 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1150 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.003f,0.02f,
1151 0.15f, 0.41f, 0.68f, 0.9f, 0.9f, 0.82f, 0.75f, 0.64f, 0.56f, 0.49f },
1152
1153// p L K+ pi- pi0 (n L K0 pi+ pi0)
1154 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1155 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1156 0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1157
1158// p L K0 pi+ pi- (n L K+ pi+ pi-)
1159 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1160 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1161 0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1162
1163// p L K0 pi0 pi0 (n L K+ pi0 pi0)
1164 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1165 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1166 0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.09f, 0.07f, 0.07f },
1167
1168// p S0 K0 pi+ pi- (n S0 K+ pi+ pi-)
1169 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1170 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1171 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1172
1173// p S0 K0 pi0 pi0 (n S0 K+ pi0 pi0)
1174 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1175 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1176 0.0f, 0.0f, 0.001f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1177
1178// p S0 K+ pi- pi0 (n S0 K0 pi+ pi0)
1179 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1180 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1181 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1182
1183// p S+ K+ pi- pi- (n S- K0 pi+ pi+)
1184 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1185 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1186 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1187
1188// p S+ K0 pi- pi0 (n S- K+ pi+ pi0)
1189 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1190 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1191 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1192
1193// p S- K0 pi+ pi0 (n S+ K+ pi- pi0)
1194 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1195 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1196 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1197
1198// p S- K+ pi+ pi- (n S+ K0 pi+ pi-)
1199 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1200 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1201 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1202
1203// p S- K+ pi0 pi0 (n S+ K0 pi0 pi0)
1204 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1205 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1206 0.0f, 0.0f, 0.002f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1207
1208// n L K+ pi+ pi- (p L K0 pi+ pi-)
1209 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1210 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1211 0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1212
1213// n L K+ pi0 pi0 (p L K0 pi0 pi0)
1214 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1215 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1216 0.0f, 0.0f, 0.001f,0.004f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1217
1218// n L K0 pi+ pi0 (p L K+ pi- pi0)
1219 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1220 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1221 0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1222
1223// n S0 K+ pi+ pi- (p S0 K0 pi+ pi-)
1224 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1225 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1226 0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1227
1228// n S0 K+ pi0 pi0 (p S0 K0 pi0 pi0)
1229 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1230 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1231 0.0f, 0.0f, 0.0f, 0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.003f },
1232
1233// n S0 K0 pi+ pi0 (p S0 K+ pi- pi0)
1234 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1235 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1236 0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1237
1238// n S+ K0 pi+ pi- (p S- K+ pi+ pi-)
1239 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1240 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1241 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1242
1243// n S+ K0 pi0 pi0 (p S- K+ pi0 pi0)
1244 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1245 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1246 0.0f, 0.0f, 0.001f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1247
1248// n S+ K+ pi- pi0 (p S- K0 pi+ pi0)
1249 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1250 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1251 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1252
1253// n S- K+ pi+ pi0 (p S+ K0 pi- pi0)
1254 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1255 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1256 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1257
1258// n S- K0 pi+ pi+ (p S+ K+ pi- pi-)
1259 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1260 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1261 0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1262
1263// p n K+ K- pi0 (p n K0 K0bar pi0)
1264 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1265 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1266 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1267
1268// p n K0 K0bar pi0 (p n K+ K- pi0)
1269 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1270 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1271 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1272
1273// p n K0 K- pi+ (p n K+ K0bar pi-)
1274 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1275 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1276 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1277
1278// p n K+ K0bar pi- (p n K0 K- pi+)
1279 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1280 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1281 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1282
1283// p p K0 K0bar pi- (n n K+ K- pi+)
1284 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1285 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1286 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1287
1288// p p K+ K- pi- (n n K0 K0bar pi+)
1289 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1290 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1291 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1292
1293// p p K0 K- pi0 (n n K+ K0bar pi0)
1294 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1295 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1296 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1297
1298// n n K+ K- pi+ (p p K0 K0bar pi-)
1299 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1300 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1301 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1302
1303// n n K0 K0bar pi+ (p p K+ K- pi-)
1304 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1305 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1306 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1307
1308// n n K+ K0bar pi0 (p p K0 K- pi0)
1309 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1310 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1311 0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1312//
1313// multiplicity 6 (7 channels)
1314//
1315// p n pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi-)
1316 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1317 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1318 0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1319
1320// p n pi+ pi- pi0 pi0 (p n pi+ pi- pi0 pi0)
1321 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1322 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1323 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1324
1325// p n pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0)
1326 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1327 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1328 0.01f, 0.02f, 0.05f, 0.1f, 0.13f, 0.12f, 0.11f, 0.1f, 0.1f, 0.09f },
1329
1330// p p pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi0)
1331 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1332 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1333 0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1334
1335// p p pi- pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0)
1336 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1337 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1338 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1339
1340// n n pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi0)
1341 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1342 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1343 0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1344
1345// n n pi+ pi0 pi0 pi0 (p p pi- pi0 pi0 pi0)
1346 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1347 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1348 0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1349//
1350// multiplicity 7 (9 channels)
1351//
1352// p n pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi0)
1353 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1354 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1355 0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
1356
1357// p n pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0)
1358 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1359 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1360 0.0f, 0.0f, 0.03f, 0.08f, 0.25f, 0.35f, 0.35f, 0.35f, 0.33f, 0.31f },
1361
1362// p n pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0)
1363 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1364 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1365 0.0f, 0.0f, 0.015f,0.04f, 0.12f, 0.17f, 0.18f, 0.17f, 0.16f, 0.15f },
1366
1367// p p pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi-)
1368 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1369 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1370 0.0f, 0.0f, 0.06f, 0.19f, 0.31f, 0.41f, 0.44f, 0.47f, 0.45f, 0.45f },
1371
1372// p p pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi0 pi0)
1373 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1374 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1375 0.0f, 0.0f, 0.03f, 0.1f, 0.15f, 0.2f, 0.22f, 0.23f, 0.22f, 0.22f },
1376
1377// p p pi- pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0)
1378 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1379 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1380 0.0f, 0.0f, 0.02f, 0.05f, 0.07f, 0.1f, 0.11f, 0.12f, 0.11f, 0.11f },
1381
1382// n n pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi-)
1383 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1384 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1385 0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
1386
1387// n n pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi0 pi0)
1388 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1389 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1390 0.0f, 0.0f, 0.03f, 0.08f, 0.25f, 0.35f, 0.35f, 0.34f, 0.33f, 0.31f },
1391
1392// n n pi+ pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0)
1393 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1394 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1395 0.0f, 0.0f, 0.02f, 0.05f, 0.07f, 0.1f, 0.11f, 0.12f, 0.11f, 0.11f },
1396//
1397// multiplicity 8 (10 channels)
1398//
1399// p n pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi-)
1400 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1401 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1402 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1403
1404// p n pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0)
1405 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1406 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1407 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1408
1409// p n pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0)
1410 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1411 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1412 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1413
1414// p n pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0)
1415 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1416 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1417 0.0f, 0.0f, 0.003f,0.01f, 0.024f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
1418
1419// p p pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi0)
1420 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1421 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1422 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1423
1424// p p pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0)
1425 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1426 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1427 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1428
1429// p p pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0)
1430 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1431 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1432 0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
1433
1434// n n pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi0)
1435 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1436 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1437 0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1438
1439// n n pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0)
1440 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1441 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1442 0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1443
1444// n n pi+ pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0)
1445 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1446 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1447 0.0f, 0.0f, 0.003f,0.009f,0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
1448//
1449// multiplicity 9 (12 channels)
1450//
1451// p n pi+ pi+ pi+ pi- pi- pi- pi0 (p n pi+ pi+ pi+ pi- pi- pi- pi0)
1452 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1453 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1454 0.0f, 0.002f,0.008f,0.025f,0.074f,0.11f, 0.14f, 0.15f, 0.15f, 0.15f },
1455
1456// p n pi+ pi+ pi- pi- pi0 pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0 pi0)
1457 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1458 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1459 0.0f, 0.001f,0.005f,0.015f,0.045f,0.07f, 0.084f,0.09f, 0.09f, 0.09f },
1460
1461// p n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0 pi0)
1462 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1463 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1464 0.0f, 0.001f,0.003f,0.009f,0.027f,0.042f,0.05f, 0.054f,0.054f,0.054f },
1465
1466// p n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
1467 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1468 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1469 0.0f, 0.0f, 0.0f, 0.005f,0.016f,0.025f,0.03f, 0.032f,0.032f,0.032f },
1470
1471// p p pi+ pi+ pi+ pi- pi- pi- pi- (n n pi+ pi+ pi+ pi+ pi- pi- pi-)
1472 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1473 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1474 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1475
1476// p p pi+ pi+ pi- pi- pi- pi0 pi0 (n n pi+ pi+ pi+ pi- pi- pi0 pi0)
1477 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1478 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1479 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1480
1481// p p pi+ pi- pi- pi0 pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0 pi0)
1482 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1483 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1484 0.0f, 0.001f,0.004f,0.012f,0.036f,0.09f, 0.12f, 0.13f, 0.13f, 0.13f },
1485
1486// p p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
1487 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1488 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1489 0.0f, 0.001f,0.002f,0.007f,0.022f,0.054f,0.072f,0.078f,0.078f,0.078f },
1490
1491// n n pi+ pi+ pi+ pi+ pi- pi- pi- (p p pi+ pi+ pi+ pi- pi- pi- pi-)
1492 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1493 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1494 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1495
1496// n n pi+ pi+ pi+ pi- pi- pi0 pi0 (p p pi+ pi+ pi- pi- pi- pi0 pi0)
1497 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1498 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1499 0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1500
1501// n n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0 pi0)
1502 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1503 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1504 0.0f, 0.001f,0.004f,0.012f,0.036f,0.09f, 0.12f, 0.13f, 0.13f, 0.13f },
1505
1506// n n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0 pi0)
1507 { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1508 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1509 0.0f, 0.001f,0.002f,0.007f,0.022f,0.054f,0.072f,0.078f,0.078f,0.078f }};
G4float
float G4float
Definition: G4Types.hh:84
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:89
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:102
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGNucleonInelastic::GetMultiplicityT0
G4int GetMultiplicityT0(G4double KE) const
Definition: G4RPGNucleonInelastic.cc:98
G4RPGNucleonInelastic::T0_2bfs
static const G4int T0_2bfs[1][2]
Definition: G4RPGNucleonInelastic.hh:86
G4RPGNucleonInelastic::GetMultiplicityT1
G4int GetMultiplicityT1(G4double KE) const
Definition: G4RPGNucleonInelastic.cc:117
G4RPGNucleonInelastic::G4RPGNucleonInelastic
G4RPGNucleonInelastic(const G4String &modelName="RPGNucleonInelastic")
Definition: G4RPGNucleonInelastic.cc:32
G4RPGNucleonInelastic::T0_8bfs
static const G4int T0_8bfs[10][8]
Definition: G4RPGNucleonInelastic.hh:92
G4RPGNucleonInelastic::t1_dSigma_dMult
static G4ThreadLocal G4double t1_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:97
G4RPGNucleonInelastic::T0_9bfs
static const G4int T0_9bfs[12][9]
Definition: G4RPGNucleonInelastic.hh:93
G4RPGNucleonInelastic::T1_3bfs
static const G4int T1_3bfs[2][6][3]
Definition: G4RPGNucleonInelastic.hh:78
G4RPGNucleonInelastic::T0_4bfs
static const G4int T0_4bfs[22][4]
Definition: G4RPGNucleonInelastic.hh:88
G4RPGNucleonInelastic::T1_2bfs
static const G4int T1_2bfs[2][1][2]
Definition: G4RPGNucleonInelastic.hh:77
G4RPGNucleonInelastic::T0_5bfs
static const G4int T0_5bfs[38][5]
Definition: G4RPGNucleonInelastic.hh:89
G4RPGNucleonInelastic::T1_4bfs
static const G4int T1_4bfs[2][18][4]
Definition: G4RPGNucleonInelastic.hh:79
G4RPGNucleonInelastic::pPtot
static G4ThreadLocal G4double pPtot[30]
Definition: G4RPGNucleonInelastic.hh:95
G4RPGNucleonInelastic::GetFSPartTypesForT0
std::vector< G4int > GetFSPartTypesForT0(G4int mult, G4double KE) const
Definition: G4RPGNucleonInelastic.cc:137
G4RPGNucleonInelastic::T0_6bfs
static const G4int T0_6bfs[7][6]
Definition: G4RPGNucleonInelastic.hh:90
G4RPGNucleonInelastic::T1_5bfs
static const G4int T1_5bfs[2][32][5]
Definition: G4RPGNucleonInelastic.hh:80
G4RPGNucleonInelastic::T1_9bfs
static const G4int T1_9bfs[2][11][9]
Definition: G4RPGNucleonInelastic.hh:84
G4RPGNucleonInelastic::T1_7bfs
static const G4int T1_7bfs[2][8][7]
Definition: G4RPGNucleonInelastic.hh:82
G4RPGNucleonInelastic::T1_6bfs
static const G4int T1_6bfs[2][7][6]
Definition: G4RPGNucleonInelastic.hh:81
G4RPGNucleonInelastic::pPCrossSections
static const G4float pPCrossSections[93][30]
Definition: G4RPGNucleonInelastic.hh:100
G4RPGNucleonInelastic::t0_dSigma_dMult
static G4ThreadLocal G4double t0_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:98
G4RPGNucleonInelastic::pPindex
static const G4int pPindex[8][2]
Definition: G4RPGNucleonInelastic.hh:74
G4RPGNucleonInelastic::GetFSPartTypesForT1
std::vector< G4int > GetFSPartTypesForT1(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGNucleonInelastic.cc:185
G4RPGNucleonInelastic::T0_7bfs
static const G4int T0_7bfs[9][7]
Definition: G4RPGNucleonInelastic.hh:91
G4RPGNucleonInelastic::pNindex
static const G4int pNindex[8][2]
Definition: G4RPGNucleonInelastic.hh:75
G4RPGNucleonInelastic::T0_3bfs
static const G4int T0_3bfs[9][3]
Definition: G4RPGNucleonInelastic.hh:87
G4RPGNucleonInelastic::pNtot
static G4ThreadLocal G4double pNtot[30]
Definition: G4RPGNucleonInelastic.hh:96
G4RPGNucleonInelastic::pNCrossSections
static const G4float pNCrossSections[108][30]
Definition: G4RPGNucleonInelastic.hh:101
G4RPGNucleonInelastic::T1_8bfs
static const G4int T1_8bfs[2][10][8]
Definition: G4RPGNucleonInelastic.hh:83
G4ThreadLocal
#define G4ThreadLocal
Definition: tls.hh:77