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

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

import all except CVS

File size: 23.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 "G4CascadeXiMinusNChannel.hh"
28
29
30G4CascadeXiMinusNChannel::G4CascadeXiMinusNChannel()
31 :G4CascadeChannel()
32{
33 G4int i, k, m;
34 G4int start, stop;
35
36 // Initialize multiplicity array
37
38 for (m = 0; m < 6; m++) {
39 start = xmnindex[m][0];
40 stop = xmnindex[m][1];
41 for (k = 0; k < 31; k++) {
42 xmnMultiplicities[m][k] = 0.0;
43 for (i = start; i < stop; i++) xmnMultiplicities[m][k] += xmnCrossSections[i][k];
44 }
45 }
46
47 // Initialize total cross section array
48
49 for (k = 0; k < 31; k++) {
50 xmntot[k] = 0.0;
51 for (m = 0; m < 6; m++) xmntot[k] += xmnMultiplicities[m][k];
52 }
53
54}
55
56
57G4CascadeXiMinusNChannel::~G4CascadeXiMinusNChannel()
58{;}
59
60
61G4double G4CascadeXiMinusNChannel::getCrossSection(G4double ke) const
62{
63 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
64 G4int k = epair.first;
65 G4double fraction = epair.second;
66
67 return xmntot[k] + fraction*(xmntot[k+1] - xmntot[k]);
68}
69
70
71G4int G4CascadeXiMinusNChannel::getMultiplicity(G4double ke) const
72{
73 G4double multint(0.);
74 std::vector<G4double> sigma;
75
76 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
77 G4int k = epair.first;
78 G4double fraction = epair.second;
79
80 for(G4int m = 0; m < 6; m++) {
81 multint = xmnMultiplicities[m][k]
82 + fraction*(xmnMultiplicities[m][k+1] - xmnMultiplicities[m][k]);
83 sigma.push_back(multint);
84 }
85
86 return sampleFlat(sigma);
87}
88
89
90std::vector<G4int>
91G4CascadeXiMinusNChannel::getOutgoingParticleTypes(G4int mult, G4double ke) const
92{
93 G4int i;
94 G4double sigint(0.);
95 std::vector<G4double> sigma;
96
97 std::pair<G4int, G4double> epair = interpolateEnergy(ke);
98 G4int k = epair.first;
99 G4double fraction = epair.second;
100
101 G4int start = xmnindex[mult-2][0];
102 G4int stop = xmnindex[mult-2][1];
103
104 for(i = start; i < stop; i++) {
105 sigint = xmnCrossSections[i][k]
106 + fraction*(xmnCrossSections[i][k+1] - xmnCrossSections[i][k]);
107 sigma.push_back(sigint);
108 }
109
110 G4int channel = sampleFlat(sigma);
111
112 std::vector<G4int> kinds;
113
114 if (mult == 2) {
115 for(i = 0; i < mult; i++) kinds.push_back(xmn2bfs[channel][i]);
116 } else if (mult == 3) {
117 for(i = 0; i < mult; i++) kinds.push_back(xmn3bfs[channel][i]);
118 } else if (mult == 4) {
119 for(i = 0; i < mult; i++) kinds.push_back(xmn4bfs[channel][i]);
120 } else if (mult == 5) {
121 for(i = 0; i < mult; i++) kinds.push_back(xmn5bfs[channel][i]);
122 } else if (mult == 6) {
123 for(i = 0; i < mult; i++) kinds.push_back(xmn6bfs[channel][i]);
124 } else if (mult == 7) {
125 for(i = 0; i < mult; i++) kinds.push_back(xmn7bfs[channel][i]);
126 } else {
127 G4cout << " Illegal multiplicity " << G4endl;
128 }
129
130 return kinds;
131}
132
133// Total cross section as a function of kinetic energy
134G4double G4CascadeXiMinusNChannel::xmntot[31];
135
136// Multiplicities as a function of kinetic energy
137G4double G4CascadeXiMinusNChannel::xmnMultiplicities[6][31];
138
139
140const G4int G4CascadeXiMinusNChannel::xmnindex[6][2] =
141 {{0, 3}, {3, 21}, {21,74}, {74,76}, {76,78}, {78,80}};
142
143// Outgoing particle types of a given multiplicity
144
145const G4int G4CascadeXiMinusNChannel::xmn2bfs[3][2] =
146 {{2, 31}, {21,27}, {25,27}};
147
148const G4int G4CascadeXiMinusNChannel::xmn3bfs[18][3] =
149 {{2,13,21}, {5,21,21}, {1,13,27}, {2,13,25}, {7,21,27}, {5,23,27},
150 {7,25,27}, {3,27,27}, {2,17,27}, {5,21,25}, {5,25,25}, {2,5,29},
151 {2,7,31}, {15,21,31},{11,27,31},{15,27,29},{1,5,31}, {15,25,31}};
152
153const G4int G4CascadeXiMinusNChannel::xmn4bfs[53][4] =
154 {{1,2,13,13}, {1,5,13,21}, {1,5,13,25}, {1,7,13,27}, {1,5,17,27},
155 {1,5,5,29}, {1,5,7,31}, {1,13,15,31}, {2,2,13,17}, {2,7,13,21},
156 {2,5,17,21}, {2,7,13,25}, {2,5,17,25}, {2,5,13,23}, {2,7,17,27},
157 {2,3,13,27}, {2,7,7,31}, {2,3,5,31}, {2,15,17,31}, {2,11,13,31},
158 {2,5,7,29}, {2,13,15,29}, {5,7,21,21}, {13,15,21,21},{5,7,21,25},
159 {13,15,21,25},{5,5,21,23}, {15,17,21,27},{11,13,21,27},{7,7,21,27},
160 {3,5,21,27}, {5,15,21,29}, {7,15,21,31}, {5,11,21,31}, {7,7,25,27},
161 {3,5,25,27}, {11,13,25,27},{15,17,25,27},{5,7,25,25}, {13,15,25,25},
162 {5,5,23,25}, {5,7,23,27}, {3,7,27,27}, {13,15,23,27},{11,17,27,27},
163 {5,15,25,29}, {7,15,25,31}, {5,11,25,31}, {5,15,23,31}, {5,11,27,29},
164 {3,15,27,31}, {7,11,27,31}, {7,15,27,29}};
165
166const G4int G4CascadeXiMinusNChannel::xmn5bfs[2][5] =
167 {{2,7,7,7,31}, {2,3,5,7,31}};
168
169const G4int G4CascadeXiMinusNChannel::xmn6bfs[2][6] =
170 {{2,7,7,7,7,31}, {2,3,5,7,7,31}};
171
172const G4int G4CascadeXiMinusNChannel::xmn7bfs[2][7] =
173 {{2,7,7,7,7,7,31}, {2,3,5,7,7,7,31}};
174
175// Cross sections for X- n -> 2-7 body final states
176//
177// first index: 0-2: channels for mult = 2
178// 3-20: channels for mult = 3
179// 21-73: channels for mult = 4
180// 74-75: channels for mult = 5
181// 76-77: channels for mult = 6
182// 78-79: channels for mult = 7
183//
184// second index: kinetic energy
185//
186const G4float G4CascadeXiMinusNChannel::xmnCrossSections[80][31] = {
187 //
188 // multiplicity 2 (3 channels)
189 //
190 // X- n
191 {16.00, 9.60, 7.00, 6.80, 6.60, 6.40, 6.20, 6.00, 5.82, 5.64,
192 5.45, 4.80, 4.58, 4.20, 3.80, 3.50, 3.20, 3.10, 2.95, 2.80,
193 2.65, 2.55, 2.45, 2.35, 2.25, 2.15, 2.10, 2.05, 2.04, 2.03, 1.95},
194
195 // L S-
196 { 0.0, 0.0, 0.0, 1.80, 3.60, 5.40, 7.20, 9.00, 9.82,10.64,
197 11.45,14.40,12.97,11.31, 8.38, 5.45, 2.50, 2.10, 1.90, 1.70,
198 1.50, 1.40, 1.30, 1.20, 1.10, 1.05, 1.00, 0.95, 0.94, 0.93, 0.85},
199
200 // S0 S-
201 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.55, 1.09,
202 1.64, 3.60, 4.32, 3.77, 2.80, 2.00, 1.50, 1.40, 1.30, 1.20,
203 1.10, 1.00, 0.95, 0.90, 0.85, 0.80, 0.75, 0.70, 0.69, 0.68, 0.60},
204
205 //
206 // multiplicity 3 (18 channels)
207 //
208 // L n K-
209 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09, 0.18,
210 0.27, 0.60, 0.71, 0.90, 1.10, 1.00, 0.90, 0.80, 0.68, 0.56,
211 0.40, 0.30, 0.20, 0.10, 0.07, 0.05, 0.03, 0.01, 0.01, 0.01, 0.0},
212
213 // L L pi-
214 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.46, 0.91,
215 1.36, 3.00, 1.93, 1.80, 1.70, 1.50, 1.30, 1.10, 1.00, 0.90,
216 0.75, 0.60, 0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.14, 0.13, 0.10},
217
218 // S- p K-
219 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
220 0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
221 0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
222
223 // S0 n K-
224 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
225 0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
226 0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
227
228 // S- L pi0
229 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.55, 1.09,
230 1.64, 3.60, 3.36, 3.10, 3.00, 2.80, 2.60, 2.30, 2.00, 1.80,
231 1.60, 1.50, 1.40, 1.30, 1.20, 1.15, 1.10, 1.05, 1.04, 1.03, 0.90},
232
233 // S+ S- pi-
234 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
235 0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
236 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
237
238 // S- S0 pi0
239 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
240 0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
241 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
242
243 // S- S- pi+
244 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
245 0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
246 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
247
248 // S- n K0bar
249 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
250 0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
251 0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
252
253 // S0 L pi-
254 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
255 0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
256 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
257
258 // S0 S0 pi-
259 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
260 0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
261 0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
262
263 // X0 n pi-
264 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.14, 0.27,
265 0.41, 0.90, 1.02, 1.50, 1.40, 1.10, 1.00, 0.90, 0.80, 0.65,
266 0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.12, 0.10, 0.09, 0.09, 0.05},
267
268 // X- n pi0
269 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.14, 0.27,
270 0.41, 0.90, 1.02, 1.50, 1.40, 1.10, 1.00, 0.90, 0.80, 0.65,
271 0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.12, 0.10, 0.09, 0.09, 0.05},
272
273 // X- L K0
274 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
275 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
276 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
277
278 // X- S- K+
279 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
280 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
281 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
282
283 // X0 S- K0
284 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
285 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
286 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
287
288 // X- p pi-
289 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.46, 0.91,
290 1.36, 3.00, 3.56, 3.32, 3.10, 3.10, 3.00, 2.80, 2.60, 2.40,
291 2.20, 2.00, 1.90, 1.80, 1.70, 1.65, 1.60, 1.55, 1.54, 1.53, 1.40},
292
293 // X- S0 K0
294 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
295 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
296 0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
297
298 //
299 // multiplicity 4 (52 channels)
300 //
301 // p n K- K-
302 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
303 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
304 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
305
306 // p L K- pi-
307 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
308 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
309 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
310
311 // p S0 K- pi-
312 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
313 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
314 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
315
316 // p S- K- pi0
317 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
318 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
319 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
320
321 // p S- K0bar pi-
322 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
323 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
324 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
325
326 // p X0 pi- pi-
327 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
328 0.0, 0.0, 0.0, 0.06, 0.23, 0.40, 0.57, 0.75, 0.70, 0.60,
329 0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.10, 0.08, 0.07, 0.07, 0.04},
330
331 // p X- pi- pi0
332 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
333 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
334 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
335
336 // p X- K- K0
337 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
338 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
339 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
340
341 // n n K0bar K-
342 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
343 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
344 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
345
346 // n L K- pi0
347 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
348 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
349 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
350
351 // n L K0bar pi-
352 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
353 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
354 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
355
356 // n S0 K- pi0
357 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
358 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
359 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
360
361 // n S0 K0bar pi-
362 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
363 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
364 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
365
366 // n S+ K- pi-
367 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
368 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
369 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
370
371 // n S- K0bar pi0
372 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
373 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
374 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
375
376 // n S- K- pi+
377 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
378 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
379 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
380
381 // n X- 2pi0
382 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
383 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.13, 0.12, 0.11,
384 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
385
386 // n X- pi+ pi-
387 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
388 0.0, 0.0, 0.0, 0.04, 0.09, 0.14, 0.20, 0.25, 0.20, 0.17,
389 0.15, 0.12, 0.10, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.02, 0.01},
390
391 // n X- K0 K0bar
392 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
393 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
394 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
395
396 // n X- K+ K-
397 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
398 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
399 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
400
401 // n X0 pi- pi0
402 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
403 0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.15, 0.14, 0.13,
404 0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.03},
405
406 // n X0 K- K0
407 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
408 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
409 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
410
411 // L L pi- pi0
412 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
413 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
414 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
415
416 // L L K- K0
417 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
418 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
419 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
420
421 // L S0 pi- pi0
422 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
423 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
424 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
425
426 // L S0 K- K0
427 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
428 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
429 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
430
431 // L S+ 2pi-
432 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
433 0.0, 0.0, 0.0, 0.07, 0.34, 0.60, 0.87, 1.14, 1.50, 1.40,
434 1.30, 1.20, 1.10, 1.00, 0.85, 0.75, 0.60, 0.50, 0.48, 0.46, 0.30},
435
436 // L S- K0bar K0
437 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
438 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
439 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
440
441 // L S- K+ K-
442 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
443 0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
444 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
445
446 // L S- 2pi0
447 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
448 0.0, 0.0, 0.0, 0.12, 0.56, 1.00, 1.45, 1.90, 2.50, 2.20,
449 2.05, 1.90, 1.75, 1.60, 1.45, 1.30, 1.10, 1.00, 0.88, 0.86, 0.70},
450
451 // L S- pi+ pi-
452 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
453 0.0, 0.0, 0.0, 0.12, 0.56, 1.00, 1.45, 1.90, 2.50, 2.40,
454 2.30, 2.10, 1.95, 1.80, 1.65, 1.50, 1.35, 1.20, 1.18, 1.16, 0.90},
455
456 // L X0 K0 pi-
457 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
458 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
459 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
460
461 // L X- K0 pi0
462 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
463 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
464 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
465
466 // L X- K+ pi-
467 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
468 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
469 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
470
471 // S0 S- 2pi0
472 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
473 0.0, 0.0, 0.0, 0.04, 0.12, 0.20, 0.28, 0.36, 0.24, 0.20,
474 0.17, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
475
476 // S0 S- pi+ pi-
477 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
478 0.0, 0.0, 0.0, 0.06, 0.19, 0.33, 0.46, 0.60, 0.40, 0.25,
479 0.15, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
480
481 // S0 S- K+ K-
482 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
483 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
484 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
485
486 // S0 S- K0 K0bar
487 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
488 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
489 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
490
491 // S0 S0 pi- pi0
492 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
493 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
494 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
495
496 // S0 S0 K0 K-
497 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
498 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
499 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
500
501 // S+ S0 2pi-
502 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
503 0.0, 0.0, 0.0, 0.06, 0.19, 0.33, 0.46, 0.60, 0.40, 0.25,
504 0.15, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
505
506 // S+ S- pi- pi0
507 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
508 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
509 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
510
511 // S- S- pi+ pi0
512 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
513 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
514 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
515
516 // S+ S- K0 K-
517 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
518 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
519 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
520
521 // S- S- K+ K0bar
522 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
523 0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
524 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
525
526 // S0 X0 K0 pi-
527 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
528 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
529 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
530
531 // S0 X- K0 pi0
532 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
533 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
534 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
535
536 // S0 X- K+ pi-
537 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
538 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
539 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
540
541 // S+ X- K0 pi-
542 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
543 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
544 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
545
546 // S- X0 K+ pi-
547 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
548 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
549 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
550
551 // S- X- K0 pi+
552 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
553 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
554 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
555
556 // S- X- K+ pi0
557 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
558 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
559 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
560
561 // S- X0 K0 pi0
562 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
563 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
564 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
565
566 //
567 // multiplicity 5 (2 channels)
568 //
569 // n X- 3pi0
570 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
571 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
572 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
573
574 // n X- pi+ 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.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
577 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
578
579 //
580 // multiplicity 6 (2 channels)
581 //
582 // n X- 4pi0
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.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
585 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
586
587 // n X- pi+ pi- 2pi0
588 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
589 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
590 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
591
592 //
593 // multiplicity 7 (2 channels)
594 //
595 // n X- 5pi0
596 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
597 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01,
598 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
599
600 // n X- pi+ pi- 3pi0
601 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
602 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01,
603 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01} };
Note: See TracBrowser for help on using the repository browser.