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

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

update ti head
File size: 58.6 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.2 2009/05/25 19:05:47 dennis Exp $
27// GEANT4 tag $Name: geant4-09-03-ref-09 $
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  //  printCrossSections();
74
75}
76
77
78/*
79void G4RPGPionInelastic::printCrossSections() const
80{
81  G4cout << " pi+ p total cross sections (mb) " << G4endl;
82  for (G4int i = 0; i < 5; i++) {
83    G4int istart = i*6;
84    G4int istop = istart + 6;
85    for (G4int t = istart; t < istop; t++) G4cout << pipPtot[t] << "  " ;
86    G4cout << G4endl;
87  }
88  G4cout << G4endl;
89
90  G4cout << " pi- p total cross sections (mb) " << G4endl;
91  for (G4int i = 0; i < 5; i++) {
92    G4int istart = i*6;
93    G4int istop = istart + 6;
94    for (G4int t = istart; t < istop; t++) G4cout << pimPtot[t] << "  " ;
95    G4cout << G4endl;
96  }
97}
98*/
99
100
101G4int G4RPGPionInelastic::GetMultiplicityT12(G4double KE) const
102{
103  G4double multint(0.);
104  std::vector<G4double> sigma;
105
106  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
107  G4int k = epair.first;
108  G4double fraction = epair.second;
109
110  for(G4int m = 0; m < 8; m++) {
111    multint = t12_dSigma_dMult[m][k]
112         + fraction*(t12_dSigma_dMult[m][k+1] - t12_dSigma_dMult[m][k]);
113      sigma.push_back(multint);
114  }
115
116  return sampleFlat(sigma) + 2;
117}
118
119
120G4int G4RPGPionInelastic::GetMultiplicityT32(G4double KE) const
121{
122  G4double multint(0.);
123  std::vector<G4double> sigma;
124
125  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
126  G4int k = epair.first;
127  G4double fraction = epair.second;
128
129  for(G4int m = 0; m < 8; m++) {
130    multint = t32_dSigma_dMult[m][k]
131         + fraction*(t32_dSigma_dMult[m][k+1] - t32_dSigma_dMult[m][k]);
132      sigma.push_back(multint);
133  }
134
135  return sampleFlat(sigma) + 2;
136}
137
138
139std::vector<G4int> 
140G4RPGPionInelastic::GetFSPartTypesForT32(G4int mult, G4double KE, G4int tzindex) const
141{
142  G4int i;
143  G4double sigint(0.);
144  std::vector<G4double> sigma;
145
146  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
147  G4int k = epair.first;
148  G4double fraction = epair.second;
149
150  G4int start = pipPindex[mult-2][0];
151  G4int stop = pipPindex[mult-2][1];
152
153  for(i = start; i < stop; i++) {
154      sigint = pipPCrossSections[i][k]
155          + fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[i][k]);
156      sigma.push_back(sigint);
157  }
158
159  G4int channel = sampleFlat(sigma);
160
161  std::vector<G4int> kinds;
162
163  if (mult == 2) {
164    for(i = 0; i < mult; i++) kinds.push_back(T32_2bfs[tzindex][channel][i]);
165  } else if (mult == 3) {
166    for(i = 0; i < mult; i++) kinds.push_back(T32_3bfs[tzindex][channel][i]);
167  } else if (mult == 4) {
168    for(i = 0; i < mult; i++) kinds.push_back(T32_4bfs[tzindex][channel][i]);
169  } else if (mult == 5) {
170    for(i = 0; i < mult; i++) kinds.push_back(T32_5bfs[tzindex][channel][i]);
171  } else if (mult == 6) {
172    for(i = 0; i < mult; i++) kinds.push_back(T32_6bfs[tzindex][channel][i]);
173  } else if (mult == 7) {
174    for(i = 0; i < mult; i++) kinds.push_back(T32_7bfs[tzindex][channel][i]);
175  } else if (mult == 8) {
176    for(i = 0; i < mult; i++) kinds.push_back(T32_8bfs[tzindex][channel][i]);
177  } else if (mult == 9) {
178    for(i = 0; i < mult; i++) kinds.push_back(T32_9bfs[tzindex][channel][i]);
179  } else {
180    G4cout << " Illegal multiplicity " << G4endl;
181  }
182
183  return kinds;
184}
185
186
187std::vector<G4int> 
188G4RPGPionInelastic::GetFSPartTypesForT12(G4int mult, G4double KE, G4int tzindex) const
189{
190  G4int i;
191  G4double sigint(0.);
192  std::vector<G4double> sigma;
193
194  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
195  G4int k = epair.first;
196  G4double fraction = epair.second;
197
198  G4int start = pimPindex[mult-2][0];
199  G4int stop = pimPindex[mult-2][1];
200
201  for(i = start; i < stop; i++) {
202      sigint = pimPCrossSections[i][k]
203          + fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
204      sigma.push_back(sigint);
205  }
206
207  G4int channel = sampleFlat(sigma);
208
209  std::vector<G4int> kinds;
210
211  if (mult == 2) {
212    for(i = 0; i < mult; i++) kinds.push_back(T12_2bfs[tzindex][channel][i]);
213  } else if (mult == 3) {
214    for(i = 0; i < mult; i++) kinds.push_back(T12_3bfs[tzindex][channel][i]);
215  } else if (mult == 4) {
216    for(i = 0; i < mult; i++) kinds.push_back(T12_4bfs[tzindex][channel][i]);
217  } else if (mult == 5) {
218    for(i = 0; i < mult; i++) kinds.push_back(T12_5bfs[tzindex][channel][i]);
219  } else if (mult == 6) {
220    for(i = 0; i < mult; i++) kinds.push_back(T12_6bfs[tzindex][channel][i]);
221  } else if (mult == 7) {
222    for(i = 0; i < mult; i++) kinds.push_back(T12_7bfs[tzindex][channel][i]);
223  } else if (mult == 8) {
224    for(i = 0; i < mult; i++) kinds.push_back(T12_8bfs[tzindex][channel][i]);
225  } else if (mult == 9) {
226    for(i = 0; i < mult; i++) kinds.push_back(T12_9bfs[tzindex][channel][i]);
227  } else {
228    G4cout << " Illegal multiplicity " << G4endl;
229  }
230
231  return kinds;
232}
233
234
235///////////////////////////////////////////////////////////////////////////////
236//                                                                           //
237//   pi+ p and pi- n (|Tz| = 3/2) cross sections                             //
238//   and final state particle types                                          //
239//                                                                           //
240///////////////////////////////////////////////////////////////////////////////
241
242// Total pi+ cross section as a function of kinetic energy
243G4double G4RPGPionInelastic::pipPtot[30];
244
245// pi+ multiplicities as a function of kinetic energy
246G4double G4RPGPionInelastic::t32_dSigma_dMult[8][30];
247
248const G4int G4RPGPionInelastic::pipPindex[8][2] =
249 {{0, 1}, {2, 8}, {9,23}, {24,47}, {48,52}, {53,58}, {59,65}, {66,73}}; 
250
251// Outgoing particle types of a given multiplicity
252// T32_nbfs = final state types for pi+ p and pi- n
253
254const G4int G4RPGPionInelastic::T32_2bfs[2][2][2] =
255  {{{pro,pip}, {sp,kp}},
256
257   {{neu,pim}, {sm,k0}}};
258
259const G4int G4RPGPionInelastic::T32_3bfs[2][7][3] =
260  {{{pro,pip,pi0}, {neu,pip,pip}, {sp,kp,pi0}, {sp,k0,pip}, 
261    {s0,kp,pip},   {lam,kp,pip},  {pro,kp,k0b}},
262
263   {{neu,pim,pi0}, {pro,pim,pim}, {sm,k0,pi0}, {sm,kp,pim},
264    {s0,k0,pim},   {lam,k0,pim},  {neu,k0,km}}};
265
266const G4int G4RPGPionInelastic::T32_4bfs[2][15][4] =
267  {{{pro,pip,pip,pim},{pro,pip,pi0,pi0},{neu,pip,pip,pi0},
268    {sp,kp,pip,pim},  {sp,kp,pi0,pi0},  {sp,k0,pip,pi0},
269    {s0,k0,pip,pip},  {s0,kp,pip,pi0},  {lam,kp,pip,pi0},
270    {lam,k0,pip,pip}, {sm,kp,pip,pip},  {pro,pip,kp,km},
271    {pro,pip,k0,k0b}, {pro,pi0,kp,k0b}, {neu,pip,kp,k0b}},
272
273   {{neu,pip,pim,pim},{neu,pim,pi0,pi0},{pro,pim,pim,pi0},
274    {sm,k0,pip,pim},  {sm,k0,pi0,pi0},  {sm,kp,pim,pi0},
275    {s0,kp,pim,pim},  {s0,k0,pim,pi0},  {lam,k0,pim,pi0},
276    {lam,kp,pim,pim}, {sp,k0,pim,pim},  {neu,pim,k0,k0b},
277    {neu,pim,kp,km},  {neu,pi0,k0,km},  {pro,pim,k0,km}}};
278
279const G4int G4RPGPionInelastic::T32_5bfs[2][24][5] =
280  {{{pro,pip,pip,pim,pi0}, {pro,pip,pi0,pi0,pi0}, {neu,pip,pip,pip,pim},
281    {neu,pip,pip,pi0,pi0}, {sp,kp,pip,pim,pi0},   {sp,kp,pi0,pi0,pi0},
282    {sp,k0,pip,pip,pim},   {sp,k0,pip,pi0,pi0},   {lam,k0,pip,pip,pi0},
283    {lam,kp,pip,pip,pim},  {lam,kp,pip,pi0,pi0},  {s0,kp,pip,pip,pim},
284    {s0,kp,pip,pi0,pi0},   {s0,k0,pip,pip,pi0},   {sm,kp,pip,pip,pi0},
285    {sm,k0,pip,pip,pip},   {pro,pip,pim,kp,k0b},  {pro,pip,pip,k0,km},
286    {pro,pip,pi0,kp,km},   {pro,pip,pi0,k0,k0b},  {pro,pi0,pi0,kp,k0b},
287    {neu,pip,pip,kp,km},   {neu,pip,pip,k0,k0b},  {neu,pip,pi0,kp,k0b}},
288
289   {{neu,pip,pim,pim,pi0}, {neu,pim,pi0,pi0,pi0}, {pro,pip,pim,pim,pim},
290    {pro,pim,pim,pi0,pi0}, {sm,k0,pip,pim,pi0},   {sm,k0,pi0,pi0,pi0},
291    {sm,kp,pip,pim,pim},   {sm,kp,pim,pi0,pi0},   {lam,kp,pim,pim,pi0},
292    {lam,k0,pip,pim,pim},  {lam,k0,pim,pi0,pi0},  {s0,k0,pip,pim,pim},
293    {s0,k0,pim,pi0,pi0},   {s0,kp,pim,pim,pi0},   {sp,k0,pim,pim,pi0},
294    {sp,kp,pim,pim,pim},   {neu,pip,pim,k0,km},   {neu,pim,pim,kp,k0b},
295    {neu,pim,pi0,k0,k0b},  {neu,pim,pi0,kp,km},   {neu,pi0,pi0,k0,km},
296    {pro,pim,pim,k0,k0b},  {pro,pim,pim,kp,km},   {pro,pim,pi0,k0,km}}};
297
298const G4int G4RPGPionInelastic::T32_6bfs[2][5][6] =
299{{{pro,pip,pip,pip,pim,pim}, {pro,pip,pip,pim,pi0,pi0},
300  {pro,pip,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0},
301  {neu,pip,pip,pip,pim,pi0}},
302
303 {{neu,pip,pip,pim,pim,pim}, {neu,pip,pim,pim,pi0,pi0},
304  {neu,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0},
305  {pro,pip,pim,pim,pim,pi0}}};
306
307const G4int G4RPGPionInelastic::T32_7bfs[2][6][7] =
308{{{pro,pip,pip,pip,pim,pim,pi0}, {pro,pip,pip,pim,pi0,pi0,pi0},
309  {pro,pip,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim},
310  {neu,pip,pip,pip,pim,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0}},
311
312 {{neu,pip,pip,pim,pim,pim,pi0}, {neu,pip,pim,pim,pi0,pi0,pi0},
313  {neu,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pim},
314  {pro,pip,pim,pim,pim,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0}}};
315
316const G4int G4RPGPionInelastic::T32_8bfs[2][7][8] =
317{{{pro,pip,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pip,pim,pim,pi0,pi0},
318  {pro,pip,pip,pim,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0},
319  {neu,pip,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pip,pim,pi0,pi0,pi0},
320  {neu,pip,pip,pi0,pi0,pi0,pi0,pi0}},
321
322 {{neu,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pip,pim,pim,pim,pi0,pi0},
323  {neu,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
324  {pro,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pim,pim,pim,pi0,pi0,pi0},
325  {pro,pim,pim,pi0,pi0,pi0,pi0,pi0}}};
326
327const G4int G4RPGPionInelastic::T32_9bfs[2][8][9] =
328{{{pro,pip,pip,pip,pip,pim,pim,pim,pi0}, {pro,pip,pip,pip,pim,pim,pi0,pi0,pi0},
329  {pro,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
330  {neu,pip,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pip,pim,pim,pi0,pi0},
331  {neu,pip,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
332
333 {{neu,pip,pip,pip,pim,pim,pim,pim,pi0}, {neu,pip,pip,pim,pim,pim,pi0,pi0,pi0},
334  {neu,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
335  {pro,pip,pip,pip,pim,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pim,pi0,pi0},
336  {pro,pip,pim,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
337
338//
339// Cross sections (in mb) for pi+ p -> 2-9 body final states
340//
341// first index:    0-1: channels for mult = 2
342//                 2-8: channels for mult = 3
343//                9-23: channels for mult = 4
344//               24-47: channels for mult = 5
345//               48-52: channels for mult = 6
346//               53-58: channels for mult = 7
347//               59-65: channels for mult = 8
348//               66-73: channels for mult = 9
349//
350// second index: kinetic energy
351//
352
353const G4float G4RPGPionInelastic::pipPCrossSections[74][30] = {
354//
355// multiplicity 2 (2 channels)
356//
357//  p pi+ (n pi-)
358 { 0.00, 1.20, 2.50, 3.80, 5.00, 7.00, 9.00, 15.0, 30.0, 64.0,
359  130.0,190.0,130.0, 55.7, 27.2, 14.0, 8.50, 13.0, 18.0, 11.0,
360   8.50, 7.00, 6.20, 5.60, 5.00, 4.50, 4.20, 4.00, 3.80, 3.60 },
361
362//  S+ K+ (S- K0)
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.0,  0.0,  0.0, 0.16, 0.60, 0.32,
365   0.19,  0.1, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
366//
367// multiplicity 3 (7 channels)
368//
369//  p pi+ pi0 (n pi- pi0)
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.2,  0.6,  2.4,  8.8, 10.0, 12.0,  6.2,
372   4.00, 2.40, 1.69, 1.10, 0.73, 0.49, 0.41, 0.31, 0.24, 0.15 },
373
374//  n pi+ pi+ (p pi- pi-)
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.10, 0.20, 0.60, 1.50, 2.30, 3.60, 3.00,
377   2.30, 1.70, 1.30, 0.95, 0.69, 0.46, 0.38, 0.27, 0.20, 0.15 },
378
379//  S+ K+ pi0 (S- K0 pi0)
380 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
381    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
382   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
383
384//  S+ K0 pi+
385 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
386    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
387   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
388
389//  S0 K+ 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.0, 0.005,0.12,
392   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
393
394//  L K+ pi+
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.0, 0.005,0.12,
397   0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
398
399//  p K+ K0bar
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.0,  0.0, 0.04,
402   0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
403//
404// multiplicity 4 (15 channels)
405//
406//  p pi+ pi+ pi-
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.06, 0.20, 0.78, 2.20, 3.20,
409   3.50, 3.10, 2.70, 2.30, 2.00, 1.50, 1.40, 1.20, 1.00, 0.90 },
410
411//  p pi+ pi0 pi0
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.04, 0.13, 0.52, 1.50, 2.20,
414   2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
415
416//  n pi+ 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.04, 0.13, 0.52, 1.50, 2.20,
419   2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
420
421//  S+ K+ pi+ pi-
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.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
425
426//  S+ K+ pi0 pi0
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.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
430
431//  S+ K0 pi+ pi0
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.015,0.06, 0.06, 0.05, 0.04,0.032,0.028, 0.02,0.017,0.014 },
435
436//  S0 K0 pi+ pi+
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.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
440
441//  S0 K+ pi+ pi0
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.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
445
446//  L K+ pi+ pi0
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.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
450
451//  L K0 pi+ pi+
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.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
455
456//  S- K+ pi+ pi+
457 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
458    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,
459   0.02, 0.16, 0.18, 0.13, 0.09, 0.06,0.045, 0.03,0.025, 0.02 },
460
461//  p pi+ K+ K-
462 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
463    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
464   0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05,0.045, 0.04 },
465
466//  p pi+ K0 K0bar
467 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
468    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
469  0.005, 0.05, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05,0.045, 0.04 },
470
471//  p pi0 K+ K0bar
472 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
473    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
474   0.01, 0.05,0.075, 0.08,0.075, 0.07, 0.06, 0.05,0.045, 0.04 },
475
476//  n pi+ K+ K0bar
477 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
478    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
479   0.01, 0.05, 0.07,0.065,0.055,0.048, 0.04, 0.03,0.027,0.023 },
480//
481// multiplicity 5 (24 channels)
482//
483//  p pi+ pi+ pi- 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.001, 0.05, 0.30, 2.00,
486   3.20, 3.70, 3.00, 2.50, 2.10, 1.60, 1.40, 1.10, 0.89, 0.70 },
487
488//  p pi+ pi0 pi0 pi0
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.001, 0.03, 0.21, 1.40,
491   2.20, 2.60, 2.10, 1.80, 1.50, 1.10, 1.00, 0.80, 0.62, 0.50 },
492
493//  n pi+ pi+ pi+ pi-
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.007, 0.02, 0.05, 0.19,
496   0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
497
498//  n pi+ pi+ pi0 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.007, 0.02, 0.05, 0.19,
501   0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
502
503//  S+ K+ pi+ pi- pi0
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.003,0.013, 0.05, 0.09, 0.08, 0.06, 0.05,  0.04,0.036,0.03 },
507
508//  S+ K+ pi0 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.006, 0.02, 0.04, 0.04, 0.03,0.025, 0.02, 0.02, 0.01 },
512
513//  S+ K0 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.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
517
518//  S+ K0 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.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
522
523//  L 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.002,0.013, 0.04,0.052,0.059,0.053, 0.05,0.043,0.037, 0.03 },
527
528//  L K+ pi+ pi+ pi-
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.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
532
533//  L K+ pi+ pi0 pi0
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.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
537
538//  S0 K+ pi+ pi+ pi-
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.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
542
543//  S0 K+ pi+ pi0 pi0
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.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
547
548//  S0 K0 pi+ pi+ pi0
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.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
552
553//  S- K+ pi+ pi+ pi0
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.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
557
558//  S- K0 pi+ pi+ pi+
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.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
562
563//  p pi+ pi- K+ K0bar
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.001,0.005, 0.02,0.065, 0.08, 0.07, 0.06,0.054,0.048, 0.04 },
567
568//  p pi+ pi+ K0 K-
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.001,0.005, 0.02,0.045, 0.05,0.047, 0.04,0.033, 0.03,0.026 },
572
573//  p pi+ pi0 K+ K-
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.001, 0.02, 0.06, 0.06,0.054,0.048,0.042,0.038,0.035, 0.03 },
577
578//  p pi+ pi0 K0 K0bar
579 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
580    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
581  0.001, 0.01, 0.04, 0.05, 0.05, 0.04, 0.04,0.035,0.032, 0.03 },
582
583//  p pi0 pi0 K+ K0bar
584 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
585    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
586  0.001,0.005, 0.01, 0.03, 0.04,0.035, 0.03,0.025,0.025, 0.02 },
587
588//  n pi+ pi+ K+ K-
589 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
590    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
591    0.0, 0.02, 0.06, 0.05,0.042,0.038,0.035, 0.03,0.027,0.022 },
592
593//  n pi+ pi+ K0 K0bar
594 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
595    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
596    0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
597
598//  n pi+ pi0 K+ K0bar
599 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
600    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
601    0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
602//
603// multiplicity 6 (5 channels)
604//
605//  p pi+ pi+ pi+ pi- pi-
606 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
607    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
608   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
609
610//  p pi+ pi+ pi- pi0 pi0
611 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
612    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
613   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
614
615//  p pi+ pi0 pi0 pi0 pi0
616 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
617    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
618   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
619
620//  n pi+ pi+ pi0 pi0 pi0
621 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
622    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
623   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
624
625//  n pi+ pi+ pi+ pi- pi0
626 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
627    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
628   0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
629//
630// multiplicity 7 (6 channels)
631//
632//  p pi+ pi+ pi+ pi- pi- 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.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
636
637//  p pi+ pi+ pi- pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
640  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
641
642//  p pi+ pi0 pi0 pi0 pi0 pi0
643 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
644    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
645  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
646
647//  n pi+ pi+ pi+ pi+ pi- pi-
648 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
649    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
650   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
651
652//  n pi+ pi+ pi+ pi- pi0 pi0
653 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
654    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
655   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
656
657//  n pi+ pi+ pi0 pi0 pi0 pi0
658 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
659    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
660   0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
661//
662// multiplicity 8 (7 channels)
663//
664//  p pi+ pi+ pi+ pi+ pi- pi- pi-
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.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
668
669//  p pi+ pi+ pi+ pi- pi- 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.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
673
674//  p pi+ pi+ pi- pi0 pi0 pi0 pi0
675 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
676    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
677  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
678
679//  p pi+ pi0 pi0 pi0 pi0 pi0 pi0
680 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
681    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
682  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
683
684//  n pi+ pi+ pi+ pi+ pi- pi- pi0
685 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
686    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
687  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
688
689//  n pi+ pi+ pi+ pi- pi0 pi0 pi0
690 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
691    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
692  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
693
694//  n pi+ pi+ pi0 pi0 pi0 pi0 pi0
695 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
696    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
697  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
698//
699// multiplicity 9 (8 channels)
700//
701//  p pi+ pi+ pi+ pi+ pi- pi- pi- 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.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
705
706//  p pi+ pi+ pi+ pi- pi- 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.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
710
711//  p pi+ pi+ pi- 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.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
715
716//  p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
719  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
720
721//  n pi+ pi+ pi+ pi+ pi+ pi- pi- pi-
722 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
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.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
725
726//  n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
727 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
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.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
730
731//  n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
732 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
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.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
735
736//  n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
737 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
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.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 } };
740
741///////////////////////////////////////////////////////////////////////////////
742//                                                                           //
743//   pi- p and pi+ n (|Tz| = 1/2) cross sections                             //
744//   and final state particle types                                          //
745//                                                                           //
746///////////////////////////////////////////////////////////////////////////////
747
748// Total pi- cross section as a function of kinetic energy
749G4double G4RPGPionInelastic::pimPtot[30];
750
751// pi- multiplicities as a function of kinetic energy
752G4double G4RPGPionInelastic::t12_dSigma_dMult[8][30];
753
754const G4int G4RPGPionInelastic::pimPindex[8][2] =
755 {{0, 4}, {5,17}, {18,39}, {40,70}, {71,76}, {77,83}, {84,91}, {92,100}}; 
756
757// Outgoing particle types of a given multiplicity
758// T12_nbfs = final state types for pi- p and pi+ n
759
760const G4int G4RPGPionInelastic::T12_2bfs[2][5][2] =
761  {{{pro,pim}, {neu,pi0}, {lam,k0}, {s0,k0}, {sm,kp}},
762
763   {{neu,pip}, {pro,pi0}, {lam,kp}, {s0,kp}, {sp,k0}}};
764
765const G4int G4RPGPionInelastic::T12_3bfs[2][13][3] =
766  {{{pro,pim,pi0}, {neu,pip,pim}, {neu,pi0,pi0}, {lam,k0,pi0}, 
767    {lam,kp,pim},  {sm,k0,pip},   {sm,kp,pi0},   {sp,k0,pim},
768    {s0,kp,pim},   {s0,k0,pi0},   {pro,k0,km},   {neu,kp,km},
769    {neu,k0,k0b}},
770
771   {{neu,pip,pi0}, {pro,pip,pim}, {pro,pi0,pi0}, {lam,kp,pi0},
772    {lam,k0,pip},  {sp,kp,pim},   {sp,k0,pi0},   {sm,kp,pip},
773    {s0,k0,pip},   {s0,kp,pi0},   {neu,kp,k0b},  {pro,k0,k0b},
774    {pro,kp,km}}};
775
776const G4int G4RPGPionInelastic::T12_4bfs[2][22][4] =
777  {{{pro,pip,pim,pim}, {pro,pim,pi0,pi0}, {neu,pip,pim,pi0},
778    {neu,pi0,pi0,pi0}, {lam,k0,pip,pim},  {lam,k0,pi0,pi0},
779    {lam,kp,pim,pi0},  {s0,k0,pip,pim},   {s0,k0,pi0,pi0},
780    {s0,kp,pim,pi0},   {sp,kp,pim,pim},   {sp,k0,pim,pi0},
781    {sm,kp,pip,pim},   {sm,kp,pi0,pi0},   {sm,k0,pip,pi0},
782    {pro,pim,kp,km},   {pro,pim,k0,k0b},  {pro,pi0,k0,km},
783    {neu,pip,k0,km},   {neu,pi0,k0,k0b},  {neu,pi0,kp,km},
784    {neu,pim,kp,k0b}},
785
786   {{neu,pip,pip,pim},  {neu,pip,pi0,pi0}, {pro,pip,pim,pi0},
787    {pro,pi0,pi0,pi0},  {lam,kp,pip,pim},  {lam,kp,pi0,pi0},
788    {lam,k0,pip,pi0},   {s0,kp,pip,pim},   {s0,kp,pi0,pi0},
789    {s0,k0,pip,pi0},    {sm,k0,pip,pip},   {sm,kp,pip,pi0},
790    {sp,k0,pip,pim},    {sp,k0,pi0,pi0},   {sp,kp,pim,pi0},
791    {neu,pip,k0,k0b},   {neu,pip,kp,km},   {neu,pi0,kp,k0b},
792    {pro,pim,kp,k0b},   {pro,pi0,kp,km},   {pro,pi0,k0,k0b},
793    {pro,pip,k0,km}}};
794
795const G4int G4RPGPionInelastic::T12_5bfs[2][31][5] =
796  {{{pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}, {neu,pip,pip,pim,pim},
797    {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0}, {lam,k0,pip,pim,pi0},
798    {lam,kp,pim,pi0,pi0},  {lam,kp,pip,pim,pim},  {lam,k0,pi0,pi0,pi0},
799    {s0,kp,pip,pim,pim},   {s0,kp,pim,pi0,pi0},   {s0,k0,pip,pim,pi0},
800    {s0,k0,pi0,pi0,pi0},   {sp,k0,pip,pim,pim},   {sp,k0,pim,pi0,pi0},
801    {sp,kp,pim,pim,pi0},   {sm,k0,pip,pip,pim},   {sm,k0,pip,pi0,pi0},
802    {sm,kp,pip,pim,pi0},   {sm,kp,pi0,pi0,pi0},   {pro,pim,pi0,kp,km},
803    {pro,pim,pi0,k0,k0b},  {pro,pip,pim,k0,km},   {pro,pi0,pi0,k0,km},
804    {pro,pim,pim,kp,k0b},  {neu,pip,pim,kp,km},   {neu,pip,pim,k0,k0b},
805    {neu,pip,pi0,k0,km},   {neu,pim,pi0,kp,k0b},  {neu,pi0,pi0,k0,k0b},
806    {neu,pi0,pi0,kp,km}}, 
807
808   {{neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}, {pro,pip,pip,pim,pim},
809    {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0}, {lam,kp,pip,pim,pi0},
810    {lam,k0,pip,pi0,pi0},  {lam,k0,pip,pip,pim},  {lam,kp,pi0,pi0,pi0},
811    {s0,k0,pip,pip,pim},   {s0,k0,pip,pi0,pi0},   {s0,kp,pip,pim,pi0},
812    {s0,kp,pi0,pi0,pi0},   {sm,kp,pip,pip,pim},   {sm,kp,pip,pi0,pi0},
813    {sm,k0,pip,pip,pi0},   {sp,kp,pip,pim,pim},   {sp,kp,pim,pi0,pi0},
814    {sp,k0,pip,pim,pi0},   {sp,k0,pi0,pi0,pi0},   {neu,pip,pi0,k0,k0b},
815    {neu,pip,pi0,kp,km},   {neu,pip,pim,kp,k0b},  {neu,pi0,pi0,kp,k0b},
816    {neu,pip,pip,k0,km},   {pro,pip,pim,k0,k0b},  {pro,pip,pim,kp,km},
817    {pro,pim,pi0,kp,k0b},  {pro,pip,pi0,k0,km},   {pro,pi0,pi0,kp,km},
818    {pro,pi0,pi0,k0,k0b}}};
819
820const G4int G4RPGPionInelastic::T12_6bfs[2][6][6] =
821{{{pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0},
822  {pro,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pim,pim,pi0},
823  {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0}},
824
825 {{neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0},
826  {neu,pip,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pi0},
827  {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0}}};
828
829const G4int G4RPGPionInelastic::T12_7bfs[2][7][7] =
830{{{pro,pip,pip,pim,pim,pim,pi0}, {pro,pip,pim,pim,pi0,pi0,pi0},
831  {pro,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim},
832  {neu,pip,pip,pim,pim,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0},
833  {neu,pi0,pi0,pi0,pi0,pi0,pi0}},
834
835 {{neu,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pim,pi0,pi0,pi0},
836  {neu,pip,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim},
837  {pro,pip,pip,pim,pim,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0},
838  {pro,pi0,pi0,pi0,pi0,pi0,pi0}}};
839
840const G4int G4RPGPionInelastic::T12_8bfs[2][8][8] =
841{{{pro,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
842  {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
843  {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
844  {neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0}},
845
846 {{neu,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
847  {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
848  {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
849  {pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0}}};
850
851const G4int G4RPGPionInelastic::T12_9bfs[2][9][9] =
852{{{pro,pip,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0},
853  {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
854  {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
855  {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0,pi0},
856  {neu,pip,pip,pip,pip,pim,pim,pim,pim}},
857
858 {{neu,pip,pip,pip,pip,pim,pim,pim,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0},
859  {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
860  {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
861  {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0,pi0},
862  {pro,pip,pip,pip,pip,pim,pim,pim,pim}}};
863
864//
865// Cross sections (in mb) for pi- p -> 2-9 body final states
866//
867// first index:    0-4: channels for mult = 2
868//                5-17: channels for mult = 3
869//               18-39: channels for mult = 4
870//               40-70: channels for mult = 5
871//               71-76: channels for mult = 6
872//               77-83: channels for mult = 7
873//               84-91: channels for mult = 8
874//              92-100: channels for mult = 9
875//
876// second index: kinetic energy
877//
878const G4float G4RPGPionInelastic::pimPCrossSections[101][30] = {
879//
880// multiplicity 2 (5 channels)
881//
882//  pi- p (pi+ n)
883//
884 {  0.0,  1.1,  1.2,  1.4,  1.5,  1.8,  2.0,  3.0,  3.4,  7.0,
885   14.0, 24.0, 14.7, 10.5, 11.8, 20.0, 14.0, 25.0, 12.0,  9.5,
886    8.0,  7.0,  6.0,  5.7,  5.0,  4.6,  4.3,  4.0,  3.8,  3.7 },
887
888//  n pi0  (p pi0)
889 {  0.0,  2.4,  2.8,  3.3,  4.5,  5.7,  6.3,  9.0, 11.0, 17.0,
890   30.0, 43.0, 30.0, 16.5, 11.0,  7.0,  4.3,  5.0,  2.0,  0.9,
891    0.5, 0.24, 0.15,0.094,0.061,0.048,0.035,0.023,0.018,0.014 },
892
893//  L K0  (L K+)
894 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
895    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.65, 0.28, 0.16,
896   0.13,0.075, 0.05,0.032,0.022,0.015,0.011,0.008,0.006,0.004 },
897
898//  S0 K0  (S0 K+)
899 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
900    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.28, 0.18, 0.11,
901  0.091,0.055,0.037,0.025,0.018,0.012,0.008,0.005,0.004,0.003 },
902
903//  S- K+  (S+ K0)
904 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
905    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.2, 0.24, 0.09,
906   0.04,0.012,0.004,0.002,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0 },
907//
908// multiplicity 3 (13 channels)
909//
910//  p pi- pi0  (n pi+ pi0)
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.05, 0.18, 0.86,  4.4,  5.2,  6.6,  5.4,  4.4,
913    3.5,  2.5,  2.0,  1.4, 0.97, 0.68, 0.55, 0.36,  0.3, 0.22 },
914
915//  n pi+ pi-  (p pi+ pi-)
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.09,  0.9,  3.3,  6.6,  8.9,  8.2,  8.8,  6.6,
918    5.2,  3.8,  2.9,  1.9,  1.3,  0.9, 0.75, 0.38, 0.31, 0.22 },
919
920//  n pi0 pi0  (p pi0 pi0)
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.22, 0.74,  1.8,  2.7,  3.0,  2.5,  1.3, 0.64,
923   0.41, 0.24, 0.15, 0.10,0.065,0.042,  0.0,  0.0,  0.0,  0.0 },
924
925//  L K0 pi0  (L K+ pi0)
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.049, 0.18,
928   0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01,0.007,0.004 },
929
930//  L K+ pi-  (L K0 pi+)
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.035, 0.12,
933   0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01,0.007,0.004 },
934
935//  S- K0 pi+  (S+ K+ pi-)
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.10,
938   0.13, 0.06, 0.03, 0.02, 0.01,0.003,0.001,  0.0,  0.0,  0.0 },
939
940//  S- K+ pi0  (S+ K0 pi0)
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.05,
943   0.05, 0.03, 0.02, 0.01,0.005,0.002,0.001,  0.0,  0.0,  0.0 },
944
945//  S+ K0 pi-  (S- K+ pi+)
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.04,
948   0.06, 0.04, 0.02, 0.01, 0.01,0.005,0.003,0.001,  0.0,  0.0 },
949
950//  S0 K+ pi-  (S0 K0 pi+)
951 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
952    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
953   0.07, 0.04, 0.03, 0.02, 0.02, 0.01,0.004,0.002,0.001,  0.0 },
954
955//  S0 K0 pi0  (S0 K+ pi0)
956 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
957    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.02, 0.09,
958   0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01,0.005,0.003,0.002 },
959
960//  p K0 K-  (n K+ K0bar)
961 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
962    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.03,
963   0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
964
965//  n K+ K-  (p K0 K0bar)
966 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
967    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
968   0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01,  0.0,  0.0,  0.0 },
969
970//  n K0 K0bar  (p K+ K-)
971 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
972    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.06,
973   0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01,  0.0,  0.0,  0.0 },
974//
975// multiplicity 4 (22 channels)
976//
977//  p pi+ pi- pi-  (n pi+ pi+ pi-)
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.01, 0.05,  0.2, 0.48, 0.88,  1.8,
980    2.0,  2.2,  2.0,  1.9,  1.6,  1.4,  1.3,  1.1, 1.05,  1.0 },
981
982//  p pi- pi0 pi0  (n pi+ pi0 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.01, 0.06, 0.18, 0.42,  1.0,  1.9,
985    2.4,  2.4,  2.0,  1.8,  1.6,  1.4,  1.3,  1.2,  1.1,  1.0 },
986
987//  n pi+ pi- pi0  (p pi+ pi- pi0)
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.1,  0.3,  0.5,  1.0,  2.0,  2.0,
990    1.0,  0.8,  0.6,  0.5,  0.4,  0.3, 0.25,  0.2, 0.15,  0.1 },
991
992//  n pi0 pi0 pi0  (p pi0 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.1, 0.55,  1.0, 0.87,  0.7,  0.6,
995    0.5,  0.4, 0.35, 0.29, 0.21, 0.14,  0.1, 0.06, 0.03,  0.0 },
996
997//  L K0 pi+ pi-  (L K+ pi+ pi-)
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.014,
1000   0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
1001
1002//  L K0 pi0 pi0  (L K+ pi0 pi0)
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.014,
1005   0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
1006
1007//  L K+ pi- pi0  (L K0 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.009,
1010   0.06,  0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02 },
1011
1012//  S0 K0 pi+ pi-  (S0 K+ 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.007,
1015   0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
1016
1017//  S0 K0 pi0 pi0  (S0 K+ 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.007,
1020    0.04,0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
1021
1022//  S0 K+ pi- pi0  (S0 K0 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//  S+ K+ pi- pi-  (S- K0 pi+ pi+)
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.005,
1030   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1031
1032//  S+ K0 pi- pi0  (S- K+ pi+ pi0)
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.005,
1035   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1036
1037//  S- K+ pi+ pi-  (S+ K0 pi+ pi-)
1038 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1039    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1040   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1041
1042//  S- K+ pi0 pi0  (S+ K0 pi0 pi0)
1043 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1044    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1045   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1046
1047//  S- K0 pi+ pi0  (S+ K+ pi- pi0)
1048 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1049    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
1050   0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1051
1052//  p pi- K+ K-  (n pi+ 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.003,0.06,  0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05 },
1056
1057//  p pi- K0 K0bar  (n pi+ K+ 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//  p pi0 K0 K-  (n pi0 K+ K0bar)
1063 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1064    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1065   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1066
1067//  n pi+ K0 K-  (p pi- K+ K0bar)
1068 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1069    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1070   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1071
1072//  n pi0 K0 K0bar  (p pi0 K+ K-)
1073 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1074    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1075   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1076
1077//  n pi0 K+ K-  (p pi0 K0 K0bar)
1078 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1079    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1080   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1081
1082//  n pi- K+ K0bar  (p pi+ K0 K-)
1083 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1084    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1085   0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1086//
1087// multiplicity 5 (31 channels)
1088//
1089//  p pi+ pi- pi- pi0  (n pi+ 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.007,0.022, 0.10, 0.70,
1092    1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
1093
1094//  p pi- pi0 pi0 pi0  (n pi+ pi0 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.007,0.022, 0.10, 0.70,
1097    1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
1098
1099//  n pi+ pi+ pi- pi-  (p pi+ 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.014,0.029, 0.10, 0.30,
1102   0.56, 0.93,  1.2,  1.2,  1.2, 0.94, 0.74, 0.53, 0.40, 0.30 },
1103
1104//  n pi+ pi- pi0 pi0  (p pi+ pi- 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.006, 0.02, 0.09, 0.63,
1107    1.2,  1.7,  2.0,  1.8,  1.5,  1.3,  1.1, 0.80, 0.70, 0.60 },
1108
1109//  n pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0)
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.006, 0.02, 0.08, 0.60,
1112    1.0,  1.5,  2.0,  1.6,  1.4,  1.1,  1.0, 0.70, 0.60, 0.50 },
1113
1114//  L K0 pi+ pi- pi0  (L K+ pi+ pi- 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.007, 0.05, 0.07, 0.08, 0.08, 0.07,0.063, 0.06,0.053,0.048 },
1118
1119//  L K+ pi- pi0 pi0  (L K0 pi+ pi0 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.025,0.035, 0.04, 0.04,0.035, 0.03, 0.03,0.025,0.024 },
1123
1124//  L K+ pi+ pi- pi-  (L K0 pi+ pi+ pi-)
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.002, 0.01,0.036, 0.04, 0.04,0.033,0.031,0.029,0.025,0.022 },
1128
1129//  L K0 pi0 pi0 pi0  (L K+ pi0 pi0 pi0)
1130 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1131    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1132  0.001,0.005,0.018, 0.02, 0.02,0.017,0.016,0.015,0.012,0.011 },
1133
1134//  S0 K+ pi+ pi- pi-  (S0 K0 pi+ pi+ pi-)
1135 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1136    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1137    0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1138
1139//  S0 K+ pi- pi0 pi0  (S0 K0 pi+ pi0 pi0)
1140 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1141    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1142    0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1143
1144//  S0 K0 pi+ pi- pi0  (S0 K+ pi+ pi- pi0)
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.003, 0.02, 0.03, 0.04, 0.04,0.035, 0.03, 0.03,0.025, 0.02 },
1148
1149//  S0 K0 pi0 pi0 pi0  (S0 K+ pi0 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.01,0.015, 0.02, 0.02,0.017,0.015,0.015,0.012, 0.01 },
1153
1154//  S+ K0 pi+ pi- pi-  (S- K+ pi+ pi+ pi-)
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.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1158
1159//  S+ K0 pi- pi0 pi0  (S- K+ pi+ 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.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1163
1164//  S+ K+ pi- pi- pi0  (S- K0 pi+ pi+ pi0)
1165 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1166    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1167  0.001,0.002,0.006,0.009,0.008,0.007,0.006,0.005,0.005,0.004 },
1168
1169//  S- K0 pi+ pi+ pi-  (S+ K+ pi+ pi- pi-)
1170 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1171    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1172  0.002,0.005,0.015,0.014,0.012,0.009,0.008,0.006,0.005,0.004 },
1173
1174//  S- K0 pi+ pi0 pi0  (S+ K+ pi- pi0 pi0)
1175 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1176    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1177  0.001,0.002,0.008,0.007,0.006,0.005,0.004,0.003,0.003,0.002 },
1178
1179//  S- K+ pi+ pi- pi0  (S+ K0 pi+ pi- pi0)
1180 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1181    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1182  0.002,0.007,0.023,0.025,0.021,0.018,0.015,0.013,0.010,0.009 },
1183
1184//  S- K+ pi0 pi0 pi0  (S+ K0 pi0 pi0 pi0)
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.001,0.004,0.011,0.013,0.011,0.009,0.008,0.007,0.006,0.005 },
1188
1189//  p pi- pi0 K+ K-  (n pi+ pi0 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.07, 0.06,0.055, 0.05,0.047,0.042, 0.04 },
1193
1194//  p pi- pi0 K0 K0bar  (n pi+ pi0 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.003, 0.02, 0.02,0.018,0.017,0.015,0.012,0.011, 0.01 },
1198
1199//  p pi+ pi- K0 K-  (n pi+ pi- 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.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1203
1204//  p pi0 pi0 K0 K-  (n pi0 pi0 K+ K0bar)
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.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1208
1209//  p pi- pi- K+ K0bar  (n pi+ pi+ K0 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.007, 0.04, 0.05, 0.04,0.035, 0.03, 0.03,0.027,0.025 },
1213
1214//  n pi+ pi- K+ K-  (p pi+ pi- 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.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1218
1219//  n pi+ pi- K0 K0bar  (p pi+ pi- K+ K-)
1220 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1221    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1222    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1223
1224//  n pi+ pi0 K0 K-  (p pi- pi0 K+ K0bar)
1225 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1226    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1227    0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1228
1229//  n pi- pi0 K+ K0bar  (p pi+ pi0 K0 K-)
1230 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1231    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1232    0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1233
1234//  n pi0 pi0 K0 K0bar  (p pi0 pi0 K+ K-)
1235 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1236    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1237    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1238
1239//  n pi0 pi0 K+ K-  (p pi0 pi0 K0 K0bar)
1240 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1241    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1242    0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1243//
1244// multiplicity 6 (6 channels)
1245//
1246//  p pi+ pi+ pi- pi- pi-  (n pi+ pi+ pi+ pi- pi-)
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.03, 0.04,
1249   0.06, 0.11, 0.16, 0.22, 0.31, 0.34,  0.3, 0.24, 0.19, 0.16 },
1250
1251//  p pi+ pi- pi- pi0 pi0  (n pi+ pi+ pi- pi0 pi0)
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.01, 0.04,  0.1,
1254   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1255
1256//  p pi- pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 pi0)
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.01, 0.04,  0.1,
1259   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1260
1261//  n pi+ pi+ pi- pi- pi0  (p pi+ pi+ pi- pi- pi0)
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.01, 0.04,  0.1,
1264   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38, 0.4,  0.33, 0.32 },
1265
1266//  n pi+ pi- pi0 pi0 pi0  (p pi+ pi- pi0 pi0 pi0)
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.01, 0.04,  0.1,
1269   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1270
1271//  n pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 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.01, 0.04,  0.1,
1274   0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
1275//
1276// multiplicity 7 (7 channels)
1277//
1278//  p pi+ pi+ pi- pi- pi- 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.014, 0.03,
1281   0.05, 0.11, 0.19, 0.37, 0.67, 0.75,  0.7, 0.59, 0.52, 0.47 },
1282
1283//  p pi+ pi- pi- 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//  p pi- pi0 pi0 pi0 pi0 pi0
1289 {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
1290    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
1291   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1292
1293//  n pi+ pi+ 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.002,0.006,
1296  0.012,0.029,0.062, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13 },
1297
1298//  n pi+ pi+ 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.008, 0.02,
1301   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1302
1303//  n pi+ pi- 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.008, 0.02,
1306   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1307
1308//  n pi0 pi0 pi0 pi0 pi0 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.008, 0.02,
1311   0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
1312//
1313// multiplicity 8 (8 channels)
1314//
1315//  p pi+ pi+ pi+ pi- pi- pi- pi-  (n pi+ pi+ pi+ pi+ pi- pi- pi-)
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//  p pi+ pi+ pi- pi- pi- pi0 pi0  (n pi+ pi+ pi+ pi- pi- 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//  p pi+ pi- pi- pi0 pi0 pi0 pi0  (n pi+ pi+ pi- pi0 pi0 pi0 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//  p pi- pi0 pi0 pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 pi0 pi0 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1334
1335//  n pi0 pi0 pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1339
1340//  n pi+ pi- pi0 pi0 pi0 pi0 pi0  (p pi+ pi- 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1344
1345//  n pi+ pi+ pi- pi- pi0 pi0 pi0  (p pi+ pi+ pi- pi- pi0 pi0 pi0)
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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1349
1350//  n pi+ pi+ pi+ pi- pi- pi- pi0  (p pi+ pi+ pi+ pi- pi- pi- 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.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1354//
1355// multiplicity 9 (9 channels)
1356//
1357//  p pi+ pi+ pi+ pi- pi- pi- pi- 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//  p pi+ pi+ pi- pi- pi- 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//  p pi+ pi- pi- pi0 pi0 pi0 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//  p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0
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//  n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
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.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1381
1382//  n pi+ pi- pi0 pi0 pi0 pi0 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.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1386
1387//  n pi+ pi+ 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.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1391
1392//  n pi+ pi+ pi+ pi- pi- pi- 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.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1396
1397//  n pi+ pi+ pi+ pi+ pi- pi- pi- pi-
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.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 } };
1401
1402 /* end of file */
1403 
Note: See TracBrowser for help on using the repository browser.