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

Last change on this file since 1007 was 1007, checked in by garnier, 15 years ago

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