source: trunk/source/processes/hadronic/models/cascade/cascade/src/G4CascadeSigmaPlusPChannel.cc@ 819

Last change on this file since 819 was 819, checked in by garnier, 17 years ago

import all except CVS

File size: 31.3 KB
Line 
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#include "G4CascadeSigmaPlusPChannel.hh"
28
29
30G4CascadeSigmaPlusPChannel::G4CascadeSigmaPlusPChannel()
31 :G4CascadeChannel()
32{
33 G4int i, k, m;
34 G4int start, stop;
35
36 // Initialize multiplicity array
37
38 for (m = 0; m < 6; m++) {
39 start = sppindex[m][0];
40 stop = sppindex[m][1];
41 for (k = 0; k < 31; k++) {
42 sppMultiplicities[m][k] = 0.0;
43 for (i = start; i < stop; i++) sppMultiplicities[m][k] += sppCrossSections[i][k];
44 }
45 }
46
47 // Initialize total cross section array
48
49 for (k = 0; k < 31; k++) {
50 spptot[k] = 0.0;
51 for (m = 0; m < 6; m++) spptot[k] += sppMultiplicities[m][k];
52 }
53
54}
55
56
57G4CascadeSigmaPlusPChannel::~G4CascadeSigmaPlusPChannel()
58{;}
59
60
61G4double G4CascadeSigmaPlusPChannel::getCrossSection(G4double ke) const
62{
63 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
64 G4int k = epair.first;
65 G4double fraction = epair.second;
66
67 return spptot[k] + fraction*(spptot[k+1] - spptot[k]);
68}
69
70
71G4int G4CascadeSigmaPlusPChannel::getMultiplicity(G4double ke) const
72{
73 G4double multint(0.);
74 std::vector<G4double> sigma;
75
76 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
77 G4int k = epair.first;
78 G4double fraction = epair.second;
79
80 for(G4int m = 0; m < 6; m++) {
81 multint = sppMultiplicities[m][k]
82 + fraction*(sppMultiplicities[m][k+1] - sppMultiplicities[m][k]);
83 sigma.push_back(multint);
84 }
85
86 return sampleFlat(sigma);
87}
88
89
90std::vector<G4int>
91G4CascadeSigmaPlusPChannel::getOutgoingParticleTypes(G4int mult, G4double ke) const
92{
93 G4int i;
94 G4double sigint(0.);
95 std::vector<G4double> sigma;
96
97 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
98 G4int k = epair.first;
99 G4double fraction = epair.second;
100
101 G4int start = sppindex[mult-2][0];
102 G4int stop = sppindex[mult-2][1];
103
104 for(i = start; i < stop; i++) {
105 sigint = sppCrossSections[i][k]
106 + fraction*(sppCrossSections[i][k+1] - sppCrossSections[i][k]);
107 sigma.push_back(sigint);
108 }
109
110 G4int channel = sampleFlat(sigma);
111
112 std::vector<G4int> kinds;
113
114 if (mult == 2) {
115 for(i = 0; i < mult; i++) kinds.push_back(spp2bfs[channel][i]);
116 } else if (mult == 3) {
117 for(i = 0; i < mult; i++) kinds.push_back(spp3bfs[channel][i]);
118 } else if (mult == 4) {
119 for(i = 0; i < mult; i++) kinds.push_back(spp4bfs[channel][i]);
120 } else if (mult == 5) {
121 for(i = 0; i < mult; i++) kinds.push_back(spp5bfs[channel][i]);
122 } else if (mult == 6) {
123 for(i = 0; i < mult; i++) kinds.push_back(spp6bfs[channel][i]);
124 } else if (mult == 7) {
125 for(i = 0; i < mult; i++) kinds.push_back(spp7bfs[channel][i]);
126 } else {
127 G4cout << " Illegal multiplicity " << G4endl;
128 }
129
130 return kinds;
131}
132
133// Total cross section as a function of kinetic energy
134G4double G4CascadeSigmaPlusPChannel::spptot[31];
135
136// Multiplicities as a function of kinetic energy
137G4double G4CascadeSigmaPlusPChannel::sppMultiplicities[6][31];
138
139
140const G4int G4CascadeSigmaPlusPChannel::sppindex[6][2] =
141 {{0, 1}, {1, 7}, {7,27}, {27,69}, {69,94}, {94,111}};
142
143// Outgoing particle types of a given multiplicity
144
145const G4int G4CascadeSigmaPlusPChannel::spp2bfs[1][2] =
146 {{1,23}};
147
148const G4int G4CascadeSigmaPlusPChannel::spp3bfs[6][3] =
149 {{1,1,17}, {1,3,21}, {1,7,23}, {2,3,23}, {1,3,25}, {1,11,29}};
150
151const G4int G4CascadeSigmaPlusPChannel::spp4bfs[20][4] =
152 {{1,1,7,17}, {1,1,3,13}, {1,2,3,17}, {1,3,7,21}, {2,3,3,21},
153 {1,11,17,21}, {1,7,7,23}, {1,3,5,23}, {2,3,7,23}, {1,15,17,23},
154 {1,11,13,23}, {2,11,17,23}, {1,3,7,25}, {2,3,3,25}, {1,11,17,25},
155 {1,3,3,27}, {1,7,11,29}, {1,3,15,29},{2,3,11,29},{1,3,11,31}};
156
157const G4int G4CascadeSigmaPlusPChannel::spp5bfs[42][5] =
158 {{1,1,7,7,17}, {1,1,3,5,17}, {1,1,3,7,13}, {1,2,3,7,17},
159 {1,2,3,3,13}, {2,2,3,3,17}, {1,3,7,7,21}, {1,3,3,5,21},
160 {2,3,3,7,21}, {1,7,11,17,21}, {1,3,11,13,21}, {1,3,15,17,21},
161 {2,3,11,17,21}, {1,7,7,7,23}, {1,3,5,7,23}, {2,3,7,7,23},
162 {2,3,3,5,23}, {1,7,15,17,23}, {1,3,13,15,23}, {1,5,11,17,23},
163 {1,7,11,13,23}, {2,7,11,17,23}, {2,3,11,13,23}, {2,3,15,17,23},
164 {1,3,7,7,25}, {1,3,3,5,25}, {2,3,3,7,25}, {1,7,11,17,25},
165 {1,3,11,13,25}, {1,3,15,17,25}, {2,3,11,17,25}, {1,3,3,7,27},
166 {2,3,3,3,27}, {1,3,11,17,27}, {1,7,7,11,29}, {1,3,5,11,29},
167 {1,3,7,15,29}, {2,3,7,11,29}, {2,3,3,15,29}, {1,3,7,11,31},
168 {1,3,3,15,31}, {2,3,3,11,31}};
169
170const G4int G4CascadeSigmaPlusPChannel::spp6bfs[25][6] =
171 {{1,1,7,7,7,17}, {1,1,3,5,7,17}, {1,1,3,7,7,13}, {1,1,3,3,5,13},
172 {1,2,3,7,7,17}, {1,2,3,3,5,17}, {1,2,3,3,7,13}, {2,2,3,3,7,17},
173 {2,2,3,3,3,13}, {1,3,7,7,7,21}, {1,3,3,5,7,21}, {2,3,3,7,7,21},
174 {2,3,3,3,5,21}, {1,7,7,7,7,23}, {1,3,5,7,7,23}, {1,3,3,5,5,23},
175 {2,3,7,7,7,23}, {2,3,3,5,7,23}, {1,3,7,7,7,25}, {1,3,3,5,7,25},
176 {2,3,3,7,7,25}, {2,3,3,3,5,25}, {1,3,3,7,7,27}, {1,3,3,3,5,27},
177 {2,3,3,3,7,27}};
178
179const G4int G4CascadeSigmaPlusPChannel::spp7bfs[17][7] =
180 {{1,1,7,7,7,7,17}, {1,1,3,5,7,7,17}, {1,1,3,3,5,5,17},
181 {1,1,3,7,7,7,13}, {1,1,3,3,5,7,13}, {1,2,3,7,7,7,17},
182 {1,2,3,3,5,7,17}, {1,2,3,3,7,7,13}, {1,2,3,3,3,5,13},
183 {2,2,3,3,7,7,17}, {2,2,3,3,3,5,17}, {2,2,3,3,3,7,13},
184 {1,3,7,7,7,7,21}, {1,3,3,5,7,7,21}, {1,3,3,3,5,5,21},
185 {2,3,3,7,7,7,21}, {2,3,3,3,5,7,21}};
186
187// Cross sections for S+ p -> 2-7 body final states
188//
189// first index: 0: channel for mult = 2
190// 1-6: channels for mult = 3
191// 7-26: channels for mult = 4
192// 27-68: channels for mult = 5
193// 69-93: channels for mult = 6
194// 94-110: channels for mult = 7
195//
196// second index: kinetic energy
197//
198const G4float G4CascadeSigmaPlusPChannel::sppCrossSections[111][31] = {
199 //
200 // multiplicity 2 (1 channel)
201 //
202 // S+ p
203 {150.0,130.0,110.0, 90.0, 80.0, 70.0, 60.0, 50.0, 40.0, 45.0,
204 40.0, 23.30,20.14,16.90,11.80, 9.80, 7.80, 7.00, 6.60, 6.40,
205 6.20, 6.10, 5.90, 5.70, 5.50, 5.30, 5.10, 5.00, 4.98, 4.96, 4.70},
206
207 //
208 // multiplicity 3 (6 channels)
209 //
210 // p p K0bar
211 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
212 0.0, 1.00, 1.50, 1.80, 1.70, 1.50, 1.26, 1.02, 0.93, 0.84,
213 0.76, 0.68, 0.61, 0.53, 0.45, 0.38, 0.31, 0.23, 0.21, 0.19, 0.09},
214
215 // L p pi+
216 { 0.0, 0.0, 0.0, 0.10, 0.20, 0.30, 0.50, 0.70, 1.00, 1.20,
217 1.30, 2.50, 4.00, 4.50, 4.30, 4.00, 3.32, 2.68, 2.45, 2.20,
218 1.98, 1.78, 1.60, 1.38, 1.18, 0.96, 0.77, 0.58, 0.53, 0.49, 0.22},
219
220 // S+ p pi0
221 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
222 0.10, 0.90, 1.50, 1.80, 1.70, 1.50, 1.26, 1.02, 0.93, 0.84,
223 0.76, 0.68, 0.61, 0.53, 0.45, 0.38, 0.31, 0.23, 0.21, 0.19, 0.09},
224
225 // S+ n pi+
226 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
227 0.10, 0.90, 1.50, 1.80, 1.70, 1.50, 1.26, 1.02, 0.93, 0.84,
228 0.76, 0.68, 0.61, 0.53, 0.45, 0.38, 0.31, 0.23, 0.21, 0.19, 0.09},
229
230 // S0 p pi+
231 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
232 0.10, 0.90, 1.50, 1.80, 1.70, 1.50, 1.26, 1.02, 0.93, 0.84,
233 0.76, 0.68, 0.61, 0.53, 0.45, 0.38, 0.31, 0.23, 0.21, 0.19, 0.09},
234
235 // X0 p K+
236 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
237 0.0, 0.0, 0.0, 0.02, 0.03, 0.03, 0.04, 0.04, 0.03, 0.03,
238 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
239
240 //
241 // multiplicity 4 (20 channels)
242 //
243 // p p K0bar pi0
244 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
245 0.0, 0.10, 0.25, 0.40, 0.60, 0.80, 1.00, 1.20, 1.40, 1.54,
246 1.49, 1.41, 1.36, 1.29, 1.24, 1.21, 1.15, 1.11, 1.11, 1.10, 0.98},
247
248 // p p K- pi+
249 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
250 0.0, 0.10, 0.25, 0.40, 0.60, 0.80, 1.00, 1.20, 1.40, 1.54,
251 1.49, 1.41, 1.36, 1.29, 1.24, 1.21, 1.15, 1.11, 1.11, 1.10, 0.98},
252
253 // n p K0bar pi+
254 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
255 0.0, 0.10, 0.25, 0.40, 0.60, 0.80, 1.00, 1.20, 1.40, 1.54,
256 1.49, 1.41, 1.36, 1.29, 1.24, 1.21, 1.15, 1.11, 1.11, 1.10, 0.98},
257
258 // L p pi+ pi0
259 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.06,
260 0.10, 0.25, 0.40, 0.50, 0.60, 0.75, 0.90, 1.05, 1.20, 1.32,
261 1.28, 1.22, 1.18, 1.12, 1.08, 1.05, 1.00, 0.97, 0.97, 0.96, 0.85},
262
263 // L n 2pi+
264 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.06,
265 0.10, 0.25, 0.40, 0.50, 0.60, 0.75, 0.90, 1.05, 1.20, 1.32,
266 1.28, 1.22, 1.18, 1.12, 1.08, 1.05, 1.00, 0.97, 0.97, 0.96, 0.85},
267
268 // L p K+ K0bar
269 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
270 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
271 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
272
273 // S+ p 2pi0
274 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
275 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
276 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
277
278 // S+ p pi+ pi-
279 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
280 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
281 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
282
283 // S+ n pi+ pi0
284 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
285 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
286 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
287
288 // S+ p K0 K0bar
289 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
290 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
291 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
292
293 // S+ p K+ K-
294 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
295 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
296 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
297
298 // S+ n K+ K0bar
299 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
300 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
301 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
302
303 // S0 p pi+ pi0
304 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
305 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
306 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
307
308 // S0 n 2pi+
309 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
310 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
311 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
312
313 // S0 p K+ K0bar
314 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
315 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
316 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
317
318 // S- p 2pi+
319 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
320 0.0, 0.03, 0.07, 0.13, 0.25, 0.40, 0.60, 0.80, 1.00, 1.10,
321 1.07, 1.02, 0.99, 0.94, 0.90, 0.88, 0.84, 0.81, 0.81, 0.80, 0.71},
322
323 // X0 p K+ pi0
324 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
325 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
326 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
327
328 // X0 p K0 pi+
329 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
330 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
331 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
332
333 // X0 n K+ pi+
334 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
335 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
336 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
337
338 // X- p K+ pi+
339 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
340 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
341 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
342
343 //
344 // multiplicity 5 (42 channels)
345 //
346 //
347 // p p K0bar 2pi0
348 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
349 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.06, 0.07, 0.09, 0.13,
350 0.15, 0.17, 0.18, 0.22, 0.21, 0.21, 0.21, 0.21, 0.21, 0.21, 0.20},
351
352 // p p K0bar pi+ pi-
353 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
354 0.0, 0.0, 0.0, 0.0, 0.03, 0.06, 0.10, 0.11, 0.13, 0.20,
355 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
356
357 // p p K- pi+ pi0
358 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
359 0.0, 0.0, 0.0, 0.0, 0.03, 0.06, 0.10, 0.11, 0.13, 0.20,
360 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
361
362 // n p K0bar pi+ pi0
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.03, 0.06, 0.10, 0.11, 0.13, 0.20,
365 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
366
367 // n p K- 2pi+
368 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
369 0.0, 0.0, 0.0, 0.0, 0.03, 0.06, 0.10, 0.11, 0.13, 0.20,
370 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
371
372 // n n K0bar 2pi+
373 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
374 0.0, 0.0, 0.0, 0.0, 0.03, 0.06, 0.10, 0.11, 0.13, 0.20,
375 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
376
377 // L p pi+ 2pi0
378 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
379 0.0, 0.0, 0.04, 0.04, 0.08, 0.12, 0.16, 0.17, 0.18, 0.23,
380 0.26, 0.30, 0.32, 0.37, 0.36, 0.36, 0.36, 0.36, 0.35, 0.35, 0.33},
381
382 // L p 2pi+ pi-
383 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
384 0.0, 0.0, 0.04, 0.04, 0.14, 0.18, 0.24, 0.26, 0.28, 0.35,
385 0.40, 0.45, 0.48, 0.56, 0.54, 0.54, 0.54, 0.53, 0.52, 0.51, 0.47},
386
387 // L n 2pi+ pi0
388 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
389 0.0, 0.0, 0.04, 0.04, 0.14, 0.18, 0.24, 0.26, 0.28, 0.35,
390 0.40, 0.45, 0.48, 0.56, 0.54, 0.54, 0.54, 0.53, 0.52, 0.51, 0.47},
391
392 // L p K+ K0bar pi0
393 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
394 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
395 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
396
397 // L p K+ K- pi+
398 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
399 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
400 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
401
402 // L p K0 K0bar pi+
403 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
404 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
405 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
406
407 // L n K+ K0bar pi+
408 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
409 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
410 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
411
412 // S+ p 3pi0
413 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
414 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.06, 0.07, 0.09, 0.13,
415 0.15, 0.17, 0.18, 0.22, 0.21, 0.21, 0.21, 0.21, 0.21, 0.21, 0.20},
416
417 // S+ p pi+ pi- pi0
418 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
419 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
420 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
421
422 // S+ n pi+ 2pi0
423 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
424 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.06, 0.07, 0.09, 0.13,
425 0.15, 0.17, 0.18, 0.22, 0.21, 0.21, 0.21, 0.21, 0.21, 0.21, 0.20},
426
427 // S+ n 2pi+ pi-
428 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
429 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
430 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
431
432 // S+ p K0 K0bar pi0
433 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
434 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
435 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
436
437 // S+ p K0 K- pi+
438 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
439 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
440 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
441
442 // S+ p K+ K0bar pi-
443 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
444 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
445 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
446
447 // S+ p K+ K- pi0
448 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
449 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
450 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
451
452 // S+ n K+ K0bar pi0
453 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
454 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
455 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
456
457 // S+ n K+ K- pi+
458 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
459 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
460 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
461
462 // S+ n K0 K0bar pi+
463 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
464 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
465 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
466
467 // S0 p pi+ 2pi0
468 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
469 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.06, 0.07, 0.09, 0.13,
470 0.15, 0.17, 0.18, 0.22, 0.21, 0.21, 0.21, 0.21, 0.21, 0.21, 0.20},
471
472 // S0 p 2pi+ pi-
473 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
474 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
475 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
476
477 // S0 n 2pi+ pi0
478 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
479 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
480 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
481
482 // S0 p K+ K0bar pi0
483 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
484 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
485 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
486
487 // S0 p K+ K- pi+
488 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
489 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
490 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
491
492 // S0 p K0 K0bar pi+
493 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
494 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
495 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
496
497 // S0 n K+ K0bar pi+
498 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
499 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
500 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
501
502 // S- p 2pi+ pi0
503 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
504 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
505 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
506
507 // S- n 3pi+
508 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
509 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.10, 0.11, 0.13, 0.20,
510 0.23, 0.26, 0.28, 0.33, 0.32, 0.32, 0.32, 0.32, 0.32, 0.32, 0.30},
511
512 // S- p K+ K0bar pi+
513 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
514 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
515 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
516
517 // X0 p K+ 2pi0
518 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
519 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
520 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
521
522 // X0 p K+ pi+ pi-
523 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
524 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
525 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
526
527 // X0 p K0 pi+ pi0
528 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
529 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
530 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
531
532 // X0 n K+ pi+ pi0
533 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
534 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
535 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
536
537 // X0 n K0 2pi+
538 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
539 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
540 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
541
542 // X- p K+ pi+ pi0
543 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
544 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
545 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
546
547 // X- p K0 2pi+
548 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
549 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
550 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
551
552 // X- n K+ 2pi+
553 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
554 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02,
555 0.02, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01},
556
557 //
558 // multiplicity 6 (25 channels)
559 //
560 // p p K0bar 3pi0
561 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
562 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
563 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
564
565 // p p K0bar pi+ pi- pi0
566 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
567 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
568 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
569
570 // p p K- pi+ 2pi0
571 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
572 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
573 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
574
575 // p p K- 2pi+ pi-
576 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
577 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
578 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
579
580 // n p K0bar pi+ 2pi0
581 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
582 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
583 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
584
585 // n p K0bar 2pi+ pi-
586 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
587 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
588 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
589
590 // n p K- 2pi+ pi0
591 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
592 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
593 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
594
595 // n n K0bar 2pi+ pi0
596 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
597 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
598 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
599
600 // n n K- 3pi+
601 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
602 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
603 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
604
605 // L p pi+ 3pi0
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.03, 0.06, 0.08, 0.10, 0.12, 0.14,
608 0.16, 0.20, 0.24, 0.28, 0.33, 0.34, 0.36, 0.38, 0.38, 0.38, 0.34},
609
610 // L p 2pi+ pi- 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.05, 0.10, 0.13, 0.15, 0.18, 0.22,
613 0.25, 0.31, 0.36, 0.42, 0.50, 0.51, 0.54, 0.57, 0.58, 0.58, 0.52},
614
615 // L n 2pi+ 2pi0
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.03, 0.06, 0.08, 0.10, 0.12, 0.14,
618 0.16, 0.20, 0.24, 0.28, 0.33, 0.34, 0.36, 0.38, 0.38, 0.38, 0.34},
619
620 // L n 3pi+ pi-
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.05, 0.10, 0.13, 0.15, 0.18, 0.22,
623 0.25, 0.31, 0.36, 0.42, 0.50, 0.51, 0.54, 0.57, 0.58, 0.58, 0.52},
624
625 // S+ p 4pi0
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.01, 0.02, 0.03, 0.04, 0.05, 0.06,
628 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
629
630 // S+ p pi+ pi- 2pi0
631 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
632 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.05, 0.06, 0.07, 0.08,
633 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
634
635 // S+ p 2pi+ 2pi-
636 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
637 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
638 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
639
640 // S+ n pi+ 3pi0
641 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
642 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
643 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
644
645 // S+ n 2pi+ pi- pi0
646 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
647 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
648 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
649
650 // S0 p pi+ 3pi0
651 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
652 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
653 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
654
655 // S0 p 2pi+ pi- pi0
656 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
657 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
658 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
659
660 // S0 n 2pi+ 2pi0
661 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
662 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
663 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
664
665 // S0 n 3pi+ pi-
666 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
667 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
668 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
669
670 // S- p 2pi+ 2pi0
671 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
672 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
673 0.06, 0.07, 0.08, 0.10, 0.12, 0.12, 0.13, 0.13, 0.13, 0.13, 0.12},
674
675 // S- p 3pi+ pi-
676 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
677 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
678 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
679
680 // S- n 3pi+ pi0
681 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
682 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.06, 0.07, 0.08,
683 0.09, 0.11, 0.13, 0.15, 0.18, 0.18, 0.19, 0.20, 0.20, 0.20, 0.18},
684
685 //
686 // multiplicity 7 (17 channels)
687 //
688 // p p K0bar 4pi0
689 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
690 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
691 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
692
693 // p p K0bar pi+ pi- 2pi0
694 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
695 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
696 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
697
698 // p p K0bar 2pi+ 2pi-
699 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
700 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
701 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
702
703 // p p K- pi+ 3pi0
704 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
705 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
706 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
707
708 // p p K- 2pi+ pi- pi0
709 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
710 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
711 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
712
713 // n p K0bar pi+ 3pi0
714 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
715 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
716 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
717
718 // n p K0bar 2pi+ pi- pi0
719 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
720 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
721 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
722
723 // n p K- 2pi+ 2pi0
724 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
725 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
726 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
727
728 // n p K- 3pi+ pi-
729 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
730 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
731 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
732
733 // n n K0bar 2pi+ 2pi0
734 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
735 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03,
736 0.04, 0.04, 0.05, 0.07, 0.08, 0.09, 0.10, 0.12, 0.12, 0.12, 0.12},
737
738 // n n K0bar 3pi+ pi-
739 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
740 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
741 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
742
743 // n n K- 3pi+ pi0
744 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
745 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.04,
746 0.05, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.19, 0.19, 0.19, 0.19},
747
748 // L p pi+ 4pi0
749 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
750 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.04, 0.05, 0.06,
751 0.08, 0.10, 0.13, 0.14, 0.17, 0.19, 0.22, 0.26, 0.26, 0.26, 0.26},
752
753 // L p 2pi+ pi- 2pi0
754 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
755 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.07, 0.08,
756 0.11, 0.13, 0.17, 0.22, 0.26, 0.29, 0.33, 0.40, 0.40, 0.40, 0.40},
757
758 // L p 3pi+ 2pi-
759 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
760 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.06, 0.09, 0.11, 0.12,
761 0.16, 0.20, 0.26, 0.33, 0.40, 0.47, 0.55, 0.60, 0.60, 0.60, 0.60},
762
763 // L n 2pi+ 3pi0
764 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
765 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.07, 0.08,
766 0.11, 0.13, 0.17, 0.22, 0.26, 0.29, 0.33, 0.40, 0.40, 0.40, 0.40},
767
768 // L n 3pi+ pi- pi0
769 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
770 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.06, 0.09, 0.11, 0.12,
771 0.16, 0.20, 0.26, 0.33, 0.40, 0.47, 0.55, 0.60, 0.60, 0.60, 0.60}};
Note: See TracBrowser for help on using the repository browser.