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