source: trunk/source/processes/hadronic/models/rpg/src/G4RPGPionInelastic.cc @ 982

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

fichier ajoutes
File size: 58.1 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: G4RPGPionInelastic.cc,v 1.1 2008/01/05 00:15:12 dennis Exp $
27// GEANT4 tag $Name: geant4-09-02-ref-02 $
28//
29 
30#include "G4RPGPionInelastic.hh"
31#include "Randomize.hh"
32
33G4RPGPionInelastic::G4RPGPionInelastic(const G4String& modelName)
34 :G4RPGInelastic(modelName)
35{
36  SetMinEnergy( 0.0 );
37  SetMaxEnergy( 30.*GeV );
38
39  // Initialize t32_dSigma_dMult, t12_dSigma_dMult,
40  //   pi-nucleon inelastic cross sections for a given multiplicity
41  //   for |T_z| = 3/2 and 1/2, respectively
42
43  G4int i, k, m;
44  G4int start, stop;
45
46  for (m = 0; m < 8; m++) {
47    start = pipPindex[m][0];
48    stop = pipPindex[m][1] + 1;
49    for (k = 0; k < 30; k++) {
50      t32_dSigma_dMult[m][k] = 0.0;
51      for (i = start; i < stop; i++) t32_dSigma_dMult[m][k] += pipPCrossSections[i][k];
52    }
53
54    start = pimPindex[m][0];
55    stop = pimPindex[m][1] + 1;
56    for (k = 0; k < 30; k++) {
57      t12_dSigma_dMult[m][k] = 0.0;
58      for (i = start; i < stop; i++) t12_dSigma_dMult[m][k] += pimPCrossSections[i][k];
59    }
60  }
61
62  // Initialize total cross section array
63
64  for (k = 0; k < 30; k++) {
65    pipPtot[k] = 0.0;
66    pimPtot[k] = 0.0;
67    for (m = 0; m < 8; m++) {
68      pipPtot[k] += t32_dSigma_dMult[m][k];
69      pimPtot[k] += t12_dSigma_dMult[m][k];
70    }
71  }
72
73}
74
75
76G4int G4RPGPionInelastic::GetMultiplicityT12(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 = t12_dSigma_dMult[m][k]
87         + fraction*(t12_dSigma_dMult[m][k+1] - t12_dSigma_dMult[m][k]);
88      sigma.push_back(multint);
89  }
90
91  return sampleFlat(sigma) + 2;
92}
93
94
95G4int G4RPGPionInelastic::GetMultiplicityT32(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 = t32_dSigma_dMult[m][k]
106         + fraction*(t32_dSigma_dMult[m][k+1] - t32_dSigma_dMult[m][k]);
107      sigma.push_back(multint);
108  }
109
110  return sampleFlat(sigma) + 2;
111}
112
113
114std::vector<G4int> 
115G4RPGPionInelastic::GetFSPartTypesForT32(G4int mult, G4double KE, G4int tzindex) 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 = pipPindex[mult-2][0];
126  G4int stop = pipPindex[mult-2][1];
127
128  for(i = start; i < stop; i++) {
129      sigint = pipPCrossSections[i][k]
130          + fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[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(T32_2bfs[tzindex][channel][i]);
140  } else if (mult == 3) {
141    for(i = 0; i < mult; i++) kinds.push_back(T32_3bfs[tzindex][channel][i]);
142  } else if (mult == 4) {
143    for(i = 0; i < mult; i++) kinds.push_back(T32_4bfs[tzindex][channel][i]);
144  } else if (mult == 5) {
145    for(i = 0; i < mult; i++) kinds.push_back(T32_5bfs[tzindex][channel][i]);
146  } else if (mult == 6) {
147    for(i = 0; i < mult; i++) kinds.push_back(T32_6bfs[tzindex][channel][i]);
148  } else if (mult == 7) {
149    for(i = 0; i < mult; i++) kinds.push_back(T32_7bfs[tzindex][channel][i]);
150  } else if (mult == 8) {
151    for(i = 0; i < mult; i++) kinds.push_back(T32_8bfs[tzindex][channel][i]);
152  } else if (mult == 9) {
153    for(i = 0; i < mult; i++) kinds.push_back(T32_9bfs[tzindex][channel][i]);
154  } else {
155    G4cout << " Illegal multiplicity " << G4endl;
156  }
157
158  return kinds;
159}
160
161
162std::vector<G4int> 
163G4RPGPionInelastic::GetFSPartTypesForT12(G4int mult, G4double KE, G4int tzindex) const
164{
165  G4int i;
166  G4double sigint(0.);
167  std::vector<G4double> sigma;
168
169  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
170  G4int k = epair.first;
171  G4double fraction = epair.second;
172
173  G4int start = pimPindex[mult-2][0];
174  G4int stop = pimPindex[mult-2][1];
175
176  for(i = start; i < stop; i++) {
177      sigint = pimPCrossSections[i][k]
178          + fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
179      sigma.push_back(sigint);
180  }
181
182  G4int channel = sampleFlat(sigma);
183
184  std::vector<G4int> kinds;
185
186  if (mult == 2) {
187    for(i = 0; i < mult; i++) kinds.push_back(T12_2bfs[tzindex][channel][i]);
188  } else if (mult == 3) {
189    for(i = 0; i < mult; i++) kinds.push_back(T12_3bfs[tzindex][channel][i]);
190  } else if (mult == 4) {
191    for(i = 0; i < mult; i++) kinds.push_back(T12_4bfs[tzindex][channel][i]);
192  } else if (mult == 5) {
193    for(i = 0; i < mult; i++) kinds.push_back(T12_5bfs[tzindex][channel][i]);
194  } else if (mult == 6) {
195    for(i = 0; i < mult; i++) kinds.push_back(T12_6bfs[tzindex][channel][i]);
196  } else if (mult == 7) {
197    for(i = 0; i < mult; i++) kinds.push_back(T12_7bfs[tzindex][channel][i]);
198  } else if (mult == 8) {
199    for(i = 0; i < mult; i++) kinds.push_back(T12_8bfs[tzindex][channel][i]);
200  } else if (mult == 9) {
201    for(i = 0; i < mult; i++) kinds.push_back(T12_9bfs[tzindex][channel][i]);
202  } else {
203    G4cout << " Illegal multiplicity " << G4endl;
204  }
205
206  return kinds;
207}
208
209
210///////////////////////////////////////////////////////////////////////////////
211//                                                                           //
212//   pi+ p and pi- n (|Tz| = 3/2) cross sections                             //
213//   and final state particle types                                          //
214//                                                                           //
215///////////////////////////////////////////////////////////////////////////////
216
217// Total pi+ cross section as a function of kinetic energy
218G4double G4RPGPionInelastic::pipPtot[30];
219
220// pi+ multiplicities as a function of kinetic energy
221G4double G4RPGPionInelastic::t32_dSigma_dMult[8][30];
222
223const G4int G4RPGPionInelastic::pipPindex[8][2] =
224 {{0, 1}, {2, 8}, {9,23}, {24,47}, {48,52}, {53,58}, {59,65}, {66,73}}; 
225
226// Outgoing particle types of a given multiplicity
227// T32_nbfs = final state types for pi+ p and pi- n
228
229const G4int G4RPGPionInelastic::T32_2bfs[2][2][2] =
230  {{{pro,pip}, {sp,kp}},
231
232   {{neu,pim}, {sm,k0}}};
233
234const G4int G4RPGPionInelastic::T32_3bfs[2][7][3] =
235  {{{pro,pip,pi0}, {neu,pip,pip}, {sp,kp,pi0}, {sp,k0,pip}, 
236    {s0,kp,pip},   {lam,kp,pip},  {pro,kp,k0b}},
237
238   {{neu,pim,pi0}, {pro,pim,pim}, {sm,k0,pi0}, {sm,kp,pim},
239    {s0,k0,pim},   {lam,k0,pim},  {neu,k0,km}}};
240
241const G4int G4RPGPionInelastic::T32_4bfs[2][15][4] =
242  {{{pro,pip,pip,pim},{pro,pip,pi0,pi0},{neu,pip,pip,pi0},
243    {sp,kp,pip,pim},  {sp,kp,pi0,pi0},  {sp,k0,pip,pi0},
244    {s0,k0,pip,pip},  {s0,kp,pip,pi0},  {lam,kp,pip,pi0},
245    {lam,k0,pip,pip}, {sm,kp,pip,pip},  {pro,pip,kp,km},
246    {pro,pip,k0,k0b}, {pro,pi0,kp,k0b}, {neu,pip,kp,k0b}},
247
248   {{neu,pip,pim,pim},{neu,pim,pi0,pi0},{pro,pim,pim,pi0},
249    {sm,k0,pip,pim},  {sm,k0,pi0,pi0},  {sm,kp,pim,pi0},
250    {s0,kp,pim,pim},  {s0,k0,pim,pi0},  {lam,k0,pim,pi0},
251    {lam,kp,pim,pim}, {sp,k0,pim,pim},  {neu,pim,k0,k0b},
252    {neu,pim,kp,km},  {neu,pi0,k0,km},  {pro,pim,k0,km}}};
253
254const G4int G4RPGPionInelastic::T32_5bfs[2][24][5] =
255  {{{pro,pip,pip,pim,pi0}, {pro,pip,pi0,pi0,pi0}, {neu,pip,pip,pip,pim},
256    {neu,pip,pip,pi0,pi0}, {sp,kp,pip,pim,pi0},   {sp,kp,pi0,pi0,pi0},
257    {sp,k0,pip,pip,pim},   {sp,k0,pip,pi0,pi0},   {lam,k0,pip,pip,pi0},
258    {lam,kp,pip,pip,pim},  {lam,kp,pip,pi0,pi0},  {s0,kp,pip,pip,pim},
259    {s0,kp,pip,pi0,pi0},   {s0,k0,pip,pip,pi0},   {sm,kp,pip,pip,pi0},
260    {sm,k0,pip,pip,pip},   {pro,pip,pim,kp,k0b},  {pro,pip,pip,k0,km},
261    {pro,pip,pi0,kp,km},   {pro,pip,pi0,k0,k0b},  {pro,pi0,pi0,kp,k0b},
262    {neu,pip,pip,kp,km},   {neu,pip,pip,k0,k0b},  {neu,pip,pi0,kp,k0b}},
263
264   {{neu,pip,pim,pim,pi0}, {neu,pim,pi0,pi0,pi0}, {pro,pip,pim,pim,pim},
265    {pro,pim,pim,pi0,pi0}, {sm,k0,pip,pim,pi0},   {sm,k0,pi0,pi0,pi0},
266    {sm,kp,pip,pim,pim},   {sm,kp,pim,pi0,pi0},   {lam,kp,pim,pim,pi0},
267    {lam,k0,pip,pim,pim},  {lam,k0,pim,pi0,pi0},  {s0,k0,pip,pim,pim},
268    {s0,k0,pim,pi0,pi0},   {s0,kp,pim,pim,pi0},   {sp,k0,pim,pim,pi0},
269    {sp,kp,pim,pim,pim},   {neu,pip,pim,k0,km},   {neu,pim,pim,kp,k0b},
270    {neu,pim,pi0,k0,k0b},  {neu,pim,pi0,kp,km},   {neu,pi0,pi0,k0,km},
271    {pro,pim,pim,k0,k0b},  {pro,pim,pim,kp,km},   {pro,pim,pi0,k0,km}}};
272
273const G4int G4RPGPionInelastic::T32_6bfs[2][5][6] =
274{{{pro,pip,pip,pip,pim,pim}, {pro,pip,pip,pim,pi0,pi0},
275  {pro,pip,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0},
276  {neu,pip,pip,pip,pim,pi0}},
277
278 {{neu,pip,pip,pim,pim,pim}, {neu,pip,pim,pim,pi0,pi0},
279  {neu,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0},
280  {pro,pip,pim,pim,pim,pi0}}};
281
282const G4int G4RPGPionInelastic::T32_7bfs[2][6][7] =
283{{{pro,pip,pip,pip,pim,pim,pi0}, {pro,pip,pip,pim,pi0,pi0,pi0},
284  {pro,pip,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim},
285  {neu,pip,pip,pip,pim,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0}},
286
287 {{neu,pip,pip,pim,pim,pim,pi0}, {neu,pip,pim,pim,pi0,pi0,pi0},
288  {neu,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pim},
289  {pro,pip,pim,pim,pim,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0}}};
290
291const G4int G4RPGPionInelastic::T32_8bfs[2][7][8] =
292{{{pro,pip,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pip,pim,pim,pi0,pi0},
293  {pro,pip,pip,pim,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0},
294  {neu,pip,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pip,pim,pi0,pi0,pi0},
295  {neu,pip,pip,pi0,pi0,pi0,pi0,pi0}},
296
297 {{neu,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pip,pim,pim,pim,pi0,pi0},
298  {neu,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
299  {pro,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pim,pim,pim,pi0,pi0,pi0},
300  {pro,pim,pim,pi0,pi0,pi0,pi0,pi0}}};
301
302const G4int G4RPGPionInelastic::T32_9bfs[2][8][9] =
303{{{pro,pip,pip,pip,pip,pim,pim,pim,pi0}, {pro,pip,pip,pip,pim,pim,pi0,pi0,pi0},
304  {pro,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
305  {neu,pip,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pip,pim,pim,pi0,pi0},
306  {neu,pip,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
307
308 {{neu,pip,pip,pip,pim,pim,pim,pim,pi0}, {neu,pip,pip,pim,pim,pim,pi0,pi0,pi0},
309  {neu,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
310  {pro,pip,pip,pip,pim,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pim,pi0,pi0},
311  {pro,pip,pim,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
312
313//
314// Cross sections (in mb) for pi+ p -> 2-9 body final states
315//
316// first index:    0-1: channels for mult = 2
317//                 2-8: channels for mult = 3
318//                9-23: channels for mult = 4
319//               24-47: channels for mult = 5
320//               48-52: channels for mult = 6
321//               53-58: channels for mult = 7
322//               59-65: channels for mult = 8
323//               66-73: channels for mult = 9
324//
325// second index: kinetic energy
326//
327
328const G4float G4RPGPionInelastic::pipPCrossSections[74][30] = {
329//
330// multiplicity 2 (2 channels)
331//
332//  p pi+ (n pi-)
333 { 0.00, 1.20, 2.50, 3.80, 5.00, 7.00, 9.00, 15.0, 30.0, 64.0,
334  130.0,190.0,130.0, 55.7, 27.2, 14.0, 8.50, 13.0, 18.0, 11.0,
335   8.50, 7.00, 6.20, 5.60, 5.00, 4.50, 4.20, 4.00, 3.80, 3.60 },
336
337//  S+ K+ (S- K0)
338 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
339    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.16, 0.60, 0.32,
340   0.19,  0.1, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
341//
342// multiplicity 3 (7 channels)
343//
344//  p pi+ pi0 (n pi- pi0)
345 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
346    0.0,  0.0,  0.0,  0.2,  0.6,  2.4,  8.8, 10.0, 12.0,  6.2,
347   4.00, 2.40, 1.69, 1.10, 0.73, 0.49, 0.41, 0.31, 0.24, 0.15 },
348
349//  n pi+ pi+ (p pi- pi-)
350 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
351    0.0,  0.0,  0.0, 0.10, 0.20, 0.60, 1.50, 2.30, 3.60, 3.00,
352   2.30, 1.70, 1.30, 0.95, 0.69, 0.46, 0.38, 0.27, 0.20, 0.15 },
353
354//  S+ K+ pi0 (S- K0 pi0)
355 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
356    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
357   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
358
359//  S+ K0 pi+
360 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
361    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
362   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
363
364//  S0 K+ pi+
365 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
366    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
367   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
368
369//  L K+ pi+
370 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
371    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
372   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
373
374//  p K+ K0bar
375 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
376    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
377   0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
378//
379// multiplicity 4 (15 channels)
380//
381//  p pi+ pi+ pi-
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, 0.06, 0.20, 0.78, 2.20, 3.20,
384   3.50, 3.10, 2.70, 2.30, 2.00, 1.50, 1.40, 1.20, 1.00, 0.90 },
385
386//  p pi+ pi0 pi0
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.0, 0.04, 0.13, 0.52, 1.50, 2.20,
389   2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
390
391//  n pi+ pi+ pi0
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.04, 0.13, 0.52, 1.50, 2.20,
394   2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
395
396//  S+ K+ pi+ pi-
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.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
400
401//  S+ K+ pi0 pi0
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.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
405
406//  S+ K0 pi+ pi0
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.002,
409   0.015,0.06, 0.06, 0.05, 0.04,0.032,0.028, 0.02,0.017,0.014 },
410
411//  S0 K0 pi+ pi+
412 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
413    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
414   0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
415
416//  S0 K+ pi+ pi0
417 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
418    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
419   0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
420
421//  L K+ pi+ pi0
422 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
423    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
424   0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
425
426//  L K0 pi+ pi+
427 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
428    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
429   0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
430
431//  S- K+ pi+ pi+
432 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
433    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,
434   0.02, 0.16, 0.18, 0.13, 0.09, 0.06,0.045, 0.03,0.025, 0.02 },
435
436//  p pi+ K+ K-
437 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
438    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
439   0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05,0.045, 0.04 },
440
441//  p pi+ K0 K0bar
442 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
443    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
444  0.005, 0.05, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05,0.045, 0.04 },
445
446//  p pi0 K+ K0bar
447 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
448    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
449   0.01, 0.05,0.075, 0.08,0.075, 0.07, 0.06, 0.05,0.045, 0.04 },
450
451//  n pi+ K+ K0bar
452 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
453    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
454   0.01, 0.05, 0.07,0.065,0.055,0.048, 0.04, 0.03,0.027,0.023 },
455//
456// multiplicity 5 (24 channels)
457//
458//  p pi+ pi+ pi- pi0
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.001, 0.05, 0.30, 2.00,
461   3.20, 3.70, 3.00, 2.50, 2.10, 1.60, 1.40, 1.10, 0.89, 0.70 },
462
463//  p pi+ pi0 pi0 pi0
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.001, 0.03, 0.21, 1.40,
466   2.20, 2.60, 2.10, 1.80, 1.50, 1.10, 1.00, 0.80, 0.62, 0.50 },
467
468//  n pi+ pi+ pi+ pi-
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.007, 0.02, 0.05, 0.19,
471   0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
472
473//  n pi+ pi+ pi0 pi0
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.007, 0.02, 0.05, 0.19,
476   0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
477
478//  S+ K+ pi+ pi- pi0
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.003,0.013, 0.05, 0.09, 0.08, 0.06, 0.05,  0.04,0.036,0.03 },
482
483//  S+ K+ pi0 pi0 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.001,0.006, 0.02, 0.04, 0.04, 0.03,0.025, 0.02, 0.02, 0.01 },
487
488//  S+ K0 pi+ pi+ pi-
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.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
492
493//  S+ K0 pi+ pi0 pi0
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.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
497
498//  L K0 pi+ pi+ pi0
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.002,0.013, 0.04,0.052,0.059,0.053, 0.05,0.043,0.037, 0.03 },
502
503//  L K+ pi+ pi+ pi-
504 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
505    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
506  0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
507
508//  L K+ pi+ pi0 pi0
509 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
510    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
511  0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
512
513//  S0 K+ pi+ pi+ pi-
514 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
515    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
516  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
517
518//  S0 K+ pi+ pi0 pi0
519 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
520    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
521  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
522
523//  S0 K0 pi+ pi+ pi0
524 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
525    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
526  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
527
528//  S- K+ pi+ pi+ pi0
529 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
530    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
531  0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
532
533//  S- K0 pi+ pi+ pi+
534 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
535    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
536  0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
537
538//  p pi+ pi- K+ K0bar
539 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
540    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
541  0.001,0.005, 0.02,0.065, 0.08, 0.07, 0.06,0.054,0.048, 0.04 },
542
543//  p pi+ pi+ K0 K-
544 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
545    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
546  0.001,0.005, 0.02,0.045, 0.05,0.047, 0.04,0.033, 0.03,0.026 },
547
548//  p pi+ pi0 K+ K-
549 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
550    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
551  0.001, 0.02, 0.06, 0.06,0.054,0.048,0.042,0.038,0.035, 0.03 },
552
553//  p pi+ pi0 K0 K0bar
554 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
555    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
556  0.001, 0.01, 0.04, 0.05, 0.05, 0.04, 0.04,0.035,0.032, 0.03 },
557
558//  p pi0 pi0 K+ K0bar
559 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
560    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
561  0.001,0.005, 0.01, 0.03, 0.04,0.035, 0.03,0.025,0.025, 0.02 },
562
563//  n pi+ pi+ K+ K-
564 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
565    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
566    0.0, 0.02, 0.06, 0.05,0.042,0.038,0.035, 0.03,0.027,0.022 },
567
568//  n pi+ pi+ K0 K0bar
569 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
570    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
571    0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
572
573//  n pi+ pi0 K+ K0bar
574 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
575    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
576    0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
577//
578// multiplicity 6 (5 channels)
579//
580//  p pi+ pi+ pi+ 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.001, 0.02,
583   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
584
585//  p pi+ pi+ pi- pi0 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.001, 0.02,
588   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
589
590//  p pi+ pi0 pi0 pi0 pi0
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.001, 0.02,
593   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
594
595//  n pi+ pi+ pi0 pi0 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.001, 0.02,
598   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
599
600//  n pi+ pi+ pi+ pi- pi0
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.001, 0.02,
603   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
604//
605// multiplicity 7 (6 channels)
606//
607//  p pi+ pi+ pi+ pi- pi- pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
610  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
611
612//  p pi+ pi+ pi- pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
615  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
616
617//  p pi+ pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
620  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
621
622//  n pi+ pi+ pi+ pi+ pi- pi-
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
625   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
626
627//  n pi+ pi+ pi+ pi- pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
630   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
631
632//  n pi+ pi+ pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
635   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
636//
637// multiplicity 8 (7 channels)
638//
639//  p pi+ pi+ pi+ pi+ pi- pi- pi-
640 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
641    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
642  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
643
644//  p pi+ pi+ pi+ pi- pi- pi0 pi0
645 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
646    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
647  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
648
649//  p pi+ pi+ pi- pi0 pi0 pi0 pi0
650 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
651    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
652  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
653
654//  p pi+ pi0 pi0 pi0 pi0 pi0 pi0
655 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
656    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
657  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
658
659//  n pi+ pi+ pi+ pi+ pi- pi- pi0
660 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
661    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
662  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
663
664//  n pi+ pi+ pi+ pi- pi0 pi0 pi0
665 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
666    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
667  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
668
669//  n pi+ pi+ pi0 pi0 pi0 pi0 pi0
670 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
671    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
672  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
673//
674// multiplicity 9 (8 channels)
675//
676//  p pi+ pi+ pi+ pi+ pi- pi- pi- pi0
677 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
678    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
679  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
680
681//  p pi+ pi+ pi+ pi- pi- pi0 pi0 pi0
682 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
683    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
684  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
685
686//  p pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0
687 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
688    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
689  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
690
691//  p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
692 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
693    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
694  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
695
696//  n pi+ pi+ pi+ pi+ pi+ pi- pi- pi-
697 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
698    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
699    0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
700
701//  n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
704    0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
705
706//  n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
709    0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
710
711//  n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
714    0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 } };
715
716///////////////////////////////////////////////////////////////////////////////
717//                                                                           //
718//   pi- p and pi+ n (|Tz| = 1/2) cross sections                             //
719//   and final state particle types                                          //
720//                                                                           //
721///////////////////////////////////////////////////////////////////////////////
722
723// Total pi- cross section as a function of kinetic energy
724G4double G4RPGPionInelastic::pimPtot[30];
725
726// pi- multiplicities as a function of kinetic energy
727G4double G4RPGPionInelastic::t12_dSigma_dMult[8][30];
728
729const G4int G4RPGPionInelastic::pimPindex[8][2] =
730 {{0, 4}, {5,17}, {18,39}, {40,70}, {71,76}, {77,83}, {84,91}, {92,100}}; 
731
732// Outgoing particle types of a given multiplicity
733// T12_nbfs = final state types for pi- p and pi+ n
734
735const G4int G4RPGPionInelastic::T12_2bfs[2][5][2] =
736  {{{pro,pim}, {neu,pi0}, {lam,k0}, {s0,k0}, {sm,kp}},
737
738   {{neu,pip}, {pro,pi0}, {lam,kp}, {s0,kp}, {sp,k0}}};
739
740const G4int G4RPGPionInelastic::T12_3bfs[2][13][3] =
741  {{{pro,pim,pi0}, {neu,pip,pim}, {neu,pi0,pi0}, {lam,k0,pi0}, 
742    {lam,kp,pim},  {sm,k0,pip},   {sm,kp,pi0},   {sp,k0,pim},
743    {s0,kp,pim},   {s0,k0,pi0},   {pro,k0,km},   {neu,kp,km},
744    {neu,k0,k0b}},
745
746   {{neu,pip,pi0}, {pro,pip,pim}, {pro,pi0,pi0}, {lam,kp,pi0},
747    {lam,k0,pip},  {sp,kp,pim},   {sp,k0,pi0},   {sm,kp,pip},
748    {s0,k0,pip},   {s0,kp,pi0},   {neu,kp,k0b},  {pro,k0,k0b},
749    {pro,kp,km}}};
750
751const G4int G4RPGPionInelastic::T12_4bfs[2][22][4] =
752  {{{pro,pip,pim,pim}, {pro,pim,pi0,pi0}, {neu,pip,pim,pi0},
753    {neu,pi0,pi0,pi0}, {lam,k0,pip,pim},  {lam,k0,pi0,pi0},
754    {lam,kp,pim,pi0},  {s0,k0,pip,pim},   {s0,k0,pi0,pi0},
755    {s0,kp,pim,pi0},   {sp,kp,pim,pim},   {sp,k0,pim,pi0},
756    {sm,kp,pip,pim},   {sm,kp,pi0,pi0},   {sm,k0,pip,pi0},
757    {pro,pim,kp,km},   {pro,pim,k0,k0b},  {pro,pi0,k0,km},
758    {neu,pip,k0,km},   {neu,pi0,k0,k0b},  {neu,pi0,kp,km},
759    {neu,pim,kp,k0b}},
760
761   {{neu,pip,pip,pim},  {neu,pip,pi0,pi0}, {pro,pip,pim,pi0},
762    {pro,pi0,pi0,pi0},  {lam,kp,pip,pim},  {lam,kp,pi0,pi0},
763    {lam,k0,pip,pi0},   {s0,kp,pip,pim},   {s0,kp,pi0,pi0},
764    {s0,k0,pip,pi0},    {sm,k0,pip,pip},   {sm,kp,pip,pi0},
765    {sp,k0,pip,pim},    {sp,k0,pi0,pi0},   {sp,kp,pim,pi0},
766    {neu,pip,k0,k0b},   {neu,pip,kp,km},   {neu,pi0,kp,k0b},
767    {pro,pim,kp,k0b},   {pro,pi0,kp,km},   {pro,pi0,k0,k0b},
768    {pro,pip,k0,km}}};
769
770const G4int G4RPGPionInelastic::T12_5bfs[2][31][5] =
771  {{{pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}, {neu,pip,pip,pim,pim},
772    {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0}, {lam,k0,pip,pim,pi0},
773    {lam,kp,pim,pi0,pi0},  {lam,kp,pip,pim,pim},  {lam,k0,pi0,pi0,pi0},
774    {s0,kp,pip,pim,pim},   {s0,kp,pim,pi0,pi0},   {s0,k0,pip,pim,pi0},
775    {s0,k0,pi0,pi0,pi0},   {sp,k0,pip,pim,pim},   {sp,k0,pim,pi0,pi0},
776    {sp,kp,pim,pim,pi0},   {sm,k0,pip,pip,pim},   {sm,k0,pip,pi0,pi0},
777    {sm,kp,pip,pim,pi0},   {sm,kp,pi0,pi0,pi0},   {pro,pim,pi0,kp,km},
778    {pro,pim,pi0,k0,k0b},  {pro,pip,pim,k0,km},   {pro,pi0,pi0,k0,km},
779    {pro,pim,pim,kp,k0b},  {neu,pip,pim,kp,km},   {neu,pip,pim,k0,k0b},
780    {neu,pip,pi0,k0,km},   {neu,pim,pi0,kp,k0b},  {neu,pi0,pi0,k0,k0b},
781    {neu,pi0,pi0,kp,km}}, 
782
783   {{neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}, {pro,pip,pip,pim,pim},
784    {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0}, {lam,kp,pip,pim,pi0},
785    {lam,k0,pip,pi0,pi0},  {lam,k0,pip,pip,pim},  {lam,kp,pi0,pi0,pi0},
786    {s0,k0,pip,pip,pim},   {s0,k0,pip,pi0,pi0},   {s0,kp,pip,pim,pi0},
787    {s0,kp,pi0,pi0,pi0},   {sm,kp,pip,pip,pim},   {sm,kp,pip,pi0,pi0},
788    {sm,k0,pip,pip,pi0},   {sp,kp,pip,pim,pim},   {sp,kp,pim,pi0,pi0},
789    {sp,k0,pip,pim,pi0},   {sp,k0,pi0,pi0,pi0},   {neu,pip,pi0,k0,k0b},
790    {neu,pip,pi0,kp,km},   {neu,pip,pim,kp,k0b},  {neu,pi0,pi0,kp,k0b},
791    {neu,pip,pip,k0,km},   {pro,pip,pim,k0,k0b},  {pro,pip,pim,kp,km},
792    {pro,pim,pi0,kp,k0b},  {pro,pip,pi0,k0,km},   {pro,pi0,pi0,kp,km},
793    {pro,pi0,pi0,k0,k0b}}};
794
795const G4int G4RPGPionInelastic::T12_6bfs[2][6][6] =
796{{{pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0},
797  {pro,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pim,pim,pi0},
798  {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0}},
799
800 {{neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0},
801  {neu,pip,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pi0},
802  {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0}}};
803
804const G4int G4RPGPionInelastic::T12_7bfs[2][7][7] =
805{{{pro,pip,pip,pim,pim,pim,pi0}, {pro,pip,pim,pim,pi0,pi0,pi0},
806  {pro,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim},
807  {neu,pip,pip,pim,pim,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0},
808  {neu,pi0,pi0,pi0,pi0,pi0,pi0}},
809
810 {{neu,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pim,pi0,pi0,pi0},
811  {neu,pip,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim},
812  {pro,pip,pip,pim,pim,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0},
813  {pro,pi0,pi0,pi0,pi0,pi0,pi0}}};
814
815const G4int G4RPGPionInelastic::T12_8bfs[2][8][8] =
816{{{pro,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
817  {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
818  {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
819  {neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0}},
820
821 {{neu,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
822  {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
823  {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
824  {pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0}}};
825
826const G4int G4RPGPionInelastic::T12_9bfs[2][9][9] =
827{{{pro,pip,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0},
828  {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
829  {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
830  {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0,pi0},
831  {neu,pip,pip,pip,pip,pim,pim,pim,pim}},
832
833 {{neu,pip,pip,pip,pip,pim,pim,pim,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0},
834  {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
835  {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
836  {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0,pi0},
837  {pro,pip,pip,pip,pip,pim,pim,pim,pim}}};
838
839//
840// Cross sections (in mb) for pi- p -> 2-9 body final states
841//
842// first index:    0-4: channels for mult = 2
843//                5-17: channels for mult = 3
844//               18-39: channels for mult = 4
845//               40-70: channels for mult = 5
846//               71-76: channels for mult = 6
847//               77-83: channels for mult = 7
848//               84-91: channels for mult = 8
849//              92-100: channels for mult = 9
850//
851// second index: kinetic energy
852//
853const G4float G4RPGPionInelastic::pimPCrossSections[101][30] = {
854//
855// multiplicity 2 (5 channels)
856//
857//  pi- p (pi+ n)
858//
859 {  0.0,  1.1,  1.2,  1.4,  1.5,  1.8,  2.0,  3.0,  3.4,  7.0,
860   14.0, 24.0, 14.7, 10.5, 11.8, 20.0, 14.0, 25.0, 12.0,  9.5,
861    8.0,  7.0,  6.0,  5.7,  5.0,  4.6,  4.3,  4.0,  3.8,  3.7 },
862
863//  n pi0  (p pi0)
864 {  0.0,  2.4,  2.8,  3.3,  4.5,  5.7,  6.3,  9.0, 11.0, 17.0,
865   30.0, 43.0, 30.0, 16.5, 11.0,  7.0,  4.3,  5.0,  2.0,  0.9,
866    0.5, 0.24, 0.15,0.094,0.061,0.048,0.035,0.023,0.018,0.014 },
867
868//  L K0  (L K+)
869 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
870    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.65, 0.28, 0.16,
871   0.13,0.075, 0.05,0.032,0.022,0.015,0.011,0.008,0.006,0.004 },
872
873//  S0 K0  (S0 K+)
874 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
875    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.28, 0.18, 0.11,
876  0.091,0.055,0.037,0.025,0.018,0.012,0.008,0.005,0.004,0.003 },
877
878//  S- K+  (S+ K0)
879 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
880    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.2, 0.24, 0.09,
881   0.04,0.012,0.004,0.002,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0 },
882//
883// multiplicity 3 (13 channels)
884//
885//  p pi- pi0  (n pi+ pi0)
886 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
887    0.0,  0.0, 0.05, 0.18, 0.86,  4.4,  5.2,  6.6,  5.4,  4.4,
888    3.5,  2.5,  2.0,  1.4, 0.97, 0.68, 0.55, 0.36,  0.3, 0.22 },
889
890//  n pi+ pi-  (p pi+ pi-)
891 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
892    0.0,  0.0, 0.09,  0.9,  3.3,  6.6,  8.9,  8.2,  8.8,  6.6,
893    5.2,  3.8,  2.9,  1.9,  1.3,  0.9, 0.75, 0.38, 0.31, 0.22 },
894
895//  n pi0 pi0  (p pi0 pi0)
896 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
897    0.0,  0.0, 0.22, 0.74,  1.8,  2.7,  3.0,  2.5,  1.3, 0.64,
898   0.41, 0.24, 0.15, 0.10,0.065,0.042,  0.0,  0.0,  0.0,  0.0 },
899
900//  L K0 pi0  (L K+ pi0)
901 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
902    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.049, 0.18,
903   0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01,0.007,0.004 },
904
905//  L K+ pi-  (L K0 pi+)
906 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
907    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.035, 0.12,
908   0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01,0.007,0.004 },
909
910//  S- K0 pi+  (S+ K+ pi-)
911 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
912    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.10,
913   0.13, 0.06, 0.03, 0.02, 0.01,0.003,0.001,  0.0,  0.0,  0.0 },
914
915//  S- K+ pi0  (S+ K0 pi0)
916 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
917    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.05,
918   0.05, 0.03, 0.02, 0.01,0.005,0.002,0.001,  0.0,  0.0,  0.0 },
919
920//  S+ K0 pi-  (S- K+ pi+)
921 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
922    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
923   0.06, 0.04, 0.02, 0.01, 0.01,0.005,0.003,0.001,  0.0,  0.0 },
924
925//  S0 K+ pi-  (S0 K0 pi+)
926 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
927    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
928   0.07, 0.04, 0.03, 0.02, 0.02, 0.01,0.004,0.002,0.001,  0.0 },
929
930//  S0 K0 pi0  (S0 K+ pi0)
931 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
932    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.02, 0.09,
933   0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01,0.005,0.003,0.002 },
934
935//  p K0 K-  (n K+ K0bar)
936 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
937    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.03,
938   0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
939
940//  n K+ K-  (p K0 K0bar)
941 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
942    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
943   0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01,  0.0,  0.0,  0.0 },
944
945//  n K0 K0bar  (p K+ K-)
946 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
947    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.06,
948   0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01,  0.0,  0.0,  0.0 },
949//
950// multiplicity 4 (22 channels)
951//
952//  p pi+ pi- pi-  (n pi+ pi+ 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.01, 0.05,  0.2, 0.48, 0.88,  1.8,
955    2.0,  2.2,  2.0,  1.9,  1.6,  1.4,  1.3,  1.1, 1.05,  1.0 },
956
957//  p pi- pi0 pi0  (n pi+ pi0 pi0)
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.01, 0.06, 0.18, 0.42,  1.0,  1.9,
960    2.4,  2.4,  2.0,  1.8,  1.6,  1.4,  1.3,  1.2,  1.1,  1.0 },
961
962//  n pi+ pi- pi0  (p pi+ pi- pi0)
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.1,  0.3,  0.5,  1.0,  2.0,  2.0,
965    1.0,  0.8,  0.6,  0.5,  0.4,  0.3, 0.25,  0.2, 0.15,  0.1 },
966
967//  n pi0 pi0 pi0  (p pi0 pi0 pi0)
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.1, 0.55,  1.0, 0.87,  0.7,  0.6,
970    0.5,  0.4, 0.35, 0.29, 0.21, 0.14,  0.1, 0.06, 0.03,  0.0 },
971
972//  L K0 pi+ pi-  (L K+ pi+ pi-)
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.014,
975   0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
976
977//  L K0 pi0 pi0  (L K+ pi0 pi0)
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.014,
980   0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
981
982//  L K+ pi- pi0  (L K0 pi+ pi0)
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.009,
985   0.06,  0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02 },
986
987//  S0 K0 pi+ pi-  (S0 K+ pi+ pi-)
988 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
989    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007,
990   0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
991
992//  S0 K0 pi0 pi0  (S0 K+ pi0 pi0)
993 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
994    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007,
995    0.04,0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
996
997//  S0 K+ pi- pi0  (S0 K0 pi+ pi0)
998 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
999    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1000   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1001
1002//  S+ K+ pi- pi-  (S- K0 pi+ pi+)
1003 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1004    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1005   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1006
1007//  S+ K0 pi- pi0  (S- K+ pi+ pi0)
1008 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1009    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1010   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1011
1012//  S- K+ pi+ pi-  (S+ K0 pi+ pi-)
1013 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1014    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1015   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1016
1017//  S- K+ pi0 pi0  (S+ K0 pi0 pi0)
1018 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1019    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1020   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1021
1022//  S- K0 pi+ pi0  (S+ K+ pi- pi0)
1023 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1024    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1025   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1026
1027//  p pi- K+ K-  (n pi+ K0 K0bar)
1028 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1029    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1030   0.003,0.06,  0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05 },
1031
1032//  p pi- K0 K0bar  (n pi+ K+ K-)
1033 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1034    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1035   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1036
1037//  p pi0 K0 K-  (n pi0 K+ K0bar)
1038 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1039    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1040   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1041
1042//  n pi+ K0 K-  (p pi- K+ K0bar)
1043 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1044    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1045   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1046
1047//  n pi0 K0 K0bar  (p pi0 K+ K-)
1048 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1049    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1050   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1051
1052//  n pi0 K+ K-  (p pi0 K0 K0bar)
1053 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1054    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1055   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1056
1057//  n pi- K+ K0bar  (p pi+ K0 K-)
1058 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1059    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1060   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1061//
1062// multiplicity 5 (31 channels)
1063//
1064//  p pi+ pi- pi- pi0  (n pi+ pi+ pi- 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.007,0.022, 0.10, 0.70,
1067    1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
1068
1069//  p pi- pi0 pi0 pi0  (n pi+ pi0 pi0 pi0)
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.007,0.022, 0.10, 0.70,
1072    1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
1073
1074//  n pi+ pi+ pi- pi-  (p pi+ pi+ pi- 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.014,0.029, 0.10, 0.30,
1077   0.56, 0.93,  1.2,  1.2,  1.2, 0.94, 0.74, 0.53, 0.40, 0.30 },
1078
1079//  n pi+ pi- pi0 pi0  (p pi+ pi- pi0 pi0)
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.006, 0.02, 0.09, 0.63,
1082    1.2,  1.7,  2.0,  1.8,  1.5,  1.3,  1.1, 0.80, 0.70, 0.60 },
1083
1084//  n pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0)
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.006, 0.02, 0.08, 0.60,
1087    1.0,  1.5,  2.0,  1.6,  1.4,  1.1,  1.0, 0.70, 0.60, 0.50 },
1088
1089//  L K0 pi+ pi- pi0  (L K+ pi+ pi- pi0)
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.007, 0.05, 0.07, 0.08, 0.08, 0.07,0.063, 0.06,0.053,0.048 },
1093
1094//  L K+ pi- pi0 pi0  (L K0 pi+ pi0 pi0)
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.003,0.025,0.035, 0.04, 0.04,0.035, 0.03, 0.03,0.025,0.024 },
1098
1099//  L K+ pi+ pi- pi-  (L K0 pi+ pi+ pi-)
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1102  0.002, 0.01,0.036, 0.04, 0.04,0.033,0.031,0.029,0.025,0.022 },
1103
1104//  L K0 pi0 pi0 pi0  (L K+ pi0 pi0 pi0)
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1107  0.001,0.005,0.018, 0.02, 0.02,0.017,0.016,0.015,0.012,0.011 },
1108
1109//  S0 K+ pi+ pi- pi-  (S0 K0 pi+ pi+ pi-)
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1112    0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1113
1114//  S0 K+ pi- pi0 pi0  (S0 K0 pi+ pi0 pi0)
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1117    0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1118
1119//  S0 K0 pi+ pi- pi0  (S0 K+ pi+ pi- pi0)
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1122  0.003, 0.02, 0.03, 0.04, 0.04,0.035, 0.03, 0.03,0.025, 0.02 },
1123
1124//  S0 K0 pi0 pi0 pi0  (S0 K+ pi0 pi0 pi0)
1125 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1126    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1127  0.001, 0.01,0.015, 0.02, 0.02,0.017,0.015,0.015,0.012, 0.01 },
1128
1129//  S+ K0 pi+ pi- pi-  (S- K+ pi+ pi+ pi-)
1130 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1131    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1132  0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1133
1134//  S+ K0 pi- pi0 pi0  (S- K+ pi+ pi0 pi0)
1135 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1136    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1137  0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1138
1139//  S+ K+ pi- pi- pi0  (S- K0 pi+ pi+ pi0)
1140 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1141    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1142  0.001,0.002,0.006,0.009,0.008,0.007,0.006,0.005,0.005,0.004 },
1143
1144//  S- K0 pi+ pi+ pi-  (S+ K+ pi+ pi- pi-)
1145 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1146    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1147  0.002,0.005,0.015,0.014,0.012,0.009,0.008,0.006,0.005,0.004 },
1148
1149//  S- K0 pi+ pi0 pi0  (S+ K+ pi- pi0 pi0)
1150 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1151    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1152  0.001,0.002,0.008,0.007,0.006,0.005,0.004,0.003,0.003,0.002 },
1153
1154//  S- K+ pi+ pi- pi0  (S+ K0 pi+ pi- pi0)
1155 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1156    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1157  0.002,0.007,0.023,0.025,0.021,0.018,0.015,0.013,0.010,0.009 },
1158
1159//  S- K+ pi0 pi0 pi0  (S+ K0 pi0 pi0 pi0)
1160 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1161    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1162  0.001,0.004,0.011,0.013,0.011,0.009,0.008,0.007,0.006,0.005 },
1163
1164//  p pi- pi0 K+ K-  (n pi+ pi0 K0 K0bar)
1165 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1166    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1167    0.0, 0.01, 0.07, 0.07, 0.06,0.055, 0.05,0.047,0.042, 0.04 },
1168
1169//  p pi- pi0 K0 K0bar  (n pi+ pi0 K+ K-)
1170 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1171    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1172    0.0,0.003, 0.02, 0.02,0.018,0.017,0.015,0.012,0.011, 0.01 },
1173
1174//  p pi+ pi- K0 K-  (n pi+ pi- K+ K0bar)
1175 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1176    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1177    0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1178
1179//  p pi0 pi0 K0 K-  (n pi0 pi0 K+ K0bar)
1180 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1181    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1182    0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1183
1184//  p pi- pi- K+ K0bar  (n pi+ pi+ K0 K-)
1185 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1186    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1187    0.0,0.007, 0.04, 0.05, 0.04,0.035, 0.03, 0.03,0.027,0.025 },
1188
1189//  n pi+ pi- K+ K-  (p pi+ pi- K0 K0bar)
1190 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1191    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1192    0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1193
1194//  n pi+ pi- K0 K0bar  (p pi+ pi- K+ K-)
1195 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1196    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1197    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1198
1199//  n pi+ pi0 K0 K-  (p pi- pi0 K+ K0bar)
1200 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1201    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1202    0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1203
1204//  n pi- pi0 K+ K0bar  (p pi+ pi0 K0 K-)
1205 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1206    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1207    0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1208
1209//  n pi0 pi0 K0 K0bar  (p pi0 pi0 K+ K-)
1210 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1211    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1212    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1213
1214//  n pi0 pi0 K+ K-  (p pi0 pi0 K0 K0bar)
1215 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1216    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1217    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1218//
1219// multiplicity 6 (6 channels)
1220//
1221//  p pi+ pi+ pi- pi- pi-  (n pi+ pi+ pi+ pi- pi-)
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.01, 0.03, 0.04,
1224   0.06, 0.11, 0.16, 0.22, 0.31, 0.34,  0.3, 0.24, 0.19, 0.16 },
1225
1226//  p pi+ pi- pi- pi0 pi0  (n pi+ pi+ pi- pi0 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.01, 0.04,  0.1,
1229   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1230
1231//  p pi- pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 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.01, 0.04,  0.1,
1234   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1235
1236//  n pi+ pi+ pi- pi- pi0  (p pi+ pi+ pi- pi- pi0)
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.01, 0.04,  0.1,
1239   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38, 0.4,  0.33, 0.32 },
1240
1241//  n pi+ pi- pi0 pi0 pi0  (p pi+ pi- pi0 pi0 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.01, 0.04,  0.1,
1244   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1245
1246//  n pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 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.01, 0.04,  0.1,
1249   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1250//
1251// multiplicity 7 (7 channels)
1252//
1253//  p pi+ pi+ pi- pi- pi- pi0
1254 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1255    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.014, 0.03,
1256   0.05, 0.11, 0.19, 0.37, 0.67, 0.75,  0.7, 0.59, 0.52, 0.47 },
1257
1258//  p pi+ pi- pi- pi0 pi0 pi0
1259 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1260    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1261   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1262
1263//  p pi- pi0 pi0 pi0 pi0 pi0
1264 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1265    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1266   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1267
1268//  n pi+ pi+ pi+ pi- pi- pi-
1269 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1270    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,0.006,
1271  0.012,0.029,0.062, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13 },
1272
1273//  n pi+ pi+ pi- pi- pi0 pi0
1274 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1275    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1276   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1277
1278//  n pi+ pi- pi0 pi0 pi0 pi0
1279 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1280    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1281   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1282
1283//  n pi0 pi0 pi0 pi0 pi0 pi0
1284 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1285    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1286   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1287//
1288// multiplicity 8 (8 channels)
1289//
1290//  p pi+ pi+ pi+ pi- pi- pi- pi-  (n pi+ pi+ pi+ pi+ pi- pi- 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1294
1295//  p pi+ pi+ pi- pi- pi- pi0 pi0  (n pi+ pi+ pi+ pi- pi- pi0 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1299
1300//  p pi+ pi- pi- pi0 pi0 pi0 pi0  (n pi+ pi+ pi- pi0 pi0 pi0 pi0)
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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1304
1305//  p pi- pi0 pi0 pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1309
1310//  n pi0 pi0 pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 pi0 pi0 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1314
1315//  n pi+ pi- pi0 pi0 pi0 pi0 pi0  (p pi+ pi- pi0 pi0 pi0 pi0 pi0)
1316 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1317    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1318  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1319
1320//  n pi+ pi+ pi- pi- pi0 pi0 pi0  (p pi+ pi+ pi- pi- pi0 pi0 pi0)
1321 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1322    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1323  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1324
1325//  n pi+ pi+ pi+ pi- pi- pi- pi0  (p pi+ pi+ pi+ pi- pi- pi- pi0)
1326 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1327    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1328  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1329//
1330// multiplicity 9 (9 channels)
1331//
1332//  p pi+ pi+ pi+ pi- 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.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1336
1337//  p pi+ pi+ pi- pi- 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.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1341
1342//  p pi+ pi- pi- pi0 pi0 pi0 pi0 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.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1346
1347//  p pi- pi0 pi0 pi0 pi0 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.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1351
1352//  n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
1353 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1354    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1355    0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1356
1357//  n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0
1358 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1359    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1360    0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1361
1362//  n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0
1363 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1364    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1365    0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1366
1367//  n pi+ pi+ pi+ pi- pi- pi- pi0 pi0
1368 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1369    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1370    0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1371
1372//  n pi+ pi+ pi+ pi+ pi- pi- pi- pi-
1373 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1374    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1375    0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 } };
1376
1377 /* end of file */
1378 
Note: See TracBrowser for help on using the repository browser.