Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
G4CascadeT31piNChannel.cc
Go to the documentation of this file.
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$
27//
28// 20100804 M. Kelsey -- Add name string to ctor
29// 20110719 M. Kelsey -- Add initial state code to ctor
30// 20110725 M. Kelsey -- Instantiate cross-section object for self-registration
31// 20110916 M. Kelsey -- Drop self-registration due to platform inconsistencies
32
33#include "G4CascadePiMinusPChannel.hh"
34#include "G4CascadePiPlusNChannel.hh"
35#include "G4InuclParticleNames.hh"
36using namespace G4InuclParticleNames;
37
38namespace {
39 // pi- p : Outgoing particle types of a given multiplicity
40 static const G4int pimP2bfs[5][2] =
41 {{pim,pro}, {pi0,neu}, {k0,lam}, {k0,s0}, {kpl,sm}};
42
43 static const G4int pimP3bfs[13][3] =
44 {{pim,pro,pi0}, {pim,neu,pip}, {pi0,neu,pi0}, {pi0,lam,k0},
45 {pim,lam,kpl}, {pip,sm,k0}, {pi0,sm,kpl}, {pim,sp,k0},
46 {pim,s0,kpl}, {pi0,s0,k0}, {kmi,pro,k0}, {kmi,neu,kpl},
47 {k0b,neu,k0}};
48
49 static const G4int pimP4bfs[22][4] =
50 {{pim,pro,pip,pim}, {pim,pro,pi0,pi0}, {pim,neu,pip,pi0},
51 {pi0,neu,pi0,pi0}, {pim,lam,k0,pip}, {pi0,lam,k0,pi0},
52 {pim,lam,kpl,pi0}, {pim,s0,k0,pip}, {pi0,s0,k0,pi0},
53 {pim,s0,kpl,pi0}, {pim,sp,kpl,pim}, {pim,sp,k0,pi0},
54 {pim,sm,kpl,pip}, {pi0,sm,kpl,pi0}, {pip,sm,k0,pi0},
55 {pim,pro,kpl,kmi}, {pim,pro,k0,k0b}, {pi0,pro,k0,kmi},
56 {pip,neu,k0,kmi}, {pi0,neu,k0,k0b}, {pi0,neu,kpl,kmi},
57 {pim,neu,kpl,k0b}};
58
59 static const G4int pimP5bfs[31][5] =
60 {{pim,pro,pip,pim,pi0}, {pim,pro,pi0,pi0,pi0}, {pim,neu,pip,pip,pim},
61 {pim,neu,pip,pi0,pi0}, {pi0,neu,pi0,pi0,pi0}, {pim,lam,k0,pip,pi0},
62 {pim,lam,kpl,pi0,pi0}, {pim,lam,kpl,pip,pim}, {pi0,lam,k0,pi0,pi0},
63 {pim,s0,kpl,pip,pim}, {pim,s0,kpl,pi0,pi0}, {pim,s0,k0,pip,pi0},
64 {pi0,s0,k0,pi0,pi0}, {pim,sp,k0,pip,pim}, {pim,sp,k0,pi0,pi0},
65 {pim,sp,kpl,pim,pi0}, {pim,sm,k0,pip,pip}, {pip,sm,k0,pi0,pi0},
66 {pim,sm,kpl,pip,pi0}, {pi0,sm,kpl,pi0,pi0}, {pim,pro,pi0,kpl,kmi},
67 {pim,pro,pi0,k0,k0b}, {pim,pro,pip,k0,kmi}, {pi0,pro,pi0,k0,kmi},
68 {pim,pro,pim,kpl,k0b}, {pim,neu,pip,kpl,kmi}, {pim,neu,pip,k0,k0b},
69 {pip,neu,pi0,k0,kmi}, {pim,neu,pi0,kpl,k0b}, {pi0,neu,pi0,k0,k0b},
70 {pi0,neu,pi0,kpl,kmi}};
71
72 static const G4int pimP6bfs[6][6] =
73 {{pim,pro,pip,pip,pim,pim}, {pim,pro,pip,pim,pi0,pi0},
74 {pim,pro,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pim,pi0},
75 {pim,neu,pip,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0}};
76
77 static const G4int pimP7bfs[7][7] =
78 {{pim,pro,pip,pip,pim,pim,pi0}, {pim,pro,pip,pim,pi0,pi0,pi0},
79 {pim,pro,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim},
80 {pim,neu,pip,pip,pim,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0},
81 {pi0,neu,pi0,pi0,pi0,pi0,pi0}};
82
83 static const G4int pimP8bfs[8][8] =
84 {{pim,pro,pip,pip,pip,pim,pim,pim}, {pim,pro,pip,pip,pim,pim,pi0,pi0},
85 {pim,pro,pip,pim,pi0,pi0,pi0,pi0}, {pim,pro,pi0,pi0,pi0,pi0,pi0,pi0},
86 {pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0,pi0},
87 {pim,neu,pip,pip,pim,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim,pi0}};
88
89 static const G4int pimP9bfs[9][9] =
90 {{pim,pro,pip,pip,pip,pim,pim,pim,pi0}, {pim,pro,pip,pip,pim,pim,pi0,pi0,pi0},
91 {pim,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pim,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
92 {pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
93 {pim,neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {pim,neu,pip,pip,pip,pim,pim,pi0,pi0},
94 {pim,neu,pip,pip,pip,pip,pim,pim,pim}};
95}
96
97namespace {
98 // pi+ n : Outgoing particle types of a given multiplicity
99 static const G4int pipN2bfs[5][2] =
100 {{pip,neu}, {pi0,pro}, {kpl,lam}, {kpl,s0}, {k0,sp}};
101
102 static const G4int pipN3bfs[13][3] =
103 {{pip,neu,pi0}, {pip,pro,pim}, {pi0,pro,pi0}, {pi0,lam,kpl},
104 {pip,lam,k0}, {pim,sp,kpl}, {pi0,sp,k0}, {pip,sm,kpl},
105 {pip,s0,k0}, {pi0,s0,kpl}, {k0b,neu,kpl}, {k0b,pro,k0},
106 {kmi,pro,kpl}};
107
108 static const G4int pipN4bfs[22][4] =
109 {{pip,neu,pip,pim}, {pip,neu,pi0,pi0}, {pip,pro,pim,pi0},
110 {pi0,pro,pi0,pi0}, {pip,lam,kpl,pim}, {pi0,lam,kpl,pi0},
111 {pip,lam,k0,pi0}, {pip,s0,kpl,pim}, {pi0,s0,kpl,pi0},
112 {pip,s0,k0,pi0}, {pip,sm,k0,pip}, {pip,sm,kpl,pi0},
113 {pip,sp,k0,pim}, {pi0,sp,k0,pi0}, {pim,sp,kpl,pi0},
114 {pip,neu,k0,k0b}, {pip,neu,kpl,kmi}, {pi0,neu,kpl,k0b},
115 {pim,pro,kpl,k0b}, {pi0,pro,kpl,kmi}, {pi0,pro,k0,k0b},
116 {pip,pro,k0,kmi}};
117
118 static const G4int pipN5bfs[31][5] =
119 {{pip,neu,pip,pim,pi0}, {pip,neu,pi0,pi0,pi0}, {pip,pro,pip,pim,pim},
120 {pip,pro,pim,pi0,pi0}, {pi0,pro,pi0,pi0,pi0}, {pip,lam,kpl,pim,pi0},
121 {pip,lam,k0,pi0,pi0}, {pip,lam,k0,pip,pim}, {pi0,lam,kpl,pi0,pi0},
122 {pip,s0,k0,pip,pim}, {pip,s0,k0,pi0,pi0}, {pip,s0,kpl,pim,pi0},
123 {pi0,s0,kpl,pi0,pi0}, {pip,sm,kpl,pip,pim}, {pip,sm,kpl,pi0,pi0},
124 {pip,sm,k0,pip,pi0}, {pip,sp,kpl,pim,pim}, {pim,sp,kpl,pi0,pi0},
125 {pip,sp,k0,pim,pi0}, {pi0,sp,k0,pi0,pi0}, {pip,neu,pi0,k0,k0b},
126 {pip,neu,pi0,kpl,kmi}, {pip,neu,pim,kpl,k0b}, {pi0,neu,pi0,kpl,k0b},
127 {pip,neu,pip,k0,kmi}, {pip,pro,pim,k0,k0b}, {pip,pro,pim,kpl,kmi},
128 {pim,pro,pi0,kpl,k0b}, {pip,pro,pi0,k0,kmi}, {pi0,pro,pi0,kpl,kmi},
129 {pi0,pro,pi0,k0,k0b}};
130
131 static const G4int pipN6bfs[6][6] =
132 {{pip,neu,pip,pip,pim,pim}, {pip,neu,pip,pim,pi0,pi0},
133 {pip,neu,pi0,pi0,pi0,pi0}, {pip,pro,pip,pim,pim,pi0},
134 {pip,pro,pim,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0}};
135
136 static const G4int pipN7bfs[7][7] =
137 {{pip,neu,pip,pip,pim,pim,pi0}, {pip,neu,pip,pim,pi0,pi0,pi0},
138 {pip,neu,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim},
139 {pip,pro,pip,pim,pim,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0},
140 {pi0,pro,pi0,pi0,pi0,pi0,pi0}};
141
142 static const G4int pipN8bfs[8][8] =
143 {{pip,neu,pip,pip,pip,pim,pim,pim}, {pip,neu,pip,pip,pim,pim,pi0,pi0},
144 {pip,neu,pip,pim,pi0,pi0,pi0,pi0}, {pip,neu,pi0,pi0,pi0,pi0,pi0,pi0},
145 {pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0,pi0},
146 {pip,pro,pip,pim,pim,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim,pi0}};
147
148 static const G4int pipN9bfs[9][9] =
149 {{pip,neu,pip,pip,pip,pim,pim,pim,pi0}, {pip,neu,pip,pip,pim,pim,pi0,pi0,pi0},
150 {pip,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pip,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
151 {pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
152 {pip,pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pip,pro,pip,pip,pim,pim,pim,pi0,pi0},
153 {pip,pro,pip,pip,pip,pim,pim,pim,pim}};
154}
155
156namespace {
157 // Total pi- p cross section as a function of kinetic energy
158 // New cs values after 9B tuning (27 July 09)
159 static const G4double pimPtotXSec[30] =
160 { 6.13, 6.4, 6.67, 6.94, 7.22, 7.5, 8.3, 12.0, 14.4, 24.0,
161 46.0, 72.04, 43.02, 27.19, 27.32, 43.8, 37.08, 51.37, 34.21, 34.79,
162 32.08, 31.19, 30.32, 28.5, 27.0, 25.9, 25.5, 25.2, 25.0, 24.8};
163
164 // pi- p cross sections as functions of kinetic energy and multiplicity
165 const G4double pimPCrossSections[101][30] = {
166 //
167 // multiplicity 2 (5 channels)
168 //
169 // pi- p (pi+ n)
170 { 1.43, 1.5, 1.57, 1.64, 1.72, 1.8, 2.0, 3.0, 3.4, 7.0,
171 14.0, 24.0, 14.7, 10.5, 10.84, 21.0, 13.97, 25.1, 10.46, 9.88,
172 8.0, 7.1, 6.0, 5.7, 5.0, 4.6, 4.3, 4.0, 3.8, 3.7},
173
174 // n pi0 (p pi0)
175 { 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 6.3, 9.0, 11.0, 17.0,
176 32.0, 48.0, 28.0, 14.5, 11.04, 8.99, 4.79, 5.02, 2.08, 0.94,
177 0.5, 0.25, 0.15, 0.09, 0.06, 0.05, 0.04, 0.02, 0.02, 0.01},
178
179 // L K0 (L K+)
180 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
181 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.65, 0.29, 0.17,
182 0.13, 0.08, 0.05, 0.03, 0.02, 0.01, 0.01, 0.01, 0.01, 0.0},
183
184 // S0 K0 (S0 K+)
185 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
186 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.28, 0.19, 0.12,
187 0.09, 0.06, 0.04, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0},
188
189 // S- K+ (S+ K0)
190 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
191 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.20, 0.25, 0.09,
192 0.04, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
193 //
194 // multiplicity 3 (13 channels)
195 //
196 // p pi- pi0 (n pi+ pi0)
197 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
198 0.0, 0.0, 0.05, 0.18, 0.86, 3.89, 5.19, 6.63, 5.11, 4.03,
199 3.2, 2.5, 2.0, 1.4, 0.97, 0.68, 0.55, 0.36, 0.3, 0.22},
200
201 // n pi+ pi- (p pi+ pi-)
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.03, 0.19, 1.41, 3.12, 6.2, 8.14, 8.23, 7.84, 5.91,
204 4.58, 3.58, 2.71, 1.73, 1.06, 0.85, 0.62, 0.38, 0.31, 0.22},
205
206 // n pi0 pi0 (p pi0 pi0)
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.01, 0.08, 0.6, 1.34, 2.7, 2.99, 2.51, 1.35, 0.87,
209 0.61, 0.24, 0.15, 0.10, 0.06, 0.04, 0.0, 0.0, 0.0, 0.0},
210
211 // L K0 pi0 (L K+ pi0)
212 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
213 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.19,
214 0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.0},
215
216 // L K+ pi- (L K0 pi+)
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.0, 0.0, 0.0, 0.0, 0.0, 0.04, 0.12,
219 0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01, 0.01, 0.0},
220
221 // S- K0 pi+ (S+ K+ pi-)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.1,
224 0.13, 0.07, 0.03, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
225
226 // S- K+ pi0 (S+ K0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
229 0.05, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
230
231 // S+ K0 pi- (S- K+ 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
234 0.06, 0.04, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
235
236 // S0 K+ pi- (S0 K0 pi+)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
239 0.07, 0.04, 0.03, 0.02, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0},
240
241 // S0 K0 pi0 (S0 K+ pi0)
242 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
243 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02, 0.09,
244 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.0, 0.0, 0.0},
245
246 // p K0 K- (n K+ K0bar)
247 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
248 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03,
249 0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
250
251 // n K+ K- (p K0 K0bar)
252 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
253 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
254 0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01, 0.0, 0.0, 0.0},
255
256 // n K0 K0bar (p K+ K-)
257 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
258 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06,
259 0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0},
260 //
261 // multiplicity 4 (22 channels)
262 //
263 // p pi+ pi- pi- (n pi+ pi+ pi-)
264 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
265 0.0, 0.0, 0.0, 0.0, 0.01, 0.04, 0.18, 0.43, 0.82, 1.69,
266 1.8, 1.99, 1.8, 1.71, 1.44, 1.26, 1.17, 0.99, 1.04, 0.9},
267
268 // p pi- pi0 pi0 (n pi+ pi0 pi0)
269 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
270 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.09, 0.21, 1.5, 2.5,
271 1.2, 1.2, 1.0, 0.9, 0.8, 0.7, 0.65, 0.6, 0.55, 0.5},
272
273 // n pi+ pi- pi0 (p pi+ pi- pi0)
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.4, 0.7, 1.11, 2.6, 3.9,
276 4.07, 3.8, 2.76, 1.38, 1.16, 0.97, 0.85, 0.7, 0.55, 0.4},
277
278 // n pi0 pi0 pi0 (p pi0 pi0 pi0)
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.1, 0.55, 1.0, 0.87, 0.73, 0.71,
281 0.7, 0.56, 0.35, 0.29, 0.21, 0.14, 0.1, 0.06, 0.03, 0.0},
282
283 // L K0 pi+ pi- (L K+ pi+ pi-)
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.0, 0.0, 0.0, 0.0, 0.01,
286 0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03},
287
288 // L K0 pi0 pi0 (L K+ pi0 pi0)
289 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
290 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
291 0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03},
292
293 // L K+ pi- pi0 (L K0 pi+ pi0)
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.0, 0.0, 0.0, 0.0, 0.01,
296 0.06, 0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02},
297
298 // S0 K0 pi+ pi- (S0 K+ pi+ pi-)
299 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
300 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
301 0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
302
303 // S0 K0 pi0 pi0 (S0 K+ pi0 pi0)
304 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
305 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
306 0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
307
308 // S0 K+ pi- pi0 (S0 K0 pi+ pi0)
309 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
310 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
311 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
312
313 // S+ K+ pi- pi- (S- K0 pi+ pi+)
314 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
315 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
316 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
317
318 // S+ K0 pi- pi0 (S- K+ pi+ pi0)
319 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
320 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
321 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
322
323 // S- K+ pi+ pi- (S+ K0 pi+ pi-)
324 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
325 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
326 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
327
328 // S- K+ pi0 pi0 (S+ K0 pi0 pi0)
329 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
330 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
331 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
332
333 // S- K0 pi+ pi0 (S+ K+ pi- pi0)
334 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
335 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
336 0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01},
337
338 // p pi- K+ K- (n pi+ K0 K0bar)
339 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.06, 0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05},
342
343 // p pi- K0 K0bar (n pi+ K+ K-)
344 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
345 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
346 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
347
348 // p pi0 K0 K- (n pi0 K+ K0bar)
349 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
350 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
351 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
352
353 // n pi+ K0 K- (p pi- K+ K0bar)
354 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
355 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
356 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
357
358 // n pi0 K0 K0bar (p pi0 K+ K-)
359 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
360 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
361 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
362
363 // n pi0 K+ K- (p pi0 K0 K0bar)
364 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
365 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
366 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
367
368 // n pi- K+ K0bar (p pi+ K0 K-)
369 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
370 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
371 0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05},
372 //
373 // multiplicity 5 (31 channels)
374 //
375 // p pi+ pi- pi- pi0 (n pi+ pi+ pi- pi0)
376 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
377 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.10, 0.93,
378 1.5, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62},
379
380 // p pi- pi0 pi0 pi0 (n pi+ pi0 pi0 pi0)
381 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
382 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02, 0.10, 0.73,
383 1.4, 1.9, 2.2, 2.0, 1.7, 1.4, 1.2, 0.90, 0.76, 0.62},
384
385 // n pi+ pi+ pi- pi- (p pi+ pi+ pi- pi-)
386 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
387 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03, 0.10, 0.31,
388 0.66, 0.93, 1.2, 1.2, 1.2, 0.94, 0.74, 0.53, 0.40, 0.30},
389
390 // n pi+ pi- pi0 pi0 (p pi+ pi- pi0 pi0)
391 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
392 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.22, 0.4,
393 0.8, 0.95, 1.0, 0.9, 0.75, 0.65, 0.55, 0.40, 0.35, 0.30},
394
395 // n pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0)
396 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
397 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.22, 0.39,
398 0.66, 0.9, 1.0, 0.8, 0.7, 0.55, 0.5, 0.35, 0.30, 0.25},
399
400 // L K0 pi+ pi- pi0 (L K+ pi+ pi- pi0)
401 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
402 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
403 0.01, 0.05, 0.07, 0.08, 0.08, 0.07, 0.06, 0.06, 0.05, 0.05},
404
405 // L K+ pi- pi0 pi0 (L K0 pi+ pi0 pi0)
406 { 0.0, 0.0, 0.0, 0.0, 0.0, 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,
408 0.0, 0.02, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
409
410 // L K+ pi+ pi- pi- (L K0 pi+ pi+ pi-)
411 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
414
415 // L K0 pi0 pi0 pi0 (L K+ pi0 pi0 pi0)
416 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
419
420 // S0 K+ pi+ pi- pi- (S0 K0 pi+ pi+ pi-)
421 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
424
425 // S0 K+ pi- pi0 pi0 (S0 K0 pi+ pi0 pi0)
426 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
427 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
428 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
429
430 // S0 K0 pi+ pi- pi0 (S0 K+ pi+ pi- pi0)
431 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
432 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
433 0.0, 0.02, 0.03, 0.04, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02},
434
435 // S0 K0 pi0 pi0 pi0 (S0 K+ pi0 pi0 pi0)
436 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
437 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
438 0.0, 0.01, 0.01, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
439
440 // S+ K0 pi+ pi- pi- (S- K+ pi+ pi+ pi-)
441 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
442 0.0, 0.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.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
444
445 // S+ K0 pi- pi0 pi0 (S- K+ pi+ pi0 pi0)
446 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
447 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
448 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
449
450 // S+ K+ pi- pi- pi0 (S- K0 pi+ pi+ pi0)
451 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
452 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
453 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
454
455 // S- K0 pi+ pi+ pi- (S+ K+ pi+ pi- pi-)
456 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
457 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
458 0.0, 0.01, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0},
459
460 // S- K0 pi+ pi0 pi0 (S+ K+ pi- pi0 pi0)
461 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
462 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
463 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.0, 0.0, 0.0, 0.0},
464
465 // S- K+ pi+ pi- pi0 (S+ K0 pi+ pi- pi0)
466 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
467 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
468 0.0, 0.01, 0.02, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
469
470 // S- K+ pi0 pi0 pi0 (S+ K0 pi0 pi0 pi0)
471 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
472 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
473 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.0},
474
475 // p pi- pi0 K+ K- (n pi+ pi0 K0 K0bar)
476 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
477 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
478 0.0, 0.01, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04},
479
480 // p pi- pi0 K0 K0bar (n pi+ pi0 K+ K-)
481 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
484
485 // p pi+ pi- K0 K- (n pi+ pi- K+ K0bar)
486 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
489
490 // p pi0 pi0 K0 K- (n pi0 pi0 K+ K0bar)
491 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.02, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
494
495 // p pi- pi- K+ K0bar (n pi+ pi+ K0 K-)
496 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02},
499
500 // n pi+ pi- K+ K- (p pi+ pi- K0 K0bar)
501 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03},
504
505 // n pi+ pi- K0 K0bar (p pi+ pi- K+ K-)
506 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
507 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
508 0.0, 0.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
509
510 // n pi+ pi0 K0 K- (p pi- pi0 K+ K0bar)
511 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
512 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
513 0.0, 0.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03},
514
515 // n pi- pi0 K+ K0bar (p pi+ pi0 K0 K-)
516 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
517 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
518 0.0, 0.01, 0.07, 0.06, 0.05, 0.04, 0.04, 0.04, 0.03, 0.03},
519
520 // n pi0 pi0 K0 K0bar (p pi0 pi0 K+ K-)
521 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
522 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
523 0.0, 0.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
524
525 // n pi0 pi0 K+ K- (p pi0 pi0 K0 K0bar)
526 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
527 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
528 0.0, 0.01, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
529 //
530 // multiplicity 6 (6 channels)
531 //
532 // p pi+ pi+ pi- pi- pi- (n pi+ pi+ pi+ pi- pi-)
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.01, 0.03, 0.04,
535 0.06, 0.12, 0.16, 0.22, 0.31, 0.34, 0.3, 0.24, 0.19, 0.16},
536
537 // p pi+ pi- pi- pi0 pi0 (n pi+ pi+ pi- pi0 pi0)
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.01, 0.02, 0.05,
540 0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
541
542 // p pi- pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0)
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.01, 0.02, 0.05,
545 0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
546
547 // n pi+ pi+ pi- pi- pi0 (p pi+ pi+ pi- pi- pi0)
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.01, 0.02, 0.05,
550 0.07, 0.1, 0.12, 0.15, 0.18, 0.2, 0.19, 0.2, 0.17, 0.16},
551
552 // n pi+ pi- pi0 pi0 pi0 (p pi+ pi- pi0 pi0 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.01, 0.03,
555 0.05, 0.06, 0.08, 0.1, 0.12, 0.13, 0.14, 0.13, 0.11, 0.11},
556
557 // n pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 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.01, 0.03,
560 0.04, 0.05, 0.06, 0.08, 0.09, 0.1, 0.1, 0.1, 0.08, 0.07},
561 //
562 // multiplicity 7 (7 channels)
563 //
564 // p pi+ pi+ pi- pi- pi- pi0
565 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
566 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.03,
567 0.05, 0.12, 0.19, 0.37, 0.67, 0.70, 0.65, 0.60, 0.50, 0.45},
568
569 // p pi+ pi- pi- pi0 pi0 pi0
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.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
572 0.03, 0.07, 0.11, 0.22, 0.4, 0.4, 0.35, 0.30, 0.25, 0.2},
573
574 // p pi- pi0 pi0 pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
577 0.02, 0.03, 0.05, 0.11, 0.20, 0.2, 0.15, 0.12, 0.10, 0.08},
578
579 // n pi+ pi+ pi+ pi- pi- pi- (measured R 501)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
582 0.01, 0.03, 0.06, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13},
583
584 // n pi+ pi+ pi- pi- pi0 pi0 (measured R 496)
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.02,
587 0.03, 0.07, 0.11, 0.22, 0.4, 0.45, 0.4, 0.4, 0.3, 0.3},
588
589 // n pi+ pi- pi0 pi0 pi0 pi0
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.0, 0.01,
592 0.02, 0.04, 0.07, 0.14, 0.26, 0.26, 0.22, 0.20, 0.18, 0.15},
593
594 // n pi0 pi0 pi0 pi0 pi0 pi0
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.0, 0.0,
597 0.01, 0.02, 0.03, 0.05, 0.1, 0.11, 0.1, 0.1, 0.08, 0.7},
598 //
599 // multiplicity 8 (8 channels)
600 //
601 // p pi+ pi+ pi+ pi- pi- pi- pi- (n pi+ pi+ pi+ pi+ pi- pi- pi-) (measured R 420)
602 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
603 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
604 0.0, 0.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09},
605
606 // p pi+ pi+ pi- pi- pi- pi0 pi0 (n pi+ pi+ pi+ pi- pi- pi0 pi0)
607 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
608 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
609 0.0, 0.01, 0.03, 0.09, 0.13, 0.13, 0.11, 0.10, 0.09, 0.08},
610
611 // p pi+ pi- pi- pi0 pi0 pi0 pi0 (n pi+ pi+ pi- pi0 pi0 pi0 pi0)
612 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
613 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
614 0.0, 0.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09},
615
616 // p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n pi+ pi0 pi0 pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0,
619 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.03, 0.03, 0.03},
620
621 // n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p pi0 pi0 pi0 pi0 pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0,
624 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.03, 0.03, 0.03, 0.02},
625
626 // n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p pi+ pi- pi0 pi0 pi0 pi0 pi0)
627 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
628 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
629 0.0, 0.0, 0.01, 0.03, 0.04, 0.05, 0.05, 0.05, 0.05, 0.04},
630
631 // n pi+ pi+ pi- pi- pi0 pi0 pi0 (p pi+ pi+ pi- pi- pi0 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.01, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
635
636 // n pi+ pi+ pi+ pi- pi- pi- pi0 (p pi+ pi+ pi+ pi- pi- 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.0, 0.01, 0.03, 0.09, 0.12, 0.13, 0.11, 0.09, 0.08, 0.07},
640 //
641 // multiplicity 9 (9 channels)
642 //
643 // p pi+ pi+ pi+ pi- pi- pi- pi- pi0 (measured R 418)
644 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
645 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
646 0.0, 0.0, 0.01, 0.04, 0.08, 0.14, 0.18, 0.18, 0.18, 0.17},
647
648 // p pi+ pi+ pi- pi- pi- pi0 pi0 pi0
649 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
650 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
651 0.0, 0.0, 0.01, 0.02, 0.05, 0.08, 0.11, 0.11, 0.11, 0.1},
652
653 // p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0
654 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
655 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
656 0.0, 0.0, 0.0, 0.02, 0.03, 0.06, 0.07, 0.07, 0.07, 0.07},
657
658 // p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0
659 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
660 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
661 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
662
663 // n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
664 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
665 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
666 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
667
668 // n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0
669 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
670 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
671 0.0, 0.0, 0.0, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04},
672
673 // n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0
674 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
675 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
676 0.0, 0.0, 0.0, 0.02, 0.03, 0.06, 0.07, 0.07, 0.07, 0.07},
677
678 // n pi+ pi+ pi+ pi- pi- pi- pi0 pi0
679 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
680 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
681 0.0, 0.0, 0.01, 0.02, 0.05, 0.08, 0.11, 0.11, 0.11, 0.1},
682
683 // n pi+ pi+ pi+ pi+ pi- pi- pi- pi- (measured R 503)
684 { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
685 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
686 0.0, 0.0, 0.01, 0.04, 0.08, 0.14, 0.18, 0.18, 0.18, 0.17} };
687 }
688
689// Initialize both |T Tz> = |3/2 1/2> channels, using pipP cross-section table
690
691G4CascadePiMinusPChannelData::data_t
692G4CascadePiMinusPChannelData::data(pimP2bfs, pimP3bfs, pimP4bfs,
693 pimP5bfs, pimP6bfs, pimP7bfs,
694 pimP8bfs, pimP9bfs, pimPCrossSections,
695 pimPtotXSec, pim*pro, "PiMinusP");
696
697G4CascadePiPlusNChannelData::data_t
698G4CascadePiPlusNChannelData::data(pipN2bfs, pipN3bfs, pipN4bfs,
699 pipN5bfs, pipN6bfs, pipN7bfs,
700 pipN8bfs, pipN9bfs, pimPCrossSections,
701 pimPtotXSec, pip*neu, "PiPlusN");
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66