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

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

update ti head
File size: 65.8 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: G4RPGNucleonInelastic.cc,v 1.2 2009/05/25 19:06:03 dennis Exp $
[1340]27// GEANT4 tag $Name: geant4-09-03-ref-09 $
[968]28//
29
30#include "G4RPGNucleonInelastic.hh"
31#include "Randomize.hh"
32
33G4RPGNucleonInelastic::G4RPGNucleonInelastic(const G4String& modelName)
34 :G4RPGInelastic(modelName)
35{
36 SetMinEnergy( 0.0 );
37 SetMaxEnergy( 30.*GeV );
38
39 // Initialize t1_dSigma_dMult, t0_dSigma_dMult,
40 // nucleon-nucleon inelastic cross sections for a given multiplicity
41 // for |T_z| = 1 and 0, respectively
42
43 G4int i, k, m;
44 G4int start, stop;
45
46 for (m = 0; m < 8; m++) {
47 start = pPindex[m][0];
48 stop = pPindex[m][1] + 1;
49 for (k = 0; k < 30; k++) {
50 t1_dSigma_dMult[m][k] = 0.0;
51 for (i = start; i < stop; i++) t1_dSigma_dMult[m][k] += pPCrossSections[i][k];
52 }
53
54 start = pNindex[m][0];
55 stop = pNindex[m][1] + 1;
56 for (k = 0; k < 30; k++) {
57 t0_dSigma_dMult[m][k] = 0.0;
58 for (i = start; i < stop; i++) t0_dSigma_dMult[m][k] += pNCrossSections[i][k];
59 }
60 }
61
62 // Initialize total cross section array
63
64 for (k = 0; k < 30; k++) {
65 pPtot[k] = 0.0;
66 pNtot[k] = 0.0;
67 for (m = 0; m < 8; m++) {
68 pPtot[k] += t1_dSigma_dMult[m][k];
69 pNtot[k] += t0_dSigma_dMult[m][k];
70 }
71 }
[1055]72
73 // printCrossSections();
[968]74
75}
76
[1055]77/*
78void G4RPGNucleonInelastic::printCrossSections() const
79{
80 G4cout << " pp total cross sections (mb) " << G4endl;
81 for (G4int i = 0; i < 5; i++) {
82 G4int istart = i*6;
83 G4int istop = istart + 6;
84 for (G4int t = istart; t < istop; t++) G4cout << pPtot[t] << " " ;
85 G4cout << G4endl;
86 }
87 G4cout << G4endl;
[968]88
[1055]89 G4cout << " pn total cross sections (mb) " << G4endl;
90 for (G4int i = 0; i < 5; i++) {
91 G4int istart = i*6;
92 G4int istop = istart + 6;
93 for (G4int t = istart; t < istop; t++) G4cout << pNtot[t] << " " ;
94 G4cout << G4endl;
95 }
96}
97*/
98
99
[968]100G4int G4RPGNucleonInelastic::GetMultiplicityT0(G4double KE) const
101{
102 G4double multint(0.);
103 std::vector<G4double> sigma;
104
105 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
106 G4int k = epair.first;
107 G4double fraction = epair.second;
108
109 for(G4int m = 0; m < 8; m++) {
110 multint = t0_dSigma_dMult[m][k]
111 + fraction*(t0_dSigma_dMult[m][k+1] - t0_dSigma_dMult[m][k]);
112 sigma.push_back(multint);
113 }
114
115 return sampleFlat(sigma) + 2;
116}
117
118
119G4int G4RPGNucleonInelastic::GetMultiplicityT1(G4double KE) const
120{
121 G4double multint(0.);
122 std::vector<G4double> sigma;
123
124 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
125 G4int k = epair.first;
126 G4double fraction = epair.second;
127
128 for(G4int m = 0; m < 8; m++) {
129 multint = t1_dSigma_dMult[m][k]
130 + fraction*(t1_dSigma_dMult[m][k+1] - t1_dSigma_dMult[m][k]);
131 sigma.push_back(multint);
132 }
133
134 return sampleFlat(sigma) + 2;
135}
136
137
138std::vector<G4int>
139G4RPGNucleonInelastic::GetFSPartTypesForT0(G4int mult, G4double KE) const
140{
141 G4int i;
142 G4double sigint(0.);
143 std::vector<G4double> sigma;
144
145 std::pair<G4int, G4double> epair = interpolateEnergy(KE);
146 G4int k = epair.first;
147 G4double fraction = epair.second;
148
149 G4int start = pNindex[mult-2][0];
150 G4int stop = pNindex[mult-2][1];
151
152 for(i = start; i < stop; i++) {
153 sigint = pNCrossSections[i][k]
154 + fraction*(pNCrossSections[i][k+1] - pNCrossSections[i][k]);
155 sigma.push_back(sigint);
156 }
157
158 G4int channel = sampleFlat(sigma);
159
160 std::vector<G4int> kinds;
161
162 if (mult == 2) {
163 for(i = 0; i < mult; i++) kinds.push_back(T0_2bfs[channel][i]);
164 } else if (mult == 3) {
165 for(i = 0; i < mult; i++) kinds.push_back(T0_3bfs[channel][i]);
166 } else if (mult == 4) {
167 for(i = 0; i < mult; i++) kinds.push_back(T0_4bfs[channel][i]);
168 } else if (mult == 5) {
169 for(i = 0; i < mult; i++) kinds.push_back(T0_5bfs[channel][i]);
170 } else if (mult == 6) {
171 for(i = 0; i < mult; i++) kinds.push_back(T0_6bfs[channel][i]);
172 } else if (mult == 7) {
173 for(i = 0; i < mult; i++) kinds.push_back(T0_7bfs[channel][i]);
174 } else if (mult == 8) {
175 for(i = 0; i < mult; i++) kinds.push_back(T0_8bfs[channel][i]);
176 } else if (mult == 9) {
177 for(i = 0; i < mult; i++) kinds.push_back(T0_9bfs[channel][i]);
178 } else {
179 G4cout << " Illegal multiplicity " << G4endl;
180 }
181
182 return kinds;
183}
184
185
186std::vector<G4int>
187G4RPGNucleonInelastic::GetFSPartTypesForT1(G4int mult, G4double KE,
188 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 = pPindex[mult-2][0];
199 G4int stop = pPindex[mult-2][1];
200
201 for(i = start; i < stop; i++) {
202 sigint = pPCrossSections[i][k]
203 + fraction*(pPCrossSections[i][k+1] - pPCrossSections[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(T1_2bfs[tzindex][channel][i]);
213 } else if (mult == 3) {
214 for(i = 0; i < mult; i++) kinds.push_back(T1_3bfs[tzindex][channel][i]);
215 } else if (mult == 4) {
216 for(i = 0; i < mult; i++) kinds.push_back(T1_4bfs[tzindex][channel][i]);
217 } else if (mult == 5) {
218 for(i = 0; i < mult; i++) kinds.push_back(T1_5bfs[tzindex][channel][i]);
219 } else if (mult == 6) {
220 for(i = 0; i < mult; i++) kinds.push_back(T1_6bfs[tzindex][channel][i]);
221 } else if (mult == 7) {
222 for(i = 0; i < mult; i++) kinds.push_back(T1_7bfs[tzindex][channel][i]);
223 } else if (mult == 8) {
224 for(i = 0; i < mult; i++) kinds.push_back(T1_8bfs[tzindex][channel][i]);
225 } else if (mult == 9) {
226 for(i = 0; i < mult; i++) kinds.push_back(T1_9bfs[tzindex][channel][i]);
227 } else {
228 G4cout << " Illegal multiplicity " << G4endl;
229 }
230
231 return kinds;
232}
233
234///////////////////////////////////////////////////////////////////////////////
235// //
236// p p and n n (|Tz| = 1) cross sections //
237// and final state particle types //
238// //
239///////////////////////////////////////////////////////////////////////////////
240
241// Total p p cross section as a function of kinetic energy
242G4double G4RPGNucleonInelastic::pPtot[30];
243
244// p p multiplicities as a function of kinetic energy
245G4double G4RPGNucleonInelastic::t1_dSigma_dMult[8][30];
246
247const G4int G4RPGNucleonInelastic::pPindex[8][2] =
248 {{0, 0}, {1, 6}, {7,24}, {25,56}, {57,63}, {64,71}, {72,81}, {82,92}};
249
250// Outgoing particle types of a given multiplicity
251// T1_nbfs = final state types for p p and n n
252
253const G4int G4RPGNucleonInelastic::T1_2bfs[2][1][2] =
254 {{{pro,pro}},
255
256 {{neu,neu}}};
257
258const G4int G4RPGNucleonInelastic::T1_3bfs[2][6][3] =
259 {{{pro,pro,pi0}, {pro,neu,pip}, {pro,lam,kp},
260 {pro,s0,kp}, {pro,sp,k0}, {neu,sp,kp}},
261
262 {{neu,neu,pi0}, {pro,neu,pim}, {neu,lam,k0},
263 {neu,s0,k0}, {neu,sm,kp}, {pro,sm,k0}}};
264
265const G4int G4RPGNucleonInelastic::T1_4bfs[2][18][4] =
266 {{{pro,pro,pip,pim},{pro,neu,pip,pi0},{pro,pro,pi0,pi0},
267 {neu,neu,pip,pip},{pro,lam,kp,pi0}, {pro,lam,k0,pip},
268 {neu,lam,kp,pip}, {neu,s0,kp,pip}, {pro,s0,kp,pi0},
269 {pro,s0,k0,pip}, {pro,sm,kp,pip}, {pro,sp,k0,pi0},
270 {neu,sp,k0,pip}, {pro,sp,kp,pim}, {neu,sp,kp,pi0},
271 {pro,pro,k0,k0b}, {pro,pro,kp,km}, {pro,neu,kp,k0b}},
272
273 {{neu,neu,pip,pim},{pro,neu,pim,pi0},{neu,neu,pi0,pi0},
274 {pro,pro,pim,pim},{neu,lam,k0,pi0}, {neu,lam,kp,pim},
275 {pro,lam,k0,pim}, {pro,s0,k0,pim}, {neu,s0,k0,pi0},
276 {neu,s0,kp,pim}, {neu,sp,k0,pim}, {neu,sm,kp,pi0},
277 {pro,sm,kp,pim}, {neu,sm,k0,pip}, {pro,sm,k0,pi0},
278 {neu,neu,kp,km}, {neu,neu,k0,k0b}, {pro,neu,k0,km}}};
279
280const G4int G4RPGNucleonInelastic::T1_5bfs[2][32][5] =
281 {{{pro,pro,pip,pim,pi0},{pro,pro,pi0,pi0,pi0},{pro,neu,pip,pip,pim},
282 {pro,neu,pip,pi0,pi0},{neu,neu,pip,pip,pi0},{pro,lam,kp,pip,pim},
283 {pro,lam,kp,pi0,pi0}, {pro,lam,k0,pip,pi0}, {pro,s0,kp,pip,pim},
284 {pro,s0,kp,pi0,pi0}, {pro,s0,k0,pip,pi0}, {pro,sp,k0,pip,pim},
285 {pro,sp,k0,pi0,pi0}, {pro,sp,kp,pim,pi0}, {pro,sm,kp,pip,pi0},
286 {pro,sm,k0,pip,pip}, {neu,lam,kp,pip,pi0}, {neu,lam,k0,pip,pip},
287 {neu,s0,kp,pip,pi0}, {neu,s0,k0,pip,pip}, {neu,sp,k0,pip,pi0},
288 {neu,sp,kp,pip,pim}, {neu,sp,kp,pi0,pi0}, {neu,sm,kp,pip,pip},
289 {pro,pro,pip,k0,km}, {pro,pro,pim,kp,k0b}, {pro,pro,pi0,k0,k0b},
290 {pro,pro,pi0,kp,km}, {pro,neu,pip,k0,k0b}, {pro,neu,pip,kp,km},
291 {pro,neu,pi0,kp,k0b}, {neu,neu,pip,kp,k0b}},
292
293 {{neu,neu,pip,pim,pi0},{neu,neu,pi0,pi0,pi0},{pro,neu,pip,pim,pim},
294 {pro,neu,pim,pi0,pi0},{pro,pro,pim,pim,pi0},{neu,lam,k0,pip,pim},
295 {neu,lam,k0,pi0,pi0}, {neu,lam,kp,pim,pi0}, {neu,s0,k0,pip,pim},
296 {neu,s0,k0,pi0,pi0}, {neu,s0,kp,pim,pi0}, {neu,sm,kp,pip,pim},
297 {neu,sm,kp,pi0,pi0}, {neu,sm,k0,pip,pi0}, {neu,sp,k0,pim,pi0},
298 {neu,sp,kp,pim,pim}, {pro,lam,k0,pim,pi0}, {pro,lam,kp,pim,pim},
299 {pro,s0,k0,pim,pi0}, {pro,s0,kp,pim,pim}, {pro,sm,kp,pim,pi0},
300 {pro,sm,k0,pip,pim}, {pro,sm,k0,pi0,pi0}, {pro,sp,k0,pim,pim},
301 {neu,neu,pim,kp,k0b}, {neu,neu,pip,k0,km}, {neu,neu,pi0,kp,km},
302 {neu,neu,pi0,k0,k0b}, {pro,neu,pim,kp,km}, {pro,neu,pim,k0,k0b},
303 {pro,neu,pi0,k0,km}, {pro,pro,pim,k0,km}}};
304
305const G4int G4RPGNucleonInelastic::T1_6bfs[2][7][6] =
306 {{{pro,pro,pip,pip,pim,pim},{pro,pro,pip,pim,pi0,pi0},
307 {pro,pro,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pi0},
308 {pro,neu,pip,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim},
309 {neu,neu,pip,pip,pi0,pi0}},
310
311 {{neu,neu,pip,pip,pim,pim},{neu,neu,pip,pim,pi0,pi0},
312 {neu,neu,pi0,pi0,pi0,pi0},{pro,neu,pip,pim,pim,pi0},
313 {pro,neu,pim,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pim},
314 {pro,pro,pim,pim,pi0,pi0}}};
315
316const G4int G4RPGNucleonInelastic::T1_7bfs[2][8][7] =
317 {{{pro,pro,pip,pip,pim,pim,pi0},{pro,pro,pip,pim,pi0,pi0,pi0},
318 {pro,pro,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pip,pim,pim},
319 {pro,neu,pip,pip,pim,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0},
320 {neu,neu,pip,pip,pip,pim,pi0},{neu,neu,pip,pip,pi0,pi0,pi0}},
321
322 {{neu,neu,pip,pip,pim,pim,pi0},{neu,neu,pip,pim,pi0,pi0,pi0},
323 {neu,neu,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pim,pim},
324 {pro,neu,pip,pim,pim,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0},
325 {pro,pro,pip,pim,pim,pim,pi0},{pro,pro,pim,pim,pi0,pi0,pi0}}};
326
327const G4int G4RPGNucleonInelastic::T1_8bfs[2][10][8] =
328 {{{pro,pro,pip,pip,pip,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pi0,pi0},
329 {pro,pro,pip,pim,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0},
330 {pro,neu,pip,pip,pip,pim,pim,pi0},{pro,neu,pip,pip,pim,pi0,pi0,pi0},
331 {pro,neu,pip,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pip,pim,pim},
332 {neu,neu,pip,pip,pip,pim,pi0,pi0},{neu,neu,pip,pip,pi0,pi0,pi0,pi0}},
333
334 {{neu,neu,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pim,pim,pi0,pi0},
335 {neu,neu,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0},
336 {pro,neu,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pim,pim,pi0,pi0,pi0},
337 {pro,neu,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim,pim},
338 {pro,pro,pip,pim,pim,pim,pi0,pi0},{pro,pro,pim,pim,pi0,pi0,pi0,pi0}}};
339
340const G4int G4RPGNucleonInelastic::T1_9bfs[2][11][9] =
341{{{pro,pro,pip,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pip,pim,pim,pi0,pi0,pi0},
342 {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
343 {pro,neu,pip,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pip,pim,pim,pi0,pi0},
344 {pro,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
345 {neu,neu,pip,pip,pip,pip,pim,pim,pi0},{neu,neu,pip,pip,pip,pim,pi0,pi0,pi0},
346 {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}},
347
348 {{neu,neu,pip,pip,pip,pim,pim,pim,pi0},{neu,neu,pip,pip,pim,pim,pi0,pi0,pi0},
349 {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
350 {pro,neu,pip,pip,pip,pim,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pim,pi0,pi0},
351 {pro,neu,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
352 {pro,pro,pip,pip,pim,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pim,pi0,pi0,pi0},
353 {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}}};
354
355//
356// Cross sections (in mb) for p p -> 2-9 body final states
357//
358// first index: 0: channels for mult = 2
359// 1-6: channels for mult = 3
360// 7-24: channels for mult = 4
361// 25-56: channels for mult = 5
362// 57-63: channels for mult = 6
363// 64-71: channels for mult = 7
364// 72-81: channels for mult = 8
365// 82-92: channels for mult = 9
366//
367// second index: kinetic energy
368//
369
370const G4float G4RPGNucleonInelastic::pPCrossSections[93][30] = {
371//
372// multiplicity 2 (1 channel)
373//
374// p p (n n)
375 { 0.0,330.0,240.0,160.0,110.0, 85.0, 63.0, 44.0, 33.0, 28.0,
376 25.0, 24.0, 23.0, 23.0, 26.3, 26.1, 25.0, 23.5, 21.0, 18.0,
377 16.0, 14.3, 12.5, 11.2, 10.3, 9.6, 9.0, 8.5, 8.0, 7.7 },
378//
379// multiplicity 3 (6 channels)
380//
381// p p pi0 (n n pi0)
382 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
383 0.0, 0.0, 0.0, 0.0, 0.0, 1.4, 4.0, 4.3, 4.0, 4.0,
384 3.6, 3.0, 2.8, 2.5, 1.7, 1.3, 1.1, 1.0, 0.9, 0.85 },
385
386// p n pi+ (p n pi-)
387 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
388 0.0, 0.0, 0.0, 0.0, 0.7, 4.5, 15.0, 19.1, 18.0, 16.0,
389 13.0, 10.0, 8.2, 6.0, 4.3, 3.3, 2.6, 2.0, 1.65, 1.4 },
390
391// p L K+ (n L K0)
392 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
393 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.012,
394 0.03, 0.06, 0.06, 0.055,0.05, 0.047,0.043,0.04, 0.037,0.033 },
395
396// p S0 K+ (n S0 K0)
397 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
398 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
399 0.006,0.02, 0.027,0.026,0.021,0.018,0.015,0.011,0.009,0.007 },
400
401// p S+ K0 (n S- K+)
402 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
403 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
404 0.013,0.025,0.03, 0.029,0.027,0.026,0.024,0.022,0.021,0.019 },
405
406// n S+ K+ (p S- K0)
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.0, 0.0, 0.0, 0.0, 0.0,
409 0.015,0.06, 0.07, 0.065,0.05, 0.04, 0.033,0.026,0.02, 0.015 },
410//
411// multiplicity 4 (18 channels)
412//
413// p p pi+ pi- (n n pi+ pi-)
414 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
415 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.6, 1.9,
416 2.8, 3.0, 3.0, 2.8, 2.5, 2.1, 1.9, 1.6, 1.4, 1.2 },
417
418// p n pi+ pi0 (p n pi- pi0)
419 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
420 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.6, 3.5,
421 4.0, 3.9, 3.5, 3.1, 2.8, 2.4, 2.2, 1.9, 1.7, 1.5 },
422
423// p p pi0 pi0 (n n pi0 pi0)
424 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
425 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.24, 0.76,
426 1.1, 1.2, 1.2, 1.1, 1.0, 0.84, 0.76, 0.64, 0.56, 0.48 },
427
428// n n pi+ pi+ (p p pi- pi-)
429 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
430 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.24, 1.4,
431 1.6, 1.6, 1.4, 1.2, 1.1, 1.0, 0.88, 0.76, 0.68, 0.6 },
432
433// L K+ p pi0 (L K0 n pi0)
434 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
435 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
436 0.002,0.015,0.05, 0.06, 0.052,0.042,0.037,0.029,0.025,0.020 },
437
438// L K0 p pi+ (L K+ n pi-)
439 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
440 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
441 0.002,0.015,0.06, 0.086,0.09, 0.082,0.072,0.06, 0.051,0.043 },
442
443// L K+ n pi+ (L K0 p pi-)
444 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
445 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
446 0.002,0.01, 0.04, 0.08, 0.066,0.058,0.05, 0.04, 0.035,0.03 },
447
448// S0 K+ n pi+ (S0 K0 p pi-)
449 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
450 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
451 0.0, 0.003,0.02, 0.03, 0.03, 0.025,0.02, 0.015,0.011,0.01 },
452
453// S0 K+ p pi0 (S0 K0 n pi0)
454 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
455 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
456 0.0, 0.005,0.02, 0.025,0.022,0.02, 0.015,0.01, 0.008,0.007 },
457
458// S0 K0 p pi+ (S0 K+ n pi-)
459 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
460 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
461 0.0, 0.003,0.012,0.04, 0.037,0.03, 0.027,0.022,0.019,0.016 },
462
463// S- K+ p pi+ (S+ K0 n pi-)
464 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
465 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
466 0.0, 0.004,0.016,0.037,0.031,0.028,0.023,0.02, 0.017,0.014 },
467
468// S+ K0 p pi0 (S- K+ n pi0)
469 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
470 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
471 0.0, 0.003,0.007,0.021,0.021,0.019,0.017,0.014,0.012,0.01 },
472
473// S+ K0 n pi+ (S- K+ p pi-)
474 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
475 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
476 0.0, 0.003,0.009,0.036,0.052,0.043,0.038,0.03, 0.026,0.02 },
477
478// S+ K+ p pi- (S- K0 n pi+)
479 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
480 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
481 0.0, 0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
482
483// S+ K+ n pi0 (S- K0 p 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.0, 0.0, 0.0, 0.0,
486 0.0, 0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
487
488// p p K0 K0bar (n n K+ K-)
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.0, 0.0, 0.0, 0.0,
491 0.0, 0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
492
493// p p K+ K- (n n K0 K0bar)
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.0, 0.0, 0.0, 0.0,
496 0.0, 0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
497
498// p n K+ K0bar (p n K0 K-)
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.0, 0.0, 0.0, 0.0,
501 0.0, 0.002,0.007,0.029,0.024,0.02, 0.017,0.014,0.012,0.009 },
502//
503// multiplicity 5 (32 channels)
504//
505// p p pi+ pi- pi0 (n n pi+ pi- pi0)
506 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
507 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.06,
508 0.4, 1.1, 1.8, 2.4, 2.4, 2.2, 2.0, 1.7, 1.5, 1.3 },
509
510// p p pi0 pi0 pi0 (n n pi0 pi0 pi0)
511 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
512 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.003,0.018,
513 0.12, 0.33, 0.54, 0.72, 0.72, 0.66, 0.6, 0.51, 0.45, 0.39 },
514
515// p n pi+ pi+ pi- (p n pi+ pi- pi-)
516 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
517 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12, 0.26,
518 0.7, 1.6, 2.4, 2.6, 2.3, 2.0, 1.8, 1.6, 1.4, 1.2 },
519
520// p n pi+ pi0 pi0 (p n pi- pi0 pi0)
521 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
522 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.036,0.078,
523 0.21, 0.48, 0.72, 0.78, 0.69, 0.6, 0.54, 0.48, 0.42, 0.36 },
524
525// n n pi+ pi+ pi0 (p p pi- pi- pi0)
526 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
527 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.036,
528 0.24, 0.66, 1.08, 1.44, 1.44, 1.32, 1.2, 1.0, 0.9, 0.78 },
529
530// p L K+ pi+ pi- (n L K0 pi+ pi-)
531 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
532 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
533 0.0, 0.0, 0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
534
535// p L K+ pi0 pi0 (n L K0 pi0 pi0)
536 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
537 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
538 0.0, 0.0, 0.001,0.005,0.020,0.022,0.020,0.017,0.015,0.010 },
539
540// p L K0 pi+ pi0 (n L K+ pi- pi0)
541 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
542 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
543 0.0, 0.0, 0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
544
545// p S0 K+ pi+ pi- (n S0 K0 pi+ pi-)
546 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
547 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
548 0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
549
550// p S0 K+ pi0 pi0 (n S0 K0 pi0 pi0)
551 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
552 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
553 0.0, 0.0, 0.001,0.004,0.015,0.017,0.015,0.014,0.011,0.009 },
554
555// p S0 K0 pi+ pi0 (n S0 K+ pi- pi0)
556 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
557 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
558 0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
559
560// p S+ K0 pi+ pi- (n S- K+ pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
563 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
564
565// p S+ K0 pi0 pi0 (n S- K+ pi0 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
568 0.0, 0.0, 0.001,0.002,0.009,0.020,0.013,0.014,0.011,0.009 },
569
570// p S+ K+ pi- pi0 (n S- K0 pi+ pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
573 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
574
575// p S- K+ pi+ pi0 (n S+ K0 pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
578 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
579
580// p S- K0 pi+ pi+ (n S+ K+ pi- pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
583 0.0, 0.0, 0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
584
585// n L K+ pi+ pi0 (p L K0 pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
588 0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
589
590// n L K0 pi+ pi+ (p L K+ pi- pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
593 0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
594
595// n S0 K+ pi+ pi0 (p S0 K0 pi- 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
598 0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
599
600// n S0 K0 pi+ pi+ (p S0 K+ pi- pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
603 0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
604
605// n S+ K0 pi+ pi0 (p S- K+ pi- pi0)
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.0, 0.0,
608 0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
609
610// n S+ K+ pi+ pi- (p S- K0 pi+ pi-)
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.0, 0.0,
613 0.0, 0.0, 0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
614
615// n S+ K+ pi0 pi0 (p S- K0 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.0, 0.0,
618 0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
619
620// n S- K+ pi+ pi+ (p S+ K0 pi- 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
623 0.0, 0.0, 0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
624
625// p p pi+ K0 K- (n n pi- K+ K0bar)
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.0, 0.0,
628 0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
629
630// p p pi- K+ K0bar (n n pi+ K0 K-)
631 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
632 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
633 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
634
635// p p pi0 K0 K0bar (n n pi0 K+ K-)
636 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
637 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
638 0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
639
640// p p pi0 K+ K- (n n pi0 K0 K0bar)
641 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
642 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
643 0.0, 0.0, 0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
644
645// p n pi+ K0 K0bar (p n pi- K+ K-)
646 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
647 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
648 0.0, 0.0, 0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
649
650// p n pi+ K+ K- (p n pi- K0 K0bar)
651 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
652 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
653 0.0, 0.0, 0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
654
655// p n pi0 K+ K0bar (p n pi0 K0 K-)
656 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
657 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
658 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
659
660// n n pi+ K+ K0bar (p p pi- K0 K-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
663 0.0, 0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
664//
665// multiplicity 6 (7 channels)
666//
667// p p pi+ pi+ pi- pi- (n n pi+ pi+ pi- pi-)
668 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
669 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
670 0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
671
672// p p pi+ pi- pi0 pi0 (n n pi+ pi- pi0 pi0)
673 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
674 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
675 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
676
677// p p pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0)
678 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
679 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
680 0.01, 0.02, 0.05, 0.1, 0.13, 0.12, 0.11, 0.1, 0.1, 0.09 },
681
682// p n pi+ pi+ pi- pi0 (p n pi+ pi- pi- pi0)
683 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
684 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
685 0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
686
687// p n pi+ pi0 pi0 pi0 (p n pi- pi0 pi0 pi0)
688 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
689 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
690 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
691
692// n n pi+ pi+ pi+ pi- (p p pi+ pi- pi- pi-)
693 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
694 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
695 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
696
697// n n pi+ pi+ pi0 pi0 (p p pi- pi- pi0 pi0)
698 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
699 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
700 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
701//
702// multiplicity 7 (8 channels)
703//
704// p p pi+ pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi- pi0)
705 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
706 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
707 0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
708
709// p p pi+ pi- pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0)
710 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
711 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
712 0.0, 0.0, 0.036,0.096,0.30, 0.42, 0.42, 0.42, 0.40, 0.37 },
713
714// p p pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0)
715 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
716 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
717 0.0, 0.0, 0.018,0.048,0.14, 0.20, 0.22, 0.20, 0.19, 0.18 },
718
719// p n pi+ pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi- pi-)
720 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
721 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
722 0.0, 0.0, 0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
723
724// p n pi+ pi+ pi- pi0 pi0 (p n pi+ pi- pi- pi0 pi0)
725 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
726 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
727 0.0, 0.0, 0.036,0.12, 0.18, 0.24, 0.26, 0.23, 0.28, 0.26 },
728
729// p n pi+ pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0)
730 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
731 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
732 0.0, 0.0, 0.024,0.06, 0.08, 0.12, 0.13, 0.14, 0.13, 0.13 },
733
734// n n pi+ pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi- pi0)
735 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
736 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
737 0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
738
739// n n pi+ pi+ pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0)
740 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
741 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
742 0.0, 0.0, 0.036,0.096,0.30, 0.42, 0.42, 0.41, 0.40, 0.37 },
743//
744// multiplicity 8 (10 channels)
745//
746// p p pi+ pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi- pi-)
747 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
748 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
749 0.0, 0.003,0.01, 0.024,0.075,0.18,0.27, 0.30, 0.27, 0.24 },
750
751// p p pi+ pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0)
752 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
753 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
754 0.0, 0.003,0.01, 0.024,0.075,0.18,0.27, 0.30, 0.27, 0.24 },
755
756// p p pi+ pi- pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0)
757 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
758 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
759 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
760
761// p p pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0)
762 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
763 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
764 0.0, 0.0, 0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
765
766// p n pi+ pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi- pi0)
767 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
768 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
769 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
770
771// p n pi+ pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0)
772 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
773 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
774 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
775
776// p n pi+ pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0)
777 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
778 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
779 0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
780
781// n n pi+ pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi- pi-)
782 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
783 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
784 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
785
786// n n pi+ pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0)
787 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
788 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
789 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
790
791// n n pi+ pi+ pi0 pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0 pi0)
792 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
793 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
794 0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
795//
796// multiplicity 9 (11 channels)
797//
798// p p pi+ pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi- pi0)
799 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
800 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
801 0.0, 0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
802
803// p p pi+ pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0 pi0)
804 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
805 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
806 0.0, 0.001,0.006,0.020,0.058,0.091,0.11, 0.12, 0.12, 0.12 },
807
808// p p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0 pi0)
809 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
810 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
811 0.0, 0.001,0.004,0.012,0.035,0.055,0.065,0.07, 0.07, 0.07 },
812
813// p p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
814 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
815 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
816 0.0, 0.0, 0.0, 0.006,0.027,0.032,0.04, 0.042,0.042,0.042 },
817
818// p n pi+ pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi- pi-)
819 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
820 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
821 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
822
823// p n pi+ pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi- pi0 pi0)
824 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
825 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
826 0.0, 0.003,0.008,0.026,0.078,0.20, 0.25, 0.29, 0.29, 0.29 },
827
828// p n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0 pi0)
829 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
830 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
831 0.0, 0.001,0.005,0.016,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
832
833// p n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0 pi0)
834 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
835 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
836 0.0, 0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 },
837
838// n n pi+ pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi- pi0)
839 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
840 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
841 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
842
843// n n pi+ pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0 pi0)
844 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
845 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
846 0.0, 0.001,0.005,0.015,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
847
848// n n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p p pi+ pi- pi0 pi0 pi0 pi0 pi0)
849 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
850 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
851 0.0, 0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 }};
852
853///////////////////////////////////////////////////////////////////////////////
854// //
855// p n and n p (|Tz| = 0) cross sections //
856// and final state particle types //
857// //
858///////////////////////////////////////////////////////////////////////////////
859
860// Total p n cross section as a function of kinetic energy
861G4double G4RPGNucleonInelastic::pNtot[30];
862
863// p n multiplicities as a function of kinetic energy
864G4double G4RPGNucleonInelastic::t0_dSigma_dMult[8][30];
865
866const G4int G4RPGNucleonInelastic::pNindex[8][2] =
867 {{0, 0}, {1,9}, {10,31}, {32,69}, {70,76}, {77,85}, {86,95}, {96,107}};
868// first index: 0: channels for mult = 2
869// 1-9: channels for mult = 3
870// 10-31: channels for mult = 4
871// 32-69: channels for mult = 5
872// 70-76: channels for mult = 6
873// 77-85: channels for mult = 7
874// 86-95: channels for mult = 8
875// 96-107: channels for mult = 9
876
877// Outgoing particle types of a given multiplicity
878// T0_nbfs = final state types for p n and n p
879
880const G4int G4RPGNucleonInelastic::T0_2bfs[1][2] =
881 {{pro,neu}};
882
883const G4int G4RPGNucleonInelastic::T0_3bfs[9][3] =
884 {{pro,pro,pim},{pro,neu,pi0},{neu,neu,pip},{pro,lam,k0},
885 {pro,s0,k0}, {pro,sm,kp}, {neu,lam,kp}, {neu,s0,kp},
886 {neu,sp,k0}};
887
888const G4int G4RPGNucleonInelastic::T0_4bfs[22][4] =
889 {{pro,neu,pip,pim},{pro,pro,pim,pi0},{pro,neu,pi0,pi0},
890 {neu,neu,pip,pi0},{pro,lam,kp,pim}, {pro,s0,kp,pim},
891 {pro,lam,k0,pi0}, {pro,s0,k0,pi0}, {pro,sp,k0,pim},
892 {pro,sm,kp,pi0}, {pro,sm,k0,pip}, {neu,lam,kp,pi0},
893 {neu,lam,k0,pip}, {neu,sp,kp,pim}, {neu,sp,k0,pi0},
894 {neu,s0,kp,pi0}, {neu,s0,k0,pip}, {neu,sm,kp,pip},
895 {pro,neu,kp,km}, {pro,neu,k0,k0b}, {pro,pro,k0,km},
896 {neu,neu,kp,k0b}};
897
898const G4int G4RPGNucleonInelastic::T0_5bfs[38][5] =
899 {{pro,neu,pip,pim,pi0},{pro,neu,pi0,pi0,pi0},{pro,pro,pip,pim,pim},
900 {pro,pro,pim,pi0,pi0},{neu,neu,pip,pip,pim},{neu,neu,pip,pi0,pi0},
901 {pro,lam,kp,pim,pi0}, {pro,lam,k0,pip,pim}, {pro,lam,k0,pi0,pi0},
902 {pro,s0,k0,pip,pim}, {pro,s0,k0,pi0,pi0}, {pro,s0,kp,pim,pi0},
903 {pro,sp,kp,pim,pim}, {pro,sp,k0,pim,pi0}, {pro,sm,k0,pip,pi0},
904 {pro,sm,kp,pip,pim}, {pro,sm,kp,pi0,pi0}, {neu,lam,kp,pip,pim},
905 {neu,lam,kp,pi0,pi0}, {neu,lam,k0,pip,pi0}, {neu,s0,kp,pip,pim},
906 {neu,s0,kp,pi0,pi0}, {neu,s0,k0,pip,pi0}, {neu,sp,k0,pip,pim},
907 {neu,sp,k0,pi0,pi0}, {neu,sp,kp,pim,pi0}, {neu,sm,kp,pip,pi0},
908 {neu,sm,k0,pip,pip}, {pro,neu,kp,km,pi0}, {pro,neu,k0,k0b,pi0},
909 {pro,neu,k0,km,pip}, {pro,neu,kp,k0b,pim}, {pro,pro,k0,k0b,pim},
910 {pro,pro,kp,km,pim}, {pro,pro,k0,km,pi0}, {neu,neu,kp,km,pip},
911 {neu,neu,k0,k0b,pip}, {neu,neu,kp,k0b,pi0}};
912
913const G4int G4RPGNucleonInelastic::T0_6bfs[7][6] =
914 {{pro,neu,pip,pip,pim,pim},{pro,neu,pip,pim,pi0,pi0},
915 {pro,neu,pi0,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pi0},
916 {pro,pro,pim,pi0,pi0,pi0},{neu,neu,pip,pip,pim,pi0},
917 {neu,neu,pip,pi0,pi0,pi0}};
918
919const G4int G4RPGNucleonInelastic::T0_7bfs[9][7] =
920 {{pro,neu,pip,pip,pim,pim,pi0},{pro,neu,pip,pim,pi0,pi0,pi0},
921 {pro,neu,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim},
922 {pro,pro,pip,pim,pim,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0},
923 {neu,neu,pip,pip,pip,pim,pim},{neu,neu,pip,pip,pim,pi0,pi0},
924 {neu,neu,pip,pi0,pi0,pi0,pi0}};
925
926const G4int G4RPGNucleonInelastic::T0_8bfs[10][8] =
927{{pro,neu,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pi0,pi0},
928 {pro,neu,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0},
929 {pro,pro,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pi0,pi0,pi0},
930 {pro,pro,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim,pim,pi0},
931 {neu,neu,pip,pip,pim,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0}};
932
933const G4int G4RPGNucleonInelastic::T0_9bfs[12][9] =
934{{pro,neu,pip,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pip,pim,pim,pi0,pi0,pi0},
935 {pro,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
936 {pro,pro,pip,pip,pip,pim,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pim,pi0,pi0},
937 {pro,pro,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
938 {neu,neu,pip,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pip,pim,pim,pi0,pi0},
939 {neu,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}};
940
941//
942// Cross sections (in mb) for p n -> 2-9 body final states
943//
944// first index: 0: channels for mult = 2
945// 1-9: channels for mult = 3
946// 10-31: channels for mult = 4
947// 32-69: channels for mult = 5
948// 70-76: channels for mult = 6
949// 77-85: channels for mult = 7
950// 86-95: channels for mult = 8
951// 96-107: channels for mult = 9
952//
953// second index: kinetic energy
954//
955const G4float G4RPGNucleonInelastic::pNCrossSections[108][30] = {
956//
957// multiplicity 2 (1 channel)
958//
959// p n (p n)
960 { 0.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0,
961 44.0, 42.0, 40.0, 35.0, 31.0, 27.0, 23.0, 19.0, 17.0, 15.5,
962 14.0, 13.0, 12.0, 11.0, 10.0, 9.5, 9.0, 8.5, 8.0, 7.7 },
963//
964// multiplicity 3 (9 channels)
965//
966// p p pi- (n n pi+)
967 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
968 0.0, 0.0, 0.0, 0.0, 0.25, 0.9, 1.75, 2.3, 2.8, 2.8,
969 2.2, 1.9, 1.6, 1.35, 1.1, 0.95, 0.8, 0.7, 0.6, 0.53 },
970
971// p n pi0 (p n pi0)
972 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
973 0.0, 0.0, 0.0, 0.0, 1.8, 4.7, 8.3, 11.3, 12.0, 10.2,
974 8.2, 6.0, 4.9, 3.6, 2.5, 2.0, 1.6, 1.2, 1.0, 0.08 },
975
976// n n pi+ (p p pi-)
977 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
978 0.0, 0.0, 0.0, 0.0, 0.95, 2.4, 4.2, 5.6, 6.1, 5.1,
979 4.1, 3.0, 2.5, 1.8, 1.2, 1.0, 0.8, 0.6, 0.5, 0.41},
980
981// p L K0 (n L K+)
982 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
983 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
984 0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
985
986// p S0 K0 (n S0 K+)
987 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
988 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
989 0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
990
991// p S- K+ (n S+ K0)
992 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
993 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
994 0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
995
996// n L K+ (p L K0)
997 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
998 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
999 0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
1000
1001// n S0 K+ (p S0 K0)
1002 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1003 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1004 0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
1005
1006// n S+ K0 (p S- K+)
1007 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1008 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1009 0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
1010//
1011// multiplicity 4 (22 channels)
1012//
1013// p n pi+ pi- (p n pi+ pi-)
1014 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1015 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12, 0.38, 1.1, 3.5,
1016 5.9, 5.9, 5.1, 4.2, 3.7, 3.0, 2.6, 2.1, 1.8, 1.4 },
1017
1018// p p pi- pi0 (n n pi+ pi0)
1019 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1020 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.1, 0.24, 0.55,
1021 1.2, 1.5, 1.45, 1.25, 1.0, 0.9, 0.8, 0.7, 0.6, 0.53 },
1022
1023// p n pi0 pi0 (p n pi0 pi0)
1024 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1025 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.07, 0.24, 0.66, 2.1,
1026 3.6, 3.6, 3.1, 2.5, 2.2, 1.8, 1.5, 1.2, 1.1, 0.84 },
1027
1028// n n pi+ pi0 (p p pi- pi0)
1029 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1030 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.1, 0.24, 0.55,
1031 1.2, 1.5, 1.45, 1.25, 1.0, 0.9, 0.8, 0.7, 0.6, 0.53 },
1032
1033// p L K+ pi- (n L K0 pi+)
1034 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1035 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1036 0.0, 0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
1037
1038// p S0 K+ pi- (n S0 K0 pi+)
1039 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1040 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1041 0.0, 0.0, 0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
1042
1043// p L K0 pi0 (n L K+ pi0)
1044 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1045 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1046 0.0, 0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
1047
1048// p S0 K0 pi0 (n S0 K+ pi0)
1049 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1050 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1051 0.0, 0.0, 0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
1052
1053// p S+ K0 pi- (n S- K+ pi+)
1054 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1055 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1056 0.0, 0.0015,0.004,0.013,0.02,0.016,0.013,0.01, 0.009,0.007 },
1057
1058// p S- K+ pi0 (n S+ K0 pi0)
1059 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1060 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1061 0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1062
1063// p S- K0 pi+ (n S+ K+ pi-)
1064 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1065 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1066 0.0, 0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1067
1068// n L K+ pi0 (p L K0 pi0)
1069 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1070 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1071 0.0, 0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
1072
1073// n L K0 pi+ (p L K+ pi-)
1074 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1075 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1076 0.0, 0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
1077
1078// n S+ K+ pi- (p S- K0 pi+)
1079 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1080 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1081 0.0, 0.0, 0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
1082
1083// n S+ K0 pi0 (p S- K+ pi0)
1084 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1085 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1086 0.0, 0.0, 0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
1087
1088// n S0 K+ pi0 (p S0 K0 pi0)
1089 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1092
1093// n S0 K0 pi+ (p S0 K+ pi-)
1094 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
1097
1098// n S- K+ pi+ (p S+ K0 pi-)
1099 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.0015,0.004,0.013,0.02, 0.016,0.013,0.01, 0.009,0.007 },
1102
1103// p n K+ K- (p n K0 K0bar)
1104 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1107
1108// p n K0 K0bar (p n K+ K-)
1109 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1112
1113// p p K0 K- (n n K+ K0bar)
1114 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1117
1118// n n K+ K0bar (p p K0 K-)
1119 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
1122//
1123// multiplicity 5 (38 channels)
1124//
1125// p n pi+ pi- pi0 (p n pi+ pi- pi0)
1126 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1127 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
1128 0.3, 0.82, 1.35, 1.8, 1.8, 1.65, 1.5, 1.28, 1.12, 0.98 },
1129
1130// p n pi0 pi0 pi0 (p n pi0 pi0 pi0)
1131 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1132 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.004,0.022,
1133 0.15, 0.41, 0.68, 0.9, 0.9, 0.82, 0.75, 0.64, 0.55, 0.49 },
1134
1135// p p pi+ pi- pi- (n n pi+ pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09, 0.2,
1138 0.52, 1.2, 1.8, 2.0, 1.7, 1.5, 1.35, 1.2, 1.05, 0.9 },
1139
1140// p p pi- pi0 pi0 (n n pi+ pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.1,
1143 0.26, 0.6, 0.9, 0.98, 0.86, 0.75, 0.68, 0.6, 0.52, 0.45 },
1144
1145// n n pi+ pi+ pi- (p p pi+ pi- pi-)
1146 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1147 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
1148 0.3, 0.82, 1.35, 1.8, 1.8, 1.65, 1.5, 1.28, 1.12, 0.98 },
1149
1150// n n pi+ pi0 pi0 (p p pi- pi0 pi0)
1151 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1152 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.003,0.02,
1153 0.15, 0.41, 0.68, 0.9, 0.9, 0.82, 0.75, 0.64, 0.56, 0.49 },
1154
1155// p L K+ pi- pi0 (n L K0 pi+ pi0)
1156 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1157 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1158 0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1159
1160// p L K0 pi+ pi- (n L K+ pi+ pi-)
1161 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1162 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1163 0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1164
1165// p L K0 pi0 pi0 (n L K+ pi0 pi0)
1166 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1167 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1168 0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.09, 0.07, 0.07 },
1169
1170// p S0 K0 pi+ pi- (n S0 K+ pi+ pi-)
1171 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1172 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1173 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1174
1175// p S0 K0 pi0 pi0 (n S0 K+ pi0 pi0)
1176 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1177 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1178 0.0, 0.0, 0.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1179
1180// p S0 K+ pi- pi0 (n S0 K0 pi+ pi0)
1181 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1182 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1183 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1184
1185// p S+ K+ pi- pi- (n S- K0 pi+ pi+)
1186 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1187 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1188 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1189
1190// p S+ K0 pi- pi0 (n S- K+ pi+ pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1193 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1194
1195// p S- K0 pi+ pi0 (n S+ K+ pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1198 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1199
1200// p S- K+ pi+ pi- (n S+ K0 pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1203 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1204
1205// p S- K+ pi0 pi0 (n S+ K0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1208 0.0, 0.0, 0.002,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1209
1210// n L K+ pi+ pi- (p L K0 pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1213 0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1214
1215// n L K+ pi0 pi0 (p L K0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1218 0.0, 0.0, 0.001,0.004,0.011,0.012,0.011,0.009,0.007,0.006 },
1219
1220// n L K0 pi+ pi0 (p L K+ pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1223 0.0, 0.0, 0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
1224
1225// n S0 K+ pi+ pi- (p S0 K0 pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1228 0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
1229
1230// n S0 K+ pi0 pi0 (p S0 K0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1233 0.0, 0.0, 0.0, 0.002,0.005,0.006,0.005,0.005,0.004,0.003 },
1234
1235// n S0 K0 pi+ pi0 (p S0 K+ pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1238 0.0, 0.0, 0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
1239
1240// n S+ K0 pi+ pi- (p S- K+ pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1243 0.0, 0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1244
1245// n S+ K0 pi0 pi0 (p S- K+ pi0 pi0)
1246 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
1249
1250// n S+ K+ pi- pi0 (p S- K0 pi+ pi0)
1251 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1254
1255// n S- K+ pi+ pi0 (p S+ K0 pi- pi0)
1256 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1259
1260// n S- K0 pi+ pi+ (p S+ K+ pi- pi-)
1261 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
1264
1265// p n K+ K- pi0 (p n K0 K0bar pi0)
1266 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1269
1270// p n K0 K0bar pi0 (p n K+ K- pi0)
1271 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1274
1275// p n K0 K- pi+ (p n K+ K0bar pi-)
1276 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1277 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1278 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1279
1280// p n K+ K0bar pi- (p n K0 K- pi+)
1281 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1282 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1283 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1284
1285// p p K0 K0bar pi- (n n K+ K- pi+)
1286 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1287 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1288 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1289
1290// p p K+ K- pi- (n n K0 K0bar pi+)
1291 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1292 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1293 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1294
1295// p p K0 K- pi0 (n n K+ K0bar pi0)
1296 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1297 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1298 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1299
1300// n n K+ K- pi+ (p p K0 K0bar pi-)
1301 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1302 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1303 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1304
1305// n n K0 K0bar pi+ (p p K+ K- pi-)
1306 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1307 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1308 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1309
1310// n n K+ K0bar pi0 (p p K0 K- pi0)
1311 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1312 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1313 0.0, 0.0, 0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
1314//
1315// multiplicity 6 (7 channels)
1316//
1317// p n pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi-)
1318 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1319 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1320 0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1321
1322// p n pi+ pi- pi0 pi0 (p n pi+ pi- pi0 pi0)
1323 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1324 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1325 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1326
1327// p n pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0)
1328 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1329 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1330 0.01, 0.02, 0.05, 0.1, 0.13, 0.12, 0.11, 0.1, 0.1, 0.09 },
1331
1332// p p pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi0)
1333 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1334 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1335 0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1336
1337// p p pi- pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0)
1338 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1339 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1340 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1341
1342// n n pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi0)
1343 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1344 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1345 0.06, 0.1, 0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
1346
1347// n n pi+ pi0 pi0 pi0 (p p pi- pi0 pi0 pi0)
1348 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1349 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1350 0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2, 0.19, 0.18 },
1351//
1352// multiplicity 7 (9 channels)
1353//
1354// p n pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi0)
1355 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1356 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1357 0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
1358
1359// p n pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0)
1360 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1361 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1362 0.0, 0.0, 0.03, 0.08, 0.25, 0.35, 0.35, 0.35, 0.33, 0.31 },
1363
1364// p n pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0)
1365 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1366 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1367 0.0, 0.0, 0.015,0.04, 0.12, 0.17, 0.18, 0.17, 0.16, 0.15 },
1368
1369// p p pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi-)
1370 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1371 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1372 0.0, 0.0, 0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
1373
1374// p p pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi0 pi0)
1375 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1376 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1377 0.0, 0.0, 0.03, 0.1, 0.15, 0.2, 0.22, 0.23, 0.22, 0.22 },
1378
1379// p p pi- pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0)
1380 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1381 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1382 0.0, 0.0, 0.02, 0.05, 0.07, 0.1, 0.11, 0.12, 0.11, 0.11 },
1383
1384// n n pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi-)
1385 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1386 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1387 0.0, 0.0, 0.06, 0.17, 0.5, 0.7, 0.7, 0.69, 0.66, 0.62 },
1388
1389// n n pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi0 pi0)
1390 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1391 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1392 0.0, 0.0, 0.03, 0.08, 0.25, 0.35, 0.35, 0.34, 0.33, 0.31 },
1393
1394// n n pi+ pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0)
1395 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1396 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1397 0.0, 0.0, 0.02, 0.05, 0.07, 0.1, 0.11, 0.12, 0.11, 0.11 },
1398//
1399// multiplicity 8 (10 channels)
1400//
1401// p n pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi-)
1402 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1403 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1404 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1405
1406// p n pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0)
1407 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1408 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1409 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1410
1411// p n pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0)
1412 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1413 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1414 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1415
1416// p n pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0)
1417 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1418 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1419 0.0, 0.0, 0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
1420
1421// p p pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi0)
1422 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1423 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1424 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1425
1426// p p pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0)
1427 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1428 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1429 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1430
1431// p p pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0)
1432 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1433 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1434 0.0, 0.0, 0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
1435
1436// n n pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi0)
1437 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1438 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1439 0.0, 0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
1440
1441// n n pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0)
1442 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1443 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1444 0.0, 0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
1445
1446// n n pi+ pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0)
1447 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1448 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1449 0.0, 0.0, 0.003,0.009,0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
1450//
1451// multiplicity 9 (12 channels)
1452//
1453// p n pi+ pi+ pi+ pi- pi- pi- pi0 (p n pi+ pi+ pi+ pi- pi- pi- pi0)
1454 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1455 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1456 0.0, 0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
1457
1458// p n pi+ pi+ pi- pi- pi0 pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0 pi0)
1459 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1460 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1461 0.0, 0.001,0.005,0.015,0.045,0.07, 0.084,0.09, 0.09, 0.09 },
1462
1463// p n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0 pi0)
1464 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1465 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1466 0.0, 0.001,0.003,0.009,0.027,0.042,0.05, 0.054,0.054,0.054 },
1467
1468// p n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
1469 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1470 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1471 0.0, 0.0, 0.0, 0.005,0.016,0.025,0.03, 0.032,0.032,0.032 },
1472
1473// p p pi+ pi+ pi+ pi- pi- pi- pi- (n n pi+ pi+ pi+ pi+ pi- pi- pi-)
1474 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1475 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1476 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1477
1478// p p pi+ pi+ pi- pi- pi- pi0 pi0 (n n pi+ pi+ pi+ pi- pi- pi0 pi0)
1479 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1480 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1481 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1482
1483// p p pi+ pi- pi- pi0 pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0 pi0)
1484 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1485 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1486 0.0, 0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
1487
1488// p p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
1489 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1490 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1491 0.0, 0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 },
1492
1493// n n pi+ pi+ pi+ pi+ pi- pi- pi- (p p pi+ pi+ pi+ pi- pi- pi- pi-)
1494 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1495 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1496 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1497
1498// n n pi+ pi+ pi+ pi- pi- pi0 pi0 (p p pi+ pi+ pi- pi- pi- pi0 pi0)
1499 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1500 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1501 0.0, 0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
1502
1503// n n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0 pi0)
1504 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1505 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1506 0.0, 0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
1507
1508// n n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0 pi0)
1509 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1510 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1511 0.0, 0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 }};
Note: See TracBrowser for help on using the repository browser.