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