source: trunk/source/processes/hadronic/models/rpg/src/G4RPGNucleonInelastic.cc @ 1340

Last change on this file since 1340 was 1340, checked in by garnier, 14 years ago

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