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