Context Navigation

source: trunk/source/processes/hadronic/models/cascade/cascade/src/G4PionSampler.cc @ 1228

Last change on this file since 1228 was 1228, checked in by garnier, 14 years ago
update geant4.9.3 tag
File size: 98.9 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: G4PionSampler.cc,v 1.5 2009/12/02 17:35:10 dennis Exp $
27	// GEANT4 tag $Name: geant4-09-03 $
28	//
29
30	#include "G4PionSampler.hh"
31	#include "Randomize.hh"
32
33	G4PionSampler::G4PionSampler()
34	:G4FinalStateSampler()
35	{
36	initCrossSections();
37
38	// Initialize t33_dSigma_dMult, t31_dSigma_dMult, t11_dSigma_dMult:
39	// pi+, pi- and pi0 - nucleon inelastic cross sections as a function
40	// of multiplicity
41	// \|T, Tz> = \|3/2,3/2> , \|3/2,-1/2> , \|1/2, 1/2> respectively
42
43	// First set up indeces to arrays
44	const G4int pipPChanNums[8] = {2, 7, 15, 24, 5, 6, 7, 8};
45	const G4int pimPChanNums[8] = {5, 13, 22, 31, 6, 7, 8, 9};
46	const G4int pizPChanNums[8] = {5, 13, 21, 30, 6, 7, 8, 9};
47	G4int pipPTotChans = -1;
48	G4int pimPTotChans = -1;
49	G4int pizPTotChans = -1;
50	for (G4int i = 0; i < 8; i++) {
51	pipPTotChans += pipPChanNums[i];
52	pimPTotChans += pimPChanNums[i];
53	pizPTotChans += pizPChanNums[i];
54	pipPindex[i][1] = pipPTotChans;
55	pipPindex[i][0] = pipPTotChans - pipPChanNums[i] + 1;
56	pimPindex[i][1] = pimPTotChans;
57	pimPindex[i][0] = pimPTotChans - pimPChanNums[i] + 1;
58	pizPindex[i][1] = pizPTotChans;
59	pizPindex[i][0] = pizPTotChans - pizPChanNums[i] + 1;
60	}
61
62	G4int j, k, m;
63	G4int start, stop;
64
65	for (m = 0; m < 8; m++) {
66	start = pipPindex[m][0];
67	stop = pipPindex[m][1] + 1;
68	for (k = 0; k < 30; k++) {
69	t33_dSigma_dMult[m][k] = 0.0;
70	for (j = start; j < stop; j++) t33_dSigma_dMult[m][k] += pipPCrossSections[j][k];
71	}
72
73	start = pimPindex[m][0];
74	stop = pimPindex[m][1] + 1;
75	for (k = 0; k < 30; k++) {
76	t31_dSigma_dMult[m][k] = 0.0;
77	for (j = start; j < stop; j++) t31_dSigma_dMult[m][k] += pimPCrossSections[j][k];
78	}
79
80	start = pizPindex[m][0];
81	stop = pizPindex[m][1] + 1;
82	for (k = 0; k < 30; k++) {
83	t11_dSigma_dMult[m][k] = 0.0;
84	for (j = start; j < stop; j++) t11_dSigma_dMult[m][k] += pizPCrossSections[j][k];
85	}
86	}
87
88	// Initialize total cross section array
89
90	for (k = 0; k < 30; k++) {
91	pipPsummed[k] = 0.0;
92	pimPsummed[k] = 0.0;
93	pizPsummed[k] = 0.0;
94	for (m = 0; m < 8; m++) {
95	pipPsummed[k] += t33_dSigma_dMult[m][k];
96	pimPsummed[k] += t31_dSigma_dMult[m][k];
97	pizPsummed[k] += t11_dSigma_dMult[m][k];
98	}
99	}
100
101	// printCrossSections();
102
103	}
104
105
106	void G4PionSampler::printCrossSections() const
107	{
108	G4cout << " pi+ p total cross sections (mb) " << G4endl;
109	for (G4int i = 0; i < 5; i++) {
110	G4int istart = i*6;
111	G4int istop = istart + 6;
112	for (G4int t = istart; t < istop; t++) G4cout << pipPtot[t] << " " ;
113	G4cout << G4endl;
114	}
115	G4cout << G4endl;
116
117	G4cout << " pi+ p summed partial cross sections (mb) " << G4endl;
118	for (G4int i = 0; i < 5; i++) {
119	G4int istart = i*6;
120	G4int istop = istart + 6;
121	for (G4int t = istart; t < istop; t++) G4cout << pipPsummed[t] << " " ;
122	G4cout << G4endl;
123	}
124	G4cout << G4endl;
125
126	G4cout << " pi- p total cross sections (mb) " << G4endl;
127	for (G4int i = 0; i < 5; i++) {
128	G4int istart = i*6;
129	G4int istop = istart + 6;
130	for (G4int t = istart; t < istop; t++) G4cout << pimPtot[t] << " " ;
131	G4cout << G4endl;
132	}
133	G4cout << G4endl;
134
135	G4cout << " pi- p summed partial cross sections (mb) " << G4endl;
136	for (G4int i = 0; i < 5; i++) {
137	G4int istart = i*6;
138	G4int istop = istart + 6;
139	for (G4int t = istart; t < istop; t++) G4cout << pimPsummed[t] << " " ;
140	G4cout << G4endl;
141	}
142	G4cout << G4endl;
143
144	G4cout << " pi0 p total cross sections (mb) " << G4endl;
145	for (G4int i = 0; i < 5; i++) {
146	G4int istart = i*6;
147	G4int istop = istart + 6;
148	for (G4int t = istart; t < istop; t++) G4cout << pizPtot[t] << " " ;
149	G4cout << G4endl;
150	}
151	G4cout << G4endl;
152
153	G4cout << " pi0 p summed partial cross sections (mb) " << G4endl;
154	for (G4int i = 0; i < 5; i++) {
155	G4int istart = i*6;
156	G4int istop = istart + 6;
157	for (G4int t = istart; t < istop; t++) G4cout << pizPsummed[t] << " " ;
158	G4cout << G4endl;
159	}
160	G4cout << G4endl;
161	}
162
163
164	G4int G4PionSampler::GetMultiplicityT31(G4double KE) const
165	{
166	G4double multint(0.);
167	std::vector<G4double> sigma;
168
169	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
170	G4int k = epair.first;
171	G4double fraction = epair.second;
172
173	// Compare summed partial cross section with total cross section
174	// Truncate multiplicity at 9 if summed < total
175
176	G4double summed = pimPsummed[k] + fraction*(pimPsummed[k+1] - pimPsummed[k]);
177	G4double total = pimPtot[k] + fraction*(pimPtot[k+1] - pimPtot[k]);
178
179	if (G4UniformRand() > summed/total) {
180	// G4cout << " T31: partial sum = " << summed << " , total = " << total << " : truncating to 9 " << G4endl;
181	return 9;
182	} else {
183	for(G4int m = 0; m < 8; m++) {
184	multint = t31_dSigma_dMult[m][k]
185	+ fraction*(t31_dSigma_dMult[m][k+1] - t31_dSigma_dMult[m][k]);
186	sigma.push_back(multint);
187	}
188	return sampleFlat(sigma) + 2;
189	}
190	}
191
192
193	G4int G4PionSampler::GetMultiplicityT33(G4double KE) const
194	{
195	G4double multint(0.);
196	std::vector<G4double> sigma;
197
198	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
199	G4int k = epair.first;
200	G4double fraction = epair.second;
201
202	// Compare summed partial cross section with total cross section
203	// Truncate multiplicity at 9 if summed < total
204
205	G4double summed = pipPsummed[k] + fraction*(pipPsummed[k+1] - pipPsummed[k]);
206	G4double total = pipPtot[k] + fraction*(pipPtot[k+1] - pipPtot[k]);
207
208	if (G4UniformRand() > summed/total) {
209	// G4cout << " T33: partial sum = " << summed << " , total = " << total << " : truncating to 9 " << G4endl;
210	return 9;
211	} else {
212	for(G4int m = 0; m < 8; m++) {
213	multint = t33_dSigma_dMult[m][k]
214	+ fraction*(t33_dSigma_dMult[m][k+1] - t33_dSigma_dMult[m][k]);
215	sigma.push_back(multint);
216	}
217	return sampleFlat(sigma) + 2;
218	}
219	}
220
221
222	G4int G4PionSampler::GetMultiplicityT11(G4double KE) const
223	{
224	G4double multint(0.);
225	std::vector<G4double> sigma;
226
227	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
228	G4int k = epair.first;
229	G4double fraction = epair.second;
230
231	// Compare summed partial cross section with total cross section
232	// Truncate multiplicity at 8 if summed < total
233
234	G4double summed = pizPsummed[k] + fraction*(pizPsummed[k+1] - pizPsummed[k]);
235	G4double total = pizPtot[k] + fraction*(pizPtot[k+1] - pizPtot[k]);
236
237	if (G4UniformRand() > summed/total) {
238	// G4cout << " T11: partial sum = " << summed << " , total = " << total << " : truncating to 8 " << G4endl;
239	return 9;
240	} else {
241	for(G4int m = 0; m < 8; m++) {
242	multint = t11_dSigma_dMult[m][k]
243	+ fraction*(t11_dSigma_dMult[m][k+1] - t11_dSigma_dMult[m][k]);
244	sigma.push_back(multint);
245	}
246	return sampleFlat(sigma) + 2;
247	}
248	}
249
250
251	std::vector<G4int>
252	G4PionSampler::GetFSPartTypesForT33(G4int mult, G4double KE, G4int tzindex) const
253	{
254	G4int i;
255	G4double sigint(0.);
256	std::vector<G4double> sigma;
257
258	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
259	G4int k = epair.first;
260	G4double fraction = epair.second;
261
262	G4int start = pipPindex[mult-2][0];
263	G4int stop = pipPindex[mult-2][1];
264
265	for(i = start; i < stop; i++) {
266	sigint = pipPCrossSections[i][k]
267	+ fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[i][k]);
268	sigma.push_back(sigint);
269	}
270
271	G4int channel = sampleFlat(sigma);
272
273	std::vector<G4int> kinds;
274
275	if (mult == 2) {
276	for(i = 0; i < mult; i++) kinds.push_back(T33_2bfs[tzindex][channel][i]);
277	} else if (mult == 3) {
278	for(i = 0; i < mult; i++) kinds.push_back(T33_3bfs[tzindex][channel][i]);
279	} else if (mult == 4) {
280	for(i = 0; i < mult; i++) kinds.push_back(T33_4bfs[tzindex][channel][i]);
281	} else if (mult == 5) {
282	for(i = 0; i < mult; i++) kinds.push_back(T33_5bfs[tzindex][channel][i]);
283	} else if (mult == 6) {
284	for(i = 0; i < mult; i++) kinds.push_back(T33_6bfs[tzindex][channel][i]);
285	} else if (mult == 7) {
286	for(i = 0; i < mult; i++) kinds.push_back(T33_7bfs[tzindex][channel][i]);
287	} else if (mult == 8) {
288	for(i = 0; i < mult; i++) kinds.push_back(T33_8bfs[tzindex][channel][i]);
289	} else if (mult == 9) {
290	for(i = 0; i < mult; i++) kinds.push_back(T33_9bfs[tzindex][channel][i]);
291	} else {
292	G4cout << " Illegal multiplicity " << G4endl;
293	}
294
295	return kinds;
296	}
297
298
299	std::vector<G4int>
300	G4PionSampler::GetFSPartTypesForT31(G4int mult, G4double KE, G4int tzindex) const
301	{
302	G4int i;
303	G4double sigint(0.);
304	std::vector<G4double> sigma;
305
306	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
307	G4int k = epair.first;
308	G4double fraction = epair.second;
309
310	G4int start = pimPindex[mult-2][0];
311	G4int stop = pimPindex[mult-2][1];
312
313	for(i = start; i < stop; i++) {
314	sigint = pimPCrossSections[i][k]
315	+ fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
316	sigma.push_back(sigint);
317	}
318
319	G4int channel = sampleFlat(sigma);
320
321	std::vector<G4int> kinds;
322
323	if (mult == 2) {
324	for(i = 0; i < mult; i++) kinds.push_back(T31_2bfs[tzindex][channel][i]);
325	} else if (mult == 3) {
326	for(i = 0; i < mult; i++) kinds.push_back(T31_3bfs[tzindex][channel][i]);
327	} else if (mult == 4) {
328	for(i = 0; i < mult; i++) kinds.push_back(T31_4bfs[tzindex][channel][i]);
329	} else if (mult == 5) {
330	for(i = 0; i < mult; i++) kinds.push_back(T31_5bfs[tzindex][channel][i]);
331	} else if (mult == 6) {
332	for(i = 0; i < mult; i++) kinds.push_back(T31_6bfs[tzindex][channel][i]);
333	} else if (mult == 7) {
334	for(i = 0; i < mult; i++) kinds.push_back(T31_7bfs[tzindex][channel][i]);
335	} else if (mult == 8) {
336	for(i = 0; i < mult; i++) kinds.push_back(T31_8bfs[tzindex][channel][i]);
337	} else if (mult == 9) {
338	for(i = 0; i < mult; i++) kinds.push_back(T31_9bfs[tzindex][channel][i]);
339	} else {
340	G4cout << " Illegal multiplicity " << G4endl;
341	}
342
343	return kinds;
344	}
345
346
347	std::vector<G4int>
348	G4PionSampler::GetFSPartTypesForT11(G4int mult, G4double KE, G4int tzindex) const
349	{
350	G4int i;
351	G4double sigint(0.);
352	std::vector<G4double> sigma;
353
354	std::pair<G4int, G4double> epair = interpolateEnergy(KE);
355	G4int k = epair.first;
356	G4double fraction = epair.second;
357
358	G4int start = pizPindex[mult-2][0];
359	G4int stop = pizPindex[mult-2][1];
360
361	for(i = start; i < stop; i++) {
362	sigint = pizPCrossSections[i][k]
363	+ fraction*(pizPCrossSections[i][k+1] - pizPCrossSections[i][k]);
364	sigma.push_back(sigint);
365	}
366
367	G4int channel = sampleFlat(sigma);
368
369	std::vector<G4int> kinds;
370
371	if (mult == 2) {
372	for(i = 0; i < mult; i++) kinds.push_back(T11_2bfs[tzindex][channel][i]);
373	} else if (mult == 3) {
374	for(i = 0; i < mult; i++) kinds.push_back(T11_3bfs[tzindex][channel][i]);
375	} else if (mult == 4) {
376	for(i = 0; i < mult; i++) kinds.push_back(T11_4bfs[tzindex][channel][i]);
377	} else if (mult == 5) {
378	for(i = 0; i < mult; i++) kinds.push_back(T11_5bfs[tzindex][channel][i]);
379	} else if (mult == 6) {
380	for(i = 0; i < mult; i++) kinds.push_back(T11_6bfs[tzindex][channel][i]);
381	} else if (mult == 7) {
382	for(i = 0; i < mult; i++) kinds.push_back(T11_7bfs[tzindex][channel][i]);
383	} else if (mult == 8) {
384	for(i = 0; i < mult; i++) kinds.push_back(T11_8bfs[tzindex][channel][i]);
385	} else if (mult == 9) {
386	for(i = 0; i < mult; i++) kinds.push_back(T11_9bfs[tzindex][channel][i]);
387	} else {
388	G4cout << " Illegal multiplicity " << G4endl;
389	}
390
391	return kinds;
392	}
393
394
395	void G4PionSampler::initCrossSections()
396	{
397	/////////////////////////////////////////////////////////////////////////////
398	// //
399	// pi+ p and pi- n (\|Tz\| = 3/2) cross sections //
400	// and final state particle types //
401	// //
402	/////////////////////////////////////////////////////////////////////////////
403
404	// Total pi+ p cross sections as a function of kinetic energy
405	// New cs after 9-body tuning (27 July 09)
406	const G4double pipPtotData[30] =
407	{ 0.0, 1.2, 2.5, 3.8, 5.0, 7.0, 9.0, 15.0, 30.0, 64.0,
408	130.0, 190.0, 130.0, 56.0, 28.0, 17.14, 19.28, 27.4, 40.05, 32.52,
409	30.46, 29.0, 27.26, 25.84, 25.5, 24.5, 24.0, 23.5, 23.0, 23.0};
410
411	for (G4int i = 0; i < 30; i++) pipPtot[i] = pipPtotData[i];
412
413	// Outgoing particle types of a given multiplicity
414	// T33_nbfs = final state types for pi+ p and pi- n
415
416	const G4int T33_2bfsData[2][2][2] =
417	{{{pip,pro}, {kp,sp}},
418
419	{{pim,neu}, {k0,sm}}};
420
421	const G4int T33_3bfsData[2][7][3] =
422	{{{pip,pro,pi0}, {pip,neu,pip}, {pi0,sp,kp}, {pip,sp,k0},
423	{pip,s0,kp}, {pip,lam,kp}, {kp,pro,k0b}},
424
425	{{pim,neu,pi0}, {pim,pro,pim}, {pi0,sm,k0}, {pim,sm,kp},
426	{pim,s0,k0}, {pim,lam,k0}, {k0,neu,km}}};
427
428	const G4int T33_4bfsData[2][15][4] =
429	{{{pip,pro,pip,pim},{pip,pro,pi0,pi0},{pip,neu,pip,pi0},
430	{pip,sp,kp,pim}, {pi0,sp,kp,pi0}, {pip,sp,k0,pi0},
431	{pip,s0,k0,pip}, {pip,s0,kp,pi0}, {pip,lam,kp,pi0},
432	{pip,lam,k0,pip}, {pip,sm,kp,pip}, {pip,pro,kp,km},
433	{pip,pro,k0,k0b}, {pi0,pro,kp,k0b}, {pip,neu,kp,k0b}},
434
435	{{pim,neu,pip,pim},{pim,neu,pi0,pi0},{pim,pro,pim,pi0},
436	{pim,sm,k0,pip}, {pi0,sm,k0,pi0}, {pim,sm,kp,pi0},
437	{pim,s0,kp,pim}, {pim,s0,k0,pi0}, {pim,lam,k0,pi0},
438	{pim,lam,kp,pim}, {pim,sp,k0,pim}, {pim,neu,k0,k0b},
439	{pim,neu,kp,km}, {pi0,neu,k0,km}, {pim,pro,k0,km}}};
440
441	const G4int T33_5bfsData[2][24][5] =
442	{{{pip,pro,pip,pim,pi0}, {pip,pro,pi0,pi0,pi0}, {pip,neu,pip,pip,pim},
443	{pip,neu,pip,pi0,pi0}, {pip,sp,kp,pim,pi0}, {pi0,sp,kp,pi0,pi0},
444	{pip,sp,k0,pip,pim}, {pip,sp,k0,pi0,pi0}, {pip,lam,k0,pip,pi0},
445	{pip,lam,kp,pip,pim}, {pip,lam,kp,pi0,pi0}, {pip,s0,kp,pip,pim},
446	{pip,s0,kp,pi0,pi0}, {pip,s0,k0,pip,pi0}, {pip,sm,kp,pip,pi0},
447	{pip,sm,k0,pip,pip}, {pip,pro,pim,kp,k0b}, {pip,pro,pip,k0,km},
448	{pip,pro,pi0,kp,km}, {pip,pro,pi0,k0,k0b}, {pi0,pro,pi0,kp,k0b},
449	{pip,neu,pip,kp,km}, {pip,neu,pip,k0,k0b}, {pip,neu,pi0,kp,k0b}},
450
451	{{pim,neu,pip,pim,pi0}, {pim,neu,pi0,pi0,pi0}, {pim,pro,pip,pim,pim},
452	{pim,pro,pim,pi0,pi0}, {pim,sm,k0,pip,pi0}, {pi0,sm,k0,pi0,pi0},
453	{pim,sm,kp,pip,pim}, {pim,sm,kp,pi0,pi0}, {pim,lam,kp,pim,pi0},
454	{pim,lam,k0,pip,pim}, {pim,lam,k0,pi0,pi0}, {pim,s0,k0,pip,pim},
455	{pim,s0,k0,pi0,pi0}, {pim,s0,kp,pim,pi0}, {pim,sp,k0,pim,pi0},
456	{pim,sp,kp,pim,pim}, {pim,neu,pip,k0,km}, {pim,neu,pim,kp,k0b},
457	{pim,neu,pi0,k0,k0b}, {pim,neu,pi0,kp,km}, {pi0,neu,pi0,k0,km},
458	{pim,pro,pim,k0,k0b}, {pim,pro,pim,kp,km}, {pim,pro,pi0,k0,km}}};
459
460	const G4int T33_6bfsData[2][5][6] =
461	{{{pip,pro,pip,pip,pim,pim}, {pip,pro,pip,pim,pi0,pi0},
462	{pip,pro,pi0,pi0,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0},
463	{pip,neu,pip,pip,pim,pi0}},
464
465	{{pim,neu,pip,pip,pim,pim}, {pim,neu,pip,pim,pi0,pi0},
466	{pim,neu,pi0,pi0,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0},
467	{pim,pro,pip,pim,pim,pi0}}};
468
469	const G4int T33_7bfsData[2][6][7] =
470	{{{pip,pro,pip,pip,pim,pim,pi0}, {pip,pro,pip,pim,pi0,pi0,pi0},
471	{pip,pro,pi0,pi0,pi0,pi0,pi0}, {pip,neu,pip,pip,pip,pim,pim},
472	{pip,neu,pip,pip,pim,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0,pi0}},
473
474	{{pim,neu,pip,pip,pim,pim,pi0}, {pim,neu,pip,pim,pi0,pi0,pi0},
475	{pim,neu,pi0,pi0,pi0,pi0,pi0}, {pim,pro,pip,pip,pim,pim,pim},
476	{pim,pro,pip,pim,pim,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0,pi0}}};
477
478	const G4int T33_8bfsData[2][7][8] =
479	{{{pip,pro,pip,pip,pip,pim,pim,pim}, {pip,pro,pip,pip,pim,pim,pi0,pi0},
480	{pip,pro,pip,pim,pi0,pi0,pi0,pi0}, {pip,pro,pi0,pi0,pi0,pi0,pi0,pi0},
481	{pip,neu,pip,pip,pip,pim,pim,pi0}, {pip,neu,pip,pip,pim,pi0,pi0,pi0},
482	{pip,neu,pip,pi0,pi0,pi0,pi0,pi0}},
483
484	{{pim,neu,pip,pip,pip,pim,pim,pim}, {pim,neu,pip,pip,pim,pim,pi0,pi0},
485	{pim,neu,pip,pim,pi0,pi0,pi0,pi0}, {pim,neu,pi0,pi0,pi0,pi0,pi0,pi0},
486	{pim,pro,pip,pip,pim,pim,pim,pi0}, {pim,pro,pip,pim,pim,pi0,pi0,pi0},
487	{pim,pro,pim,pi0,pi0,pi0,pi0,pi0}}};
488
489	const G4int T33_9bfsData[2][8][9] =
490	{{{pip,pro,pip,pip,pip,pim,pim,pim,pi0}, {pip,pro,pip,pip,pim,pim,pi0,pi0,pi0},
491	{pip,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
492	{pip,neu,pip,pip,pip,pip,pim,pim,pim}, {pip,neu,pip,pip,pip,pim,pim,pi0,pi0},
493	{pip,neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
494
495	{{pim,neu,pip,pip,pip,pim,pim,pim,pi0}, {pim,neu,pip,pip,pim,pim,pi0,pi0,pi0},
496	{pim,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
497	{pim,pro,pip,pip,pip,pim,pim,pim,pim}, {pim,pro,pip,pip,pim,pim,pim,pi0,pi0},
498	{pim,pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
499
500	for (G4int i = 0; i < 2; i++) {
501	for (G4int j = 0; j < 2; j++) {
502	for (G4int k = 0; k < 2; k++) {
503	T33_2bfs[i][j][k] = T33_2bfsData[i][j][k];
504	}
505	}
506	for (G4int j = 0; j < 7; j++) {
507	for (G4int k = 0; k < 3; k++) {
508	T33_3bfs[i][j][k] = T33_3bfsData[i][j][k];
509	}
510	}
511	for (G4int j = 0; j < 15; j++) {
512	for (G4int k = 0; k < 4; k++) {
513	T33_4bfs[i][j][k] = T33_4bfsData[i][j][k];
514	}
515	}
516	for (G4int j = 0; j < 24; j++) {
517	for (G4int k = 0; k < 5; k++) {
518	T33_5bfs[i][j][k] = T33_5bfsData[i][j][k];
519	}
520	}
521	for (G4int j = 0; j < 5; j++) {
522	for (G4int k = 0; k < 6; k++) {
523	T33_6bfs[i][j][k] = T33_6bfsData[i][j][k];
524	}
525	}
526	for (G4int j = 0; j < 6; j++) {
527	for (G4int k = 0; k < 7; k++) {
528	T33_7bfs[i][j][k] = T33_7bfsData[i][j][k];
529	}
530	}
531	for (G4int j = 0; j < 7; j++) {
532	for (G4int k = 0; k < 8; k++) {
533	T33_8bfs[i][j][k] = T33_8bfsData[i][j][k];
534	}
535	}
536	for (G4int j = 0; j < 8; j++) {
537	for (G4int k = 0; k < 9; k++) {
538	T33_9bfs[i][j][k] = T33_9bfsData[i][j][k];
539	}
540	}
541	}
542
543	//
544	// Cross sections (in mb) for pi+ p -> 2-9 body final states
545	//
546	// first index: 0-1: channels for mult = 2
547	// 2-8: channels for mult = 3
548	// 9-23: channels for mult = 4
549	// 24-47: channels for mult = 5
550	// 48-52: channels for mult = 6
551	// 53-58: channels for mult = 7
552	// 59-65: channels for mult = 8
553	// 66-73: channels for mult = 9
554	//
555	// second index: kinetic energy
556	//
557
558	const G4float pipPCrossSectionsData[74][30] = {
559	//
560	// multiplicity 2 (2 channels)
561	//
562	// p pi+ (n pi-)
563	{ 0.0, 1.2, 2.5, 3.8, 5.0, 7.00, 9.00, 15.0, 30.0, 64.0,
564	130.0, 190.0, 130.0, 55.7, 27.2, 13.95, 8.38, 12.98, 18.53, 11.81,
565	9.4, 7.7, 6.8, 5.9, 5.3, 4.8, 4.2, 4.0, 3.80, 3.6},
566
567	// S+ K+ (S- K0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.16, 0.62, 0.34,
570	0.25, 0.2, 0.1, 0.05, 0.04, 0.02, 0.02, 0.01, 0.01, 0.01},
571	//
572	// multiplicity 3 (7 channels)
573	//
574	// p pi+ pi0 (n pi- pi0)
575	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
576	0.0, 0.0, 0.0, 0.2, 0.6, 2.39, 8.67, 9.99, 12.36, 6.66,
577	4.4, 2.6, 1.95, 1.15, 0.8, 0.6, 0.41, 0.31, 0.24, 0.15},
578
579	// n pi+ pi+ (p pi- pi-)
580	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
581	0.0, 0.0, 0.0, 0.1, 0.20, 0.6, 1.48, 2.30, 3.71, 3.22,
582	2.7, 2.1, 1.5, 1.0, 0.74, 0.5, 0.38, 0.27, 0.20, 0.15},
583
584	// S+ K+ pi0 (S- K0 pi0)
585	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
586	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
587	0.15, 0.12, 0.09, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
588
589	// S+ K0 pi+
590	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
591	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
592	0.15, 0.12, 0.09, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
593
594	// S0 K+ pi+
595	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
596	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
597	0.15, 0.12, 0.09, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
598
599	// L K+ pi+
600	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
601	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
602	0.15, 0.12, 0.09, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
603
604	// p K+ K0bar
605	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
606	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
607	0.07, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
608	//
609	// multiplicity 4 (15 channels)
610	//
611	// p pi+ pi+ pi-
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.06, 0.2, 0.78, 2.26, 3.43,
614	3.9, 3.3, 2.85, 2.4, 2.16, 1.7, 1.4, 1.20, 1.00, 0.90},
615
616	// p pi+ pi0 pi0
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.03, 0.08, 0.15, 1.03, 2.16,
619	2.2, 1.7, 1.15, 0.55, 0.41, 0.33, 0.24, 0.19, 0.14, 0.09},
620
621	// n pi+ pi+ pi0
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.11, 0.45, 0.95, 1.06, 1.51,
624	1.65, 1.5, 1.12, 0.7, 0.43, 0.24, 0.11, 0.05, 0.02, 0.01},
625
626	// S+ K+ pi+ pi-
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.05, 0.13, 0.09, 0.06, 0.03, 0.02, 0.01, 0.01, 0.0},
630
631	// S+ K+ pi0 pi0
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.05, 0.13, 0.09, 0.06, 0.03, 0.02, 0.01, 0.01, 0.0},
635
636	// S+ K0 pi+ pi0
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.02, 0.06, 0.07, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.01},
640
641	// S0 K0 pi+ pi+
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
644	0.0, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
645
646	// S0 K+ pi+ pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
649	0.0, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
650
651	// L K+ pi+ pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
654	0.0, 0.02, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01},
655
656	// L K0 pi+ pi+
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
659	0.0, 0.02, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01},
660
661	// S- K+ pi+ pi+
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
664	0.02, 0.16, 0.18, 0.14, 0.1, 0.06, 0.04, 0.03, 0.03, 0.02},
665
666	// p pi+ K+ K-
667	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
668	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
669	0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04},
670
671	// p pi+ K0 K0bar
672	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
673	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
674	0.01, 0.06, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05, 0.05, 0.04},
675
676	// p pi0 K+ K0bar
677	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
678	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
679	0.01, 0.05, 0.09, 0.08, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04},
680
681	// n pi+ K+ K0bar
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
684	0.01, 0.05, 0.07, 0.06, 0.06, 0.05, 0.04, 0.03, 0.03, 0.02},
685	//
686	// multiplicity 5 (24 channels)
687	//
688	// p pi+ pi+ pi- pi0
689	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
690	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.31, 2.14,
691	3.45, 4.4, 3.5, 2.62, 2.20, 1.8, 1.40, 1.10, 0.89, 0.70},
692
693	// p pi+ pi0 pi0 pi0
694	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
695	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.21,
696	0.34, 0.55, 0.4, 0.29, 0.23, 0.18, 0.14, 0.11, 0.09, 0.07},
697
698	// n pi+ pi+ pi+ pi-
699	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
700	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.05, 0.20,
701	0.4, 0.7, 0.99, 0.92, 0.75, 0.58, 0.42, 0.31, 0.24, 0.18},
702
703	// n pi+ pi+ pi0 pi0
704	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
705	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.05, 0.20,
706	0.4, 0.8, 0.99, 0.92, 0.75, 0.58, 0.42, 0.31, 0.24, 0.18},
707
708	// S+ K+ pi+ pi- pi0
709	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
710	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
711	0.0, 0.02, 0.05, 0.1, 0.08, 0.06, 0.05, 0.04, 0.04, 0.03},
712
713	// S+ K+ pi0 pi0 pi0
714	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
715	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
716	0.0, 0.01, 0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.01},
717
718	// S+ K0 pi+ pi+ pi-
719	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
720	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
721	0.0, 0.01, 0.02, 0.04, 0.04, 0.03, 0.02, 0.02, 0.02, 0.01},
722
723	// S+ K0 pi+ pi0 pi0
724	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
725	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
726	0.0, 0.01, 0.02, 0.04, 0.04, 0.03, 0.02, 0.02, 0.02, 0.01},
727
728	// L K0 pi+ pi+ pi0
729	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
730	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
731	0.0, 0.01, 0.04, 0.05, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
732
733	// L K+ pi+ pi+ pi-
734	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
735	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
736	0.0, 0.01, 0.02, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
737
738	// L K+ pi+ pi0 pi0
739	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
740	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
741	0.0, 0.01, 0.02, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
742
743	// S0 K+ pi+ pi+ pi-
744	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
745	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
746	0.0, 0.0, 0.0, 0.01, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
747
748	// S0 K+ pi+ pi0 pi0
749	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
750	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
751	0.0, 0.0, 0.0, 0.01, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
752
753	// S0 K0 pi+ pi+ pi0
754	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
755	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
756	0.0, 0.0, 0.0, 0.01, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
757
758	// S- K+ pi+ pi+ pi0
759	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
760	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
761	0.0, 0.0, 0.02, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
762
763	// S- K0 pi+ pi+ pi+
764	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
765	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
766	0.0, 0.0, 0.02, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
767
768	// p pi+ pi- K+ K0bar
769	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
770	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
771	0.0, 0.0, 0.02, 0.06, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04},
772
773	// p pi+ pi+ K0 K-
774	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.04, 0.05, 0.05, 0.04, 0.03, 0.03, 0.03},
777
778	// p pi+ pi0 K+ K-
779	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03},
782
783	// p pi+ pi0 K0 K0bar
784	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.04, 0.05, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03},
787
788	// p pi0 pi0 K+ K0bar
789	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02},
792
793	// n pi+ pi+ K+ K-
794	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.06, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02},
797
798	// n pi+ pi+ K0 K0bar
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.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
802
803	// n pi+ pi0 K+ K0bar
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.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
807	//
808	// multiplicity 6 (5 channels)
809	//
810	// p pi+ pi+ pi+ pi- pi-
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
813	0.09, 0.27, 0.35, 0.42, 0.5, 0.44, 0.4, 0.32, 0.29, 0.23},
814
815	// p pi+ pi+ pi- pi0 pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
818	0.09, 0.27, 0.35, 0.42, 0.5, 0.44, 0.4, 0.32, 0.29, 0.23},
819
820	// p pi+ pi0 pi0 pi0 pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
823	0.02, 0.04, 0.06, 0.08, 0.09, 0.08, 0.08, 0.07, 0.06, 0.05},
824
825	// n pi+ pi+ pi0 pi0 pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
828	0.09, 0.27, 0.35, 0.42, 0.5, 0.44, 0.4, 0.32, 0.29, 0.23},
829
830	// n pi+ pi+ pi+ pi- pi0
831	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
832	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
833	0.09, 0.27, 0.35, 0.42, 0.5, 0.44, 0.4, 0.32, 0.29, 0.23},
834	//
835	// multiplicity 7 (6 channels)
836	//
837	// p pi+ pi+ pi+ pi- pi- pi0 (measured R 363)
838	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
839	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
840	0.02, 0.25, 0.5, 1.00, 0.97, 0.96, 0.9, 0.84, 0.78, 0.72},
841
842	// p pi+ pi+ pi- pi0 pi0 pi0
843	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
844	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
845	0.02, 0.25, 0.5, 0.95, 1.0, 0.9, 0.8, 0.75, 0.7, 0.65},
846
847	// p pi+ pi0 pi0 pi0 pi0 pi0
848	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.04, 0.1, 0.24, 0.3, 0.23, 0.2, 0.18, 0.16, 0.14},
851
852	// n pi+ pi+ pi+ pi+ pi- pi-
853	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
854	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
855	0.01, 0.04, 0.08, 0.13, 0.24, 0.32, 0.3, 0.26, 0.24, 0.22},
856
857	// n pi+ pi+ pi+ pi- pi0 pi0
858	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
859	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
860	0.01, 0.04, 0.08, 0.13, 0.3, 0.25, 0.2, 0.18, 0.15, 0.12},
861
862	// n pi+ pi+ pi0 pi0 pi0 pi0
863	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
864	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
865	0.01, 0.04, 0.08, 0.13, 0.26, 0.23, 0.2, 0.18, 0.15, 0.12},
866	//
867	// multiplicity 8 (7 channels)
868	//
869	// p pi+ pi+ pi+ pi+ pi- pi- pi- (measured R 386)
870	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
871	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
872	0.02, 0.03, 0.05, 0.08, 0.13, 0.17, 0.16, 0.15, 0.15, 0.14},
873
874	// p pi+ pi+ pi+ pi- pi- pi0 pi0
875	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
876	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
877	0.01, 0.03, 0.04, 0.07, 0.13, 0.15, 0.14, 0.12, 0.11, 0.10},
878
879	// p pi+ pi+ pi- pi0 pi0 pi0 pi0
880	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
881	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
882	0.01, 0.03, 0.05, 0.08, 0.13, 0.17, 0.16, 0.15, 0.15, 0.14},
883
884	// p pi+ pi0 pi0 pi0 pi0 pi0 pi0
885	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
886	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
887	0.01, 0.02, 0.03, 0.04, 0.06, 0.08, 0.09, 0.08, 0.06, 0.07},
888
889	// n pi+ pi+ pi+ pi+ pi- pi- pi0
890	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
891	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
892	0.04, 0.07, 0.12, 0.19, 0.3, 0.42, 0.4, 0.36, 0.36, 0.32},
893
894	// n pi+ pi+ pi+ pi- pi0 pi0 pi0
895	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
896	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
897	0.04, 0.07, 0.12, 0.19, 0.3, 0.34, 0.32, 0.30, 0.28, 0.26},
898
899	// n pi+ pi+ pi0 pi0 pi0 pi0 pi0
900	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
901	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
902	0.02, 0.03, 0.05, 0.10, 0.15, 0.21, 0.2, 0.18, 0.18, 0.16},
903	//
904	// multiplicity 9 (8 channels)
905	//
906	// p pi+ pi+ pi+ pi+ pi- pi- pi- pi0 (measured R 383)
907	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
908	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
909	0.01, 0.02, 0.04, 0.07, 0.12, 0.19, 0.26, 0.36, 0.36, 0.36},
910
911	// p pi+ pi+ pi+ pi- pi- pi0 pi0 pi0
912	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
913	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
914	0.0, 0.01, 0.02, 0.04, 0.08, 0.13, 0.18, 0.24, 0.24, 0.24},
915
916	// p pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0
917	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
918	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
919	0.0, 0.01, 0.02, 0.03, 0.05, 0.07, 0.10, 0.14, 0.14, 0.14},
920
921	// p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
922	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
923	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
924	0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.06, 0.09, 0.09, 0.09},
925
926	// n pi+ pi+ pi+ pi+ pi+ pi- pi- pi- (measured R 446 )
927	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
928	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
929	0.0, 0.0, 0.02, 0.03, 0.05, 0.08, 0.10, 0.15, 0.15, 0.15},
930
931	// n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
932	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
933	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
934	0.0, 0.0, 0.01, 0.02, 0.03, 0.05, 0.06, 0.10, 0.10, 0.10},
935
936	// n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
937	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
938	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
939	0.0, 0.0, 0.01, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.06},
940
941	// n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
942	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04} };
945
946	// Put array in class scope
947	for (G4int i = 0; i < 74; i++) {
948	for (G4int j = 0; j < 30; j++) {
949	pipPCrossSections[i][j] = pipPCrossSectionsData[i][j];
950	}
951	}
952
953	/////////////////////////////////////////////////////////////////////////////
954	// //
955	// pi- p and pi+ n (\|Tz\| = 1/2) cross sections //
956	// and final state particle types //
957	// //
958	/////////////////////////////////////////////////////////////////////////////
959
960	// Total pi- p cross section as a function of kinetic energy
961	// New cs values after 9B tuning (27 July 09)
962	const G4double pimPtotData[30] =
963	{ 6.13, 6.4, 6.67, 6.94, 7.22, 7.5, 8.3, 12.0, 14.4, 24.0,
964	46.0, 72.04, 43.02, 27.19, 27.32, 43.8, 37.08, 51.37, 34.21, 34.79,
965	32.08, 31.19, 30.32, 28.5, 27.0, 25.9, 25.5, 25.2, 25.0, 24.8};
966
967	for (G4int i = 0; i < 30; i++) pimPtot[i] = pimPtotData[i];
968
969	// Outgoing particle types of a given multiplicity
970	// T31_nbfs = final state types for pi- p and pi+ n
971
972	const G4int T31_2bfsData[2][5][2] =
973	{{{pim,pro}, {pi0,neu}, {k0,lam}, {k0,s0}, {kp,sm}},
974
975	{{pip,neu}, {pi0,pro}, {kp,lam}, {kp,s0}, {k0,sp}}};
976
977	const G4int T31_3bfsData[2][13][3] =
978	{{{pim,pro,pi0}, {pim,neu,pip}, {pi0,neu,pi0}, {pi0,lam,k0},
979	{pim,lam,kp}, {pip,sm,k0}, {pi0,sm,kp}, {pim,sp,k0},
980	{pim,s0,kp}, {pi0,s0,k0}, {km,pro,k0}, {km,neu,kp},
981	{k0b,neu,k0}},
982
983	{{pip,neu,pi0}, {pip,pro,pim}, {pi0,pro,pi0}, {pi0,lam,kp},
984	{pip,lam,k0}, {pim,sp,kp}, {pi0,sp,k0}, {pip,sm,kp},
985	{pip,s0,k0}, {pi0,s0,kp}, {k0b,neu,kp}, {k0b,pro,k0},
986	{km,pro,kp}}};
987
988	const G4int T31_4bfsData[2][22][4] =
989	{{{pim,pro,pip,pim}, {pim,pro,pi0,pi0}, {pim,neu,pip,pi0},
990	{pi0,neu,pi0,pi0}, {pim,lam,k0,pip}, {pi0,lam,k0,pi0},
991	{pim,lam,kp,pi0}, {pim,s0,k0,pip}, {pi0,s0,k0,pi0},
992	{pim,s0,kp,pi0}, {pim,sp,kp,pim}, {pim,sp,k0,pi0},
993	{pim,sm,kp,pip}, {pi0,sm,kp,pi0}, {pip,sm,k0,pi0},
994	{pim,pro,kp,km}, {pim,pro,k0,k0b}, {pi0,pro,k0,km},
995	{pip,neu,k0,km}, {pi0,neu,k0,k0b}, {pi0,neu,kp,km},
996	{pim,neu,kp,k0b}},
997
998	{{pip,neu,pip,pim}, {pip,neu,pi0,pi0}, {pip,pro,pim,pi0},
999	{pi0,pro,pi0,pi0}, {pip,lam,kp,pim}, {pi0,lam,kp,pi0},
1000	{pip,lam,k0,pi0}, {pip,s0,kp,pim}, {pi0,s0,kp,pi0},
1001	{pip,s0,k0,pi0}, {pip,sm,k0,pip}, {pip,sm,kp,pi0},
1002	{pip,sp,k0,pim}, {pi0,sp,k0,pi0}, {pim,sp,kp,pi0},
1003	{pip,neu,k0,k0b}, {pip,neu,kp,km}, {pi0,neu,kp,k0b},
1004	{pim,pro,kp,k0b}, {pi0,pro,kp,km}, {pi0,pro,k0,k0b},
1005	{pip,pro,k0,km}}};
1006
1007	const G4int T31_5bfsData[2][31][5] =
1008	{{{pim,pro,pip,pim,pi0}, {pim,pro,pi0,pi0,pi0}, {pim,neu,pip,pip,pim},
1009	{pim,neu,pip,pi0,pi0}, {pi0,neu,pi0,pi0,pi0}, {pim,lam,k0,pip,pi0},
1010	{pim,lam,kp,pi0,pi0}, {pim,lam,kp,pip,pim}, {pi0,lam,k0,pi0,pi0},
1011	{pim,s0,kp,pip,pim}, {pim,s0,kp,pi0,pi0}, {pim,s0,k0,pip,pi0},
1012	{pi0,s0,k0,pi0,pi0}, {pim,sp,k0,pip,pim}, {pim,sp,k0,pi0,pi0},
1013	{pim,sp,kp,pim,pi0}, {pim,sm,k0,pip,pip}, {pip,sm,k0,pi0,pi0},
1014	{pim,sm,kp,pip,pi0}, {pi0,sm,kp,pi0,pi0}, {pim,pro,pi0,kp,km},
1015	{pim,pro,pi0,k0,k0b}, {pim,pro,pip,k0,km}, {pi0,pro,pi0,k0,km},
1016	{pim,pro,pim,kp,k0b}, {pim,neu,pip,kp,km}, {pim,neu,pip,k0,k0b},
1017	{pip,neu,pi0,k0,km}, {pim,neu,pi0,kp,k0b}, {pi0,neu,pi0,k0,k0b},
1018	{pi0,neu,pi0,kp,km}},
1019
1020	{{pip,neu,pip,pim,pi0}, {pip,neu,pi0,pi0,pi0}, {pip,pro,pip,pim,pim},
1021	{pip,pro,pim,pi0,pi0}, {pi0,pro,pi0,pi0,pi0}, {pip,lam,kp,pim,pi0},
1022	{pip,lam,k0,pi0,pi0}, {pip,lam,k0,pip,pim}, {pi0,lam,kp,pi0,pi0},
1023	{pip,s0,k0,pip,pim}, {pip,s0,k0,pi0,pi0}, {pip,s0,kp,pim,pi0},
1024	{pi0,s0,kp,pi0,pi0}, {pip,sm,kp,pip,pim}, {pip,sm,kp,pi0,pi0},
1025	{pip,sm,k0,pip,pi0}, {pip,sp,kp,pim,pim}, {pim,sp,kp,pi0,pi0},
1026	{pip,sp,k0,pim,pi0}, {pi0,sp,k0,pi0,pi0}, {pip,neu,pi0,k0,k0b},
1027	{pip,neu,pi0,kp,km}, {pip,neu,pim,kp,k0b}, {pi0,neu,pi0,kp,k0b},
1028	{pip,neu,pip,k0,km}, {pip,pro,pim,k0,k0b}, {pip,pro,pim,kp,km},
1029	{pim,pro,pi0,kp,k0b}, {pip,pro,pi0,k0,km}, {pi0,pro,pi0,kp,km},
1030	{pi0,pro,pi0,k0,k0b}}};
1031
1032	const G4int T31_6bfsData[2][6][6] =
1033	{{{pim,pro,pip,pip,pim,pim}, {pim,pro,pip,pim,pi0,pi0},
1034	{pim,pro,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pim,pi0},
1035	{pim,neu,pip,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0}},
1036
1037	{{pip,neu,pip,pip,pim,pim}, {pip,neu,pip,pim,pi0,pi0},
1038	{pip,neu,pi0,pi0,pi0,pi0}, {pip,pro,pip,pim,pim,pi0},
1039	{pip,pro,pim,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0}}};
1040
1041	const G4int T31_7bfsData[2][7][7] =
1042	{{{pim,pro,pip,pip,pim,pim,pi0}, {pim,pro,pip,pim,pi0,pi0,pi0},
1043	{pim,pro,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim},
1044	{pim,neu,pip,pip,pim,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0},
1045	{pi0,neu,pi0,pi0,pi0,pi0,pi0}},
1046
1047	{{pip,neu,pip,pip,pim,pim,pi0}, {pip,neu,pip,pim,pi0,pi0,pi0},
1048	{pip,neu,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim},
1049	{pip,pro,pip,pim,pim,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0},
1050	{pi0,pro,pi0,pi0,pi0,pi0,pi0}}};
1051
1052	const G4int T31_8bfsData[2][8][8] =
1053	{{{pim,pro,pip,pip,pip,pim,pim,pim}, {pim,pro,pip,pip,pim,pim,pi0,pi0},
1054	{pim,pro,pip,pim,pi0,pi0,pi0,pi0}, {pim,pro,pi0,pi0,pi0,pi0,pi0,pi0},
1055	{pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0,pi0},
1056	{pim,neu,pip,pip,pim,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim,pi0}},
1057
1058	{{pip,neu,pip,pip,pip,pim,pim,pim}, {pip,neu,pip,pip,pim,pim,pi0,pi0},
1059	{pip,neu,pip,pim,pi0,pi0,pi0,pi0}, {pip,neu,pi0,pi0,pi0,pi0,pi0,pi0},
1060	{pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0,pi0},
1061	{pip,pro,pip,pim,pim,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim,pi0}}};
1062
1063	const G4int T31_9bfsData[2][9][9] =
1064	{{{pim,pro,pip,pip,pip,pim,pim,pim,pi0}, {pim,pro,pip,pip,pim,pim,pi0,pi0,pi0},
1065	{pim,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pim,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
1066	{pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
1067	{pim,neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim,pi0,pi0},
1068	{pim,neu,pip,pip,pip,pip,pim,pim,pim}},
1069
1070	{{pip,neu,pip,pip,pip,pim,pim,pim,pi0}, {pip,neu,pip,pip,pim,pim,pi0,pi0,pi0},
1071	{pip,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pip,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
1072	{pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
1073	{pip,pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim,pi0,pi0},
1074	{pip,pro,pip,pip,pip,pim,pim,pim,pim}}};
1075
1076	for (G4int i = 0; i < 2; i++) {
1077	for (G4int j = 0; j < 5; j++) {
1078	for (G4int k = 0; k < 2; k++) {
1079	T31_2bfs[i][j][k] = T31_2bfsData[i][j][k];
1080	}
1081	}
1082	for (G4int j = 0; j < 13; j++) {
1083	for (G4int k = 0; k < 3; k++) {
1084	T31_3bfs[i][j][k] = T31_3bfsData[i][j][k];
1085	}
1086	}
1087	for (G4int j = 0; j < 22; j++) {
1088	for (G4int k = 0; k < 4; k++) {
1089	T31_4bfs[i][j][k] = T31_4bfsData[i][j][k];
1090	}
1091	}
1092	for (G4int j = 0; j < 31; j++) {
1093	for (G4int k = 0; k < 5; k++) {
1094	T31_5bfs[i][j][k] = T31_5bfsData[i][j][k];
1095	}
1096	}
1097	for (G4int j = 0; j < 6; j++) {
1098	for (G4int k = 0; k < 6; k++) {
1099	T31_6bfs[i][j][k] = T31_6bfsData[i][j][k];
1100	}
1101	}
1102	for (G4int j = 0; j < 7; j++) {
1103	for (G4int k = 0; k < 7; k++) {
1104	T31_7bfs[i][j][k] = T31_7bfsData[i][j][k];
1105	}
1106	}
1107	for (G4int j = 0; j < 8; j++) {
1108	for (G4int k = 0; k < 8; k++) {
1109	T31_8bfs[i][j][k] = T31_8bfsData[i][j][k];
1110	}
1111	}
1112	for (G4int j = 0; j < 9; j++) {
1113	for (G4int k = 0; k < 9; k++) {
1114	T31_9bfs[i][j][k] = T31_9bfsData[i][j][k];
1115	}
1116	}
1117	}
1118
1119	//
1120	// Cross sections (in mb) for pi- p -> 2-9 body final states
1121	//
1122	// first index: 0-4: channels for mult = 2
1123	// 5-17: channels for mult = 3
1124	// 18-39: channels for mult = 4
1125	// 40-70: channels for mult = 5
1126	// 71-76: channels for mult = 6
1127	// 77-83: channels for mult = 7
1128	// 84-91: channels for mult = 8
1129	// 92-100: channels for mult = 9
1130	//
1131	// second index: kinetic energy
1132	//
1133	const G4float pimPCrossSectionsData[101][30] = {
1134	//
1135	// multiplicity 2 (5 channels)
1136	//
1137	// pi- p (pi+ n)
1138	{ 1.43, 1.5, 1.57, 1.64, 1.72, 1.8, 2.0, 3.0, 3.4, 7.0,
1139	14.0, 24.0, 14.7, 10.5, 10.84, 21.0, 13.97, 25.1, 10.46, 9.88,
1140	8.0, 7.1, 6.0, 5.7, 5.0, 4.6, 4.3, 4.0, 3.8, 3.7},
1141
1142	// n pi0 (p pi0)
1143	{ 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 6.3, 9.0, 11.0, 17.0,
1144	32.0, 48.0, 28.0, 14.5, 11.04, 8.99, 4.79, 5.02, 2.08, 0.94,
1145	0.5, 0.25, 0.15, 0.09, 0.06, 0.05, 0.04, 0.02, 0.02, 0.01},
1146
1147	// L K0 (L K+)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.65, 0.29, 0.17,
1150	0.13, 0.08, 0.05, 0.03, 0.02, 0.01, 0.01, 0.01, 0.01, 0.0},
1151
1152	// S0 K0 (S0 K+)
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.0, 0.0, 0.0, 0.0, 0.0, 0.28, 0.19, 0.12,
1155	0.09, 0.06, 0.04, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0},
1156
1157	// S- K+ (S+ K0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.20, 0.25, 0.09,
1160	0.04, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
1161	//
1162	// multiplicity 3 (13 channels)
1163	//
1164	// p pi- pi0 (n pi+ pi0)
1165	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1166	0.0, 0.0, 0.05, 0.18, 0.86, 3.89, 5.19, 6.63, 5.11, 4.03,
1167	3.2, 2.5, 2.0, 1.4, 0.97, 0.68, 0.55, 0.36, 0.3, 0.22},
1168
1169	// n pi+ pi- (p pi+ pi-)
1170	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1171	0.0, 0.03, 0.19, 1.41, 3.12, 6.2, 8.14, 8.23, 7.84, 5.91,
1172	4.58, 3.58, 2.71, 1.73, 1.06, 0.85, 0.62, 0.38, 0.31, 0.22},
1173
1174	// n pi0 pi0 (p pi0 pi0)
1175	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1176	0.0, 0.01, 0.08, 0.6, 1.34, 2.7, 2.99, 2.51, 1.35, 0.87,
1177	0.61, 0.24, 0.15, 0.10, 0.06, 0.04, 0.0, 0.0, 0.0, 0.0},
1178
1179	// L K0 pi0 (L K+ pi0)
1180	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1181	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.19,
1182	0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.0},
1183
1184	// L K+ pi- (L K0 pi+)
1185	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1186	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.12,
1187	0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01, 0.01, 0.0},
1188
1189	// S- K0 pi+ (S+ K+ pi-)
1190	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1191	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.1,
1192	0.13, 0.07, 0.03, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
1193
1194	// S- K+ pi0 (S+ K0 pi0)
1195	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1196	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
1197	0.05, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
1198
1199	// S+ K0 pi- (S- K+ pi+)
1200	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1201	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
1202	0.06, 0.04, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
1203
1204	// S0 K+ pi- (S0 K0 pi+)
1205	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1206	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
1207	0.07, 0.04, 0.03, 0.02, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0},
1208
1209	// S0 K0 pi0 (S0 K+ pi0)
1210	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1211	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.09,
1212	0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0},
1213
1214	// p K0 K- (n K+ K0bar)
1215	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1216	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03,
1217	0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
1218
1219	// n K+ K- (p K0 K0bar)
1220	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1221	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
1222	0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01, 0.0, 0.0, 0.0},
1223
1224	// n K0 K0bar (p K+ K-)
1225	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1226	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06,
1227	0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0},
1228	//
1229	// multiplicity 4 (22 channels)
1230	//
1231	// p pi+ pi- pi- (n pi+ pi+ pi-)
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.01, 0.04, 0.18, 0.43, 0.82, 1.69,
1234	1.8, 1.99, 1.8, 1.71, 1.44, 1.26, 1.17, 0.99, 1.04, 0.9},
1235
1236	// p pi- pi0 pi0 (n pi+ pi0 pi0)
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.01, 0.03, 0.09, 0.21, 1.5, 2.5,
1239	1.2, 1.2, 1.0, 0.9, 0.8, 0.7, 0.65, 0.6, 0.55, 0.5},
1240
1241	// n pi+ pi- pi0 (p pi+ pi- 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.4, 0.7, 1.11, 2.6, 3.9,
1244	4.07, 3.8, 2.76, 1.38, 1.16, 0.97, 0.85, 0.7, 0.55, 0.4},
1245
1246	// n pi0 pi0 pi0 (p pi0 pi0 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.1, 0.55, 1.0, 0.87, 0.73, 0.71,
1249	0.7, 0.56, 0.35, 0.29, 0.21, 0.14, 0.1, 0.06, 0.03, 0.0},
1250
1251	// L K0 pi+ pi- (L K+ pi+ 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.01,
1254	0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03},
1255
1256	// L K0 pi0 pi0 (L K+ pi0 pi0)
1257	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1258	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1259	0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03},
1260
1261	// L K+ pi- pi0 (L K0 pi+ pi0)
1262	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1263	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1264	0.06, 0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02},
1265
1266	// S0 K0 pi+ pi- (S0 K+ pi+ 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.01,
1269	0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
1270
1271	// S0 K0 pi0 pi0 (S0 K+ pi0 pi0)
1272	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1273	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1274	0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
1275
1276	// S0 K+ pi- pi0 (S0 K0 pi+ pi0)
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.01,
1279	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1280
1281	// S+ K+ pi- pi- (S- K0 pi+ 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.01,
1284	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1285
1286	// S+ K0 pi- pi0 (S- K+ pi+ 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.01,
1289	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1290
1291	// S- K+ pi+ pi- (S+ K0 pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1294	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1295
1296	// S- K+ pi0 pi0 (S+ K0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1299	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1300
1301	// S- K0 pi+ pi0 (S+ K+ pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1304	0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
1305
1306	// p pi- K+ K- (n pi+ K0 K0bar)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1309	0.0, 0.06, 0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05},
1310
1311	// p pi- K0 K0bar (n pi+ K+ K-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1314	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1315
1316	// p pi0 K0 K- (n pi0 K+ K0bar)
1317	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1318	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1319	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1320
1321	// n pi+ K0 K- (p pi- K+ K0bar)
1322	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1323	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1324	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1325
1326	// n pi0 K0 K0bar (p pi0 K+ K-)
1327	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1328	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1329	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1330
1331	// n pi0 K+ K- (p pi0 K0 K0bar)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1334	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1335
1336	// n pi- K+ K0bar (p pi+ K0 K-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1339	0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
1340	//
1341	// multiplicity 5 (31 channels)
1342	//
1343	// p pi+ pi- pi- pi0 (n pi+ pi+ pi- pi0)
1344	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1345	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.10, 0.93,
1346	1.5, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62},
1347
1348	// p pi- pi0 pi0 pi0 (n pi+ pi0 pi0 pi0)
1349	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1350	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.10, 0.73,
1351	1.4, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62},
1352
1353	// n pi+ pi+ pi- pi- (p pi+ pi+ pi- pi-)
1354	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1355	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.10, 0.31,
1356	0.66, 0.93, 1.2, 1.2, 1.2, 0.94, 0.74, 0.53, 0.40, 0.30},
1357
1358	// n pi+ pi- pi0 pi0 (p pi+ pi- pi0 pi0)
1359	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1360	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.22, 0.4,
1361	0.8, 0.95, 1.0, 0.9, 0.75, 0.65, 0.55, 0.40, 0.35, 0.30},
1362
1363	// n pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0)
1364	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1365	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.22, 0.39,
1366	0.66, 0.9, 1.0, 0.8, 0.7, 0.55, 0.5, 0.35, 0.30, 0.25},
1367
1368	// L K0 pi+ pi- pi0 (L K+ pi+ pi- pi0)
1369	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1370	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1371	0.01, 0.05, 0.07, 0.08, 0.08, 0.07, 0.06, 0.06, 0.05, 0.05},
1372
1373	// L K+ pi- pi0 pi0 (L K0 pi+ pi0 pi0)
1374	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1375	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1376	0.0, 0.02, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
1377
1378	// L K+ pi+ pi- pi- (L K0 pi+ pi+ pi-)
1379	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1380	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1381	0.0, 0.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
1382
1383	// L K0 pi0 pi0 pi0 (L K+ pi0 pi0 pi0)
1384	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1385	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1386	0.0, 0.01, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1387
1388	// S0 K+ pi+ pi- pi- (S0 K0 pi+ pi+ pi-)
1389	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1390	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1391	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
1392
1393	// S0 K+ pi- pi0 pi0 (S0 K0 pi+ pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1396	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
1397
1398	// S0 K0 pi+ pi- pi0 (S0 K+ pi+ pi- pi0)
1399	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1400	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1401	0.0, 0.02, 0.03, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02},
1402
1403	// S0 K0 pi0 pi0 pi0 (S0 K+ pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1406	0.0, 0.01, 0.01, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1407
1408	// S+ K0 pi+ pi- pi- (S- K+ pi+ pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1411	0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
1412
1413	// S+ K0 pi- pi0 pi0 (S- K+ pi+ pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1416	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
1417
1418	// S+ K+ pi- pi- pi0 (S- K0 pi+ pi+ pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1421	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
1422
1423	// S- K0 pi+ pi+ pi- (S+ K+ pi+ pi- pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1426	0.0, 0.01, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
1427
1428	// S- K0 pi+ pi0 pi0 (S+ K+ pi- pi0 pi0)
1429	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
1432
1433	// S- K+ pi+ pi- pi0 (S+ K0 pi+ pi- pi0)
1434	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.02, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
1437
1438	// S- K+ pi0 pi0 pi0 (S+ K0 pi0 pi0 pi0)
1439	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
1442
1443	// p pi- pi0 K+ K- (n pi+ pi0 K0 K0bar)
1444	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04},
1447
1448	// p pi- pi0 K0 K0bar (n pi+ pi0 K+ K-)
1449	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1452
1453	// p pi+ pi- K0 K- (n pi+ pi- K+ K0bar)
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.0, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1457
1458	// p pi0 pi0 K0 K- (n pi0 pi0 K+ K0bar)
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.01, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1462
1463	// p pi- pi- K+ K0bar (n pi+ pi+ K0 K-)
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.01, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02},
1467
1468	// n pi+ pi- K+ K- (p pi+ pi- K0 K0bar)
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.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03},
1472
1473	// n pi+ pi- K0 K0bar (p pi+ pi- K+ K-)
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.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
1477
1478	// n pi+ pi0 K0 K- (p pi- pi0 K+ K0bar)
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.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03},
1482
1483	// n pi- pi0 K+ K0bar (p pi+ pi0 K0 K-)
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.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03},
1487
1488	// n pi0 pi0 K0 K0bar (p pi0 pi0 K+ K-)
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.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
1492
1493	// n pi0 pi0 K+ K- (p pi0 pi0 K0 K0bar)
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.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
1497	//
1498	// multiplicity 6 (6 channels)
1499	//
1500	// p pi+ pi+ pi- pi- pi- (n pi+ pi+ pi+ pi- pi-)
1501	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1502	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.04,
1503	0.06, 0.12, 0.16, 0.22, 0.31, 0.34, 0.3, 0.24, 0.19, 0.16},
1504
1505	// p pi+ pi- pi- pi0 pi0 (n pi+ pi+ pi- pi0 pi0)
1506	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1507	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.05,
1508	0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
1509
1510	// p pi- pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0)
1511	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1512	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.05,
1513	0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
1514
1515	// n pi+ pi+ pi- pi- pi0 (p pi+ pi+ pi- pi- pi0)
1516	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1517	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.05,
1518	0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
1519
1520	// n pi+ pi- pi0 pi0 pi0 (p pi+ pi- pi0 pi0 pi0)
1521	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1522	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
1523	0.05, 0.06, 0.08, 0.1, 0.12, 0.13, 0.14, 0.13, 0.11, 0.11},
1524
1525	// n pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 pi0)
1526	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1527	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
1528	0.04, 0.05, 0.06, 0.08, 0.09, 0.1, 0.1, 0.1, 0.08, 0.07},
1529	//
1530	// multiplicity 7 (7 channels)
1531	//
1532	// p pi+ pi+ pi- pi- pi- pi0
1533	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1534	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
1535	0.05, 0.12, 0.19, 0.37, 0.67, 0.70, 0.65, 0.60, 0.50, 0.45},
1536
1537	// p pi+ pi- pi- pi0 pi0 pi0
1538	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1539	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
1540	0.03, 0.07, 0.11, 0.22, 0.4, 0.4, 0.35, 0.30, 0.25, 0.2},
1541
1542	// p pi- pi0 pi0 pi0 pi0 pi0
1543	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1544	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1545	0.02, 0.03, 0.05, 0.11, 0.20, 0.2, 0.15, 0.12, 0.10, 0.08},
1546
1547	// n pi+ pi+ pi+ pi- pi- pi- (measured R 501)
1548	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1549	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1550	0.01, 0.03, 0.06, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13},
1551
1552	// n pi+ pi+ pi- pi- pi0 pi0 (measured R 496)
1553	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1554	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
1555	0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3},
1556
1557	// n pi+ pi- pi0 pi0 pi0 pi0
1558	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1559	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1560	0.02, 0.04, 0.07, 0.14, 0.26, 0.26, 0.22, 0.20, 0.18, 0.15},
1561
1562	// n pi0 pi0 pi0 pi0 pi0 pi0
1563	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1564	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1565	0.01, 0.02, 0.03, 0.05, 0.1, 0.11, 0.1, 0.1, 0.08, 0.7},
1566	//
1567	// multiplicity 8 (8 channels)
1568	//
1569	// p pi+ pi+ pi+ pi- pi- pi- pi- (n pi+ pi+ pi+ pi+ pi- pi- pi-) (measured R 420)
1570	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1571	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1572	0.0, 0.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09},
1573
1574	// p pi+ pi+ pi- pi- pi- pi0 pi0 (n pi+ pi+ pi+ pi- pi- pi0 pi0)
1575	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1576	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1577	0.0, 0.01, 0.03, 0.09, 0.13, 0.13, 0.11, 0.10, 0.09, 0.08},
1578
1579	// p pi+ pi- pi- pi0 pi0 pi0 pi0 (n pi+ pi+ pi- pi0 pi0 pi0 pi0)
1580	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1581	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1582	0.0, 0.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09},
1583
1584	// p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
1585	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1586	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1587	0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03},
1588
1589	// n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 pi0 pi0 pi0)
1590	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1591	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1592	0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.03, 0.03, 0.02},
1593
1594	// n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p pi+ pi- pi0 pi0 pi0 pi0 pi0)
1595	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1596	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1597	0.0, 0.0, 0.01, 0.03, 0.04, 0.05, 0.05, 0.05, 0.05, 0.04},
1598
1599	// n pi+ pi+ pi- pi- pi0 pi0 pi0 (p pi+ pi+ pi- pi- pi0 pi0 pi0)
1600	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1601	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1602	0.0, 0.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
1603
1604	// n pi+ pi+ pi+ pi- pi- pi- pi0 (p pi+ pi+ pi+ pi- pi- pi- pi0)
1605	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1606	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1607	0.0, 0.01, 0.03, 0.09, 0.12, 0.13, 0.11, 0.09, 0.08, 0.07},
1608	//
1609	// multiplicity 9 (9 channels)
1610	//
1611	// p pi+ pi+ pi+ pi- pi- pi- pi- pi0 (measured R 418)
1612	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1613	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1614	0.0, 0.0, 0.01, 0.04, 0.08, 0.14, 0.18, 0.18, 0.18, 0.17},
1615
1616	// p pi+ pi+ pi- pi- pi- pi0 pi0 pi0
1617	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1618	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1619	0.0, 0.0, 0.01, 0.02, 0.05, 0.08, 0.11, 0.11, 0.11, 0.1},
1620
1621	// p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0
1622	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1623	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1624	0.0, 0.0, 0.0, 0.02, 0.03, 0.06, 0.07, 0.07, 0.07, 0.07},
1625
1626	// p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0
1627	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1628	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1629	0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
1630
1631	// n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
1632	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1633	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1634	0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
1635
1636	// n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0
1637	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1638	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1639	0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
1640
1641	// n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0
1642	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1643	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1644	0.0, 0.0, 0.0, 0.02, 0.03, 0.06, 0.07, 0.07, 0.07, 0.07},
1645
1646	// n pi+ pi+ pi+ pi- pi- pi- pi0 pi0
1647	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1648	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1649	0.0, 0.0, 0.01, 0.02, 0.05, 0.08, 0.11, 0.11, 0.11, 0.1},
1650
1651	// n pi+ pi+ pi+ pi+ pi- pi- pi- pi- (measured R 503)
1652	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1653	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1654	0.0, 0.0, 0.01, 0.04, 0.08, 0.14, 0.18, 0.18, 0.18, 0.17} };
1655
1656	// Put array in class scope
1657	for (G4int i = 0; i < 101; i++) {
1658	for (G4int j = 0; j < 30; j++) {
1659	pimPCrossSections[i][j] = pimPCrossSectionsData[i][j];
1660	}
1661	}
1662
1663	/////////////////////////////////////////////////////////////////////////////
1664	// //
1665	// pi0 p and pi0 n ( \|T, Tz> = \|1/2, 1/2>, \|1/2, -1/2> ) cross sections //
1666	// and final state particle types //
1667	// //
1668	/////////////////////////////////////////////////////////////////////////////
1669
1670	// Total pi0 p cross section as a function of kinetic energy
1671	// G4double pizPtotData[30] =
1672	// { 0.0, 3.55, 4.65, 5.9, 7.75, 10.1, 11.8, 18.0, 27.7, 52.5,
1673	// 102.0, 150.0, 102.64, 51.03, 34.94, 34.52, 32.45, 44.05, 40.2, 34.93,
1674	// 32.0, 30.0, 28.29, 26.91, 26.25, 25.25, 24.75, 24.35, 24.0, 23.9};
1675
1676	// New test
1677	const G4double pizPtotData[30] =
1678	{ 6.43, 7.18, 7.54, 8.01, 8.52, 9.13, 10.22, 14.37, 20.96, 34.73,
1679	61.07, 98.23, 61.97, 32.62, 28.07, 31.37, 35.15, 40.17, 37.27, 33.49,
1680	31.06, 29.52, 28.29, 26.91, 26.25, 25.25, 24.75, 24.35, 24.0, 23.9};
1681
1682	for (G4int i = 0; i < 30; i++) pizPtot[i] = pizPtotData[i];
1683
1684	// Outgoing particle types of a given multiplicity
1685	// T11_nbfs = final state types for pi0 p and pi0 n
1686
1687	const G4int T11_2bfsData[2][5][2] =
1688	{{{pi0,pro}, {pip,neu}, {kp,lam}, {kp,s0}, {k0,sp}},
1689
1690	{{pi0,neu}, {pim,pro}, {k0,lam}, {k0,s0}, {kp,sm}}};
1691
1692	const G4int T11_3bfsData[2][13][3] =
1693	{{{pip,pro,pim}, {pi0,pro,pi0}, {pi0,neu,pip}, {pi0,lam,kp},
1694	{pip,lam,k0}, {pi0,s0,kp}, {pip,s0,k0}, {pi0,sp,k0},
1695	{pim,sp,kp}, {pip,sm,kp}, {km,pro,kp}, {k0b,pro,k0},
1696	{k0b,neu,kp}},
1697
1698	{{pim,neu,pip}, {pi0,neu,pi0}, {pi0,pro,pim}, {pi0,lam,k0},
1699	{pim,lam,kp}, {pi0,s0,k0}, {pim,s0,kp}, {pi0,sm,kp},
1700	{pip,sm,k0}, {pim,sp,k0}, {k0b,neu,k0}, {km,neu,kp},
1701	{km,pro,k0}}};
1702
1703	const G4int T11_4bfsData[2][21][4] =
1704	{{{pi0,pro,pim,pip}, {pi0,pro,pi0,pi0}, {pip,neu,pip,pim},
1705	{pi0,neu,pip,pi0}, {pip,lam,kp,pim}, {pip,s0,kp,pim},
1706	{pi0,lam,kp,pi0}, {pi0,s0,kp,pi0}, {pi0,sp,k0,pi0},
1707	{pi0,lam,k0,pip}, {pi0,s0,k0,pip}, {pip,sp,k0,pim},
1708	{pi0,sp,kp,pim}, {pi0,sm,kp,pip}, {pi0,pro,kp,km},
1709	{pi0,pro,k0,k0b}, {pip,pro,k0,km}, {pim,pro,kp,k0b},
1710	{pip,neu,kp,km}, {pip,neu,k0,k0b}, {pi0,neu,kp,k0b}},
1711
1712	{{pi0,neu,pim,pip}, {pi0,neu,pi0,pi0}, {pim,pro,pip,pim},
1713	{pi0,pro,pim,pi0}, {pip,lam,k0,pim}, {pip,s0,k0,pim},
1714	{pi0,lam,k0,pi0}, {pi0,s0,k0,pi0}, {pi0,sm,kp,pi0},
1715	{pi0,lam,kp,pim}, {pi0,s0,kp,pim}, {pip,sm,kp,pim},
1716	{pi0,sm,k0,pip}, {pi0,sp,k0,pim}, {pi0,neu,k0,k0b},
1717	{pi0,neu,kp,km}, {pim,neu,kp,k0b}, {pip,neu,k0,km},
1718	{pim,pro,k0,k0b}, {pim,pro,kp,km}, {pi0,pro,k0,km}}};
1719
1720	const G4int T11_5bfsData[2][30][5] =
1721	{{{pip,pro,pip,pim,pim}, {pi0,pro,pip,pim,pi0}, {pi0,pro,pi0,pi0,pi0},
1722	{pi0,neu,pip,pip,pim}, {pi0,neu,pip,pi0,pi0}, {pi0,lam,kp,pip,pim},
1723	{pi0,lam,kp,pi0,pi0}, {pip,lam,k0,pip,pim}, {pi0,lam,k0,pip,pi0},
1724	{pi0,s0,kp,pip,pim}, {pi0,s0,kp,pi0,pi0}, {pip,s0,k0,pip,pim},
1725	{pi0,s0,k0,pip,pi0}, {pi0,sp,k0,pip,pim}, {pi0,sp,k0,pi0,pi0},
1726	{pim,sp,kp,pip,pim}, {pi0,sp,kp,pim,pi0}, {pip,sm,kp,pip,pim},
1727	{pi0,sm,kp,pip,pi0}, {pi0,pro,kp,km,pi0}, {pi0,pro,k0,k0b,pi0},
1728	{pi0,pro,kp,k0b,pim}, {pi0,pro,k0,km,pip}, {pip,pro,kp,km,pim},
1729	{pip,pro,k0,k0b,pim}, {pi0,neu,kp,km,pip}, {pi0,neu,k0,k0b,pip},
1730	{pi0,neu,kp,k0b,pi0}, {pip,neu,k0,km,pip}, {pip,neu,kp,k0b,pim}},
1731
1732	{{pim,neu,pip,pip,pim}, {pi0,neu,pip,pim,pi0}, {pi0,neu,pi0,pi0,pi0},
1733	{pi0,pro,pip,pim,pim}, {pi0,pro,pim,pi0,pi0}, {pi0,lam,k0,pip,pim},
1734	{pi0,lam,k0,pi0,pi0}, {pim,lam,kp,pip,pim}, {pi0,lam,kp,pim,pi0},
1735	{pi0,s0,k0,pip,pim}, {pi0,s0,k0,pi0,pi0}, {pim,s0,kp,pip,pim},
1736	{pi0,s0,kp,pim,pi0}, {pi0,sm,kp,pip,pim}, {pi0,sm,kp,pi0,pi0},
1737	{pip,sm,k0,pip,pim}, {pi0,sm,k0,pip,pi0}, {pim,sp,k0,pip,pim},
1738	{pi0,sp,k0,pim,pi0}, {pi0,neu,k0,k0b,pi0}, {pi0,neu,kp,km,pi0},
1739	{pi0,neu,k0,km,pip}, {pi0,neu,kp,k0b,pim}, {pim,neu,k0,k0b,pip},
1740	{pim,neu,kp,km,pip}, {pi0,pro,k0,k0b,pim}, {pi0,pro,kp,km,pim},
1741	{pi0,pro,k0,km,pi0}, {pim,pro,kp,k0b,pim}, {pim,pro,k0,km,pip}}};
1742
1743	const G4int T11_6bfsData[2][6][6] =
1744	{{{pi0,pro,pip,pip,pim,pim}, {pi0,pro,pip,pim,pi0,pi0},
1745	{pi0,pro,pi0,pi0,pi0,pi0}, {pip,neu,pip,pip,pim,pim},
1746	{pi0,neu,pip,pip,pim,pi0}, {pi0,neu,pip,pi0,pi0,pi0}},
1747
1748	{{pi0,neu,pip,pip,pim,pim}, {pi0,neu,pip,pim,pi0,pi0},
1749	{pi0,neu,pi0,pi0,pi0,pi0}, {pim,pro,pip,pip,pim,pim},
1750	{pi0,pro,pip,pim,pim,pi0}, {pi0,pro,pim,pi0,pi0,pi0}}};
1751
1752	const G4int T11_7bfsData[2][7][7] =
1753	{{{pip,pro,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pi0},
1754	{pi0,pro,pip,pim,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0,pi0},
1755	{pi0,neu,pip,pip,pip,pim,pim}, {pi0,neu,pip,pip,pim,pi0,pi0},
1756	{pi0,neu,pip,pi0,pi0,pi0,pi0}},
1757
1758	{{pim,neu,pip,pip,pip,pim,pim}, {pi0,neu,pip,pip,pim,pim,pi0},
1759	{pi0,neu,pip,pim,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0,pi0},
1760	{pi0,pro,pip,pip,pim,pim,pim}, {pi0,pro,pip,pim,pim,pi0,pi0},
1761	{pi0,pro,pim,pi0,pi0,pi0,pi0}}};
1762
1763	const G4int T11_8bfsData[2][8][8] =
1764	{{{pi0,pro,pip,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pi0,pi0},
1765	{pi0,pro,pip,pim,pi0,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0},
1766	{pip,neu,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pip,pim,pim,pi0},
1767	{pi0,neu,pip,pip,pim,pi0,pi0,pi0}, {pi0,neu,pip,pi0,pi0,pi0,pi0,pi0}},
1768
1769	{{pi0,neu,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pim,pim,pi0,pi0},
1770	{pi0,neu,pip,pim,pi0,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0},
1771	{pim,pro,pip,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pim,pi0},
1772	{pi0,pro,pip,pim,pim,pi0,pi0,pi0}, {pi0,pro,pim,pi0,pi0,pi0,pi0,pi0}}};
1773
1774	const G4int T11_9bfsData[2][9][9] =
1775	{{{pip,pro,pip,pip,pip,pim,pim,pim,pim}, {pi0,pro,pip,pip,pip,pim,pim,pim,pi0},
1776	{pi0,pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pi0,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
1777	{pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pi0,neu,pip,pip,pip,pip,pim,pim,pim},
1778	{pi0,neu,pip,pip,pip,pim,pim,pi0,pi0}, {pi0,neu,pip,pip,pim,pi0,pi0,pi0,pi0},
1779	{pi0,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
1780
1781	{{pim,neu,pip,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pip,pim,pim,pim,pi0},
1782	{pi0,neu,pip,pip,pim,pim,pi0,pi0,pi0}, {pi0,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
1783	{pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pi0,pro,pip,pip,pip,pim,pim,pim,pim},
1784	{pi0,pro,pip,pip,pim,pim,pim,pi0,pi0}, {pi0,pro,pip,pim,pim,pi0,pi0,pi0,pi0},
1785	{pi0,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
1786
1787	for (G4int i = 0; i < 2; i++) {
1788	for (G4int j = 0; j < 5; j++) {
1789	for (G4int k = 0; k < 2; k++) {
1790	T11_2bfs[i][j][k] = T11_2bfsData[i][j][k];
1791	}
1792	}
1793	for (G4int j = 0; j < 13; j++) {
1794	for (G4int k = 0; k < 3; k++) {
1795	T11_3bfs[i][j][k] = T11_3bfsData[i][j][k];
1796	}
1797	}
1798	for (G4int j = 0; j < 21; j++) {
1799	for (G4int k = 0; k < 4; k++) {
1800	T11_4bfs[i][j][k] = T11_4bfsData[i][j][k];
1801	}
1802	}
1803	for (G4int j = 0; j < 30; j++) {
1804	for (G4int k = 0; k < 5; k++) {
1805	T11_5bfs[i][j][k] = T11_5bfsData[i][j][k];
1806	}
1807	}
1808	for (G4int j = 0; j < 6; j++) {
1809	for (G4int k = 0; k < 6; k++) {
1810	T11_6bfs[i][j][k] = T11_6bfsData[i][j][k];
1811	}
1812	}
1813	for (G4int j = 0; j < 7; j++) {
1814	for (G4int k = 0; k < 7; k++) {
1815	T11_7bfs[i][j][k] = T11_7bfsData[i][j][k];
1816	}
1817	}
1818	for (G4int j = 0; j < 8; j++) {
1819	for (G4int k = 0; k < 8; k++) {
1820	T11_8bfs[i][j][k] = T11_8bfsData[i][j][k];
1821	}
1822	}
1823	for (G4int j = 0; j < 9; j++) {
1824	for (G4int k = 0; k < 9; k++) {
1825	T11_9bfs[i][j][k] = T11_9bfsData[i][j][k];
1826	}
1827	}
1828	}
1829
1830	//
1831	// Cross sections (in mb) for pi0 p -> 2-9 body final states
1832	//
1833	// first index: 0-4: channels for mult = 2
1834	// 5-17: channels for mult = 3
1835	// 18-38: channels for mult = 4
1836	// 39-68: channels for mult = 5
1837	// 69-74: channels for mult = 6
1838	// 75-81: channels for mult = 7
1839	// 82-89: channels for mult = 8
1840	// 90-98: channels for mult = 9
1841	//
1842	// second index: kinetic energy
1843	//
1844
1845	const G4float pizPCrossSectionsData[99][30] = {
1846	//
1847	// multiplicity 2 (5 channels)
1848	//
1849	// p pi0 (n pi0)
1850	// { 0.0, 1.15, 1.85, 2.6, 3.25, 4.4, 5.5, 9.0, 16.7, 35.5,
1851	// 72.0, 107.0, 72.35, 33.1, 19.51, 17.02, 11.27, 19.11, 15.35, 11.0,
1852	// 8.83, 7.4, 6.4, 5.75, 5.2, 4.7, 4.25, 4.0, 3.8, 3.65},
1853
1854	// p pi0 (n pi0) go back to original for now
1855	{ 1.73, 2.28, 2.44, 2.71, 3.02, 3.43, 3.92, 5.37, 9.96, 17.73,
1856	29.07, 50.23, 33.68, 16.69, 12.60, 11.89, 13.49, 15.24, 12.39, 9.59,
1857	7.92, 6.97, 6.13, 5.37, 4.82, 4.68, 4.54, 4.0, 3.8, 3.65},
1858
1859	// n pi+ (p pi-)
1860	// { 0.0, 2.4, 2.8, 3.3, 4.5, 5.7, 6.3, 9.0, 11.0, 17.0,
1861	// 30.0, 43.0, 30.0, 16.5, 11.0, 7.01, 4.31, 5.03, 2.05, 0.97,
1862	// 0.53, 0.3, 0.2, 0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01},
1863
1864	// n pi+ (p pi-) taken from pi- p -> pi0 n
1865	{ 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 6.3, 9.0, 11.0, 17.0,
1866	32.0, 48.0, 28.0, 14.5, 11.04, 8.99, 4.79, 5.02, 2.08, 0.94,
1867	0.5, 0.25, 0.15, 0.09, 0.06, 0.05, 0.04, 0.02, 0.02, 0.01},
1868
1869	// L K+ (L K0)
1870	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1871	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.65, 0.29, 0.17,
1872	0.14, 0.1, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.0},
1873
1874	// S0 K+ (S0 K0)
1875	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1876	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.28, 0.18, 0.12,
1877	0.10, 0.06, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0},
1878
1879	// S+ K0 (S- K+)
1880	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1881	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.18, 0.43, 0.22,
1882	0.12, 0.06, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0},
1883	//
1884	// multiplicity 3 (13 channels)
1885	//
1886	// p pi+ pi- (n pi+ pi-)
1887	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1888	0.0, 0.0, 0.03, 0.19, 0.73, 3.4, 7.01, 8.35, 8.9, 5.69,
1889	4.01, 2.7, 2.0, 1.30, 0.9, 0.68, 0.48, 0.34, 0.27, 0.19},
1890
1891	// p pi0 pi0 (n pi0 pi0)
1892	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1893	0.0, 0.0, 0.22, 0.74, 1.8, 2.7, 3.0, 2.52, 1.33, 0.69,
1894	0.44, 0.3, 0.2, 0.12, 0.07, 0.04, 0.0, 0.0, 0.0, 0.0},
1895
1896	// n pi+ pi0 (p pi- pi0)
1897	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1898	0.0, 0.0, 0.04, 0.5, 1.75, 3.6, 5.21, 5.28, 6.34, 5.15,
1899	4.01, 3.1, 2.2, 1.42, 1.05, 0.8, 0.57, 0.33, 0.26, 0.19},
1900
1901	// L K+ pi0 (L K0 pi0)
1902	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1903	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.13,
1904	0.14, 0.1, 0.08, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
1905
1906	// L K0 pi+ (L K+ pi-)
1907	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1908	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.16,
1909	0.15, 0.1, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
1910
1911	// S0 K+ pi0 (S0 K0 pi0)
1912	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1913	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09,
1914	0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
1915
1916	// S0 K0 pi+ (S0 K+ pi-)
1917	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1918	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.11,
1919	0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
1920
1921	// S+ K0 pi0 (S- K+ pi0)
1922	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1923	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09,
1924	0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0},
1925
1926	// S+ K+ pi- (S- K0 pi+)
1927	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1928	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09,
1929	0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.01},
1930
1931	// S- K+ pi+ (S+ K0 pi-)
1932	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1933	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.08,
1934	0.10, 0.05, 0.03, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
1935
1936	// p K+ K- (n K0 K0b)
1937	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1938	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.11,
1939	0.12, 0.08, 0.07, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
1940
1941	// p K0 K0b (n K+ K-)
1942	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1943	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.11,
1944	0.12, 0.08, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
1945
1946	// n K+ K0b (p K0 K-)
1947	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1948	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
1949	0.11, 0.22, 0.2, 0.09, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0},
1950	//
1951	// multiplicity 4 (21 channels)
1952	//
1953	// p pi+ pi- pi0 (n pi+ pi- pi0)
1954	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1955	0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.14, 0.44, 1.11, 1.88,
1956	2.05, 2.07, 1.75, 1.5, 1.3, 1.1, 0.95, 0.80, 0.72, 0.66},
1957
1958	// p pi0 pi0 pi0 (n pi0 pi0 pi0)
1959	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1960	0.0, 0.0, 0.0, 0.0, 0.1, 0.55, 1.0, 0.88, 0.72, 0.64,
1961	0.53, 0.4, 0.35, 0.29, 0.22, 0.14, 0.1, 0.06, 0.03, 0.0},
1962
1963	// n pi+ pi+ pi- (p pi+ pi- pi-)
1964	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1965	0.0, 0.0, 0.0, 0.0, 0.05, 0.17, 0.32, 0.77, 1.79, 2.25,
1966	1.82, 1.5, 1.3, 1.05, 0.9, 0.72, 0.63, 0.5, 0.43, 0.35},
1967
1968	// n pi+ pi0 pi0 (p pi- pi0 pi0)
1969	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1970	0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.09, 0.28, 0.77, 1.32,
1971	1.54, 1.37, 1.2, 1.02, 0.9, 0.77, 0.69, 0.6, 0.54, 0.48},
1972
1973	// L K+ pi+ pi- (L K0 pi+ pi-)
1974	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1975	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1976	0.03, 0.06, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02},
1977
1978	// S0 K+ pi+ pi- (S0 K0 pi+ pi-)
1979	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1980	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1981	0.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1982
1983	// L K+ pi0 pi0 (L K0 pi0 pi0)
1984	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1985	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
1986	0.03, 0.06, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02},
1987
1988	// S0 K+ pi0 pi0 (S0 K0 pi0 pi0)
1989	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1990	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1991	0.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
1992
1993	// S+ K0 pi0 pi0 (S- K+ pi0 pi0)
1994	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1995	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1996	0.02, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
1997
1998	// L K0 pi+ pi0 (L K+ pi- pi0)
1999	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2000	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
2001	0.04, 0.07, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02},
2002
2003	// S0 K0 pi+ pi0 (S0 K+ pi- pi0)
2004	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2005	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2006	0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
2007
2008	// S+ K0 pi+ pi- (S- K+ pi+ pi-)
2009	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2010	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2011	0.02, 0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
2012
2013	// S+ K+ pi- pi0 (S- K0 pi+ pi0)
2014	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2015	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2016	0.02, 0.05, 0.08, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
2017
2018	// S- K+ pi+ pi0 (S+ K0 pi- pi0)
2019	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2020	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2021	0.03, 0.11, 0.11, 0.09, 0.07, 0.04, 0.03, 0.03, 0.02, 0.02},
2022
2023	// p pi0 K+ K- (n pi0 K0 K0b)
2024	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2025	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2026	0.01, 0.1, 0.11, 0.1, 0.08, 0.07, 0.07, 0.06, 0.05, 0.04},
2027
2028	// p pi0 K0 K0b (n pi0 K+ K-)
2029	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2030	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2031	0.02, 0.05, 0.11, 0.1, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04},
2032
2033	// p pi+ K0 K- (n pi- K+ K0b)
2034	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2035	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2036	0.02, 0.05, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04},
2037
2038	// p pi- K+ K0b (n pi+ K0 K-)
2039	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2040	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2041	0.02, 0.05, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04},
2042
2043	// n pi+ K+ K- (p pi- K0 K0b)
2044	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2045	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2046	0.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
2047
2048	// n pi+ K0 K0b (p pi- K+ K-)
2049	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2050	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2051	0.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
2052
2053	// n pi0 K+ K0b (p pi0 K0 K-)
2054	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2055	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2056	0.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
2057	//
2058	// multiplicity 5 (30 channels)
2059	//
2060	// p pi+ pi+ pi- pi- (n pi+ pi+ pi- pi-)
2061	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2062	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.2, 1.45,
2063	2.41, 3.1, 2.7, 2.3, 2.0, 1.54, 1.3, 1.0, 0.83, 0.66},
2064
2065	// p pi+ pi- pi0 pi0 (n pi+ pi- pi0 pi0)
2066	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2067	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.07, 0.48,
2068	0.80, 1.1, 0.9, 0.76, 0.63, 0.51, 0.43, 0.33, 0.28, 0.22},
2069
2070	// p pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0)
2071	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2072	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.21,
2073	0.35, 0.5, 0.67, 0.55, 0.47, 0.38, 0.33, 0.23, 0.2, 0.17},
2074
2075	// n pi+ pi+ pi- pi0 (p pi+ pi- pi- pi0)
2076	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2077	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.1, 0.2, 0.41, 1.16,
2078	1.75, 1.8, 1.8, 1.38, 1.0, 0.76, 0.58, 0.42, 0.32, 0.24},
2079
2080	// n pi+ pi0 pi0 pi0 (p pi- pi0 pi0 pi0)
2081	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2082	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.21,
2083	0.42, 0.6, 0.72, 0.68, 0.57, 0.47, 0.38, 0.28, 0.24, 0.2},
2084
2085	// L K+ pi+ pi- pi0 (L K0 pi+ pi- pi0)
2086	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2087	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2088	0.0, 0.01, 0.03, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
2089
2090	// L K+ pi0 pi0 pi0 (L K0 pi0 pi0 pi0)
2091	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2092	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2093	0.0, 0.01, 0.03, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
2094
2095	// L K0 pi+ pi+ pi- (L K+ pi+ pi- pi-)
2096	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2097	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2098	0.0, 0.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
2099
2100	// L K0 pi+ pi0 pi0 (L K+ pi- pi0 pi0)
2101	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2102	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2103	0.0, 0.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
2104
2105	// S0 K+ pi+ pi- pi0 (S0 K0 pi+ pi- pi0)
2106	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2107	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2108	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
2109
2110	// S0 K+ pi0 pi0 pi0 (S0 K0 pi0 pi0 pi0)
2111	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2112	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2113	0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
2114
2115	// S0 K0 pi+ pi+ pi- (S0 K+ pi+ pi- pi-)
2116	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2117	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2118	0.0, 0.01, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
2119
2120	// S0 K0 pi+ pi0 pi0 (S0 K+ pi- pi0 pi0)
2121	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2122	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2123	0.0, 0.01, 0.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
2124
2125	// S+ K0 pi+ pi- pi0 (S- K+ pi+ pi- pi0)
2126	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2127	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2128	0.0, 0.0, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
2129
2130	// S+ K0 pi0 pi0 pi0 (S- K+ pi0 pi0 pi0)
2131	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2132	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2133	0.0, 0.0, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
2134
2135	// S+ K+ pi+ pi- pi- (S- K0 pi+ pi+ pi-)
2136	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2137	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2138	0.0, 0.0, 0.01, 0.02, 0.03, 0.02, 0.02, 0.01, 0.01, 0.02},
2139
2140	// S+ K+ pi- pi0 pi0 (S- K0 pi+ pi0 pi0)
2141	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2142	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2143	0.0, 0.01, 0.03, 0.05, 0.05, 0.03, 0.03, 0.02, 0.02, 0.02},
2144
2145	// S- K+ pi+ pi+ pi- (S+ K0 pi+ pi- pi-)
2146	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2147	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2148	0.0, 0.0, 0.01, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
2149
2150	// S- K+ pi+ pi0 pi0 (S+ K0 pi- pi0 pi0)
2151	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2152	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2153	0.0, 0.01, 0.02, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
2154
2155	// p K+ K- pi0 pi0 (n K0 K0b pi0 pi0)
2156	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2157	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2158	0.0, 0.02, 0.06, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
2159
2160	// p K0 K0b pi0 pi0 (n K+ K- pi0 pi0)
2161	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2162	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2163	0.0, 0.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
2164
2165	// p K+ K0b pi- pi0 (n K0 K- pi+ pi0)
2166	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2167	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2168	0.0, 0.01, 0.03, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
2169
2170	// p K0 K- pi+ pi0 (n K+ K0b pi- pi0)
2171	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2172	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2173	0.0, 0.01, 0.02, 0.03, 0.04, 0.03, 0.03, 0.02, 0.02, 0.03},
2174
2175	// p K+ K- pi+ pi- (n K0 K0b pi+ pi-)
2176	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2177	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2178	0.0, 0.02, 0.06, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
2179
2180	// p K0 K0b pi+ pi- (n K+ K- pi+ pi-)
2181	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2182	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2183	0.0, 0.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
2184
2185	// n K+ K- pi+ pi0 (p K0 K0b pi- pi0)
2186	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2187	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2188	0.0, 0.02, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03},
2189
2190	// n K0 K0b pi+ pi0 (p K+ K- pi- pi0)
2191	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2192	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2193	0.0, 0.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
2194
2195	// n K+ K0b pi0 pi0 (p K0 K- pi0 pi0)
2196	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2197	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2198	0.0, 0.01, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
2199
2200	// n K0 K- pi+ pi+ (p K+ K0b pi- pi-)
2201	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2202	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2203	0.0, 0.01, 0.06, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
2204
2205	// n K+ K0b pi+ pi- (p K0 K- pi+ pi-)
2206	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2207	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2208	0.0, 0.01, 0.06, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
2209	//
2210	// multiplicity 6 (6 channels)
2211	//
2212	// p pi+ pi+ pi- pi- pi0 (n pi+ pi+ pi- pi- pi0)
2213	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2214	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
2215	0.04, 0.08, 0.12, 0.15, 0.2, 0.19, 0.17, 0.14, 0.12, 0.1},
2216
2217	// p pi+ pi- pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0)
2218	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2219	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
2220	0.06, 0.1, 0.14, 0.16, 0.20, 0.19, 0.19, 0.16, 0.14, 0.13},
2221
2222	// p pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0 pi0)
2223	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2224	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
2225	0.03, 0.05, 0.07, 0.09, 0.11, 0.10, 0.10, 0.09, 0.08, 0.07},
2226
2227	// n pi+ pi+ pi+ pi- pi- (p pi+ pi+ pi- pi- pi-)
2228	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2229	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.07,
2230	0.12, 0.21, 0.28, 0.36, 0.4, 0.42, 0.39, 0.36, 0.31, 0.28},
2231
2232	// n pi+ pi+ pi- pi0 pi0 (p pi+ pi- pi- pi0 pi0)
2233	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2234	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
2235	0.07, 0.15, 0.2, 0.24, 0.26, 0.29, 0.27, 0.24, 0.21, 0.18},
2236
2237	// n pi+ pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 pi0)
2238	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2239	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.04,
2240	0.09, 0.15, 0.2, 0.25, 0.30, 0.29, 0.28, 0.25, 0.21, 0.19},
2241	//
2242	// multiplicity 7 (7 channels)
2243	//
2244	// p pi+ pi+ pi+ pi- pi- pi- (n pi+ pi+ pi+ pi- pi- pi-)
2245	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2246	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
2247	0.03, 0.15, 0.4, 0.66, 0.85, 0.82, 0.75, 0.7, 0.65, 0.60},
2248
2249	// p pi+ pi+ pi- pi- pi0 pi0 (n pi+ pi+ pi- pi- pi0 pi0)
2250	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2251	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
2252	0.03, 0.15, 0.4, 0.66, 0.85, 0.82, 0.75, 0.7, 0.65, 0.60},
2253
2254	// p pi+ pi- pi0 pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0 pi0)
2255	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2256	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2257	0.01, 0.07, 0.16, 0.28, 0.31, 0.29, 0.25, 0.20, 0.15, 0.10},
2258
2259	// p pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0 pi0 pi0)
2260	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2261	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2262	0.0, 0.02, 0.05, 0.10, 0.12, 0.12, 0.11, 0.10, 0.09, 0.08},
2263
2264	// n pi+ pi+ pi+ pi- pi- pi0 (p pi+ pi+ pi- pi- pi- pi0)
2265	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2266	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2267	0.01, 0.04, 0.08, 0.14, 0.27, 0.3, 0.27, 0.22, 0.2, 0.18},
2268
2269	// n pi+ pi+ pi- pi0 pi0 pi0 (p pi+ pi- pi- pi0 pi0 pi0)
2270	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2271	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
2272	0.01, 0.04, 0.08, 0.12, 0.16, 0.15, 0.13, 0.11, 0.09, 0.07},
2273
2274	// n pi+ pi0 pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 pi0 pi0)
2275	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2276	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2277	0.01, 0.03, 0.05, 0.08, 0.16, 0.19, 0.17, 0.16, 0.14, 0.13},
2278	//
2279	// multiplicity 8 (8 channels)
2280	//
2281	// p pi+ pi+ pi+ pi- pi- pi- pi0 (n pi+ pi+ pi+ pi- pi- pi- pi0)
2282	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2283	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2284	0.01, 0.02, 0.04, 0.07, 0.11, 0.13, 0.14, 0.13, 0.12, 0.11},
2285
2286	// p pi+ pi+ pi- pi- pi0 pi0 pi0 (n pi+ pi+ pi- pi- pi0 pi0 pi0)
2287	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2288	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2289	0.01, 0.01, 0.03, 0.05, 0.08, 0.10, 0.10, 0.09, 0.09, 0.08},
2290
2291	// p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0 pi0 pi0)
2292	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2293	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2294	0.0, 0.01, 0.02, 0.04, 0.06, 0.07, 0.07, 0.06, 0.06, 0.05},
2295
2296	// p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
2297	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2298	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2299	0.0, 0.01, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04, 0.04},
2300
2301	// n pi+ pi+ pi+ pi+ pi- pi- pi- (p pi+ pi+ pi+ pi- pi- pi- pi-)
2302	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2303	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2304	0.02, 0.04, 0.07, 0.12, 0.19, 0.26, 0.26, 0.24, 0.23, 0.21},
2305
2306	// n pi+ pi+ pi+ pi- pi- pi0 pi0 (p pi+ pi+ pi- pi- pi- pi0 pi0)
2307	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2308	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2309	0.02, 0.04, 0.07, 0.12, 0.19, 0.26, 0.25, 0.23, 0.23, 0.2},
2310
2311	// n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p pi+ pi- pi- pi0 pi0 pi0 pi0)
2312	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2313	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2314	0.01, 0.03, 0.05, 0.08, 0.13, 0.13, 0.12, 0.11, 0.09, 0.08},
2315
2316	// n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 pi0 pi0 pi0)
2317	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2318	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2319	0.01, 0.02, 0.04, 0.06, 0.10, 0.13, 0.12, 0.11, 0.11, 0.1},
2320
2321	//
2322	// multiplicity 9 (9 channels)
2323	//
2324	// p pi+ pi+ pi+ pi+ pi- pi- pi- pi- (n pi+ pi+ pi+ pi+ pi- pi- pi- pi-)
2325	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2326	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2327	0.0, 0.01, 0.02, 0.05, 0.11, 0.16, 0.22, 0.27, 0.27, 0.27},
2328
2329	// p pi+ pi+ pi+ pi- pi- pi- pi0 pi0 (n pi+ pi+ pi+ pi- pi- pi- pi0 pi0)
2330	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2331	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2332	0.0, 0.01, 0.02, 0.03, 0.07, 0.10, 0.13, 0.18, 0.18, 0.18},
2333
2334	// p pi+ pi+ pi- pi- pi0 pi0 pi0 pi0 (n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0)
2335	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2336	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2337	0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.09, 0.11, 0.11, 0.11},
2338
2339	// p pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0)
2340	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2341	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2342	0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.07, 0.07, 0.07},
2343
2344	// p pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0)
2345	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2346	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2347	0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.07, 0.07, 0.07},
2348
2349	// n pi+ pi+ pi+ pi+ pi- pi- pi- pi0 (p pi+ pi+ pi+ pi- pi- pi- pi- pi0)
2350	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2351	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2352	0.0, 0.0, 0.01, 0.03, 0.06, 0.11, 0.14, 0.16, 0.16, 0.16},
2353
2354	// n pi+ pi+ pi+ pi- pi- pi0 pi0 pi0 (p pi+ pi+ pi- pi- pi- pi0 pi0 pi0)
2355	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2356	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2357	0.0, 0.0, 0.01, 0.02, 0.05, 0.07, 0.08, 0.1, 0.1, 0.1},
2358
2359	// n pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0 (p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0)
2360	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2361	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2362	0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.06, 0.06, 0.06},
2363
2364	// n pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0)
2365	{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2366	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
2367	0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04}};
2368
2369	// Put array in class scope
2370	for (G4int i = 0; i < 99; i++) {
2371	for (G4int j = 0; j < 30; j++) {
2372	pizPCrossSections[i][j] = pizPCrossSectionsData[i][j];
2373	}
2374	}
2375
2376	}
2377	/* end of file */
2378

Note: See TracBrowser for help on using the repository browser.

Download in other formats: