Context Navigation

source: trunk/source/processes/hadronic/models/cascade/cascade/src/G4CascadeKminusPChannel.cc @ 1315

Last change on this file since 1315 was 1315, checked in by garnier, 14 years ago
update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08
File size: 38.2 KB

Line
1	//
2	// ********************************************************************
3	// * License and Disclaimer *
4	// * *
5	// * The Geant4 software is copyright of the Copyright Holders of *
6	// * the Geant4 Collaboration. It is provided under the terms and *
7	// * conditions of the Geant4 Software License, included in the file *
8	// * LICENSE and available at http://cern.ch/geant4/license . These *
9	// * include a list of copyright holders. *
10	// * *
11	// * Neither the authors of this software system, nor their employing *
12	// * institutes,nor the agencies providing financial support for this *
13	// * work make any representation or warranty, express or implied, *
14	// * regarding this software system or assume any liability for its *
15	// * use. Please see the license in the file LICENSE and URL above *
16	// * for the full disclaimer and the limitation of liability. *
17	// * *
18	// * This code implementation is the result of the scientific and *
19	// * technical work of the GEANT4 collaboration. *
20	// * By using, copying, modifying or distributing the software (or *
21	// * any work based on the software) you agree to acknowledge its *
22	// * use in resulting scientific publications, and indicate your *
23	// * acceptance of all terms of the Geant4 Software license. *
24	// ********************************************************************
25	//
26
27	#include "G4CascadeKminusPChannel.hh"
28
29	namespace {
30	// Outgoing particle types of a given multiplicity
31
32	static const G4int kmp2bfs[8][2] =
33	{{1, 13}, {2, 17}, {7, 21}, {5, 23},
34	{7, 25}, {3, 27}, {15, 29}, {11, 31} };
35
36	static const G4int kmp3bfs[20][3] =
37	{{1,7,13}, {1,5,17}, {2,3,13}, {2,7,17}, {7,7,21},
38	{3,5,21}, {11,13,21}, {15,17,21}, {5,7,23}, {13,15,23},
39	{7,7,25}, {3,5,25}, {11,13,25}, {15,17,25}, {3,7,27},
40	{11,17,27}, {7,15,29}, {5,11,29}, {3,15,31}, {7,11,31} };
41
42	static const G4int kmp4bfs[34][4] =
43	{{1,7,7,13}, {1,3,5,13}, {1,5,7,17}, {2,3,7,13},
44	{2,7,7,17}, {2,3,5,17}, {7,7,7,21}, {3,5,7,21},
45	{3,13,15,21}, {5,11,17,21}, {7,11,13,21}, {7,15,17,21},
46	{5,7,7,23}, {3,5,5,23}, {7,7,7,25}, {3,5,7,25},
47	{3,7,7,27}, {3,3,5,27}, {7,13,15,23}, {5,11,13,23},
48	{5,15,17,23}, {3,13,15,25}, {5,11,17,25}, {7,11,13,25},
49	{7,15,17,25}, {3,11,13,27}, {3,15,17,27}, {7,11,17,27},
50	{7,7,15,29}, {3,5,15,29}, {5,7,11,29}, {3,7,15,31},
51	{7,7,11,31}, {3,5,11,31} };
52
53	static const G4int kmp5bfs[48][5] =
54	{{1,7,7,7,13}, {1,3,5,7,13}, {1,5,7,7,17}, {1,3,5,5,17},
55	{2,3,7,7,13}, {2,3,3,5,13}, {2,7,7,7,17}, {2,3,5,7,17},
56	{7,7,7,7,21}, {3,5,7,7,21}, {3,3,5,5,21}, {3,7,13,15,21},
57	{5,7,11,17,21}, {7,7,11,13,21}, {3,5,11,13,21}, {7,7,15,17,21},
58	{3,5,15,17,21}, {5,7,7,7,23}, {3,5,5,7,23}, {7,7,13,15,23},
59	{3,5,13,15,23}, {5,7,11,13,23}, {5,7,15,17,23}, {5,5,11,17,23},
60	{7,7,7,7,25}, {3,5,7,7,25}, {3,3,5,5,25}, {3,7,13,15,25},
61	{5,7,11,17,25}, {7,7,11,13,25}, {3,5,11,13,25}, {7,7,15,17,25},
62	{3,5,15,17,25}, {3,7,7,7,27}, {3,3,5,7,27}, {3,7,11,13,27},
63	{3,7,15,17,27}, {7,7,11,17,27}, {3,5,11,17,27}, {3,3,13,15,27},
64	{7,7,7,15,29}, {3,5,7,15,29}, {5,7,7,11,29}, {3,5,5,11,29},
65	{7,7,7,11,31}, {3,5,7,11,31}, {3,7,7,15,31}, {3,3,5,15,31} };
66
67	static const G4int kmp6bfs[22][6] =
68	{{1,7,7,7,7,13}, {1,3,5,7,7,13}, {1,3,3,5,5,13}, {1,5,7,7,7,17},
69	{1,3,5,5,7,17}, {2,3,7,7,7,13}, {2,3,3,5,7,13}, {2,7,7,7,7,17},
70	{2,3,5,7,7,17}, {2,3,3,5,5,17}, {7,7,7,7,7,21}, {3,5,7,7,7,21},
71	{3,3,5,5,7,21}, {5,7,7,7,7,23}, {3,5,5,7,7,23}, {3,3,5,5,5,23},
72	{7,7,7,7,7,25}, {3,5,7,7,7,25}, {3,3,5,5,7,25}, {3,7,7,7,7,27},
73	{3,3,5,7,7,27}, {3,3,3,5,5,27} };
74
75	static const G4int kmp7bfs[16][7] =
76	{{1,7,7,7,7,7,13}, {1,3,5,7,7,7,13}, {1,3,3,5,5,7,13},
77	{1,5,7,7,7,7,17}, {1,3,5,5,7,7,17}, {1,3,3,5,5,5,17},
78	{2,3,7,7,7,7,13}, {2,3,3,5,7,7,13}, {2,3,3,3,5,5,13},
79	{2,7,7,7,7,7,17}, {2,3,5,7,7,7,17}, {2,3,3,5,5,7,17},
80	{7,7,7,7,7,7,21}, {3,5,7,7,7,7,21}, {3,3,5,5,7,7,21},
81	{3,3,3,5,5,5,21} };
82
83	// Cross sections for K- p -> 2-7 body final states
84	//
85	// first index: 0-7: channels for mult = 2
86	// 8-27: channels for mult = 3
87	// 28-61: channels for mult = 4
88	// 62-109: channels for mult = 5
89	// 110-131: channels for mult = 6
90	// 132-147: channels for mult = 7
91	//
92	// second index: kinetic energy
93	//
94	static const G4double kmpCrossSections[148][31] = {
95	//
96	// multiplicity 2 (8 channels)
97	//
98	// K- p
99	{400.0,110.0, 37.6, 29.0, 22.29,21.0, 21.0, 21.0, 21.0, 21.0,
100	21.0, 10.98, 7.66, 6.51, 5.41, 5.02, 4.63, 4.30, 4.00, 3.90,
101	3.80, 3.70, 3.60, 3.60, 3.50, 3.40, 3.30, 3.20, 3.20, 3.18, 3.00},
102
103	// K0bar n
104	{173.0, 36.75, 8.0, 9.75, 4.81, 4.61, 4.86, 5.12, 5.37, 5.63,
105	5.88, 2.78, 1.13, 0.75, 0.50, 0.36, 0.27, 0.22, 0.18, 0.15,
106	0.13, 0.11, 0.10, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.06, 0.03},
107
108	// L pi0
109	{325.0, 62.08, 5.37, 3.88, 3.13, 3.00, 3.00, 3.00, 3.00, 3.00,
110	3.00, 1.55, 0.62, 0.29, 0.17, 0.11, 0.07, 0.05, 0.04, 0.03,
111	0.03, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
112
113	// S+ pi-
114	{241.0, 50.17,11.41, 9.52, 6.54, 5.58, 5.06, 4.55, 4.03, 3.52,
115	3.01, 1.22, 0.69, 0.42, 0.29, 0.21, 0.15, 0.12, 0.10, 0.09,
116	0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02},
117
118	// S0 pi0
119	{243.0, 47.17, 5.02, 6.78, 2.80, 2.52, 2.30, 2.08, 1.86, 1.64,
120	1.42, 0.33, 0.09, 0.03, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0,
121	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
122
123	// S- pi+
124	{615.0,112.9, 7.40, 6.45, 3.46, 2.84, 2.65, 2.45, 2.26, 2.07,
125	1.87, 0.56, 0.19, 0.07, 0.03, 0.12, 0.02, 0.01, 0.0, 0.0,
126	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
127
128	// X0 K0
129	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
130	0.0, 0.06, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0,
131	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
132
133	// X- K+
134	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
135	0.0, 0.15, 0.11, 0.07, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0,
136	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
137
138	//
139	// multiplicity 3 (20 channels)
140	//
141	// p K- pi0
142	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.10, 0.22, 0.34, 0.46, 0.58,
143	0.70, 1.44, 1.51, 1.23, 1.02, 0.88, 0.76, 0.67, 0.60, 0.54,
144	0.49, 0.45, 0.41, 0.38, 0.36, 0.34, 0.32, 0.30, 0.29, 0.28, 0.19},
145
146	// p K0bar pi-
147	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.06, 0.13, 0.20, 0.28, 0.35,
148	0.42, 1.57, 1.89, 1.54, 1.23, 0.97, 0.78, 0.65, 0.54, 0.46,
149	0.40, 0.36, 0.33, 0.30, 0.28, 0.26, 0.25, 0.23, 0.21, 0.20, 0.14},
150
151	// n K- pi+
152	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.07, 0.15, 0.24, 0.32, 0.41,
153	0.49, 2.07, 2.49, 2.10, 1.74, 1.52, 1.32, 1.15, 1.00, 0.89,
154	0.80, 0.75, 0.70, 0.65, 0.60, 0.56, 0.53, 0.50, 0.48, 0.46, 0.31},
155
156	// n K0bar pi0
157	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.08, 0.18, 0.27, 0.37, 0.47,
158	0.56, 1.13, 1.07, 0.83, 0.80, 0.80, 0.76, 0.74, 0.72, 0.67,
159	0.61, 0.58, 0.53, 0.50, 0.47, 0.44, 0.42, 0.40, 0.38, 0.37, 0.31},
160
161	// L 2pi0
162	{ 0.0, 0.0, 0.24, 1.48, 1.10, 1.24, 1.40, 1.56, 1.73, 1.89,
163	2.06, 1.59, 0.86, 0.54, 0.46, 0.41, 0.31, 0.24, 0.19, 0.15,
164	0.12, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.02},
165
166	// L pi+ pi-
167	{ 0.0, 0.0, 0.44, 1.05, 1.42, 1.76, 2.07, 2.39, 2.70, 3.02,
168	3.33, 2.65, 1.62, 0.92, 0.62, 0.43, 0.31, 0.24, 0.19, 0.15,
169	0.12, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.02},
170
171	// L K+ K-
172	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
173	0.0, 0.0, 0.05, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.04,
174	0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01},
175
176	// L K0 K0bar
177	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
178	0.0, 0.0, 0.03, 0.05, 0.06, 0.05, 0.04, 0.04, 0.04, 0.05,
179	0.05, 0.05, 0.05, 0.05, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03, 0.01},
180
181	// S+ pi- pi0
182	{ 0.0, 0.0, 0.06, 0.10, 0.17, 0.27, 0.35, 0.43, 0.51, 0.60,
183	0.68, 1.06, 0.90, 0.67, 0.46, 0.33, 0.26, 0.20, 0.16, 0.13,
184	0.11, 0.09, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04, 0.02},
185
186	// S+ K- K0
187	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
188	0.0, 0.0, 0.0, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0,
189	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
190
191	// S0 2pi0
192	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.12, 0.14, 0.17, 0.19, 0.22,
193	0.24, 0.45, 0.37, 0.25, 0.18, 0.11, 0.08, 0.06, 0.04, 0.03,
194	0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0},
195
196	// S0 pi+ pi-
197	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.14, 0.19, 0.24, 0.28, 0.33,
198	0.38, 0.74, 0.61, 0.42, 0.29, 0.20, 0.14, 0.10, 0.07, 0.05,
199	0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
200
201	// S0 K+ K-
202	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
203	0.0, 0.0, 0.01, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.0,
204	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
205
206	// S0 K0 K0bar
207	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
208	0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02,
209	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.0},
210
211	// S- pi+ pi0
212	{ 0.0, 0.0, 0.0, 0.10, 0.10, 0.16, 0.24, 0.32, 0.39, 0.47,
213	0.54, 0.79, 0.72, 0.45, 0.25, 0.15, 0.09, 0.06, 0.04, 0.03,
214	0.02, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
215
216	// S- K+ K0bar
217	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
218	0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01,
219	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
220
221	// X0 K0 pi0
222	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
223	0.0, 0.0, 0.02, 0.03, 0.04, 0.04, 0.04, 0.05, 0.05, 0.04,
224	0.04, 0.02, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
225
226	// X0 K+ pi-
227	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
228	0.0, 0.0, 0.05, 0.06, 0.02, 0.0, 0.0, 0.0, 0.0, 0.0,
229	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
230
231	// X- K0 pi+
232	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
233	0.0, 0.0, 0.08, 0.08, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01,
234	0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
235
236	// X- K+ pi0
237	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
238	0.0, 0.0, 0.03, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01,
239	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
240
241	//
242	// multiplicity 4 (34 channels)
243	//
244	// p K- 2pi0
245	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
246	0.01, 0.11, 0.44, 0.72, 0.78, 0.72, 0.69, 0.66, 0.63, 0.61,
247	0.60, 0.58, 0.57, 0.55, 0.54, 0.53, 0.52, 0.50, 0.49, 0.48, 0.43},
248
249	// p K- pi+ pi-
250	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
251	0.01, 0.17, 0.73, 1.20, 1.30, 1.21, 1.15, 1.10, 1.05, 1.02,
252	1.00, 0.97, 0.95, 0.92, 0.90, 0.88, 0.86, 0.84, 0.82, 0.81, 0.71},
253
254	// p K0bar pi- pi0
255	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
256	0.01, 0.08, 0.46, 0.99, 1.10, 1.00, 0.94, 0.89, 0.85, 0.81,
257	0.78, 0.75, 0.72, 0.70, 0.68, 0.66, 0.64, 0.62, 0.61, 0.60, 0.51},
258
259	// n K- pi+ pi0
260	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
261	0.01, 0.11, 0.28, 0.52, 0.55, 0.50, 0.47, 0.45, 0.43, 0.40,
262	0.39, 0.38, 0.36, 0.35, 0.34, 0.33, 0.32, 0.31, 0.30, 0.30, 0.25},
263
264	// n K0bar 2pi0
265	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
266	0.0, 0.06, 0.28, 0.50, 0.63, 0.71, 0.79, 0.78, 0.72, 0.65,
267	0.59, 0.54, 0.49, 0.45, 0.41, 0.38, 0.36, 0.33, 0.31, 0.30, 0.19},
268
269	// n K0bar pi+ pi-
270	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
271	0.0, 0.10, 0.45, 0.83, 1.05, 1.07, 0.94, 0.83, 0.72, 0.65,
272	0.59, 0.54, 0.49, 0.45, 0.41, 0.38, 0.36, 0.33, 0.31, 0.30, 0.19},
273
274	// L 3pi0
275	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.08, 0.11, 0.13,
276	0.16, 1.13, 1.59, 1.33, 1.12, 1.10, 0.86, 0.70, 0.34, 0.28,
277	0.38, 0.32, 0.27, 0.23, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05, 0.02},
278
279	// L pi+ pi- pi0
280	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.09, 0.13, 0.18, 0.23,
281	0.27, 1.77, 2.19, 1.84, 1.44, 1.12, 0.86, 0.70, 0.56, 0.46,
282	0.38, 0.32, 0.27, 0.23, 0.19, 0.16, 0.14, 0.12, 0.10, 0.09, 0.03},
283
284	// L K- K0 pi+
285	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
286	0.0, 0.0, 0.0, 0.03, 0.02, 0.02, 0.03, 0.03, 0.03, 0.03,
287	0.03, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
288
289	// L K+ K0bar pi-
290	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
291	0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.03, 0.02, 0.02,
292	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01},
293
294	// L K+ K- pi0
295	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
296	0.0, 0.0, 0.0, 0.0, 0.03, 0.05, 0.07, 0.08, 0.08, 0.08,
297	0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.08, 0.08, 0.08, 0.08, 0.07},
298
299	// L K0 K0bar pi0
300	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
301	0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03, 0.03,
302	0.04, 0.04, 0.05, 0.05, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
303
304	// S+ pi- 2pi0
305	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.02,
306	0.02, 0.09, 0.16, 0.16, 0.14, 0.13, 0.11, 0.10, 0.09, 0.08,
307	0.08, 0.07, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.05, 0.05, 0.01},
308
309	// S+ pi+ 2pi-
310	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.03,
311	0.04, 0.14, 0.26, 0.26, 0.24, 0.21, 0.19, 0.17, 0.15, 0.14,
312	0.13, 0.12, 0.11, 0.11, 0.10, 0.10, 0.09, 0.09, 0.08, 0.07, 0.02},
313
314	// S0 3pi0
315	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.06,
316	0.08, 0.22, 0.29, 0.35, 0.41, 0.36, 0.25, 0.20, 0.17, 0.14,
317	0.11, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.01},
318
319	// S0 pi+ pi- pi0
320	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.06, 0.09, 0.11,
321	0.13, 0.36, 0.49, 0.59, 0.69, 0.59, 0.41, 0.33, 0.28, 0.23,
322	0.19, 0.16, 0.14, 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.05, 0.02},
323
324	// S- pi+ 2pi0
325	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
326	0.01, 0.06, 0.10, 0.12, 0.11, 0.09, 0.07, 0.07, 0.06, 0.05,
327	0.05, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
328
329	// S- 2pi+ pi-
330	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01,
331	0.01, 0.10, 0.18, 0.20, 0.18, 0.15, 0.12, 0.11, 0.10, 0.08,
332	0.08, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04, 0.03, 0.01},
333
334	// S+ K- K0 pi0
335	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
336	0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.02,
337	0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
338
339	// S+ K+ K- pi-
340	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
341	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.03, 0.03,
342	0.03, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01},
343
344	// S+ K0 K0bar pi-
345	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
346	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.04, 0.04,
347	0.04, 0.04, 0.04, 0.04, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
348
349	// S0 K- K0 pi+
350	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
351	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.04, 0.04,
352	0.04, 0.02, 0.02, 0.02, 0.02, 0.02, 0.0, 0.0, 0.0, 0.0, 0.0},
353
354	// S0 K+ K0bar pi-
355	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
356	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.01, 0.02,
357	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
358
359	// S0 K+ K- pi0
360	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
361	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.02, 0.02,
362	0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.02},
363
364	// S0 K0 K0bar pi0
365	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
366	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.02,
367	0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.01},
368
369	// S- K+ K- pi+
370	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
371	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.02, 0.02,
372	0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
373
374	// S- K0 K0bar pi+
375	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
376	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02,
377	0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.0},
378
379	// S- K+ K0bar pi0
380	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
381	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
382	0.01, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
383
384	// X0 K0 2pi0
385	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
386	0.0, 0.0, 0.0, 0.0, 0.02, 0.02, 0.03, 0.03, 0.04, 0.04,
387	0.04, 0.04, 0.04, 0.04, 0.04, 0.02, 0.02, 0.02, 0.01, 0.01, 0.0},
388
389	// X0 K0 pi+ pi-
390	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
391	0.0, 0.0, 0.0, 0.01, 0.03, 0.03, 0.03, 0.03, 0.05, 0.05,
392	0.06, 0.06, 0.06, 0.06, 0.06, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02},
393
394	// X0 K+ pi- pi0
395	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
396	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.02, 0.02, 0.02,
397	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0},
398
399	// X- K0 pi+ pi0
400	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
401	0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02,
402	0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
403
404	// X- K+ 2pi0
405	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
406	0.0, 0.0, 0.0, 0.0, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0,
407	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
408
409	// X- K+ pi+ pi-
410	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
411	0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.02, 0.02, 0.02,
412	0.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
413
414	//
415	// multiplicity 5 (48 channels)
416	//
417	// p K- 3pi0
418	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
419	0.0, 0.01, 0.06, 0.16, 0.27, 0.38, 0.48, 0.54, 0.55, 0.54,
420	0.53, 0.52, 0.52, 0.52, 0.51, 0.50, 0.50, 0.50, 0.49, 0.49, 0.47},
421
422	// p K- pi+ pi- pi0
423	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
424	0.0, 0.02, 0.10, 0.26, 0.45, 0.64, 0.80, 0.90, 0.92, 0.90,
425	0.89, 0.87, 0.87, 0.86, 0.85, 0.84, 0.83, 0.83, 0.82, 0.82, 0.79},
426
427	// p K0bar pi- 2pi0
428	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
429	0.0, 0.0, 0.01, 0.04, 0.09, 0.18, 0.27, 0.32, 0.32, 0.32,
430	0.32, 0.32, 0.32, 0.32, 0.31, 0.30, 0.30, 0.30, 0.29, 0.28, 0.24},
431
432	// p K0bar pi+ pi- pi-
433	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
434	0.0, 0.0, 0.02, 0.06, 0.15, 0.27, 0.32, 0.32, 0.32, 0.32,
435	0.32, 0.32, 0.32, 0.32, 0.31, 0.30, 0.30, 0.30, 0.29, 0.28, 0.24},
436
437	// n K- pi+ 2pi0
438	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
439	0.0, 0.0, 0.02, 0.05, 0.12, 0.19, 0.24, 0.25, 0.25, 0.26,
440	0.26, 0.27, 0.28, 0.28, 0.28, 0.29, 0.28, 0.28, 0.28, 0.27, 0.20},
441
442	// n K- 2pi+ pi-
443	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
444	0.0, 0.01, 0.03, 0.08, 0.20, 0.32, 0.39, 0.41, 0.42, 0.43,
445	0.44, 0.45, 0.46, 0.46, 0.47, 0.48, 0.47, 0.47, 0.46, 0.45, 0.33},
446
447	// n K0bar 3pi0
448	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
449	0.0, 0.0, 0.02, 0.05, 0.12, 0.21, 0.34, 0.42, 0.44, 0.47,
450	0.48, 0.49, 0.51, 0.51, 0.52, 0.52, 0.53, 0.53, 0.53, 0.52, 0.47},
451
452	// n K0bar pi+ pi- pi0
453	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
454	0.0, 0.01, 0.03, 0.08, 0.20, 0.36, 0.57, 0.69, 0.74, 0.77,
455	0.79, 0.82, 0.83, 0.84, 0.86, 0.87, 0.88, 0.89, 0.88, 0.87, 0.81},
456
457	// L 4pi0
458	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
459	0.0, 0.01, 0.04, 0.06, 0.09, 0.15, 0.14, 0.13, 0.13, 0.12,
460	0.11, 0.11, 0.10, 0.10, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.01},
461
462	// L pi+ pi- 2pi0
463	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
464	0.0, 0.02, 0.05, 0.10, 0.13, 0.15, 0.14, 0.13, 0.13, 0.12,
465	0.11, 0.11, 0.10, 0.10, 0.09, 0.09, 0.08, 0.07, 0.06, 0.06, 0.03},
466
467	// L 2pi+ 2pi-
468	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
469	0.0, 0.03, 0.09, 0.16, 0.17, 0.15, 0.14, 0.13, 0.13, 0.12,
470	0.11, 0.11, 0.10, 0.10, 0.09, 0.09, 0.08, 0.07, 0.06, 0.06, 0.03},
471
472	// L K- K0 pi+ pi0
473	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
474	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.03,
475	0.04, 0.04, 0.05, 0.05, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
476
477	// L K+ K0bar pi- pi0
478	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
479	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03,
480	0.04, 0.06, 0.07, 0.07, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04},
481
482	// L K+ K- 2pi0
483	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
484	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02,
485	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
486
487	// L K+ K- pi+ pi-
488	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
489	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02,
490	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
491
492	// L K0 K0bar 2pi0
493	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
494	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.07, 0.05, 0.06,
495	0.07, 0.08, 0.08, 0.07, 0.07, 0.07, 0.06, 0.06, 0.06, 0.06, 0.03},
496
497	// L K0 K0bar pi+ pi-
498	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
499	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
500	0.07, 0.08, 0.08, 0.07, 0.07, 0.07, 0.06, 0.06, 0.06, 0.06, 0.03},
501
502	// S+ pi- 3pi0
503	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
504	0.0, 0.0, 0.04, 0.11, 0.14, 0.14, 0.14, 0.14, 0.15, 0.15,
505	0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.14, 0.14, 0.12},
506
507	// S+ pi+ 2pi- pi0
508	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
509	0.0, 0.01, 0.07, 0.19, 0.24, 0.24, 0.24, 0.24, 0.25, 0.25,
510	0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.24, 0.24, 0.20},
511
512	// S+ K- K0 2pi0
513	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
514	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
515	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
516
517	// S+ K- K0 pi+ pi-
518	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
519	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01,
520	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.02},
521
522	// S+ K+ K- pi- pi0
523	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
524	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01,
525	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
526
527	// S+ K0 K0bar pi- pi0
528	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
529	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01,
530	0.02, 0.04, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
531
532	// S+ K+ K0bar 2pi-
533	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
534	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
535	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
536
537	// S0 4pi0
538	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
539	0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.02, 0.02, 0.03, 0.02,
540	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
541
542	// S0 pi+ pi- 2pi0
543	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
544	0.0, 0.0, 0.01, 0.02, 0.07, 0.04, 0.04, 0.04, 0.04, 0.04,
545	0.04, 0.04, 0.03, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
546
547	// S0 2pi+ 2pi-
548	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
549	0.0, 0.0, 0.02, 0.03, 0.03, 0.07, 0.07, 0.07, 0.07, 0.06,
550	0.06, 0.06, 0.05, 0.05, 0.05, 0.05, 0.04, 0.03, 0.03, 0.03, 0.02},
551
552	// S0 K- K0 pi+ pi0
553	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
554	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.04,
555	0.04, 0.06, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03},
556
557	// S0 K+ K0bar pi- pi0
558	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
559	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.04,
560	0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
561
562	// S0 K+ K- 2pi0
563	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
564	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01,
565	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
566
567	// S0 K+ K- pi+ pi-
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.0, 0.01, 0.01,
570	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
571
572	// S0 K0 K0bar 2pi0
573	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
574	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.02,
575	0.04, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.03, 0.03, 0.03, 0.02},
576
577	// S0 K0 K0bar pi+ pi-
578	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
579	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.04,
580	0.06, 0.07, 0.07, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.05, 0.04},
581
582	// S- pi+ 3pi0
583	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
584	0.0, 0.0, 0.04, 0.09, 0.11, 0.11, 0.11, 0.11, 0.11, 0.11,
585	0.11, 0.11, 0.12, 0.12, 0.12, 0.12, 0.11, 0.11, 0.11, 0.11, 0.08},
586
587	// S- 2pi+ pi- pi0
588	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
589	0.0, 0.01, 0.06, 0.14, 0.19, 0.19, 0.19, 0.19, 0.19, 0.19,
590	0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.19, 0.19, 0.18, 0.18, 0.14},
591
592	// S- K+ K- pi+ pi0
593	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
594	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
595	0.02, 0.03, 0.04, 0.05, 0.06, 0.06, 0.07, 0.08, 0.09, 0.09, 0.09},
596
597	// S- K0 K0bar pi+ pi0
598	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
599	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
600	0.02, 0.03, 0.04, 0.05, 0.06, 0.06, 0.07, 0.08, 0.09, 0.09, 0.09},
601
602	// S- K+ K0bar 2pi0
603	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
604	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
605	0.01, 0.01, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
606
607	// S- K+ K0bar pi+ pi-
608	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
609	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
610	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.02},
611
612	// S- K- K0 2pi+
613	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
614	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
615	0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
616
617	// X0 K0 3pi0
618	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
619	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.04,
620	0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.03, 0.02, 0.02, 0.02, 0.02},
621
622	// X0 K0 pi+ pi- pi0
623	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
624	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.03, 0.04, 0.05,
625	0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.05, 0.03, 0.05, 0.05, 0.04},
626
627	// X0 K+ pi- 2pi0
628	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
629	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.02,
630	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
631
632	// X0 K+ pi+ 2pi- pi0
633	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
634	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
635	0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03},
636
637	// X- K+ 3pi0
638	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
639	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.02,
640	0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02, 0.02},
641
642	// X- K+ pi+ pi- pi0
643	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
644	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
645	0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03},
646
647	// X- K0 pi+ 2pi0
648	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
649	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.04,
650	0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
651
652	// X- K0 2pi+ pi-
653	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
654	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.03, 0.04, 0.05,
655	0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
656
657	//
658	// multiplicity 6 (22 channels)
659	//
660	// p K- 4pi0
661	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
662	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.03, 0.03, 0.04,
663	0.04, 0.04, 0.05, 0.05, 0.06, 0.06, 0.07, 0.07, 0.07, 0.07, 0.08},
664
665	// p K- pi+ pi- 2pi0
666	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
667	0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.05, 0.05, 0.06,
668	0.07, 0.07, 0.08, 0.08, 0.10, 0.10, 0.11, 0.11, 0.12, 0.12, 0.13},
669
670	// p K- 2pi+ 2pi-
671	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
672	0.0, 0.0, 0.0, 0.01, 0.03, 0.05, 0.06, 0.08, 0.09, 0.10,
673	0.11, 0.12, 0.13, 0.14, 0.16, 0.17, 0.18, 0.19, 0.20, 0.20, 0.21},
674
675	// p K0bar pi- 3pi0
676	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
677	0.0, 0.0, 0.0, 0.0, 0.0, 0.06, 0.15, 0.24, 0.29, 0.34,
678	0.39, 0.40, 0.41, 0.43, 0.44, 0.46, 0.48, 0.50, 0.53, 0.53, 0.57},
679
680	// p K0bar pi+ 2pi- pi0
681	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
682	0.0, 0.0, 0.0, 0.0, 0.01, 0.08, 0.17, 0.24, 0.29, 0.34,
683	0.39, 0.40, 0.41, 0.43, 0.44, 0.46, 0.48, 0.50, 0.53, 0.53, 0.57},
684
685	// n K- pi+ 3pi0
686	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
687	0.0, 0.0, 0.0, 0.0, 0.02, 0.06, 0.09, 0.12, 0.14, 0.17,
688	0.19, 0.20, 0.20, 0.22, 0.22, 0.23, 0.24, 0.25, 0.26, 0.26, 0.28},
689
690	// n K- 2pi+ pi- pi0
691	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
692	0.0, 0.0, 0.0, 0.0, 0.03, 0.09, 0.15, 0.20, 0.24, 0.28,
693	0.32, 0.33, 0.34, 0.36, 0.36, 0.38, 0.40, 0.41, 0.44, 0.44, 0.47},
694
695	// n K0bar 4pi0
696	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
697	0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.05, 0.05, 0.06,
698	0.07, 0.07, 0.09, 0.11, 0.11, 0.14, 0.16, 0.16, 0.18, 0.18, 0.17},
699
700	// n K0bar pi+ pi- 2pi0
701	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
702	0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.06, 0.09, 0.09, 0.10,
703	0.11, 0.12, 0.14, 0.16, 0.17, 0.19, 0.21, 0.23, 0.25, 0.25, 0.20},
704
705	// n K0bar 2pi+ 2pi-
706	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
707	0.0, 0.0, 0.0, 0.0, 0.0, 0.06, 0.07, 0.08, 0.09, 0.10,
708	0.11, 0.12, 0.14, 0.16, 0.17, 0.19, 0.21, 0.23, 0.25, 0.25, 0.20},
709
710	// L 5pi0
711	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
712	0.0, 0.0, 0.01, 0.04, 0.07, 0.15, 0.20, 0.23, 0.27, 0.31,
713	0.34, 0.34, 0.33, 0.33, 0.32, 0.31, 0.26, 0.21, 0.16, 0.16, 0.13},
714
715	// L pi+ pi- 3pi0
716	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
717	0.0, 0.0, 0.02, 0.05, 0.09, 0.16, 0.20, 0.23, 0.27, 0.31,
718	0.34, 0.34, 0.33, 0.33, 0.32, 0.31, 0.26, 0.21, 0.16, 0.16, 0.13},
719
720	// L 2pi+ 2pi- pi0
721	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
722	0.0, 0.0, 0.03, 0.11, 0.12, 0.16, 0.20, 0.23, 0.27, 0.31,
723	0.34, 0.34, 0.33, 0.33, 0.32, 0.31, 0.26, 0.21, 0.16, 0.16, 0.13},
724
725	// S+ pi- 4pi0
726	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
727	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.07, 0.09, 0.09,
728	0.09, 0.09, 0.09, 0.09, 0.09, 0.08, 0.08, 0.08, 0.08, 0.08, 0.07},
729
730	// S+ pi+ 2pi- 2pi0
731	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
732	0.0, 0.0, 0.0, 0.0, 0.02, 0.04, 0.06, 0.12, 0.15, 0.15,
733	0.15, 0.15, 0.15, 0.14, 0.14, 0.14, 0.14, 0.13, 0.13, 0.13, 0.11},
734
735	// S+ 2pi+ 3pi-
736	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
737	0.0, 0.0, 0.0, 0.0, 0.03, 0.06, 0.10, 0.20, 0.25, 0.25,
738	0.25, 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.22, 0.21, 0.21, 0.19},
739
740	// S0 5pi0
741	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
742	0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.04, 0.05, 0.06, 0.07,
743	0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03},
744
745	// S0 pi+ pi- 3pi0
746	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
747	0.0, 0.0, 0.0, 0.0, 0.02, 0.05, 0.07, 0.08, 0.10, 0.11,
748	0.12, 0.12, 0.12, 0.12, 0.11, 0.11, 0.10, 0.08, 0.06, 0.06, 0.05},
749
750	// S0 2pi+ 2pi- pi0
751	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
752	0.0, 0.0, 0.0, 0.0, 0.03, 0.08, 0.11, 0.14, 0.16, 0.19,
753	0.20, 0.20, 0.20, 0.20, 0.19, 0.19, 0.16, 0.13, 0.10, 0.10, 0.08},
754
755	// S- pi+ 4pi0
756	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
757	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.07, 0.07,
758	0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.06, 0.06, 0.06, 0.06, 0.06},
759
760	// S- 2pi+ pi- 2pi0
761	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
762	0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.06, 0.09, 0.11, 0.11,
763	0.12, 0.12, 0.12, 0.12, 0.11, 0.11, 0.11, 0.11, 0.10, 0.10, 0.10},
764
765	// S- 3pi+ 2pi-
766	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
767	0.0, 0.0, 0.0, 0.0, 0.02, 0.06, 0.10, 0.15, 0.19, 0.19,
768	0.20, 0.20, 0.20, 0.20, 0.19, 0.19, 0.18, 0.18, 0.17, 0.17, 0.16},
769
770	//
771	// multiplicity 7 (16 channels)
772	//
773	// p K- 5pi0
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.01, 0.01, 0.02, 0.03, 0.05,
776	0.06, 0.07, 0.08, 0.09, 0.10, 0.12, 0.13, 0.14, 0.16, 0.17, 0.21},
777
778	// p K- pi+ pi- 3pi0
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.01, 0.01, 0.02, 0.02, 0.05, 0.08,
781	0.10, 0.11, 0.14, 0.16, 0.17, 0.20, 0.22, 0.24, 0.26, 0.28, 0.35},
782
783	// p K- 2pi+ 2pi- pi0
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.01, 0.02, 0.03, 0.04, 0.09, 0.13,
786	0.16, 0.19, 0.23, 0.26, 0.29, 0.33, 0.36, 0.40, 0.43, 0.45, 0.58},
787
788	// p K0bar pi- 4pi0
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.02, 0.02, 0.02,
791	0.02, 0.02, 0.02, 0.03, 0.04, 0.04, 0.04, 0.06, 0.06, 0.06, 0.07},
792
793	// p K0bar pi+ 2pi- 2pi0
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.01, 0.01, 0.02, 0.02, 0.03,
796	0.03, 0.04, 0.04, 0.05, 0.06, 0.06, 0.07, 0.08, 0.09, 0.09, 0.09},
797
798	// p K0bar 2pi+ 3pi-
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.01, 0.01, 0.02, 0.02, 0.03,
801	0.03, 0.04, 0.04, 0.05, 0.06, 0.06, 0.07, 0.08, 0.09, 0.09, 0.09},
802
803	// n K- pi+ 4pi0
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.01, 0.02, 0.02, 0.02, 0.03,
806	0.03, 0.04, 0.04, 0.05, 0.05, 0.05, 0.06, 0.06, 0.07, 0.07, 0.07},
807
808	// n K- 2pi+ pi- 2pi0
809	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
810	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.03, 0.04, 0.05,
811	0.05, 0.06, 0.07, 0.08, 0.08, 0.09, 0.10, 0.10, 0.11, 0.11, 0.11},
812
813	// n K- 3pi+ 2pi-
814	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
815	0.0, 0.0, 0.0, 0.0, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08,
816	0.09, 0.10, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.18},
817
818	// n K0bar 5pi0
819	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
820	0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.02, 0.05, 0.09,
821	0.11, 0.13, 0.15, 0.17, 0.19, 0.23, 0.26, 0.28, 0.32, 0.34, 0.45},
822
823	// n K0bar pi+ pi- 3pi0
824	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
825	0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.03, 0.03, 0.09, 0.14,
826	0.18, 0.20, 0.26, 0.30, 0.32, 0.39, 0.45, 0.49, 0.53, 0.53, 0.55},
827
828	// n K0bar 2pi+ 2pi- pi0
829	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
830	0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.07, 0.16, 0.23,
831	0.29, 0.35, 0.42, 0.48, 0.55, 0.64, 0.73, 0.81, 0.83, 0.83, 0.87},
832
833	// L 6pi0
834	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
835	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
836	0.01, 0.01, 0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01},
837
838	// L pi+ pi- 4pi0
839	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
840	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
841	0.01, 0.01, 0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01},
842
843	// L 2pi+ 2pi- 2pi0
844	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
845	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01,
846	0.02, 0.02, 0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01},
847
848	// L 3pi+ 3pi-
849	{ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
850	0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01,
851	0.02, 0.02, 0.02, 0.02, 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01} };
852	}
853
854	G4CascadeKminusPChannelData::data_t
855	G4CascadeKminusPChannelData::data(kmp2bfs, kmp3bfs, kmp4bfs,
856	kmp5bfs, kmp6bfs, kmp7bfs,
857	kmpCrossSections);

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

Download in other formats: