Geant4.10
 All Data Structures 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$
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 
233 ///////////////////////////////////////////////////////////////////////////////
234 // //
235 // p p and n n (|Tz| = 1) cross sections //
236 // and final state particle types //
237 // //
238 ///////////////////////////////////////////////////////////////////////////////
239 
240 // Total p p cross section as a function of kinetic energy
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.0,330.0,240.0,160.0,110.0, 85.0, 63.0, 44.0, 33.0, 28.0,
375  25.0, 24.0, 23.0, 23.0, 26.3, 26.1, 25.0, 23.5, 21.0, 18.0,
376  16.0, 14.3, 12.5, 11.2, 10.3, 9.6, 9.0, 8.5, 8.0, 7.7 },
377 //
378 // multiplicity 3 (6 channels)
379 //
380 // p p pi0 (n n pi0)
381  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
382  0.0, 0.0, 0.0, 0.0, 0.0, 1.4, 4.0, 4.3, 4.0, 4.0,
383  3.6, 3.0, 2.8, 2.5, 1.7, 1.3, 1.1, 1.0, 0.9, 0.85 },
384 
385 // p n pi+ (p n pi-)
386  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
387  0.0, 0.0, 0.0, 0.0, 0.7, 4.5, 15.0, 19.1, 18.0, 16.0,
388  13.0, 10.0, 8.2, 6.0, 4.3, 3.3, 2.6, 2.0, 1.65, 1.4 },
389 
390 // p L K+ (n L K0)
391  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
392  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.012,
393  0.03, 0.06, 0.06, 0.055,0.05, 0.047,0.043,0.04, 0.037,0.033 },
394 
395 // p S0 K+ (n S0 K0)
396  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
397  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
398  0.006,0.02, 0.027,0.026,0.021,0.018,0.015,0.011,0.009,0.007 },
399 
400 // p S+ K0 (n S- K+)
401  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
402  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
403  0.013,0.025,0.03, 0.029,0.027,0.026,0.024,0.022,0.021,0.019 },
404 
405 // n S+ K+ (p S- K0)
406  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
407  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
408  0.015,0.06, 0.07, 0.065,0.05, 0.04, 0.033,0.026,0.02, 0.015 },
409 //
410 // multiplicity 4 (18 channels)
411 //
412 // p p pi+ pi- (n n pi+ pi-)
413  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
414  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.6, 1.9,
415  2.8, 3.0, 3.0, 2.8, 2.5, 2.1, 1.9, 1.6, 1.4, 1.2 },
416 
417 // p n pi+ pi0 (p n pi- pi0)
418  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
419  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.6, 3.5,
420  4.0, 3.9, 3.5, 3.1, 2.8, 2.4, 2.2, 1.9, 1.7, 1.5 },
421 
422 // p p pi0 pi0 (n n pi0 pi0)
423  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
424  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.24, 0.76,
425  1.1, 1.2, 1.2, 1.1, 1.0, 0.84, 0.76, 0.64, 0.56, 0.48 },
426 
427 // n n pi+ pi+ (p p pi- pi-)
428  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
429  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.24, 1.4,
430  1.6, 1.6, 1.4, 1.2, 1.1, 1.0, 0.88, 0.76, 0.68, 0.6 },
431 
432 // L K+ p pi0 (L K0 n pi0)
433  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
434  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
435  0.002,0.015,0.05, 0.06, 0.052,0.042,0.037,0.029,0.025,0.020 },
436 
437 // L K0 p pi+ (L K+ n pi-)
438  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
439  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
440  0.002,0.015,0.06, 0.086,0.09, 0.082,0.072,0.06, 0.051,0.043 },
441 
442 // L K+ n pi+ (L K0 p pi-)
443  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
444  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
445  0.002,0.01, 0.04, 0.08, 0.066,0.058,0.05, 0.04, 0.035,0.03 },
446 
447 // S0 K+ n pi+ (S0 K0 p pi-)
448  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
449  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
450  0.0, 0.003,0.02, 0.03, 0.03, 0.025,0.02, 0.015,0.011,0.01 },
451 
452 // S0 K+ p pi0 (S0 K0 n pi0)
453  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
454  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
455  0.0, 0.005,0.02, 0.025,0.022,0.02, 0.015,0.01, 0.008,0.007 },
456 
457 // S0 K0 p pi+ (S0 K+ n pi-)
458  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
459  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
460  0.0, 0.003,0.012,0.04, 0.037,0.03, 0.027,0.022,0.019,0.016 },
461 
462 // S- K+ p pi+ (S+ K0 n pi-)
463  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
464  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
465  0.0, 0.004,0.016,0.037,0.031,0.028,0.023,0.02, 0.017,0.014 },
466 
467 // S+ K0 p pi0 (S- K+ n pi0)
468  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
469  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
470  0.0, 0.003,0.007,0.021,0.021,0.019,0.017,0.014,0.012,0.01 },
471 
472 // S+ K0 n pi+ (S- K+ p pi-)
473  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
474  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
475  0.0, 0.003,0.009,0.036,0.052,0.043,0.038,0.03, 0.026,0.02 },
476 
477 // S+ K+ p pi- (S- K0 n pi+)
478  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
479  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
480  0.0, 0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
481 
482 // S+ K+ n pi0 (S- K0 p pi0)
483  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
484  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
485  0.0, 0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
486 
487 // p p K0 K0bar (n n K+ K-)
488  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
489  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
490  0.0, 0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
491 
492 // p p K+ K- (n n K0 K0bar)
493  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
494  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
495  0.0, 0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
496 
497 // p n K+ K0bar (p n K0 K-)
498  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
499  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
500  0.0, 0.002,0.007,0.029,0.024,0.02, 0.017,0.014,0.012,0.009 },
501 //
502 // multiplicity 5 (32 channels)
503 //
504 // p p pi+ pi- pi0 (n n pi+ pi- pi0)
505  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
506  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.06,
507  0.4, 1.1, 1.8, 2.4, 2.4, 2.2, 2.0, 1.7, 1.5, 1.3 },
508 
509 // p p pi0 pi0 pi0 (n n pi0 pi0 pi0)
510  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
511  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.003,0.018,
512  0.12, 0.33, 0.54, 0.72, 0.72, 0.66, 0.6, 0.51, 0.45, 0.39 },
513 
514 // p n pi+ pi+ pi- (p n pi+ pi- pi-)
515  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
516  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12, 0.26,
517  0.7, 1.6, 2.4, 2.6, 2.3, 2.0, 1.8, 1.6, 1.4, 1.2 },
518 
519 // p n pi+ pi0 pi0 (p n pi- pi0 pi0)
520  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
521  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.036,0.078,
522  0.21, 0.48, 0.72, 0.78, 0.69, 0.6, 0.54, 0.48, 0.42, 0.36 },
523 
524 // n n pi+ pi+ pi0 (p p pi- pi- pi0)
525  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
526  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.036,
527  0.24, 0.66, 1.08, 1.44, 1.44, 1.32, 1.2, 1.0, 0.9, 0.78 },
528 
529 // p L K+ pi+ pi- (n L K0 pi+ pi-)
530  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
531  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
532  0.0, 0.0, 0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
533 
534 // p L K+ pi0 pi0 (n L K0 pi0 pi0)
535  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
536  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
537  0.0, 0.0, 0.001,0.005,0.020,0.022,0.020,0.017,0.015,0.010 },
538 
539 // p L K0 pi+ pi0 (n L K+ pi- pi0)
540  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
541  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
542  0.0, 0.0, 0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
543 
544 // p S0 K+ pi+ pi- (n S0 K0 pi+ pi-)
545  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
546  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
547  0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
548 
549 // p S0 K+ pi0 pi0 (n S0 K0 pi0 pi0)
550  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
551  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
552  0.0, 0.0, 0.001,0.004,0.015,0.017,0.015,0.014,0.011,0.009 },
553 
554 // p S0 K0 pi+ pi0 (n S0 K+ pi- pi0)
555  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
556  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
557  0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
558 
559 // p S+ K0 pi+ pi- (n S- K+ pi+ pi-)
560  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
561  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
562  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
563 
564 // p S+ K0 pi0 pi0 (n S- K+ pi0 pi0)
565  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
566  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
567  0.0, 0.0, 0.001,0.002,0.009,0.020,0.013,0.014,0.011,0.009 },
568 
569 // p S+ K+ pi- pi0 (n S- K0 pi+ pi0)
570  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
571  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
572  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
573 
574 // p S- K+ pi+ pi0 (n S+ K0 pi- pi0)
575  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
576  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
577  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
578 
579 // p S- K0 pi+ pi+ (n S+ K+ pi- pi-)
580  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
581  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
582  0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
583 
584 // n L K+ pi+ pi0 (p L K0 pi- pi0)
585  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
586  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
587  0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
588 
589 // n L K0 pi+ pi+ (p L K+ pi- pi-)
590  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
591  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
592  0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
593 
594 // n S0 K+ pi+ pi0 (p S0 K0 pi- pi0)
595  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
596  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
597  0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
598 
599 // n S0 K0 pi+ pi+ (p S0 K+ pi- pi-)
600  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
601  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
602  0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
603 
604 // n S+ K0 pi+ pi0 (p S- K+ pi- pi0)
605  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
606  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
607  0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
608 
609 // n S+ K+ pi+ pi- (p S- K0 pi+ pi-)
610  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
611  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
612  0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
613 
614 // n S+ K+ pi0 pi0 (p S- K0 pi0 pi0)
615  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
616  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
617  0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
618 
619 // n S- K+ pi+ pi+ (p S+ K0 pi- pi-)
620  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
621  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
622  0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
623 
624 // p p pi+ K0 K- (n n pi- K+ K0bar)
625  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
626  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
627  0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
628 
629 // p p pi- K+ K0bar (n n pi+ K0 K-)
630  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
631  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
632  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
633 
634 // p p pi0 K0 K0bar (n n pi0 K+ K-)
635  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
636  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
637  0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
638 
639 // p p pi0 K+ K- (n n pi0 K0 K0bar)
640  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
641  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
642  0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
643 
644 // p n pi+ K0 K0bar (p n pi- K+ K-)
645  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
646  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
647  0.0, 0.0, 0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
648 
649 // p n pi+ K+ K- (p n pi- K0 K0bar)
650  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
651  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
652  0.0, 0.0, 0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
653 
654 // p n pi0 K+ K0bar (p n pi0 K0 K-)
655  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
656  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
657  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
658 
659 // n n pi+ K+ K0bar (p p pi- K0 K-)
660  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
661  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
662  0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
663 //
664 // multiplicity 6 (7 channels)
665 //
666 // p p pi+ pi+ pi- pi- (n n pi+ pi+ pi- pi-)
667  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
668  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
669  0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
670 
671 // p p pi+ pi- pi0 pi0 (n n pi+ pi- pi0 pi0)
672  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
673  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
674  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
675 
676 // p p pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0)
677  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
678  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
679  0.01, 0.02, 0.05, 0.1, 0.13, 0.12, 0.11, 0.1, 0.1, 0.09 },
680 
681 // p n pi+ pi+ pi- pi0 (p n pi+ pi- pi- pi0)
682  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
683  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
684  0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
685 
686 // p n pi+ pi0 pi0 pi0 (p n pi- pi0 pi0 pi0)
687  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
688  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
689  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
690 
691 // n n pi+ pi+ pi+ pi- (p p pi+ pi- pi- pi-)
692  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
693  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
694  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
695 
696 // n n pi+ pi+ pi0 pi0 (p p pi- pi- pi0 pi0)
697  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
698  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
699  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
700 //
701 // multiplicity 7 (8 channels)
702 //
703 // p p pi+ pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi- pi0)
704  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
705  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
706  0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
707 
708 // p p pi+ pi- pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0)
709  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
710  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
711  0.0, 0.0, 0.036,0.096,0.30, 0.42, 0.42, 0.42, 0.40, 0.37 },
712 
713 // p p pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0)
714  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
715  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
716  0.0, 0.0, 0.018,0.048,0.14, 0.20, 0.22, 0.20, 0.19, 0.18 },
717 
718 // p n pi+ pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi- pi-)
719  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
720  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
721  0.0, 0.0, 0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
722 
723 // p n pi+ pi+ pi- pi0 pi0 (p n pi+ pi- pi- pi0 pi0)
724  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
725  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
726  0.0, 0.0, 0.036,0.12, 0.18, 0.24, 0.26, 0.23, 0.28, 0.26 },
727 
728 // p n pi+ pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0)
729  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
730  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
731  0.0, 0.0, 0.024,0.06, 0.08, 0.12, 0.13, 0.14, 0.13, 0.13 },
732 
733 // n n pi+ pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi- pi0)
734  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
735  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
736  0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
737 
738 // n n pi+ pi+ pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0)
739  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
740  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
741  0.0, 0.0, 0.036,0.096,0.30, 0.42, 0.42, 0.41, 0.40, 0.37 },
742 //
743 // multiplicity 8 (10 channels)
744 //
745 // p p pi+ pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi- pi-)
746  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
747  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
748  0.0, 0.003,0.01, 0.024,0.075,0.18,0.27, 0.30, 0.27, 0.24 },
749 
750 // p p pi+ pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0)
751  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
752  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
753  0.0, 0.003,0.01, 0.024,0.075,0.18,0.27, 0.30, 0.27, 0.24 },
754 
755 // p p pi+ pi- pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0)
756  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
757  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
758  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
759 
760 // p p pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0)
761  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
762  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
763  0.0, 0.0, 0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
764 
765 // p n pi+ pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi- pi0)
766  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
767  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
768  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
769 
770 // p n pi+ pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0)
771  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
772  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
773  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
774 
775 // p n pi+ pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0)
776  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
777  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
778  0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
779 
780 // n n pi+ pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi- pi-)
781  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
782  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
783  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
784 
785 // n n pi+ pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0)
786  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
787  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
788  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
789 
790 // n n pi+ pi+ pi0 pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0 pi0)
791  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
792  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
793  0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
794 //
795 // multiplicity 9 (11 channels)
796 //
797 // p p pi+ pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi- pi0)
798  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
799  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
800  0.0, 0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
801 
802 // p p pi+ pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0 pi0)
803  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
804  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
805  0.0, 0.001,0.006,0.020,0.058,0.091,0.11, 0.12, 0.12, 0.12 },
806 
807 // p p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0 pi0)
808  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
809  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
810  0.0, 0.001,0.004,0.012,0.035,0.055,0.065,0.07, 0.07, 0.07 },
811 
812 // p p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
813  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
814  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
815  0.0, 0.0, 0.0, 0.006,0.027,0.032,0.04, 0.042,0.042,0.042 },
816 
817 // p n pi+ pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi- pi-)
818  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
819  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
820  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
821 
822 // p n pi+ pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi- pi0 pi0)
823  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
824  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
825  0.0, 0.003,0.008,0.026,0.078,0.20, 0.25, 0.29, 0.29, 0.29 },
826 
827 // p n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0 pi0)
828  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
829  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
830  0.0, 0.001,0.005,0.016,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
831 
832 // p n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0 pi0)
833  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
834  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
835  0.0, 0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 },
836 
837 // n n pi+ pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi- pi0)
838  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
839  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
840  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
841 
842 // n n pi+ pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0 pi0)
843  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
844  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
845  0.0, 0.001,0.005,0.015,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
846 
847 // n n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p p pi+ pi- pi0 pi0 pi0 pi0 pi0)
848  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
849  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
850  0.0, 0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 }};
851 
852 ///////////////////////////////////////////////////////////////////////////////
853 // //
854 // p n and n p (|Tz| = 0) cross sections //
855 // and final state particle types //
856 // //
857 ///////////////////////////////////////////////////////////////////////////////
858 
859 // Total p n cross section as a function of kinetic energy
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.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0,
960  44.0, 42.0, 40.0, 35.0, 31.0, 27.0, 23.0, 19.0, 17.0, 15.5,
961  14.0, 13.0, 12.0, 11.0, 10.0, 9.5, 9.0, 8.5, 8.0, 7.7 },
962 //
963 // multiplicity 3 (9 channels)
964 //
965 // p p pi- (n n pi+)
966  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
967  0.0, 0.0, 0.0, 0.0, 0.25, 0.9, 1.75, 2.3, 2.8, 2.8,
968  2.2, 1.9, 1.6, 1.35, 1.1, 0.95, 0.8, 0.7, 0.6, 0.53 },
969 
970 // p n pi0 (p n pi0)
971  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
972  0.0, 0.0, 0.0, 0.0, 1.8, 4.7, 8.3, 11.3, 12.0, 10.2,
973  8.2, 6.0, 4.9, 3.6, 2.5, 2.0, 1.6, 1.2, 1.0, 0.08 },
974 
975 // n n pi+ (p p pi-)
976  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
977  0.0, 0.0, 0.0, 0.0, 0.95, 2.4, 4.2, 5.6, 6.1, 5.1,
978  4.1, 3.0, 2.5, 1.8, 1.2, 1.0, 0.8, 0.6, 0.5, 0.41},
979 
980 // p L K0 (n L K+)
981  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
982  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
983  0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
984 
985 // p S0 K0 (n S0 K+)
986  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
987  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
988  0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
989 
990 // p S- K+ (n S+ K0)
991  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
992  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
993  0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
994 
995 // n L K+ (p L K0)
996  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
997  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
998  0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
999 
1000 // n S0 K+ (p S0 K0)
1001  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1002  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1003  0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
1004 
1005 // n S+ K0 (p S- K+)
1006  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1007  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1008  0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
1009 //
1010 // multiplicity 4 (22 channels)
1011 //
1012 // p n pi+ pi- (p n pi+ pi-)
1013  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1014  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12, 0.38, 1.1, 3.5,
1015  5.9, 5.9, 5.1, 4.2, 3.7, 3.0, 2.6, 2.1, 1.8, 1.4 },
1016 
1017 // p p pi- pi0 (n n pi+ pi0)
1018  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1019  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.1, 0.24, 0.55,
1020  1.2, 1.5, 1.45, 1.25, 1.0, 0.9, 0.8, 0.7, 0.6, 0.53 },
1021 
1022 // p n pi0 pi0 (p n pi0 pi0)
1023  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1024  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.07, 0.24, 0.66, 2.1,
1025  3.6, 3.6, 3.1, 2.5, 2.2, 1.8, 1.5, 1.2, 1.1, 0.84 },
1026 
1027 // n n pi+ pi0 (p p pi- pi0)
1028  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1029  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.1, 0.24, 0.55,
1030  1.2, 1.5, 1.45, 1.25, 1.0, 0.9, 0.8, 0.7, 0.6, 0.53 },
1031 
1032 // p L K+ pi- (n L K0 pi+)
1033  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1034  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1035  0.0, 0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
1036 
1037 // p S0 K+ pi- (n S0 K0 pi+)
1038  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1039  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1040  0.0, 0.0, 0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
1041 
1042 // p L K0 pi0 (n L K+ pi0)
1043  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1044  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1045  0.0, 0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
1046 
1047 // p S0 K0 pi0 (n S0 K+ pi0)
1048  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1049  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1050  0.0, 0.0, 0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
1051 
1052 // p S+ K0 pi- (n S- K+ pi+)
1053  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1054  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1055  0.0, 0.0015,0.004,0.013,0.02,0.016,0.013,0.01, 0.009,0.007 },
1056 
1057 // p S- K+ pi0 (n S+ K0 pi0)
1058  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1059  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1060  0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1061 
1062 // p S- K0 pi+ (n S+ K+ pi-)
1063  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1064  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1065  0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1066 
1067 // n L K+ pi0 (p L K0 pi0)
1068  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1069  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1070  0.0, 0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
1071 
1072 // n L K0 pi+ (p L K+ pi-)
1073  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1074  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1075  0.0, 0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
1076 
1077 // n S+ K+ pi- (p S- K0 pi+)
1078  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1079  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1080  0.0, 0.0, 0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
1081 
1082 // n S+ K0 pi0 (p S- K+ pi0)
1083  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1084  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1085  0.0, 0.0, 0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
1086 
1087 // n S0 K+ pi0 (p S0 K0 pi0)
1088  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1089  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1090  0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1091 
1092 // n S0 K0 pi+ (p S0 K+ pi-)
1093  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1094  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1095  0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1096 
1097 // n S- K+ pi+ (p S+ K0 pi-)
1098  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1099  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1100  0.0, 0.0015,0.004,0.013,0.02, 0.016,0.013,0.01, 0.009,0.007 },
1101 
1102 // p n K+ K- (p n K0 K0bar)
1103  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1104  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1105  0.0, 0.0, 0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1106 
1107 // p n K0 K0bar (p n K+ K-)
1108  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1109  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1110  0.0, 0.0, 0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1111 
1112 // p p K0 K- (n n K+ K0bar)
1113  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1114  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1115  0.0, 0.0, 0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1116 
1117 // n n K+ K0bar (p p K0 K-)
1118  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1119  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1120  0.0, 0.0, 0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1121 //
1122 // multiplicity 5 (38 channels)
1123 //
1124 // p n pi+ pi- pi0 (p n pi+ pi- pi0)
1125  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1126  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
1127  0.3, 0.82, 1.35, 1.8, 1.8, 1.65, 1.5, 1.28, 1.12, 0.98 },
1128 
1129 // p n pi0 pi0 pi0 (p n pi0 pi0 pi0)
1130  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1131  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.004,0.022,
1132  0.15, 0.41, 0.68, 0.9, 0.9, 0.82, 0.75, 0.64, 0.55, 0.49 },
1133 
1134 // p p pi+ pi- pi- (n n pi+ pi+ pi-)
1135  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1136  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09, 0.2,
1137  0.52, 1.2, 1.8, 2.0, 1.7, 1.5, 1.35, 1.2, 1.05, 0.9 },
1138 
1139 // p p pi- pi0 pi0 (n n pi+ pi0 pi0)
1140  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1141  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.1,
1142  0.26, 0.6, 0.9, 0.98, 0.86, 0.75, 0.68, 0.6, 0.52, 0.45 },
1143 
1144 // n n pi+ pi+ pi- (p p pi+ pi- pi-)
1145  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1146  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
1147  0.3, 0.82, 1.35, 1.8, 1.8, 1.65, 1.5, 1.28, 1.12, 0.98 },
1148 
1149 // n n pi+ pi0 pi0 (p p pi- pi0 pi0)
1150  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1151  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.003,0.02,
1152  0.15, 0.41, 0.68, 0.9, 0.9, 0.82, 0.75, 0.64, 0.56, 0.49 },
1153 
1154 // p L K+ pi- pi0 (n L K0 pi+ pi0)
1155  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1156  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1157  0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1158 
1159 // p L K0 pi+ pi- (n L K+ pi+ pi-)
1160  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1161  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1162  0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1163 
1164 // p L K0 pi0 pi0 (n L K+ pi0 pi0)
1165  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1166  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1167  0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.09, 0.07, 0.07 },
1168 
1169 // p S0 K0 pi+ pi- (n S0 K+ pi+ pi-)
1170  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1171  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1172  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1173 
1174 // p S0 K0 pi0 pi0 (n S0 K+ pi0 pi0)
1175  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1176  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1177  0.0, 0.0, 0.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1178 
1179 // p S0 K+ pi- pi0 (n S0 K0 pi+ pi0)
1180  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1181  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1182  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1183 
1184 // p S+ K+ pi- pi- (n S- K0 pi+ pi+)
1185  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1186  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1187  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1188 
1189 // p S+ K0 pi- pi0 (n S- K+ pi+ pi0)
1190  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1191  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1192  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1193 
1194 // p S- K0 pi+ pi0 (n S+ K+ pi- pi0)
1195  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1196  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1197  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1198 
1199 // p S- K+ pi+ pi- (n S+ K0 pi+ pi-)
1200  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1201  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1202  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1203 
1204 // p S- K+ pi0 pi0 (n S+ K0 pi0 pi0)
1205  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1206  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1207  0.0, 0.0, 0.002,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1208 
1209 // n L K+ pi+ pi- (p L K0 pi+ pi-)
1210  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1211  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1212  0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1213 
1214 // n L K+ pi0 pi0 (p L K0 pi0 pi0)
1215  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1216  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1217  0.0, 0.0, 0.001,0.004,0.011,0.012,0.011,0.009,0.007,0.006 },
1218 
1219 // n L K0 pi+ pi0 (p L K+ pi- pi0)
1220  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1221  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1222  0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1223 
1224 // n S0 K+ pi+ pi- (p S0 K0 pi+ pi-)
1225  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1226  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1227  0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
1228 
1229 // n S0 K+ pi0 pi0 (p S0 K0 pi0 pi0)
1230  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1231  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1232  0.0, 0.0, 0.0, 0.002,0.005,0.006,0.005,0.005,0.004,0.003 },
1233 
1234 // n S0 K0 pi+ pi0 (p S0 K+ pi- pi0)
1235  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1236  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1237  0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
1238 
1239 // n S+ K0 pi+ pi- (p S- K+ pi+ pi-)
1240  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1241  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1242  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1243 
1244 // n S+ K0 pi0 pi0 (p S- K+ pi0 pi0)
1245  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1246  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1247  0.0, 0.0, 0.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1248 
1249 // n S+ K+ pi- pi0 (p S- K0 pi+ pi0)
1250  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1251  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1252  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1253 
1254 // n S- K+ pi+ pi0 (p S+ K0 pi- pi0)
1255  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1256  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1257  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1258 
1259 // n S- K0 pi+ pi+ (p S+ K+ pi- pi-)
1260  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1261  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1262  0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1263 
1264 // p n K+ K- pi0 (p n K0 K0bar pi0)
1265  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1266  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1267  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1268 
1269 // p n K0 K0bar pi0 (p n K+ K- pi0)
1270  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1271  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1272  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1273 
1274 // p n K0 K- pi+ (p n K+ K0bar pi-)
1275  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1276  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1277  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1278 
1279 // p n K+ K0bar pi- (p n K0 K- pi+)
1280  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1281  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1282  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1283 
1284 // p p K0 K0bar pi- (n n K+ K- pi+)
1285  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1286  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1287  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1288 
1289 // p p K+ K- pi- (n n K0 K0bar pi+)
1290  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1291  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1292  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1293 
1294 // p p K0 K- pi0 (n n K+ K0bar pi0)
1295  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1296  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1297  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1298 
1299 // n n K+ K- pi+ (p p K0 K0bar pi-)
1300  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1301  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1302  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1303 
1304 // n n K0 K0bar pi+ (p p K+ K- pi-)
1305  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1306  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1307  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1308 
1309 // n n K+ K0bar pi0 (p p K0 K- pi0)
1310  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1311  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1312  0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1313 //
1314 // multiplicity 6 (7 channels)
1315 //
1316 // p n pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi-)
1317  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1318  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1319  0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1320 
1321 // p n pi+ pi- pi0 pi0 (p n pi+ pi- pi0 pi0)
1322  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1323  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1324  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1325 
1326 // p n pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0)
1327  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1328  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1329  0.01, 0.02, 0.05, 0.1, 0.13, 0.12, 0.11, 0.1, 0.1, 0.09 },
1330 
1331 // p p pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi0)
1332  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1333  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1334  0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1335 
1336 // p p pi- pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0)
1337  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1338  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1339  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1340 
1341 // n n pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi0)
1342  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1343  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1344  0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1345 
1346 // n n pi+ pi0 pi0 pi0 (p p pi- pi0 pi0 pi0)
1347  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1348  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1349  0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1350 //
1351 // multiplicity 7 (9 channels)
1352 //
1353 // p n pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi0)
1354  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1355  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1356  0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
1357 
1358 // p n pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0)
1359  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1360  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1361  0.0, 0.0, 0.03, 0.08, 0.25, 0.35, 0.35, 0.35, 0.33, 0.31 },
1362 
1363 // p n pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0)
1364  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1365  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1366  0.0, 0.0, 0.015,0.04, 0.12, 0.17, 0.18, 0.17, 0.16, 0.15 },
1367 
1368 // p p pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi-)
1369  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1370  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1371  0.0, 0.0, 0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
1372 
1373 // p p pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi0 pi0)
1374  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1375  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1376  0.0, 0.0, 0.03, 0.1, 0.15, 0.2, 0.22, 0.23, 0.22, 0.22 },
1377 
1378 // p p pi- pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0)
1379  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1380  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1381  0.0, 0.0, 0.02, 0.05, 0.07, 0.1, 0.11, 0.12, 0.11, 0.11 },
1382 
1383 // n n pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi-)
1384  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1385  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1386  0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
1387 
1388 // n n pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi0 pi0)
1389  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1390  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1391  0.0, 0.0, 0.03, 0.08, 0.25, 0.35, 0.35, 0.34, 0.33, 0.31 },
1392 
1393 // n n pi+ pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0)
1394  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1395  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1396  0.0, 0.0, 0.02, 0.05, 0.07, 0.1, 0.11, 0.12, 0.11, 0.11 },
1397 //
1398 // multiplicity 8 (10 channels)
1399 //
1400 // p n pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi-)
1401  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1402  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1403  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1404 
1405 // p n pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0)
1406  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1407  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1408  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1409 
1410 // p n pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0)
1411  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1412  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1413  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1414 
1415 // p n pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0)
1416  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1417  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1418  0.0, 0.0, 0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
1419 
1420 // p p pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi0)
1421  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1422  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1423  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1424 
1425 // p p pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0)
1426  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1427  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1428  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1429 
1430 // p p pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0)
1431  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1432  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1433  0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
1434 
1435 // n n pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi0)
1436  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1437  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1438  0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1439 
1440 // n n pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0)
1441  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1442  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1443  0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1444 
1445 // n n pi+ pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0)
1446  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1447  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1448  0.0, 0.0, 0.003,0.009,0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
1449 //
1450 // multiplicity 9 (12 channels)
1451 //
1452 // p n pi+ pi+ pi+ pi- pi- pi- pi0 (p n pi+ pi+ pi+ pi- pi- pi- pi0)
1453  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1454  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1455  0.0, 0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
1456 
1457 // p n pi+ pi+ pi- pi- pi0 pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0 pi0)
1458  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1459  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1460  0.0, 0.001,0.005,0.015,0.045,0.07, 0.084,0.09, 0.09, 0.09 },
1461 
1462 // p n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0 pi0)
1463  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1464  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1465  0.0, 0.001,0.003,0.009,0.027,0.042,0.05, 0.054,0.054,0.054 },
1466 
1467 // p n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
1468  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1469  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1470  0.0, 0.0, 0.0, 0.005,0.016,0.025,0.03, 0.032,0.032,0.032 },
1471 
1472 // p p pi+ pi+ pi+ pi- pi- pi- pi- (n n pi+ pi+ pi+ pi+ pi- pi- pi-)
1473  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1474  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1475  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1476 
1477 // p p pi+ pi+ pi- pi- pi- pi0 pi0 (n n pi+ pi+ pi+ pi- pi- pi0 pi0)
1478  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1479  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1480  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1481 
1482 // p p pi+ pi- pi- pi0 pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0 pi0)
1483  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1484  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1485  0.0, 0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
1486 
1487 // p p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
1488  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1489  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1490  0.0, 0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 },
1491 
1492 // n n pi+ pi+ pi+ pi+ pi- pi- pi- (p p pi+ pi+ pi+ pi- pi- pi- pi-)
1493  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1494  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1495  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1496 
1497 // n n pi+ pi+ pi+ pi- pi- pi0 pi0 (p p pi+ pi+ pi- pi- pi- pi0 pi0)
1498  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1499  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1500  0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1501 
1502 // n n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0 pi0)
1503  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1504  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1505  0.0, 0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
1506 
1507 // n n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0 pi0)
1508  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1509  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1510  0.0, 0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 }};
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:523
G4RPGNucleonInelastic::T0_7bfs
static const G4int T0_7bfs[9][7]
Definition: G4RPGNucleonInelastic.hh:92
G4ThreadLocal
#define G4ThreadLocal
Definition: tls.hh:52
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
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:506
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
Generated on Wed Apr 30 2014 15:55:03 for Geant4.10 by   doxygen 1.8.7