Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
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
33G4RPGPionInelastic::G4RPGPionInelastic(const G4String& modelName)
34 :G4RPGInelastic(modelName)
35{
36 SetMinEnergy( 0.0 );
37 SetMaxEnergy( 30.*GeV );
38
39 // Initialize t32_dSigma_dMult, t12_dSigma_dMult,
40 // pi-nucleon inelastic cross sections for a given multiplicity
41 // for |T_z| = 3/2 and 1/2, respectively
42
43 G4int i, k, 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/*
77void 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
99G4int 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
118G4int 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
137std::vector<G4int>
138G4RPGPionInelastic::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
185std::vector<G4int>
186G4RPGPionInelastic::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
241G4double G4RPGPionInelastic::pipPtot[30];
242
243// pi+ multiplicities as a function of kinetic energy
244G4double G4RPGPionInelastic::t32_dSigma_dMult[8][30];
245
246const 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
252const G4int G4RPGPionInelastic::T32_2bfs[2][2][2] =
253 {{{pro,pip}, {sp,kp}},
254
255 {{neu,pim}, {sm,k0}}};
256
257const 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
264const 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
277const 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
296const 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
305const 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
314const 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
325const 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
351const 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
747G4double G4RPGPionInelastic::pimPtot[30];
748
749// pi- multiplicities as a function of kinetic energy
750G4double G4RPGPionInelastic::t12_dSigma_dMult[8][30];
751
752const 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
758const 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
763const 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
774const 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
793const 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
818const 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
827const 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
838const 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
849const 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//
876const 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
G4double
double G4double
Definition: G4Types.hh:64
G4float
float G4float
Definition: G4Types.hh:65
G4int
int G4int
Definition: G4Types.hh:66
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGPionInelastic::T32_8bfs
static const G4int T32_8bfs[2][7][8]
Definition: G4RPGPionInelastic.hh:83
G4RPGPionInelastic::GetFSPartTypesForT32
std::vector< G4int > GetFSPartTypesForT32(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:138
G4RPGPionInelastic::pimPCrossSections
static const G4float pimPCrossSections[101][30]
Definition: G4RPGPionInelastic.hh:101
G4RPGPionInelastic::t12_dSigma_dMult
static G4double t12_dSigma_dMult[8][30]
Definition: G4RPGPionInelastic.hh:97
G4RPGPionInelastic::T12_4bfs
static const G4int T12_4bfs[2][22][4]
Definition: G4RPGPionInelastic.hh:88
G4RPGPionInelastic::T12_3bfs
static const G4int T12_3bfs[2][13][3]
Definition: G4RPGPionInelastic.hh:87
G4RPGPionInelastic::GetMultiplicityT32
G4int GetMultiplicityT32(G4double KE) const
Definition: G4RPGPionInelastic.cc:118
G4RPGPionInelastic::T12_9bfs
static const G4int T12_9bfs[2][9][9]
Definition: G4RPGPionInelastic.hh:93
G4RPGPionInelastic::T32_6bfs
static const G4int T32_6bfs[2][5][6]
Definition: G4RPGPionInelastic.hh:81
G4RPGPionInelastic::GetMultiplicityT12
G4int GetMultiplicityT12(G4double KE) const
Definition: G4RPGPionInelastic.cc:99
G4RPGPionInelastic::T12_2bfs
static const G4int T12_2bfs[2][5][2]
Definition: G4RPGPionInelastic.hh:86
G4RPGPionInelastic::T32_5bfs
static const G4int T32_5bfs[2][24][5]
Definition: G4RPGPionInelastic.hh:80
G4RPGPionInelastic::T12_8bfs
static const G4int T12_8bfs[2][8][8]
Definition: G4RPGPionInelastic.hh:92
G4RPGPionInelastic::pipPCrossSections
static const G4float pipPCrossSections[74][30]
Definition: G4RPGPionInelastic.hh:100
G4RPGPionInelastic::T32_3bfs
static const G4int T32_3bfs[2][7][3]
Definition: G4RPGPionInelastic.hh:78
G4RPGPionInelastic::GetFSPartTypesForT12
std::vector< G4int > GetFSPartTypesForT12(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:186
G4RPGPionInelastic::t32_dSigma_dMult
static G4double t32_dSigma_dMult[8][30]
Definition: G4RPGPionInelastic.hh:98
G4RPGPionInelastic::T32_2bfs
static const G4int T32_2bfs[2][2][2]
Definition: G4RPGPionInelastic.hh:77
G4RPGPionInelastic::pimPtot
static G4double pimPtot[30]
Definition: G4RPGPionInelastic.hh:96
G4RPGPionInelastic::G4RPGPionInelastic
G4RPGPionInelastic(const G4String &modelName="RPGPionInelastic")
Definition: G4RPGPionInelastic.cc:33
G4RPGPionInelastic::T12_6bfs
static const G4int T12_6bfs[2][6][6]
Definition: G4RPGPionInelastic.hh:90
G4RPGPionInelastic::T12_5bfs
static const G4int T12_5bfs[2][31][5]
Definition: G4RPGPionInelastic.hh:89
G4RPGPionInelastic::T32_7bfs
static const G4int T32_7bfs[2][6][7]
Definition: G4RPGPionInelastic.hh:82
G4RPGPionInelastic::T12_7bfs
static const G4int T12_7bfs[2][7][7]
Definition: G4RPGPionInelastic.hh:91
G4RPGPionInelastic::pipPtot
static G4double pipPtot[30]
Definition: G4RPGPionInelastic.hh:95
G4RPGPionInelastic::T32_9bfs
static const G4int T32_9bfs[2][8][9]
Definition: G4RPGPionInelastic.hh:84
G4RPGPionInelastic::pipPindex
static const G4int pipPindex[8][2]
Definition: G4RPGPionInelastic.hh:74
G4RPGPionInelastic::pimPindex
static const G4int pimPindex[8][2]
Definition: G4RPGPionInelastic.hh:75
G4RPGPionInelastic::T32_4bfs
static const G4int T32_4bfs[2][15][4]
Definition: G4RPGPionInelastic.hh:79