Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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: G4RPGNucleonInelastic.cc 94566 2015-11-24 10:25:06Z gcosmo $
27 //
28 
29 #include "G4RPGNucleonInelastic.hh"
30 #include "G4SystemOfUnits.hh"
31 #include "Randomize.hh"
32 
33 G4RPGNucleonInelastic::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 /*
77 void 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 
99 G4int 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 
118 G4int 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 
137 std::vector<G4int>
138 G4RPGNucleonInelastic::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 
185 std::vector<G4int>
186 G4RPGNucleonInelastic::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 
234 // //
235 // p p and n n (|Tz| = 1) cross sections //
236 // and final state particle types //
237 // //
239 
240 // Total p p cross section as a function of kinetic energy
241 G4ThreadLocal G4double G4RPGNucleonInelastic::pPtot[30];
242 
243 // p p multiplicities as a function of kinetic energy
244 G4ThreadLocal G4double G4RPGNucleonInelastic::t1_dSigma_dMult[8][30];
245 
246 const 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 
252 const G4int G4RPGNucleonInelastic::T1_2bfs[2][1][2] =
253  {{{pro,pro}},
254 
255  {{neu,neu}}};
256 
257 const 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 
264 const 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 
279 const 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 
304 const 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 
315 const 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 
326 const 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 
339 const 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 
369 const G4float G4RPGNucleonInelastic::pPCrossSections[93][30] = {
370 //
371 // multiplicity 2 (1 channel)
372 //
373 // p p (n n)
374  { 0.0f,330.0f,240.0f,160.0f,110.0f, 85.0f, 63.0f, 44.0f, 33.0f, 28.0f,
375  25.0f, 24.0f, 23.0f, 23.0f, 26.3f, 26.1f, 25.0f, 23.5f, 21.0f, 18.0f,
376  16.0f, 14.3f, 12.5f, 11.2f, 10.3f, 9.6f, 9.0f, 8.5f, 8.0f, 7.7f },
377 //
378 // multiplicity 3 (6 channels)
379 //
380 // p p pi0 (n n pi0)
381  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
382  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.4f, 4.0f, 4.3f, 4.0f, 4.0f,
383  3.6f, 3.0f, 2.8f, 2.5f, 1.7f, 1.3f, 1.1f, 1.0f, 0.9f, 0.85f },
384 
385 // p n pi+ (p n pi-)
386  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
387  0.0f, 0.0f, 0.0f, 0.0f, 0.7f, 4.5f, 15.0f, 19.1f, 18.0f, 16.0f,
388  13.0f, 10.0f, 8.2f, 6.0f, 4.3f, 3.3f, 2.6f, 2.0f, 1.65f, 1.4f },
389 
390 // p L K+ (n L K0)
391  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
392  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.012f,
393  0.03f, 0.06f, 0.06f, 0.055f,0.05f, 0.047f,0.043f,0.04f, 0.037f,0.033f },
394 
395 // p S0 K+ (n S0 K0)
396  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
397  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
398  0.006f,0.02f, 0.027f,0.026f,0.021f,0.018f,0.015f,0.011f,0.009f,0.007f },
399 
400 // p S+ K0 (n S- K+)
401  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
402  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
403  0.013f,0.025f,0.03f, 0.029f,0.027f,0.026f,0.024f,0.022f,0.021f,0.019f },
404 
405 // n S+ K+ (p S- K0)
406  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
407  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
408  0.015f,0.06f, 0.07f, 0.065f,0.05f, 0.04f, 0.033f,0.026f,0.02f, 0.015f },
409 //
410 // multiplicity 4 (18 channels)
411 //
412 // p p pi+ pi- (n n pi+ pi-)
413  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
414  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.05f, 0.6f, 1.9f,
415  2.8f, 3.0f, 3.0f, 2.8f, 2.5f, 2.1f, 1.9f, 1.6f, 1.4f, 1.2f },
416 
417 // p n pi+ pi0 (p n pi- pi0)
418  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
419  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.05f, 0.6f, 3.5f,
420  4.0f, 3.9f, 3.5f, 3.1f, 2.8f, 2.4f, 2.2f, 1.9f, 1.7f, 1.5f },
421 
422 // p p pi0 pi0 (n n pi0 pi0)
423  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
424  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.02f, 0.24f, 0.76f,
425  1.1f, 1.2f, 1.2f, 1.1f, 1.0f, 0.84f, 0.76f, 0.64f, 0.56f, 0.48f },
426 
427 // n n pi+ pi+ (p p pi- pi-)
428  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
429  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.02f, 0.24f, 1.4f,
430  1.6f, 1.6f, 1.4f, 1.2f, 1.1f, 1.0f, 0.88f, 0.76f, 0.68f, 0.6f },
431 
432 // L K+ p pi0 (L K0 n pi0)
433  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
434  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
435  0.002f,0.015f,0.05f, 0.06f, 0.052f,0.042f,0.037f,0.029f,0.025f,0.020f },
436 
437 // L K0 p pi+ (L K+ n pi-)
438  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
439  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
440  0.002f,0.015f,0.06f, 0.086f,0.09f, 0.082f,0.072f,0.06f, 0.051f,0.043f },
441 
442 // L K+ n pi+ (L K0 p pi-)
443  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
444  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
445  0.002f,0.01f, 0.04f, 0.08f, 0.066f,0.058f,0.05f, 0.04f, 0.035f,0.03f },
446 
447 // S0 K+ n pi+ (S0 K0 p pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
450  0.0f, 0.003f,0.02f, 0.03f, 0.03f, 0.025f,0.02f, 0.015f,0.011f,0.01f },
451 
452 // S0 K+ p pi0 (S0 K0 n pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
455  0.0f, 0.005f,0.02f, 0.025f,0.022f,0.02f, 0.015f,0.01f, 0.008f,0.007f },
456 
457 // S0 K0 p pi+ (S0 K+ n pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
460  0.0f, 0.003f,0.012f,0.04f, 0.037f,0.03f, 0.027f,0.022f,0.019f,0.016f },
461 
462 // S- K+ p pi+ (S+ K0 n pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
465  0.0f, 0.004f,0.016f,0.037f,0.031f,0.028f,0.023f,0.02f, 0.017f,0.014f },
466 
467 // S+ K0 p pi0 (S- K+ n pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
470  0.0f, 0.003f,0.007f,0.021f,0.021f,0.019f,0.017f,0.014f,0.012f,0.01f },
471 
472 // S+ K0 n pi+ (S- K+ p pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
475  0.0f, 0.003f,0.009f,0.036f,0.052f,0.043f,0.038f,0.03f, 0.026f,0.02f },
476 
477 // S+ K+ p pi- (S- K0 n pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
480  0.0f, 0.004f,0.012f,0.038f,0.037f,0.03f, 0.026f,0.02f, 0.017f,0.014f },
481 
482 // S+ K+ n pi0 (S- K0 p pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
485  0.0f, 0.004f,0.012f,0.038f,0.037f,0.03f, 0.026f,0.02f, 0.017f,0.014f },
486 
487 // p p K0 K0bar (n n K+ K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
490  0.0f, 0.001f,0.003f,0.009f,0.02f, 0.02f, 0.017f,0.014f,0.012f,0.009f },
491 
492 // p p K+ K- (n n K0 K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
495  0.0f, 0.001f,0.003f,0.009f,0.02f, 0.02f, 0.017f,0.014f,0.012f,0.009f },
496 
497 // p n K+ K0bar (p n K0 K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
500  0.0f, 0.002f,0.007f,0.029f,0.024f,0.02f, 0.017f,0.014f,0.012f,0.009f },
501 //
502 // multiplicity 5 (32 channels)
503 //
504 // p p pi+ pi- pi0 (n n pi+ pi- pi0)
505  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
506  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.06f,
507  0.4f, 1.1f, 1.8f, 2.4f, 2.4f, 2.2f, 2.0f, 1.7f, 1.5f, 1.3f },
508 
509 // p p pi0 pi0 pi0 (n n pi0 pi0 pi0)
510  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
511  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.003f,0.018f,
512  0.12f, 0.33f, 0.54f, 0.72f, 0.72f, 0.66f, 0.6f, 0.51f, 0.45f, 0.39f },
513 
514 // p n pi+ pi+ pi- (p n pi+ pi- pi-)
515  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
516  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.12f, 0.26f,
517  0.7f, 1.6f, 2.4f, 2.6f, 2.3f, 2.0f, 1.8f, 1.6f, 1.4f, 1.2f },
518 
519 // p n pi+ pi0 pi0 (p n pi- pi0 pi0)
520  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
521  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.036f,0.078f,
522  0.21f, 0.48f, 0.72f, 0.78f, 0.69f, 0.6f, 0.54f, 0.48f, 0.42f, 0.36f },
523 
524 // n n pi+ pi+ pi0 (p p pi- pi- pi0)
525  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
526  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.036f,
527  0.24f, 0.66f, 1.08f, 1.44f, 1.44f, 1.32f, 1.2f, 1.0f, 0.9f, 0.78f },
528 
529 // p L K+ pi+ pi- (n L K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
532  0.0f, 0.0f, 0.002f,0.010f,0.040f,0.045f,0.040f,0.035f,0.030f,0.020f },
533 
534 // p L K+ pi0 pi0 (n L K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
537  0.0f, 0.0f, 0.001f,0.005f,0.020f,0.022f,0.020f,0.017f,0.015f,0.010f },
538 
539 // p L K0 pi+ pi0 (n L K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
542  0.0f, 0.0f, 0.002f,0.010f,0.040f,0.045f,0.040f,0.035f,0.030f,0.020f },
543 
544 // p S0 K+ pi+ pi- (n S0 K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
547  0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
548 
549 // p S0 K+ pi0 pi0 (n S0 K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
552  0.0f, 0.0f, 0.001f,0.004f,0.015f,0.017f,0.015f,0.014f,0.011f,0.009f },
553 
554 // p S0 K0 pi+ pi0 (n S0 K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
557  0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
558 
559 // p S+ K0 pi+ pi- (n S- K+ pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
562  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
563 
564 // p S+ K0 pi0 pi0 (n S- K+ pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
567  0.0f, 0.0f, 0.001f,0.002f,0.009f,0.020f,0.013f,0.014f,0.011f,0.009f },
568 
569 // p S+ K+ pi- pi0 (n S- K0 pi+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
572  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
573 
574 // p S- K+ pi+ pi0 (n S+ K0 pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
577  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
578 
579 // p S- K0 pi+ pi+ (n S+ K+ pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
582  0.0f, 0.0f, 0.001f,0.007f,0.030f,0.035f,0.030f,0.028f,0.021f,0.017f },
583 
584 // n L K+ pi+ pi0 (p L K0 pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
587  0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
588 
589 // n L K0 pi+ pi+ (p L K+ pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
592  0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
593 
594 // n S0 K+ pi+ pi0 (p S0 K0 pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
597  0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
598 
599 // n S0 K0 pi+ pi+ (p S0 K+ pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
602  0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
603 
604 // n S+ K0 pi+ pi0 (p S- K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
607  0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
608 
609 // n S+ K+ pi+ pi- (p S- K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
612  0.0f, 0.0f, 0.002f,0.011f,0.042f,0.039f,0.030f,0.022f,0.018f,0.014f },
613 
614 // n S+ K+ pi0 pi0 (p S- K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
617  0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
618 
619 // n S- K+ pi+ pi+ (p S+ K0 pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
622  0.0f, 0.0f, 0.001f,0.005f,0.021f,0.020f,0.015f,0.011f,0.009f,0.007f },
623 
624 // p p pi+ K0 K- (n n pi- K+ K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
627  0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
628 
629 // p p pi- K+ K0bar (n n pi+ K0 K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
632  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
633 
634 // p p pi0 K0 K0bar (n n pi0 K+ K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
637  0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
638 
639 // p p pi0 K+ K- (n n pi0 K0 K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
642  0.0f, 0.0f, 0.002f,0.004f,0.012f,0.040f,0.055f,0.045f,0.036f,0.030f },
643 
644 // p n pi+ K0 K0bar (p n pi- K+ K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
647  0.0f, 0.0f, 0.002f,0.012f,0.035f,0.058f,0.050f,0.033f,0.023f,0.016f },
648 
649 // p n pi+ K+ K- (p n pi- K0 K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
652  0.0f, 0.0f, 0.002f,0.012f,0.035f,0.058f,0.050f,0.033f,0.023f,0.016f },
653 
654 // p n pi0 K+ K0bar (p n pi0 K0 K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
657  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
658 
659 // n n pi+ K+ K0bar (p p pi- K0 K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
662  0.0f, 0.001f,0.002f,0.004f,0.018f,0.040f,0.033f,0.028f,0.021f,0.017f },
663 //
664 // multiplicity 6 (7 channels)
665 //
666 // p p pi+ pi+ pi- pi- (n n pi+ pi+ pi- pi-)
667  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
668  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
669  0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
670 
671 // p p pi+ pi- pi0 pi0 (n n pi+ pi- pi0 pi0)
672  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
673  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
674  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
675 
676 // p p pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0)
677  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
678  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
679  0.01f, 0.02f, 0.05f, 0.1f, 0.13f, 0.12f, 0.11f, 0.1f, 0.1f, 0.09f },
680 
681 // p n pi+ pi+ pi- pi0 (p n pi+ pi- pi- pi0)
682  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
683  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
684  0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
685 
686 // p n pi+ pi0 pi0 pi0 (p n pi- pi0 pi0 pi0)
687  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
688  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
689  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
690 
691 // n n pi+ pi+ pi+ pi- (p p pi+ pi- pi- pi-)
692  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
693  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
694  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
695 
696 // n n pi+ pi+ pi0 pi0 (p p pi- pi- pi0 pi0)
697  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
698  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
699  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
700 //
701 // multiplicity 7 (8 channels)
702 //
703 // p p pi+ pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
706  0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
707 
708 // p p pi+ pi- pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
711  0.0f, 0.0f, 0.036f,0.096f,0.30f, 0.42f, 0.42f, 0.42f, 0.40f, 0.37f },
712 
713 // p p pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
716  0.0f, 0.0f, 0.018f,0.048f,0.14f, 0.20f, 0.22f, 0.20f, 0.19f, 0.18f },
717 
718 // p n pi+ pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
721  0.0f, 0.0f, 0.06f, 0.19f, 0.31f, 0.41f, 0.44f, 0.47f, 0.45f, 0.45f },
722 
723 // p n pi+ pi+ pi- pi0 pi0 (p n pi+ pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
726  0.0f, 0.0f, 0.036f,0.12f, 0.18f, 0.24f, 0.26f, 0.23f, 0.28f, 0.26f },
727 
728 // p n pi+ pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
731  0.0f, 0.0f, 0.024f,0.06f, 0.08f, 0.12f, 0.13f, 0.14f, 0.13f, 0.13f },
732 
733 // n n pi+ pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
736  0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
737 
738 // n n pi+ pi+ pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
741  0.0f, 0.0f, 0.036f,0.096f,0.30f, 0.42f, 0.42f, 0.41f, 0.40f, 0.37f },
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
747  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
748  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f,0.27f, 0.30f, 0.27f, 0.24f },
749 
750 // p p pi+ pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
753  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f,0.27f, 0.30f, 0.27f, 0.24f },
754 
755 // p p pi+ pi- pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
758  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
759 
760 // p p pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
763  0.0f, 0.0f, 0.003f,0.01f, 0.024f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
764 
765 // p n pi+ pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
768  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
769 
770 // p n pi+ pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
773  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
774 
775 // p n pi+ pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
778  0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
779 
780 // n n pi+ pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
783  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
784 
785 // n n pi+ pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
788  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
789 
790 // n n pi+ pi+ pi0 pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
793  0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
799  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
800  0.0f, 0.002f,0.008f,0.025f,0.074f,0.11f, 0.14f, 0.15f, 0.15f, 0.15f },
801 
802 // p p pi+ pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
805  0.0f, 0.001f,0.006f,0.020f,0.058f,0.091f,0.11f, 0.12f, 0.12f, 0.12f },
806 
807 // p p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
810  0.0f, 0.001f,0.004f,0.012f,0.035f,0.055f,0.065f,0.07f, 0.07f, 0.07f },
811 
812 // p p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
815  0.0f, 0.0f, 0.0f, 0.006f,0.027f,0.032f,0.04f, 0.042f,0.042f,0.042f },
816 
817 // p n pi+ pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
820  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
821 
822 // p n pi+ pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
825  0.0f, 0.003f,0.008f,0.026f,0.078f,0.20f, 0.25f, 0.29f, 0.29f, 0.29f },
826 
827 // p n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
830  0.0f, 0.001f,0.005f,0.016f,0.047f,0.12f, 0.15f, 0.17f, 0.17f, 0.17f },
831 
832 // p n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
835  0.0f, 0.001f,0.003f,0.009f,0.029f,0.07f, 0.094f,0.10f, 0.10f, 0.10f },
836 
837 // n n pi+ pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
840  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
841 
842 // n n pi+ pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
845  0.0f, 0.001f,0.005f,0.015f,0.047f,0.12f, 0.15f, 0.17f, 0.17f, 0.17f },
846 
847 // n n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p p pi+ pi- pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
850  0.0f, 0.001f,0.003f,0.009f,0.029f,0.07f, 0.094f,0.10f, 0.10f, 0.10f }};
851 
853 // //
854 // p n and n p (|Tz| = 0) cross sections //
855 // and final state particle types //
856 // //
858 
859 // Total p n cross section as a function of kinetic energy
860 G4ThreadLocal G4double G4RPGNucleonInelastic::pNtot[30];
861 
862 // p n multiplicities as a function of kinetic energy
863 G4ThreadLocal G4double G4RPGNucleonInelastic::t0_dSigma_dMult[8][30];
864 
865 const 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 
879 const G4int G4RPGNucleonInelastic::T0_2bfs[1][2] =
880  {{pro,neu}};
881 
882 const 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 
887 const 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 
897 const 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 
912 const 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 
918 const 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 
925 const 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 
932 const 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 //
954 const G4float G4RPGNucleonInelastic::pNCrossSections[108][30] = {
955 //
956 // multiplicity 2 (1 channel)
957 //
958 // p n (p n)
959  { 0.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f, 46.0f,
960  44.0f, 42.0f, 40.0f, 35.0f, 31.0f, 27.0f, 23.0f, 19.0f, 17.0f, 15.5f,
961  14.0f, 13.0f, 12.0f, 11.0f, 10.0f, 9.5f, 9.0f, 8.5f, 8.0f, 7.7f },
962 //
963 // multiplicity 3 (9 channels)
964 //
965 // p p pi- (n n pi+)
966  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
967  0.0f, 0.0f, 0.0f, 0.0f, 0.25f, 0.9f, 1.75f, 2.3f, 2.8f, 2.8f,
968  2.2f, 1.9f, 1.6f, 1.35f, 1.1f, 0.95f, 0.8f, 0.7f, 0.6f, 0.53f },
969 
970 // p n pi0 (p n pi0)
971  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
972  0.0f, 0.0f, 0.0f, 0.0f, 1.8f, 4.7f, 8.3f, 11.3f, 12.0f, 10.2f,
973  8.2f, 6.0f, 4.9f, 3.6f, 2.5f, 2.0f, 1.6f, 1.2f, 1.0f, 0.08f },
974 
975 // n n pi+ (p p pi-)
976  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
977  0.0f, 0.0f, 0.0f, 0.0f, 0.95f, 2.4f, 4.2f, 5.6f, 6.1f, 5.1f,
978  4.1f, 3.0f, 2.5f, 1.8f, 1.2f, 1.0f, 0.8f, 0.6f, 0.5f, 0.41f},
979 
980 // p L K0 (n L K+)
981  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
982  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
983  0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
984 
985 // p S0 K0 (n S0 K+)
986  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
987  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
988  0.003f,0.010f,0.016f,0.020f,0.016f,0.015f,0.014f,0.013f,0.011f,0.010f},
989 
990 // p S- K+ (n S+ K0)
991  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
992  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
993  0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
994 
995 // n L K+ (p L K0)
996  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
997  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
998  0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
999 
1000 // n S0 K+ (p S0 K0)
1001  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1002  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1003  0.003f,0.010f,0.016f,0.020f,0.016f,0.015f,0.014f,0.013f,0.011f,0.010f},
1004 
1005 // n S+ K0 (p S- K+)
1006  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1007  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1008  0.004f,0.013f,0.021f,0.025f,0.021f,0.019f,0.018f,0.016f,0.014f,0.012f},
1009 //
1010 // multiplicity 4 (22 channels)
1011 //
1012 // p n pi+ pi- (p n pi+ pi-)
1013  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1014  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.12f, 0.38f, 1.1f, 3.5f,
1015  5.9f, 5.9f, 5.1f, 4.2f, 3.7f, 3.0f, 2.6f, 2.1f, 1.8f, 1.4f },
1016 
1017 // p p pi- pi0 (n n pi+ pi0)
1018  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1019  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.03f, 0.1f, 0.24f, 0.55f,
1020  1.2f, 1.5f, 1.45f, 1.25f, 1.0f, 0.9f, 0.8f, 0.7f, 0.6f, 0.53f },
1021 
1022 // p n pi0 pi0 (p n pi0 pi0)
1023  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1024  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.07f, 0.24f, 0.66f, 2.1f,
1025  3.6f, 3.6f, 3.1f, 2.5f, 2.2f, 1.8f, 1.5f, 1.2f, 1.1f, 0.84f },
1026 
1027 // n n pi+ pi0 (p p pi- pi0)
1028  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1029  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.03f, 0.1f, 0.24f, 0.55f,
1030  1.2f, 1.5f, 1.45f, 1.25f, 1.0f, 0.9f, 0.8f, 0.7f, 0.6f, 0.53f },
1031 
1032 // p L K+ pi- (n L K0 pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1035  0.0f, 0.002f,0.006f,0.018f,0.037f,0.036f,0.033f,0.030f,0.028f,0.023f },
1036 
1037 // p S0 K+ pi- (n S0 K0 pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1040  0.0f, 0.0f, 0.002f,0.007f,0.021f,0.025f,0.022f,0.020f,0.018f,0.017f },
1041 
1042 // p L K0 pi0 (n L K+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1045  0.0f, 0.001f,0.002f,0.005f,0.009f,0.009f,0.008f,0.007f,0.007f,0.006f },
1046 
1047 // p S0 K0 pi0 (n S0 K+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1050  0.0f, 0.0f, 0.001f,0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.004f },
1051 
1052 // p S+ K0 pi- (n S- K+ pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1055  0.0f, 0.0015f,0.004f,0.013f,0.02f,0.016f,0.013f,0.01f, 0.009f,0.007f },
1056 
1057 // p S- K+ pi0 (n S+ K0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1060  0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1061 
1062 // p S- K0 pi+ (n S+ K+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1065  0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1066 
1067 // n L K+ pi0 (p L K0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1070  0.0f, 0.002f,0.006f,0.018f,0.037f,0.036f,0.033f,0.030f,0.028f,0.023f },
1071 
1072 // n L K0 pi+ (p L K+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1075  0.0f, 0.001f,0.002f,0.005f,0.009f,0.009f,0.008f,0.007f,0.007f,0.006f },
1076 
1077 // n S+ K+ pi- (p S- K0 pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1080  0.0f, 0.0f, 0.002f,0.007f,0.021f,0.025f,0.022f,0.020f,0.018f,0.017f },
1081 
1082 // n S+ K0 pi0 (p S- K+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1085  0.0f, 0.0f, 0.001f,0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.004f },
1086 
1087 // n S0 K+ pi0 (p S0 K0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1090  0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1091 
1092 // n S0 K0 pi+ (p S0 K+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1095  0.0f, 0.003f,0.008f,0.026f,0.042f,0.042f,0.035f,0.029f,0.023f,0.018f },
1096 
1097 // n S- K+ pi+ (p S+ K0 pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1100  0.0f, 0.0015f,0.004f,0.013f,0.02f, 0.016f,0.013f,0.01f, 0.009f,0.007f },
1101 
1102 // p n K+ K- (p n K0 K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1105  0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1106 
1107 // p n K0 K0bar (p n K+ K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1110  0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1111 
1112 // p p K0 K- (n n K+ K0bar)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1115  0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1116 
1117 // n n K+ K0bar (p p K0 K-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1120  0.0f, 0.0f, 0.001f,0.005f,0.015f,0.019f,0.015f,0.012f,0.009f,0.007f },
1121 //
1122 // multiplicity 5 (38 channels)
1123 //
1124 // p n pi+ pi- pi0 (p n pi+ pi- pi0)
1125  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1126  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.04f,
1127  0.3f, 0.82f, 1.35f, 1.8f, 1.8f, 1.65f, 1.5f, 1.28f, 1.12f, 0.98f },
1128 
1129 // p n pi0 pi0 pi0 (p n pi0 pi0 pi0)
1130  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1131  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.004f,0.022f,
1132  0.15f, 0.41f, 0.68f, 0.9f, 0.9f, 0.82f, 0.75f, 0.64f, 0.55f, 0.49f },
1133 
1134 // p p pi+ pi- pi- (n n pi+ pi+ pi-)
1135  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1136  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.09f, 0.2f,
1137  0.52f, 1.2f, 1.8f, 2.0f, 1.7f, 1.5f, 1.35f, 1.2f, 1.05f, 0.9f },
1138 
1139 // p p pi- pi0 pi0 (n n pi+ pi0 pi0)
1140  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1141  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.04f, 0.1f,
1142  0.26f, 0.6f, 0.9f, 0.98f, 0.86f, 0.75f, 0.68f, 0.6f, 0.52f, 0.45f },
1143 
1144 // n n pi+ pi+ pi- (p p pi+ pi- pi-)
1145  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1146  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.01f, 0.04f,
1147  0.3f, 0.82f, 1.35f, 1.8f, 1.8f, 1.65f, 1.5f, 1.28f, 1.12f, 0.98f },
1148 
1149 // n n pi+ pi0 pi0 (p p pi- pi0 pi0)
1150  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1151  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.003f,0.02f,
1152  0.15f, 0.41f, 0.68f, 0.9f, 0.9f, 0.82f, 0.75f, 0.64f, 0.56f, 0.49f },
1153 
1154 // p L K+ pi- pi0 (n L K0 pi+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1157  0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1158 
1159 // p L K0 pi+ pi- (n L K+ pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1162  0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1163 
1164 // p L K0 pi0 pi0 (n L K+ pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1167  0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.09f, 0.07f, 0.07f },
1168 
1169 // p S0 K0 pi+ pi- (n S0 K+ pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1172  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1173 
1174 // p S0 K0 pi0 pi0 (n S0 K+ pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1177  0.0f, 0.0f, 0.001f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1178 
1179 // p S0 K+ pi- pi0 (n S0 K0 pi+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1182  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1183 
1184 // p S+ K+ pi- pi- (n S- K0 pi+ pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1187  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1188 
1189 // p S+ K0 pi- pi0 (n S- K+ pi+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1192  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1193 
1194 // p S- K0 pi+ pi0 (n S+ K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1197  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1198 
1199 // p S- K+ pi+ pi- (n S+ K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1202  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1203 
1204 // p S- K+ pi0 pi0 (n S+ K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1207  0.0f, 0.0f, 0.002f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1208 
1209 // n L K+ pi+ pi- (p L K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1212  0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1213 
1214 // n L K+ pi0 pi0 (p L K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1217  0.0f, 0.0f, 0.001f,0.004f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1218 
1219 // n L K0 pi+ pi0 (p L K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1222  0.0f, 0.0f, 0.002f,0.009f,0.023f,0.025f,0.022f,0.018f,0.015f,0.013f },
1223 
1224 // n S0 K+ pi+ pi- (p S0 K0 pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1227  0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1228 
1229 // n S0 K+ pi0 pi0 (p S0 K0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1232  0.0f, 0.0f, 0.0f, 0.002f,0.005f,0.006f,0.005f,0.005f,0.004f,0.003f },
1233 
1234 // n S0 K0 pi+ pi0 (p S0 K+ pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1237  0.0f, 0.0f, 0.001f,0.005f,0.011f,0.012f,0.011f,0.009f,0.007f,0.006f },
1238 
1239 // n S+ K0 pi+ pi- (p S- K+ pi+ pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1242  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1243 
1244 // n S+ K0 pi0 pi0 (p S- K+ pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1247  0.0f, 0.0f, 0.001f,0.003f,0.007f,0.008f,0.007f,0.006f,0.005f,0.004f },
1248 
1249 // n S+ K+ pi- pi0 (p S- K0 pi+ pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1252  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1253 
1254 // n S- K+ pi+ pi0 (p S+ K0 pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1257  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1258 
1259 // n S- K0 pi+ pi+ (p S+ K+ pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1262  0.0f, 0.001f,0.003f,0.006f,0.014f,0.016f,0.015f,0.012f,0.010f,0.008f },
1263 
1264 // p n K+ K- pi0 (p n K0 K0bar pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1267  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1268 
1269 // p n K0 K0bar pi0 (p n K+ K- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1272  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1273 
1274 // p n K0 K- pi+ (p n K+ K0bar pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1277  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1278 
1279 // p n K+ K0bar pi- (p n K0 K- pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1282  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1283 
1284 // p p K0 K0bar pi- (n n K+ K- pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1287  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1288 
1289 // p p K+ K- pi- (n n K0 K0bar pi+)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1292  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1293 
1294 // p p K0 K- pi0 (n n K+ K0bar pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1297  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1298 
1299 // n n K+ K- pi+ (p p K0 K0bar pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1302  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1303 
1304 // n n K0 K0bar pi+ (p p K+ K- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1307  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1308 
1309 // n n K+ K0bar pi0 (p p K0 K- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1312  0.0f, 0.0f, 0.002f,0.004f,0.010f,0.012f,0.011f,0.010f,0.009f,0.007f },
1313 //
1314 // multiplicity 6 (7 channels)
1315 //
1316 // p n pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi-)
1317  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1318  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1319  0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1320 
1321 // p n pi+ pi- pi0 pi0 (p n pi+ pi- pi0 pi0)
1322  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1323  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1324  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1325 
1326 // p n pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0)
1327  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1328  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1329  0.01f, 0.02f, 0.05f, 0.1f, 0.13f, 0.12f, 0.11f, 0.1f, 0.1f, 0.09f },
1330 
1331 // p p pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi0)
1332  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1333  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1334  0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1335 
1336 // p p pi- pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0)
1337  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1338  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1339  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1340 
1341 // n n pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi0)
1342  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1343  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1344  0.06f, 0.1f, 0.18f, 0.38f, 0.49f, 0.46f, 0.43f, 0.40f, 0.38f, 0.36f },
1345 
1346 // n n pi+ pi0 pi0 pi0 (p p pi- pi0 pi0 pi0)
1347  { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1348  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1349  0.03f, 0.05f, 0.09f, 0.19f, 0.25f, 0.23f, 0.22f, 0.2f, 0.19f, 0.18f },
1350 //
1351 // multiplicity 7 (9 channels)
1352 //
1353 // p n pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1356  0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
1357 
1358 // p n pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1361  0.0f, 0.0f, 0.03f, 0.08f, 0.25f, 0.35f, 0.35f, 0.35f, 0.33f, 0.31f },
1362 
1363 // p n pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1366  0.0f, 0.0f, 0.015f,0.04f, 0.12f, 0.17f, 0.18f, 0.17f, 0.16f, 0.15f },
1367 
1368 // p p pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1371  0.0f, 0.0f, 0.06f, 0.19f, 0.31f, 0.41f, 0.44f, 0.47f, 0.45f, 0.45f },
1372 
1373 // p p pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1376  0.0f, 0.0f, 0.03f, 0.1f, 0.15f, 0.2f, 0.22f, 0.23f, 0.22f, 0.22f },
1377 
1378 // p p pi- pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1381  0.0f, 0.0f, 0.02f, 0.05f, 0.07f, 0.1f, 0.11f, 0.12f, 0.11f, 0.11f },
1382 
1383 // n n pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1386  0.0f, 0.0f, 0.06f, 0.17f, 0.5f, 0.7f, 0.7f, 0.69f, 0.66f, 0.62f },
1387 
1388 // n n pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1391  0.0f, 0.0f, 0.03f, 0.08f, 0.25f, 0.35f, 0.35f, 0.34f, 0.33f, 0.31f },
1392 
1393 // n n pi+ pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1396  0.0f, 0.0f, 0.02f, 0.05f, 0.07f, 0.1f, 0.11f, 0.12f, 0.11f, 0.11f },
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1402  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1403  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1404 
1405 // p n pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1408  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1409 
1410 // p n pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1413  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1414 
1415 // p n pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1418  0.0f, 0.0f, 0.003f,0.01f, 0.024f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
1419 
1420 // p p pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1423  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1424 
1425 // p p pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1428  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1429 
1430 // p p pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1433  0.0f, 0.0f, 0.003f,0.01f, 0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
1434 
1435 // n n pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1438  0.0f, 0.003f,0.01f, 0.024f,0.075f,0.18f, 0.27f, 0.30f, 0.27f, 0.24f },
1439 
1440 // n n pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1443  0.0f, 0.003f,0.006f,0.015f,0.045f,0.12f, 0.15f, 0.18f, 0.15f, 0.15f },
1444 
1445 // n n pi+ pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1448  0.0f, 0.0f, 0.003f,0.009f,0.027f,0.06f, 0.09f, 0.12f, 0.09f, 0.09f },
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1454  0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1455  0.0f, 0.002f,0.008f,0.025f,0.074f,0.11f, 0.14f, 0.15f, 0.15f, 0.15f },
1456 
1457 // p n pi+ pi+ pi- pi- pi0 pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1460  0.0f, 0.001f,0.005f,0.015f,0.045f,0.07f, 0.084f,0.09f, 0.09f, 0.09f },
1461 
1462 // p n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1465  0.0f, 0.001f,0.003f,0.009f,0.027f,0.042f,0.05f, 0.054f,0.054f,0.054f },
1466 
1467 // p n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1470  0.0f, 0.0f, 0.0f, 0.005f,0.016f,0.025f,0.03f, 0.032f,0.032f,0.032f },
1471 
1472 // p p pi+ pi+ pi+ pi- pi- pi- pi- (n n pi+ pi+ pi+ pi+ pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1475  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1476 
1477 // p p pi+ pi+ pi- pi- pi- pi0 pi0 (n n pi+ pi+ pi+ pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1480  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1481 
1482 // p p pi+ pi- pi- pi0 pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1485  0.0f, 0.001f,0.004f,0.012f,0.036f,0.09f, 0.12f, 0.13f, 0.13f, 0.13f },
1486 
1487 // p p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1490  0.0f, 0.001f,0.002f,0.007f,0.022f,0.054f,0.072f,0.078f,0.078f,0.078f },
1491 
1492 // n n pi+ pi+ pi+ pi+ pi- pi- pi- (p p pi+ pi+ pi+ pi- pi- pi- pi-)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1495  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1496 
1497 // n n pi+ pi+ pi+ pi- pi- pi0 pi0 (p p pi+ pi+ pi- pi- pi- pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1500  0.0f, 0.002f,0.006f,0.02f, 0.06f, 0.15f, 0.19f, 0.22f, 0.22f, 0.22f },
1501 
1502 // n n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1505  0.0f, 0.001f,0.004f,0.012f,0.036f,0.09f, 0.12f, 0.13f, 0.13f, 0.13f },
1506 
1507 // n n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0 pi0)
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.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,
1510  0.0f, 0.001f,0.002f,0.007f,0.022f,0.054f,0.072f,0.078f,0.078f,0.078f }};
km
static constexpr double km
Definition: G4SIunits.hh:133
G4RPGNucleonInelastic::T0_8bfs
static const G4int T0_8bfs[10][8]
Definition: G4RPGNucleonInelastic.hh:93
G4float
float G4float
Definition: G4Types.hh:77
G4RPGNucleonInelastic::T1_8bfs
static const G4int T1_8bfs[2][10][8]
Definition: G4RPGNucleonInelastic.hh:84
G4RPGNucleonInelastic::T0_2bfs
static const G4int T0_2bfs[1][2]
Definition: G4RPGNucleonInelastic.hh:87
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:524
G4RPGNucleonInelastic::T0_7bfs
static const G4int T0_7bfs[9][7]
Definition: G4RPGNucleonInelastic.hh:92
G4ThreadLocal
#define G4ThreadLocal
Definition: tls.hh:89
G4int
int G4int
Definition: G4Types.hh:78
G4RPGNucleonInelastic::pPCrossSections
static const G4float pPCrossSections[93][30]
Definition: G4RPGNucleonInelastic.hh:101
G4RPGNucleonInelastic::pNtot
static G4ThreadLocal G4double pNtot[30]
Definition: G4RPGNucleonInelastic.hh:97
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4RPGNucleonInelastic::GetMultiplicityT1
G4int GetMultiplicityT1(G4double KE) const
Definition: G4RPGNucleonInelastic.cc:118
G4cout
G4GLOB_DLL std::ostream G4cout
G4RPGNucleonInelastic::T1_4bfs
static const G4int T1_4bfs[2][18][4]
Definition: G4RPGNucleonInelastic.hh:80
G4RPGNucleonInelastic::T0_9bfs
static const G4int T0_9bfs[12][9]
Definition: G4RPGNucleonInelastic.hh:94
G4RPGNucleonInelastic::T0_4bfs
static const G4int T0_4bfs[22][4]
Definition: G4RPGNucleonInelastic.hh:89
G4RPGNucleonInelastic::GetFSPartTypesForT0
std::vector< G4int > GetFSPartTypesForT0(G4int mult, G4double KE) const
Definition: G4RPGNucleonInelastic.cc:138
Randomize.hh
G4RPGNucleonInelastic::T0_3bfs
static const G4int T0_3bfs[9][3]
Definition: G4RPGNucleonInelastic.hh:88
G4RPGNucleonInelastic::pPtot
static G4ThreadLocal G4double pPtot[30]
Definition: G4RPGNucleonInelastic.hh:96
G4RPGNucleonInelastic::T1_5bfs
static const G4int T1_5bfs[2][32][5]
Definition: G4RPGNucleonInelastic.hh:81
G4RPGNucleonInelastic::T1_2bfs
static const G4int T1_2bfs[2][1][2]
Definition: G4RPGNucleonInelastic.hh:78
G4RPGNucleonInelastic::pNindex
static const G4int pNindex[8][2]
Definition: G4RPGNucleonInelastic.hh:76
G4RPGNucleonInelastic.hh
G4RPGNucleonInelastic::T1_6bfs
static const G4int T1_6bfs[2][7][6]
Definition: G4RPGNucleonInelastic.hh:82
GeV
static constexpr double GeV
Definition: G4SIunits.hh:217
G4RPGNucleonInelastic::pPindex
static const G4int pPindex[8][2]
Definition: G4RPGNucleonInelastic.hh:75
G4RPGNucleonInelastic::G4RPGNucleonInelastic
G4RPGNucleonInelastic(const G4String &modelName="RPGNucleonInelastic")
Definition: G4RPGNucleonInelastic.cc:33
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4endl
#define G4endl
Definition: G4ios.hh:61
G4RPGNucleonInelastic::pNCrossSections
static const G4float pNCrossSections[108][30]
Definition: G4RPGNucleonInelastic.hh:102
G4double
double G4double
Definition: G4Types.hh:76
G4SystemOfUnits.hh
G4RPGNucleonInelastic::GetMultiplicityT0
G4int GetMultiplicityT0(G4double KE) const
Definition: G4RPGNucleonInelastic.cc:99
G4RPGNucleonInelastic::T1_7bfs
static const G4int T1_7bfs[2][8][7]
Definition: G4RPGNucleonInelastic.hh:83
G4RPGNucleonInelastic::T1_3bfs
static const G4int T1_3bfs[2][6][3]
Definition: G4RPGNucleonInelastic.hh:79
G4RPGNucleonInelastic::T0_6bfs
static const G4int T0_6bfs[7][6]
Definition: G4RPGNucleonInelastic.hh:91
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:507
G4RPGNucleonInelastic::t1_dSigma_dMult
static G4ThreadLocal G4double t1_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:98
G4RPGNucleonInelastic::T1_9bfs
static const G4int T1_9bfs[2][11][9]
Definition: G4RPGNucleonInelastic.hh:85
G4RPGNucleonInelastic::T0_5bfs
static const G4int T0_5bfs[38][5]
Definition: G4RPGNucleonInelastic.hh:90
G4RPGNucleonInelastic::t0_dSigma_dMult
static G4ThreadLocal G4double t0_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:99
G4RPGNucleonInelastic::GetFSPartTypesForT1
std::vector< G4int > GetFSPartTypesForT1(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGNucleonInelastic.cc:186