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