source: trunk/source/processes/hadronic/models/rpg/src/G4RPGPionInelastic.cc@ 1347

Last change on this file since 1347 was 1340, checked in by garnier, 15 years ago

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