Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
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