ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4RPGNucleonInelastic.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4RPGNucleonInelastic.cc
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 
32 G4RPGNucleonInelastic::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 /*
76 void 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 
98 G4int 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 
117 G4int 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 
136 std::vector<G4int>
137 G4RPGNucleonInelastic::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 
184 std::vector<G4int>
185 G4RPGNucleonInelastic::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 
233 // //
234 // p p and n n (|Tz| = 1) cross sections //
235 // and final state particle types //
236 // //
238 
239 // Total p p cross section as a function of kinetic energy
240 G4ThreadLocal G4double G4RPGNucleonInelastic::pPtot[30];
241 
242 // p p multiplicities as a function of kinetic energy
243 G4ThreadLocal G4double G4RPGNucleonInelastic::t1_dSigma_dMult[8][30];
244 
245 const 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 
251 const G4int G4RPGNucleonInelastic::T1_2bfs[2][1][2] =
252  {{{pro,pro}},
253 
254  {{neu,neu}}};
255 
256 const 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 
263 const 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 
278 const 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 
303 const 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 
314 const 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 
325 const 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 
338 const 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 
368 const 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 
852 // //
853 // p n and n p (|Tz| = 0) cross sections //
854 // and final state particle types //
855 // //
857 
858 // Total p n cross section as a function of kinetic energy
859 G4ThreadLocal G4double G4RPGNucleonInelastic::pNtot[30];
860 
861 // p n multiplicities as a function of kinetic energy
862 G4ThreadLocal G4double G4RPGNucleonInelastic::t0_dSigma_dMult[8][30];
863 
864 const 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 
878 const G4int G4RPGNucleonInelastic::T0_2bfs[1][2] =
879  {{pro,neu}};
880 
881 const 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 
886 const 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 
896 const 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 
911 const 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 
917 const 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 
924 const 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 
931 const 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 //
953 const 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 }};