Geant4  10.02
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: G4CascadeT31piNChannel.cc 67796 2013-03-08 06:18:39Z mkelsey $
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 
35 #include "G4InuclParticleNames.hh"
36 using namespace G4InuclParticleNames;
37 
38 namespace {
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},
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] =
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 
97 namespace {
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 
156 namespace {
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 
692 G4CascadePiMinusPChannelData::data(pimP2bfs, pimP3bfs, pimP4bfs,
693  pimP5bfs, pimP6bfs, pimP7bfs,
694  pimP8bfs, pimP9bfs, pimPCrossSections,
695  pimPtotXSec, pim*pro, "PiMinusP");
696 
698 G4CascadePiPlusNChannelData::data(pipN2bfs, pipN3bfs, pipN4bfs,
699  pipN5bfs, pipN6bfs, pipN7bfs,
700  pipN8bfs, pipN9bfs, pimPCrossSections,
701  pimPtotXSec, pip*neu, "PiPlusN");
int G4int
Definition: G4Types.hh:78
double G4double
Definition: G4Types.hh:76