Context Navigation

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

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

Download in other formats: