Geant4.10
 All Data Structures 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$
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 
233 ///////////////////////////////////////////////////////////////////////////////
234 // //
235 // pi+ p and pi- n (|Tz| = 3/2) cross sections //
236 // and final state particle types //
237 // //
238 ///////////////////////////////////////////////////////////////////////////////
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.00, 1.20, 2.50, 3.80, 5.00, 7.00, 9.00, 15.0, 30.0, 64.0,
357  130.0,190.0,130.0, 55.7, 27.2, 14.0, 8.50, 13.0, 18.0, 11.0,
358  8.50, 7.00, 6.20, 5.60, 5.00, 4.50, 4.20, 4.00, 3.80, 3.60 },
359 
360 // S+ K+ (S- K0)
361  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
362  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.16, 0.60, 0.32,
363  0.19, 0.1, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
364 //
365 // multiplicity 3 (7 channels)
366 //
367 // p pi+ pi0 (n pi- pi0)
368  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
369  0.0, 0.0, 0.0, 0.2, 0.6, 2.4, 8.8, 10.0, 12.0, 6.2,
370  4.00, 2.40, 1.69, 1.10, 0.73, 0.49, 0.41, 0.31, 0.24, 0.15 },
371 
372 // n pi+ pi+ (p pi- pi-)
373  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
374  0.0, 0.0, 0.0, 0.10, 0.20, 0.60, 1.50, 2.30, 3.60, 3.00,
375  2.30, 1.70, 1.30, 0.95, 0.69, 0.46, 0.38, 0.27, 0.20, 0.15 },
376 
377 // S+ K+ pi0 (S- K0 pi0)
378  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
379  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
380  0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
381 
382 // S+ K0 pi+
383  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
384  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
385  0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
386 
387 // S0 K+ pi+
388  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
389  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
390  0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
391 
392 // L K+ pi+
393  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
394  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005,0.12,
395  0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
396 
397 // p K+ K0bar
398  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
399  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
400  0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
401 //
402 // multiplicity 4 (15 channels)
403 //
404 // p pi+ pi+ pi-
405  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
406  0.0, 0.0, 0.0, 0.0, 0.0, 0.06, 0.20, 0.78, 2.20, 3.20,
407  3.50, 3.10, 2.70, 2.30, 2.00, 1.50, 1.40, 1.20, 1.00, 0.90 },
408 
409 // p pi+ pi0 pi0
410  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
411  0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
412  2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
413 
414 // n pi+ pi+ pi0
415  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
416  0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
417  2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
418 
419 // S+ K+ pi+ pi-
420  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
421  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
422  0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
423 
424 // S+ K+ pi0 pi0
425  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
426  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
427  0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
428 
429 // S+ K0 pi+ pi0
430  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
431  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,
432  0.015,0.06, 0.06, 0.05, 0.04,0.032,0.028, 0.02,0.017,0.014 },
433 
434 // S0 K0 pi+ pi+
435  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
436  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
437  0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
438 
439 // S0 K+ pi+ pi0
440  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
441  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
442  0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
443 
444 // L K+ pi+ pi0
445  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
446  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
447  0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
448 
449 // L K0 pi+ pi+
450  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
451  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
452  0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
453 
454 // S- K+ pi+ pi+
455  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
456  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,
457  0.02, 0.16, 0.18, 0.13, 0.09, 0.06,0.045, 0.03,0.025, 0.02 },
458 
459 // p pi+ K+ K-
460  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
461  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
462  0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05,0.045, 0.04 },
463 
464 // p pi+ K0 K0bar
465  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
466  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
467  0.005, 0.05, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05,0.045, 0.04 },
468 
469 // p pi0 K+ K0bar
470  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
471  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
472  0.01, 0.05,0.075, 0.08,0.075, 0.07, 0.06, 0.05,0.045, 0.04 },
473 
474 // n pi+ K+ K0bar
475  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
476  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
477  0.01, 0.05, 0.07,0.065,0.055,0.048, 0.04, 0.03,0.027,0.023 },
478 //
479 // multiplicity 5 (24 channels)
480 //
481 // p pi+ pi+ pi- pi0
482  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
483  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.05, 0.30, 2.00,
484  3.20, 3.70, 3.00, 2.50, 2.10, 1.60, 1.40, 1.10, 0.89, 0.70 },
485 
486 // p pi+ pi0 pi0 pi0
487  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
488  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.03, 0.21, 1.40,
489  2.20, 2.60, 2.10, 1.80, 1.50, 1.10, 1.00, 0.80, 0.62, 0.50 },
490 
491 // n pi+ pi+ pi+ pi-
492  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
493  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007, 0.02, 0.05, 0.19,
494  0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
495 
496 // n pi+ pi+ pi0 pi0
497  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
498  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007, 0.02, 0.05, 0.19,
499  0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
500 
501 // S+ K+ pi+ pi- pi0
502  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
503  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
504  0.003,0.013, 0.05, 0.09, 0.08, 0.06, 0.05, 0.04,0.036,0.03 },
505 
506 // S+ K+ pi0 pi0 pi0
507  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
508  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
509  0.001,0.006, 0.02, 0.04, 0.04, 0.03,0.025, 0.02, 0.02, 0.01 },
510 
511 // S+ K0 pi+ pi+ pi-
512  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
513  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
514  0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
515 
516 // S+ K0 pi+ pi0 pi0
517  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
518  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
519  0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
520 
521 // L K0 pi+ pi+ pi0
522  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
523  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
524  0.002,0.013, 0.04,0.052,0.059,0.053, 0.05,0.043,0.037, 0.03 },
525 
526 // L K+ pi+ pi+ pi-
527  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
528  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
529  0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
530 
531 // L K+ pi+ pi0 pi0
532  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
533  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
534  0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
535 
536 // S0 K+ pi+ pi+ pi-
537  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
538  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
539  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
540 
541 // S0 K+ pi+ pi0 pi0
542  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
543  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
544  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
545 
546 // S0 K0 pi+ pi+ pi0
547  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
548  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
549  0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
550 
551 // S- K+ pi+ pi+ pi0
552  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
553  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
554  0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
555 
556 // S- K0 pi+ pi+ pi+
557  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
558  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
559  0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
560 
561 // p pi+ pi- K+ K0bar
562  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
563  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
564  0.001,0.005, 0.02,0.065, 0.08, 0.07, 0.06,0.054,0.048, 0.04 },
565 
566 // p pi+ pi+ K0 K-
567  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
568  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
569  0.001,0.005, 0.02,0.045, 0.05,0.047, 0.04,0.033, 0.03,0.026 },
570 
571 // p pi+ pi0 K+ K-
572  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
573  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
574  0.001, 0.02, 0.06, 0.06,0.054,0.048,0.042,0.038,0.035, 0.03 },
575 
576 // p pi+ pi0 K0 K0bar
577  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
578  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
579  0.001, 0.01, 0.04, 0.05, 0.05, 0.04, 0.04,0.035,0.032, 0.03 },
580 
581 // p pi0 pi0 K+ K0bar
582  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
583  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
584  0.001,0.005, 0.01, 0.03, 0.04,0.035, 0.03,0.025,0.025, 0.02 },
585 
586 // n pi+ pi+ K+ K-
587  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
588  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
589  0.0, 0.02, 0.06, 0.05,0.042,0.038,0.035, 0.03,0.027,0.022 },
590 
591 // n pi+ pi+ K0 K0bar
592  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
593  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
594  0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
595 
596 // n pi+ pi0 K+ K0bar
597  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
598  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
599  0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
600 //
601 // multiplicity 6 (5 channels)
602 //
603 // p pi+ pi+ pi+ pi- pi-
604  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
605  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
606  0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
607 
608 // p pi+ pi+ pi- pi0 pi0
609  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
610  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
611  0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
612 
613 // p pi+ pi0 pi0 pi0 pi0
614  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
615  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
616  0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
617 
618 // n pi+ pi+ pi0 pi0 pi0
619  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
620  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
621  0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
622 
623 // n pi+ pi+ pi+ pi- pi0
624  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
625  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.001, 0.02,
626  0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
627 //
628 // multiplicity 7 (6 channels)
629 //
630 // p pi+ pi+ pi+ pi- pi- pi0
631  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
632  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
633  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
634 
635 // p pi+ pi+ pi- pi0 pi0 pi0
636  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
637  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
638  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
639 
640 // p pi+ pi0 pi0 pi0 pi0 pi0
641  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
642  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
643  0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
644 
645 // n pi+ pi+ pi+ pi+ pi- pi-
646  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
647  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
648  0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
649 
650 // n pi+ pi+ pi+ pi- pi0 pi0
651  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
652  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
653  0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
654 
655 // n pi+ pi+ pi0 pi0 pi0 pi0
656  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
657  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
658  0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
659 //
660 // multiplicity 8 (7 channels)
661 //
662 // p pi+ pi+ pi+ pi+ pi- pi- pi-
663  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
664  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
665  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
666 
667 // p pi+ pi+ pi+ pi- pi- pi0 pi0
668  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
669  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
670  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
671 
672 // p pi+ pi+ pi- pi0 pi0 pi0 pi0
673  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
674  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
675  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
676 
677 // p pi+ pi0 pi0 pi0 pi0 pi0 pi0
678  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
679  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
680  0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
681 
682 // n pi+ pi+ pi+ pi+ pi- pi- pi0
683  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
684  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
685  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
686 
687 // n pi+ pi+ pi+ pi- pi0 pi0 pi0
688  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
689  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
690  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
691 
692 // n pi+ pi+ pi0 pi0 pi0 pi0 pi0
693  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
694  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
695  0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
696 //
697 // multiplicity 9 (8 channels)
698 //
699 // p pi+ pi+ pi+ pi+ pi- pi- pi- pi0
700  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
701  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
702  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
703 
704 // p pi+ pi+ pi+ pi- pi- pi0 pi0 pi0
705  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
706  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
707  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
708 
709 // p pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0
710  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
711  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
712  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
713 
714 // p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
715  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
716  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
717  0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
718 
719 // n pi+ pi+ pi+ pi+ pi+ pi- pi- pi-
720  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
721  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
722  0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
723 
724 // n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
725  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
726  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
727  0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
728 
729 // n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
730  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
731  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
732  0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
733 
734 // n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
735  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
736  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
737  0.0, 0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 } };
738 
739 ///////////////////////////////////////////////////////////////////////////////
740 // //
741 // pi- p and pi+ n (|Tz| = 1/2) cross sections //
742 // and final state particle types //
743 // //
744 ///////////////////////////////////////////////////////////////////////////////
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.0, 1.1, 1.2, 1.4, 1.5, 1.8, 2.0, 3.0, 3.4, 7.0,
883  14.0, 24.0, 14.7, 10.5, 11.8, 20.0, 14.0, 25.0, 12.0, 9.5,
884  8.0, 7.0, 6.0, 5.7, 5.0, 4.6, 4.3, 4.0, 3.8, 3.7 },
885 
886 // n pi0 (p pi0)
887  { 0.0, 2.4, 2.8, 3.3, 4.5, 5.7, 6.3, 9.0, 11.0, 17.0,
888  30.0, 43.0, 30.0, 16.5, 11.0, 7.0, 4.3, 5.0, 2.0, 0.9,
889  0.5, 0.24, 0.15,0.094,0.061,0.048,0.035,0.023,0.018,0.014 },
890 
891 // L K0 (L K+)
892  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
893  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.65, 0.28, 0.16,
894  0.13,0.075, 0.05,0.032,0.022,0.015,0.011,0.008,0.006,0.004 },
895 
896 // S0 K0 (S0 K+)
897  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
898  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.28, 0.18, 0.11,
899  0.091,0.055,0.037,0.025,0.018,0.012,0.008,0.005,0.004,0.003 },
900 
901 // S- K+ (S+ K0)
902  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
903  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2, 0.24, 0.09,
904  0.04,0.012,0.004,0.002, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 },
905 //
906 // multiplicity 3 (13 channels)
907 //
908 // p pi- pi0 (n pi+ pi0)
909  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
910  0.0, 0.0, 0.05, 0.18, 0.86, 4.4, 5.2, 6.6, 5.4, 4.4,
911  3.5, 2.5, 2.0, 1.4, 0.97, 0.68, 0.55, 0.36, 0.3, 0.22 },
912 
913 // n pi+ pi- (p pi+ pi-)
914  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
915  0.0, 0.0, 0.09, 0.9, 3.3, 6.6, 8.9, 8.2, 8.8, 6.6,
916  5.2, 3.8, 2.9, 1.9, 1.3, 0.9, 0.75, 0.38, 0.31, 0.22 },
917 
918 // n pi0 pi0 (p pi0 pi0)
919  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
920  0.0, 0.0, 0.22, 0.74, 1.8, 2.7, 3.0, 2.5, 1.3, 0.64,
921  0.41, 0.24, 0.15, 0.10,0.065,0.042, 0.0, 0.0, 0.0, 0.0 },
922 
923 // L K0 pi0 (L K+ pi0)
924  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
925  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.049, 0.18,
926  0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01,0.007,0.004 },
927 
928 // L K+ pi- (L K0 pi+)
929  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
930  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.035, 0.12,
931  0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01,0.007,0.004 },
932 
933 // S- K0 pi+ (S+ K+ pi-)
934  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
935  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10,
936  0.13, 0.06, 0.03, 0.02, 0.01,0.003,0.001, 0.0, 0.0, 0.0 },
937 
938 // S- K+ pi0 (S+ K0 pi0)
939  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
940  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
941  0.05, 0.03, 0.02, 0.01,0.005,0.002,0.001, 0.0, 0.0, 0.0 },
942 
943 // S+ K0 pi- (S- K+ pi+)
944  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
945  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
946  0.06, 0.04, 0.02, 0.01, 0.01,0.005,0.003,0.001, 0.0, 0.0 },
947 
948 // S0 K+ pi- (S0 K0 pi+)
949  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
950  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
951  0.07, 0.04, 0.03, 0.02, 0.02, 0.01,0.004,0.002,0.001, 0.0 },
952 
953 // S0 K0 pi0 (S0 K+ pi0)
954  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
955  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.09,
956  0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01,0.005,0.003,0.002 },
957 
958 // p K0 K- (n K+ K0bar)
959  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
960  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03,
961  0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
962 
963 // n K+ K- (p K0 K0bar)
964  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
965  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
966  0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01, 0.0, 0.0, 0.0 },
967 
968 // n K0 K0bar (p K+ K-)
969  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
970  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06,
971  0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0 },
972 //
973 // multiplicity 4 (22 channels)
974 //
975 // p pi+ pi- pi- (n pi+ pi+ pi-)
976  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
977  0.0, 0.0, 0.0, 0.0, 0.01, 0.05, 0.2, 0.48, 0.88, 1.8,
978  2.0, 2.2, 2.0, 1.9, 1.6, 1.4, 1.3, 1.1, 1.05, 1.0 },
979 
980 // p pi- pi0 pi0 (n pi+ pi0 pi0)
981  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
982  0.0, 0.0, 0.0, 0.0, 0.01, 0.06, 0.18, 0.42, 1.0, 1.9,
983  2.4, 2.4, 2.0, 1.8, 1.6, 1.4, 1.3, 1.2, 1.1, 1.0 },
984 
985 // n pi+ pi- pi0 (p pi+ pi- pi0)
986  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
987  0.0, 0.0, 0.0, 0.0, 0.1, 0.3, 0.5, 1.0, 2.0, 2.0,
988  1.0, 0.8, 0.6, 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1 },
989 
990 // n pi0 pi0 pi0 (p pi0 pi0 pi0)
991  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
992  0.0, 0.0, 0.0, 0.0, 0.1, 0.55, 1.0, 0.87, 0.7, 0.6,
993  0.5, 0.4, 0.35, 0.29, 0.21, 0.14, 0.1, 0.06, 0.03, 0.0 },
994 
995 // L K0 pi+ pi- (L K+ pi+ pi-)
996  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
997  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,
998  0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
999 
1000 // L K0 pi0 pi0 (L K+ pi0 pi0)
1001  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1002  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,
1003  0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
1004 
1005 // L K+ pi- pi0 (L K0 pi+ pi0)
1006  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1007  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.009,
1008  0.06, 0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02 },
1009 
1010 // S0 K0 pi+ pi- (S0 K+ pi+ pi-)
1011  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1012  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,
1013  0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
1014 
1015 // S0 K0 pi0 pi0 (S0 K+ pi0 pi0)
1016  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1017  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,
1018  0.04,0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
1019 
1020 // S0 K+ pi- pi0 (S0 K0 pi+ pi0)
1021  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1022  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1023  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1024 
1025 // S+ K+ pi- pi- (S- K0 pi+ pi+)
1026  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1027  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1028  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1029 
1030 // S+ K0 pi- pi0 (S- K+ pi+ pi0)
1031  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1032  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1033  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1034 
1035 // S- K+ pi+ pi- (S+ K0 pi+ pi-)
1036  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1037  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1038  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1039 
1040 // S- K+ pi0 pi0 (S+ K0 pi0 pi0)
1041  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1042  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1043  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1044 
1045 // S- K0 pi+ pi0 (S+ K+ pi- pi0)
1046  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1047  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.005,
1048  0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
1049 
1050 // p pi- K+ K- (n pi+ K0 K0bar)
1051  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1052  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1053  0.003,0.06, 0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05 },
1054 
1055 // p pi- K0 K0bar (n pi+ K+ K-)
1056  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1057  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1058  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1059 
1060 // p pi0 K0 K- (n pi0 K+ K0bar)
1061  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1062  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1063  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1064 
1065 // n pi+ K0 K- (p pi- K+ K0bar)
1066  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1067  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1068  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1069 
1070 // n pi0 K0 K0bar (p pi0 K+ K-)
1071  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1072  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1073  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1074 
1075 // n pi0 K+ K- (p pi0 K0 K0bar)
1076  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1077  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1078  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1079 
1080 // n pi- K+ K0bar (p pi+ K0 K-)
1081  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1082  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1083  0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
1084 //
1085 // multiplicity 5 (31 channels)
1086 //
1087 // p pi+ pi- pi- pi0 (n pi+ pi+ pi- pi0)
1088  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1089  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,0.022, 0.10, 0.70,
1090  1.3, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62 },
1091 
1092 // p pi- pi0 pi0 pi0 (n pi+ pi0 pi0 pi0)
1093  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1094  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.007,0.022, 0.10, 0.70,
1095  1.3, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62 },
1096 
1097 // n pi+ pi+ pi- pi- (p pi+ pi+ pi- pi-)
1098  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1099  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014,0.029, 0.10, 0.30,
1100  0.56, 0.93, 1.2, 1.2, 1.2, 0.94, 0.74, 0.53, 0.40, 0.30 },
1101 
1102 // n pi+ pi- pi0 pi0 (p pi+ pi- pi0 pi0)
1103  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1104  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.006, 0.02, 0.09, 0.63,
1105  1.2, 1.7, 2.0, 1.8, 1.5, 1.3, 1.1, 0.80, 0.70, 0.60 },
1106 
1107 // n pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0)
1108  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1109  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.006, 0.02, 0.08, 0.60,
1110  1.0, 1.5, 2.0, 1.6, 1.4, 1.1, 1.0, 0.70, 0.60, 0.50 },
1111 
1112 // L K0 pi+ pi- pi0 (L K+ pi+ pi- pi0)
1113  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1114  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1115  0.007, 0.05, 0.07, 0.08, 0.08, 0.07,0.063, 0.06,0.053,0.048 },
1116 
1117 // L K+ pi- pi0 pi0 (L K0 pi+ pi0 pi0)
1118  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1119  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1120  0.003,0.025,0.035, 0.04, 0.04,0.035, 0.03, 0.03,0.025,0.024 },
1121 
1122 // L K+ pi+ pi- pi- (L K0 pi+ pi+ pi-)
1123  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1124  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1125  0.002, 0.01,0.036, 0.04, 0.04,0.033,0.031,0.029,0.025,0.022 },
1126 
1127 // L K0 pi0 pi0 pi0 (L K+ pi0 pi0 pi0)
1128  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1129  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1130  0.001,0.005,0.018, 0.02, 0.02,0.017,0.016,0.015,0.012,0.011 },
1131 
1132 // S0 K+ pi+ pi- pi- (S0 K0 pi+ pi+ pi-)
1133  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1134  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1135  0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1136 
1137 // S0 K+ pi- pi0 pi0 (S0 K0 pi+ pi0 pi0)
1138  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1139  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1140  0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
1141 
1142 // S0 K0 pi+ pi- pi0 (S0 K+ pi+ pi- pi0)
1143  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1144  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1145  0.003, 0.02, 0.03, 0.04, 0.04,0.035, 0.03, 0.03,0.025, 0.02 },
1146 
1147 // S0 K0 pi0 pi0 pi0 (S0 K+ pi0 pi0 pi0)
1148  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1149  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1150  0.001, 0.01,0.015, 0.02, 0.02,0.017,0.015,0.015,0.012, 0.01 },
1151 
1152 // S+ K0 pi+ pi- pi- (S- K+ pi+ pi+ pi-)
1153  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1154  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1155  0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1156 
1157 // S+ K0 pi- pi0 pi0 (S- K+ pi+ pi0 pi0)
1158  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1159  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1160  0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
1161 
1162 // S+ K+ pi- pi- pi0 (S- K0 pi+ pi+ pi0)
1163  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1164  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1165  0.001,0.002,0.006,0.009,0.008,0.007,0.006,0.005,0.005,0.004 },
1166 
1167 // S- K0 pi+ pi+ pi- (S+ K+ pi+ pi- pi-)
1168  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1169  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1170  0.002,0.005,0.015,0.014,0.012,0.009,0.008,0.006,0.005,0.004 },
1171 
1172 // S- K0 pi+ pi0 pi0 (S+ K+ pi- pi0 pi0)
1173  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1174  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1175  0.001,0.002,0.008,0.007,0.006,0.005,0.004,0.003,0.003,0.002 },
1176 
1177 // S- K+ pi+ pi- pi0 (S+ K0 pi+ pi- pi0)
1178  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1179  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1180  0.002,0.007,0.023,0.025,0.021,0.018,0.015,0.013,0.010,0.009 },
1181 
1182 // S- K+ pi0 pi0 pi0 (S+ K0 pi0 pi0 pi0)
1183  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1184  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1185  0.001,0.004,0.011,0.013,0.011,0.009,0.008,0.007,0.006,0.005 },
1186 
1187 // p pi- pi0 K+ K- (n pi+ pi0 K0 K0bar)
1188  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1189  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1190  0.0, 0.01, 0.07, 0.07, 0.06,0.055, 0.05,0.047,0.042, 0.04 },
1191 
1192 // p pi- pi0 K0 K0bar (n pi+ pi0 K+ K-)
1193  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1194  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1195  0.0,0.003, 0.02, 0.02,0.018,0.017,0.015,0.012,0.011, 0.01 },
1196 
1197 // p pi+ pi- K0 K- (n pi+ pi- K+ K0bar)
1198  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1199  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1200  0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1201 
1202 // p pi0 pi0 K0 K- (n pi0 pi0 K+ K0bar)
1203  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1204  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1205  0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
1206 
1207 // p pi- pi- K+ K0bar (n pi+ pi+ K0 K-)
1208  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1209  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1210  0.0,0.007, 0.04, 0.05, 0.04,0.035, 0.03, 0.03,0.027,0.025 },
1211 
1212 // n pi+ pi- K+ K- (p pi+ pi- K0 K0bar)
1213  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1214  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1215  0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1216 
1217 // n pi+ pi- K0 K0bar (p pi+ pi- K+ K-)
1218  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1219  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1220  0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1221 
1222 // n pi+ pi0 K0 K- (p pi- pi0 K+ K0bar)
1223  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1224  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1225  0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1226 
1227 // n pi- pi0 K+ K0bar (p pi+ pi0 K0 K-)
1228  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1229  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1230  0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
1231 
1232 // n pi0 pi0 K0 K0bar (p pi0 pi0 K+ K-)
1233  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1234  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1235  0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1236 
1237 // n pi0 pi0 K+ K- (p pi0 pi0 K0 K0bar)
1238  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1239  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1240  0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
1241 //
1242 // multiplicity 6 (6 channels)
1243 //
1244 // p pi+ pi+ pi- pi- pi- (n pi+ pi+ pi+ pi- pi-)
1245  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1246  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.04,
1247  0.06, 0.11, 0.16, 0.22, 0.31, 0.34, 0.3, 0.24, 0.19, 0.16 },
1248 
1249 // p pi+ pi- pi- pi0 pi0 (n pi+ pi+ pi- pi0 pi0)
1250  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1251  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
1252  0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
1253 
1254 // p pi- pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0)
1255  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1256  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
1257  0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
1258 
1259 // n pi+ pi+ pi- pi- pi0 (p pi+ pi+ pi- pi- pi0)
1260  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1261  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
1262  0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
1263 
1264 // n pi+ pi- pi0 pi0 pi0 (p pi+ pi- pi0 pi0 pi0)
1265  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1266  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
1267  0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
1268 
1269 // n pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 pi0)
1270  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1271  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.1,
1272  0.14, 0.19, 0.24, 0.31, 0.37, 0.4, 0.38, 0.4, 0.33, 0.32 },
1273 //
1274 // multiplicity 7 (7 channels)
1275 //
1276 // p pi+ pi+ pi- pi- pi- pi0
1277  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1278  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.014, 0.03,
1279  0.05, 0.11, 0.19, 0.37, 0.67, 0.75, 0.7, 0.59, 0.52, 0.47 },
1280 
1281 // p pi+ pi- pi- pi0 pi0 pi0
1282  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1283  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
1284  0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
1285 
1286 // p pi- pi0 pi0 pi0 pi0 pi0
1287  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1288  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
1289  0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
1290 
1291 // n pi+ pi+ pi+ pi- pi- pi-
1292  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1293  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.002,0.006,
1294  0.012,0.029,0.062, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13 },
1295 
1296 // n pi+ pi+ pi- pi- pi0 pi0
1297  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1298  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
1299  0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
1300 
1301 // n pi+ pi- pi0 pi0 pi0 pi0
1302  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1303  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
1304  0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
1305 
1306 // n pi0 pi0 pi0 pi0 pi0 pi0
1307  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1308  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,0.008, 0.02,
1309  0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3 },
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1315  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1316  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1317 
1318 // p pi+ pi+ pi- pi- pi- pi0 pi0 (n pi+ pi+ pi+ pi- pi- pi0 pi0)
1319  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1320  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1321  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1322 
1323 // p pi+ pi- pi- pi0 pi0 pi0 pi0 (n pi+ pi+ pi- pi0 pi0 pi0 pi0)
1324  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1325  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1326  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1327 
1328 // p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
1329  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1330  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1331  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1332 
1333 // n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 pi0 pi0 pi0)
1334  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1335  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1336  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1337 
1338 // n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p pi+ pi- pi0 pi0 pi0 pi0 pi0)
1339  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1340  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1341  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1342 
1343 // n pi+ pi+ pi- pi- pi0 pi0 pi0 (p pi+ pi+ pi- pi- pi0 pi0 pi0)
1344  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1345  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1346  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1347 
1348 // n pi+ pi+ pi+ pi- pi- pi- pi0 (p pi+ pi+ pi+ pi- pi- pi- pi0)
1349  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1350  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1351  0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1352 //
1353 // multiplicity 9 (9 channels)
1354 //
1355 // p pi+ pi+ pi+ pi- pi- pi- pi- pi0
1356  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1357  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1358  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1359 
1360 // p pi+ pi+ pi- pi- pi- pi0 pi0 pi0
1361  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1362  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1363  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1364 
1365 // p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0
1366  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1367  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1368  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1369 
1370 // p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0
1371  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1372  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1373  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1374 
1375 // n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
1376  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1377  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1378  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1379 
1380 // n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0
1381  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1382  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1383  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1384 
1385 // n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0
1386  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1387  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1388  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1389 
1390 // n pi+ pi+ pi+ pi- pi- pi- pi0 pi0
1391  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1392  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1393  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
1394 
1395 // n pi+ pi+ pi+ pi+ pi- pi- pi- pi-
1396  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1397  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1398  0.0, 0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 } };
1399 
1400  /* end of file */
1401 
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:523
G4RPGPionInelastic::T12_8bfs
static const G4int T12_8bfs[2][8][8]
Definition: G4RPGPionInelastic.hh:92
G4ThreadLocal
#define G4ThreadLocal
Definition: tls.hh:52
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
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:506
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
Generated on Wed Apr 30 2014 15:55:04 for Geant4.10 by   doxygen 1.8.7