Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4CascadeT11pizNChannel.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: G4CascadeT11pizNChannel.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 
33 #include "G4CascadePiZeroPChannel.hh"
34 #include "G4CascadePiZeroNChannel.hh"
35 #include "G4InuclParticleNames.hh"
36 using namespace G4InuclParticleNames;
37 
38 namespace {
39  // pi0 p : Outgoing particle types of a given multiplicity
40  static const G4int pizP2bfs[5][2] =
41  {{pi0,pro}, {pip,neu}, {kpl,lam}, {kpl,s0}, {k0,sp}};
42 
43  static const G4int pizP3bfs[13][3] =
44  {{pip,pro,pim}, {pi0,pro,pi0}, {pi0,neu,pip}, {pi0,lam,kpl},
45  {pip,lam,k0}, {pi0,s0,kpl}, {pip,s0,k0}, {pi0,sp,k0},
46  {pim,sp,kpl}, {pip,sm,kpl}, {kmi,pro,kpl}, {k0b,pro,k0},
47  {k0b,neu,kpl}};
48 
49  static const G4int pizP4bfs[21][4] =
50  {{pi0,pro,pim,pip}, {pi0,pro,pi0,pi0}, {pip,neu,pip,pim},
51  {pi0,neu,pip,pi0}, {pip,lam,kpl,pim}, {pip,s0,kpl,pim},
52  {pi0,lam,kpl,pi0}, {pi0,s0,kpl,pi0}, {pi0,sp,k0,pi0},
53  {pi0,lam,k0,pip}, {pi0,s0,k0,pip}, {pip,sp,k0,pim},
54  {pi0,sp,kpl,pim}, {pi0,sm,kpl,pip}, {pi0,pro,kpl,kmi},
55  {pi0,pro,k0,k0b}, {pip,pro,k0,kmi}, {pim,pro,kpl,k0b},
56  {pip,neu,kpl,kmi}, {pip,neu,k0,k0b}, {pi0,neu,kpl,k0b}};
57 
58  static const G4int pizP5bfs[30][5] =
59  {{pip,pro,pip,pim,pim}, {pi0,pro,pip,pim,pi0}, {pi0,pro,pi0,pi0,pi0},
60  {pi0,neu,pip,pip,pim}, {pi0,neu,pip,pi0,pi0}, {pi0,lam,kpl,pip,pim},
61  {pi0,lam,kpl,pi0,pi0}, {pip,lam,k0,pip,pim}, {pi0,lam,k0,pip,pi0},
62  {pi0,s0,kpl,pip,pim}, {pi0,s0,kpl,pi0,pi0}, {pip,s0,k0,pip,pim},
63  {pi0,s0,k0,pip,pi0}, {pi0,sp,k0,pip,pim}, {pi0,sp,k0,pi0,pi0},
64  {pim,sp,kpl,pip,pim}, {pi0,sp,kpl,pim,pi0}, {pip,sm,kpl,pip,pim},
65  {pi0,sm,kpl,pip,pi0}, {pi0,pro,kpl,kmi,pi0}, {pi0,pro,k0,k0b,pi0},
66  {pi0,pro,kpl,k0b,pim}, {pi0,pro,k0,kmi,pip}, {pip,pro,kpl,kmi,pim},
67  {pip,pro,k0,k0b,pim}, {pi0,neu,kpl,kmi,pip}, {pi0,neu,k0,k0b,pip},
68  {pi0,neu,kpl,k0b,pi0}, {pip,neu,k0,kmi,pip}, {pip,neu,kpl,k0b,pim}};
69 
70  static const G4int pizP6bfs[6][6] =
71  {{pi0,pro,pip,pip,pim,pim}, {pi0,pro,pip,pim,pi0,pi0},
72  {pi0,pro,pi0,pi0,pi0,pi0}, {pip,neu,pip,pip,pim,pim},
73  {pi0,neu,pip,pip,pim,pi0}, {pi0,neu,pip,pi0,pi0,pi0}};
74 
75  static const G4int pizP7bfs[7][7] =
76  {{pip,pro,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pi0},
77  {pi0,pro,pip,pim,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0,pi0},
78  {pi0,neu,pip,pip,pip,pim,pim}, {pi0,neu,pip,pip,pim,pi0,pi0},
79  {pi0,neu,pip,pi0,pi0,pi0,pi0}};
80 
81  static const G4int pizP8bfs[8][8] =
82  {{pi0,pro,pip,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pi0,pi0},
83  {pi0,pro,pip,pim,pi0,pi0,pi0,pi0}, {pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0},
84  {pip,neu,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pip,pim,pim,pi0},
85  {pi0,neu,pip,pip,pim,pi0,pi0,pi0}, {pi0,neu,pip,pi0,pi0,pi0,pi0,pi0}};
86 
87  static const G4int pizP9bfs[9][9] =
88  {{pip,pro,pip,pip,pip,pim,pim,pim,pim}, {pi0,pro,pip,pip,pip,pim,pim,pim,pi0},
89  {pi0,pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pi0,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
90  {pi0,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pi0,neu,pip,pip,pip,pip,pim,pim,pim},
91  {pi0,neu,pip,pip,pip,pim,pim,pi0,pi0}, {pi0,neu,pip,pip,pim,pi0,pi0,pi0,pi0},
92  {pi0,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}};
93 }
94 
95 namespace {
96  // pi0 n : Outgoing particle types of a given multiplicity
97  static const G4int pizN2bfs[5][2] =
98  {{pi0,neu}, {pim,pro}, {k0,lam}, {k0,s0}, {kpl,sm}};
99 
100  static const G4int pizN3bfs[13][3] =
101  {{pim,neu,pip}, {pi0,neu,pi0}, {pi0,pro,pim}, {pi0,lam,k0},
102  {pim,lam,kpl}, {pi0,s0,k0}, {pim,s0,kpl}, {pi0,sm,kpl},
103  {pip,sm,k0}, {pim,sp,k0}, {k0b,neu,k0}, {kmi,neu,kpl},
104  {kmi,pro,k0}};
105 
106  static const G4int pizN4bfs[21][4] =
107  {{pi0,neu,pim,pip}, {pi0,neu,pi0,pi0}, {pim,pro,pip,pim},
108  {pi0,pro,pim,pi0}, {pip,lam,k0,pim}, {pip,s0,k0,pim},
109  {pi0,lam,k0,pi0}, {pi0,s0,k0,pi0}, {pi0,sm,kpl,pi0},
110  {pi0,lam,kpl,pim}, {pi0,s0,kpl,pim}, {pip,sm,kpl,pim},
111  {pi0,sm,k0,pip}, {pi0,sp,k0,pim}, {pi0,neu,k0,k0b},
112  {pi0,neu,kpl,kmi}, {pim,neu,kpl,k0b}, {pip,neu,k0,kmi},
113  {pim,pro,k0,k0b}, {pim,pro,kpl,kmi}, {pi0,pro,k0,kmi}};
114 
115  static const G4int pizN5bfs[30][5] =
116  {{pim,neu,pip,pip,pim}, {pi0,neu,pip,pim,pi0}, {pi0,neu,pi0,pi0,pi0},
117  {pi0,pro,pip,pim,pim}, {pi0,pro,pim,pi0,pi0}, {pi0,lam,k0,pip,pim},
118  {pi0,lam,k0,pi0,pi0}, {pim,lam,kpl,pip,pim}, {pi0,lam,kpl,pim,pi0},
119  {pi0,s0,k0,pip,pim}, {pi0,s0,k0,pi0,pi0}, {pim,s0,kpl,pip,pim},
120  {pi0,s0,kpl,pim,pi0}, {pi0,sm,kpl,pip,pim}, {pi0,sm,kpl,pi0,pi0},
121  {pip,sm,k0,pip,pim}, {pi0,sm,k0,pip,pi0}, {pim,sp,k0,pip,pim},
122  {pi0,sp,k0,pim,pi0}, {pi0,neu,k0,k0b,pi0}, {pi0,neu,kpl,kmi,pi0},
123  {pi0,neu,k0,kmi,pip}, {pi0,neu,kpl,k0b,pim}, {pim,neu,k0,k0b,pip},
124  {pim,neu,kpl,kmi,pip}, {pi0,pro,k0,k0b,pim}, {pi0,pro,kpl,kmi,pim},
125  {pi0,pro,k0,kmi,pi0}, {pim,pro,kpl,k0b,pim}, {pim,pro,k0,kmi,pip}};
126 
127  static const G4int pizN6bfs[6][6] =
128  {{pi0,neu,pip,pip,pim,pim}, {pi0,neu,pip,pim,pi0,pi0},
129  {pi0,neu,pi0,pi0,pi0,pi0}, {pim,pro,pip,pip,pim,pim},
130  {pi0,pro,pip,pim,pim,pi0}, {pi0,pro,pim,pi0,pi0,pi0}};
131 
132  static const G4int pizN7bfs[7][7] =
133  {{pim,neu,pip,pip,pip,pim,pim}, {pi0,neu,pip,pip,pim,pim,pi0},
134  {pi0,neu,pip,pim,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0,pi0},
135  {pi0,pro,pip,pip,pim,pim,pim}, {pi0,pro,pip,pim,pim,pi0,pi0},
136  {pi0,pro,pim,pi0,pi0,pi0,pi0}};
137 
138  static const G4int pizN8bfs[8][8] =
139  {{pi0,neu,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pim,pim,pi0,pi0},
140  {pi0,neu,pip,pim,pi0,pi0,pi0,pi0}, {pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0},
141  {pim,pro,pip,pip,pip,pim,pim,pim}, {pi0,pro,pip,pip,pim,pim,pim,pi0},
142  {pi0,pro,pip,pim,pim,pi0,pi0,pi0}, {pi0,pro,pim,pi0,pi0,pi0,pi0,pi0}};
143 
144  static const G4int pizN9bfs[9][9] =
145  {{pim,neu,pip,pip,pip,pip,pim,pim,pim}, {pi0,neu,pip,pip,pip,pim,pim,pim,pi0},
146  {pi0,neu,pip,pip,pim,pim,pi0,pi0,pi0}, {pi0,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
147  {pi0,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pi0,pro,pip,pip,pip,pim,pim,pim,pim},
148  {pi0,pro,pip,pip,pim,pim,pim,pi0,pi0}, {pi0,pro,pip,pim,pim,pi0,pi0,pi0,pi0},
149  {pi0,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}};
150 }
151 
152 namespace {
153  // Total pi- p cross section as a function of kinetic energy
154  static const G4double pizPtotXSec[30] =
155  { 6.43, 7.18, 7.54, 8.01, 8.52, 9.13, 10.22, 14.37, 20.96, 34.73,
156  61.07, 98.23, 61.97, 32.62, 28.07, 31.37, 35.15, 40.17, 37.27, 33.49,
157  31.06, 29.52, 28.29, 26.91, 26.25, 25.25, 24.75, 24.35, 24.0, 23.9};
158 
159  // pi- p cross sections as functions of kinetic energy and multiplicity
160  static const G4double pizPCrossSections[99][30] = {
161  //
162  // multiplicity 2 (5 channels)
163  //
164  // p pi0 (n pi0)
165  // { 0.0, 1.15, 1.85, 2.6, 3.25, 4.4, 5.5, 9.0, 16.7, 35.5,
166  // 72.0, 107.0, 72.35, 33.1, 19.51, 17.02, 11.27, 19.11, 15.35, 11.0,
167  // 8.83, 7.4, 6.4, 5.75, 5.2, 4.7, 4.25, 4.0, 3.8, 3.65},
168 
169  // p pi0 (n pi0) go back to original for now
170  { 1.73, 2.28, 2.44, 2.71, 3.02, 3.43, 3.92, 5.37, 9.96, 17.73,
171  29.07, 50.23, 33.68, 16.69, 12.60, 11.89, 13.49, 15.24, 12.39, 9.59,
172  7.92, 6.97, 6.13, 5.37, 4.82, 4.68, 4.54, 4.0, 3.8, 3.65},
173 
174  // n pi+ (p pi-)
175  // { 0.0, 2.4, 2.8, 3.3, 4.5, 5.7, 6.3, 9.0, 11.0, 17.0,
176  // 30.0, 43.0, 30.0, 16.5, 11.0, 7.01, 4.31, 5.03, 2.05, 0.97,
177  // 0.53, 0.3, 0.2, 0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01},
178 
179  // n pi+ (p pi-) taken from pi- p -> pi0 n
180  { 4.7, 4.9, 5.1, 5.3, 5.5, 5.7, 6.3, 9.0, 11.0, 17.0,
181  32.0, 48.0, 28.0, 14.5, 11.04, 8.99, 4.79, 5.02, 2.08, 0.94,
182  0.5, 0.25, 0.15, 0.09, 0.06, 0.05, 0.04, 0.02, 0.02, 0.01},
183 
184  // L K+ (L K0)
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.65, 0.29, 0.17,
187  0.14, 0.1, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.0},
188 
189  // S0 K+ (S0 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.28, 0.18, 0.12,
192  0.10, 0.06, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0},
193 
194  // S+ K0 (S- K+)
195  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
196  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.18, 0.43, 0.22,
197  0.12, 0.06, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01, 0.0, 0.0},
198  //
199  // multiplicity 3 (13 channels)
200  //
201  // p pi+ pi- (n 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.0, 0.03, 0.19, 0.73, 3.4, 7.01, 8.35, 8.9, 5.69,
204  4.01, 2.7, 2.0, 1.30, 0.9, 0.68, 0.48, 0.34, 0.27, 0.19},
205 
206  // p pi0 pi0 (n 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.0, 0.22, 0.74, 1.8, 2.7, 3.0, 2.52, 1.33, 0.69,
209  0.44, 0.3, 0.2, 0.12, 0.07, 0.04, 0.0, 0.0, 0.0, 0.0},
210 
211  // n pi+ pi0 (p pi- 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.04, 0.5, 1.75, 3.6, 5.21, 5.28, 6.34, 5.15,
214  4.01, 3.1, 2.2, 1.42, 1.05, 0.8, 0.57, 0.33, 0.26, 0.19},
215 
216  // L K+ pi0 (L K0 pi0)
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.02, 0.13,
219  0.14, 0.1, 0.08, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
220 
221  // L K0 pi+ (L 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.03, 0.16,
224  0.15, 0.1, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
225 
226  // S0 K+ pi0 (S0 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.09,
229  0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
230 
231  // S0 K0 pi+ (S0 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.01, 0.11,
234  0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
235 
236  // S+ K0 pi0 (S- 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09,
239  0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.0, 0.0},
240 
241  // S+ K+ pi- (S- K0 pi+)
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.0, 0.09,
244  0.11, 0.07, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.01},
245 
246  // S- K+ pi+ (S+ K0 pi-)
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.08,
249  0.10, 0.05, 0.03, 0.02, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0},
250 
251  // p K+ K- (n K0 K0b)
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.02, 0.11,
254  0.12, 0.08, 0.07, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
255 
256  // p K0 K0b (n 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.02, 0.11,
259  0.12, 0.08, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
260 
261  // n K+ K0b (p K0 K-)
262  {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
263  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05,
264  0.11, 0.22, 0.2, 0.09, 0.05, 0.02, 0.01, 0.0, 0.0, 0.0},
265  //
266  // multiplicity 4 (21 channels)
267  //
268  // p pi+ pi- pi0 (n pi+ pi- 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.0, 0.04, 0.14, 0.44, 1.11, 1.88,
271  2.05, 2.07, 1.75, 1.5, 1.3, 1.1, 0.95, 0.80, 0.72, 0.66},
272 
273  // p pi0 pi0 pi0 (n pi0 pi0 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.1, 0.55, 1.0, 0.88, 0.72, 0.64,
276  0.53, 0.4, 0.35, 0.29, 0.22, 0.14, 0.1, 0.06, 0.03, 0.0},
277 
278  // n pi+ pi+ pi- (p pi+ pi- pi-)
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.05, 0.17, 0.32, 0.77, 1.79, 2.25,
281  1.82, 1.5, 1.3, 1.05, 0.9, 0.72, 0.63, 0.5, 0.43, 0.35},
282 
283  // n pi+ pi0 pi0 (p pi- pi0 pi0)
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.03, 0.09, 0.28, 0.77, 1.32,
286  1.54, 1.37, 1.2, 1.02, 0.9, 0.77, 0.69, 0.6, 0.54, 0.48},
287 
288  // L K+ pi+ pi- (L K0 pi+ pi-)
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.03, 0.06, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02},
292 
293  // S0 K+ pi+ pi- (S0 K0 pi+ pi-)
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.0,
296  0.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
297 
298  // L K+ pi0 pi0 (L K0 pi0 pi0)
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.03, 0.06, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02},
302 
303  // S0 K+ pi0 pi0 (S0 K0 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.0,
306  0.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01},
307 
308  // S+ K0 pi0 pi0 (S- K+ pi0 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.0,
311  0.02, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01},
312 
313  // L K0 pi+ pi0 (L K+ pi- pi0)
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.04, 0.07, 0.07, 0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02},
317 
318  // S0 K0 pi+ pi0 (S0 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.0,
321  0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01},
322 
323  // S+ K0 pi+ pi- (S- K+ 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.0,
326  0.02, 0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
327 
328  // S+ K+ pi- pi0 (S- K0 pi+ 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.0,
331  0.02, 0.05, 0.08, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
332 
333  // S- K+ pi+ pi0 (S+ K0 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.0,
336  0.03, 0.11, 0.11, 0.09, 0.07, 0.04, 0.03, 0.03, 0.02, 0.02},
337 
338  // p pi0 K+ K- (n pi0 K0 K0b)
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.01, 0.1, 0.11, 0.1, 0.08, 0.07, 0.07, 0.06, 0.05, 0.04},
342 
343  // p pi0 K0 K0b (n pi0 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.02, 0.05, 0.11, 0.1, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04},
347 
348  // p pi+ K0 K- (n pi- K+ K0b)
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.02, 0.05, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04},
352 
353  // p pi- K+ K0b (n pi+ K0 K-)
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.02, 0.05, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04},
357 
358  // n pi+ K+ K- (p pi- K0 K0b)
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.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
362 
363  // n pi+ K0 K0b (p pi- K+ K-)
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.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
367 
368  // n pi0 K+ K0b (p pi0 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.02, 0.05, 0.07, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
372  //
373  // multiplicity 5 (30 channels)
374  //
375  // p pi+ pi+ pi- pi- (n pi+ pi+ pi- pi-)
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.0, 0.04, 0.2, 1.45,
378  2.41, 3.1, 2.7, 2.3, 2.0, 1.54, 1.3, 1.0, 0.83, 0.66},
379 
380 // p pi+ pi- pi0 pi0 (n pi+ pi- 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.0, 0.01, 0.07, 0.48,
383  0.80, 1.1, 0.9, 0.76, 0.63, 0.51, 0.43, 0.33, 0.28, 0.22},
384 
385 // p pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0)
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.0, 0.01, 0.03, 0.21,
388  0.35, 0.5, 0.67, 0.55, 0.47, 0.38, 0.33, 0.23, 0.2, 0.17},
389 
390 // n pi+ pi+ pi- pi0 (p pi+ pi- pi- 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.1, 0.2, 0.41, 1.16,
393  1.75, 1.8, 1.8, 1.38, 1.0, 0.76, 0.58, 0.42, 0.32, 0.24},
394 
395 // n pi+ pi0 pi0 pi0 (p pi- 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.04, 0.21,
398  0.42, 0.6, 0.72, 0.68, 0.57, 0.47, 0.38, 0.28, 0.24, 0.2},
399 
400 // L K+ pi+ pi- pi0 (L K0 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.0, 0.01, 0.03, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
404 
405 // L K+ pi0 pi0 pi0 (L K0 pi0 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.01, 0.03, 0.04, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02},
409 
410 // L K0 pi+ pi+ pi- (L K+ 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 pi+ pi0 pi0 (L K+ pi- 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.03, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
419 
420 // S0 K+ pi+ pi- pi0 (S0 K0 pi+ pi- pi0)
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.01},
424 
425 // S0 K+ pi0 pi0 pi0 (S0 K0 pi0 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.01},
429 
430 // S0 K0 pi+ pi+ pi- (S0 K+ pi+ pi- pi-)
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.01, 0.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
434 
435 // S0 K0 pi+ pi0 pi0 (S0 K+ pi- 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.02, 0.03, 0.03, 0.02, 0.02, 0.02, 0.02, 0.01},
439 
440 // S+ K0 pi+ pi- pi0 (S- K+ pi+ pi- pi0)
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.01, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
444 
445 // S+ K0 pi0 pi0 pi0 (S- K+ pi0 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.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01},
449 
450 // S+ K+ pi+ pi- pi- (S- K0 pi+ pi+ pi-)
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.02, 0.03, 0.02, 0.02, 0.01, 0.01, 0.02},
454 
455 // S+ K+ pi- pi0 pi0 (S- K0 pi+ pi0 pi0)
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.03, 0.05, 0.05, 0.03, 0.03, 0.02, 0.02, 0.02},
459 
460 // S- K+ pi+ pi+ pi- (S+ K0 pi+ pi- pi-)
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.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
464 
465 // S- K+ pi+ pi0 pi0 (S+ K0 pi- pi0 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.02, 0.01, 0.01, 0.01},
469 
470 // p K+ K- pi0 pi0 (n K0 K0b 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.02, 0.06, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
474 
475 // p K0 K0b pi0 pi0 (n K+ K- pi0 pi0)
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.03, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
479 
480 // p K+ K0b pi- pi0 (n K0 K- pi+ pi0)
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.01, 0.03, 0.06, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03},
484 
485 // p K0 K- pi+ pi0 (n K+ K0b pi- pi0)
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.01, 0.02, 0.03, 0.04, 0.03, 0.03, 0.02, 0.02, 0.03},
489 
490 // p K+ K- pi+ pi- (n K0 K0b pi+ pi-)
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.02, 0.06, 0.07, 0.06, 0.05, 0.05, 0.04, 0.04, 0.04},
494 
495 // p K0 K0b pi+ pi- (n K+ K- pi+ pi-)
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.03, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02},
499 
500 // n K+ K- pi+ pi0 (p K0 K0b pi- pi0)
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.02, 0.06, 0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03},
504 
505 // n K0 K0b pi+ pi0 (p K+ K- pi- pi0)
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 K+ K0b pi0 pi0 (p K0 K- pi0 pi0)
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.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
514 
515 // n K0 K- pi+ pi+ (p K+ K0b pi- pi-)
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.06, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
519 
520 // n K+ K0b pi+ pi- (p K0 K- pi+ pi-)
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.06, 0.04, 0.04, 0.03, 0.03, 0.03, 0.02, 0.02},
524  //
525  // multiplicity 6 (6 channels)
526  //
527  // p pi+ pi+ pi- pi- pi0 (n pi+ pi+ pi- 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.0, 0.01, 0.02,
530  0.04, 0.08, 0.12, 0.15, 0.2, 0.19, 0.17, 0.14, 0.12, 0.1},
531 
532  // p pi+ pi- pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0)
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.01, 0.03,
535  0.06, 0.1, 0.14, 0.16, 0.20, 0.19, 0.19, 0.16, 0.14, 0.13},
536 
537  // p pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 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.0, 0.01, 0.02,
540  0.03, 0.05, 0.07, 0.09, 0.11, 0.10, 0.10, 0.09, 0.08, 0.07},
541 
542  // n pi+ pi+ pi+ pi- pi- (p pi+ pi+ pi- pi- pi-)
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.07,
545  0.12, 0.21, 0.28, 0.36, 0.4, 0.42, 0.39, 0.36, 0.31, 0.28},
546 
547  // n pi+ pi+ pi- pi0 pi0 (p pi+ pi- pi- pi0 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.0, 0.01, 0.04,
550  0.07, 0.15, 0.2, 0.24, 0.26, 0.29, 0.27, 0.24, 0.21, 0.18},
551 
552  // n pi+ pi0 pi0 pi0 pi0 (p pi- pi0 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.04,
555  0.09, 0.15, 0.2, 0.25, 0.30, 0.29, 0.28, 0.25, 0.21, 0.19},
556  //
557  // multiplicity 7 (7 channels)
558  //
559  // p pi+ pi+ pi+ pi- pi- pi- (n pi+ pi+ pi+ pi- pi- pi-)
560  {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
561  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.02,
562  0.03, 0.15, 0.4, 0.66, 0.85, 0.82, 0.75, 0.7, 0.65, 0.60},
563 
564  // p pi+ pi+ pi- pi- pi0 pi0 (n pi+ pi+ pi- pi- pi0 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.02,
567  0.03, 0.15, 0.4, 0.66, 0.85, 0.82, 0.75, 0.7, 0.65, 0.60},
568 
569  // p pi+ pi- pi0 pi0 pi0 pi0 (n pi+ pi- pi0 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.0, 0.0,
572  0.01, 0.07, 0.16, 0.28, 0.31, 0.29, 0.25, 0.20, 0.15, 0.10},
573 
574  // p pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 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.0,
577  0.0, 0.02, 0.05, 0.10, 0.12, 0.12, 0.11, 0.10, 0.09, 0.08},
578 
579  // n pi+ pi+ pi+ pi- pi- pi0 (p pi+ pi+ pi- pi- pi- pi0)
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.0,
582  0.01, 0.04, 0.08, 0.14, 0.27, 0.3, 0.27, 0.22, 0.2, 0.18},
583 
584  // n pi+ pi+ pi- pi0 pi0 pi0 (p pi+ pi- pi- pi0 pi0 pi0)
585  {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
586  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01,
587  0.01, 0.04, 0.08, 0.12, 0.16, 0.15, 0.13, 0.11, 0.09, 0.07},
588 
589  // n pi+ pi0 pi0 pi0 pi0 pi0 (p pi- pi0 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.0,
592  0.01, 0.03, 0.05, 0.08, 0.16, 0.19, 0.17, 0.16, 0.14, 0.13},
593  //
594  // multiplicity 8 (8 channels)
595  //
596  // p pi+ pi+ pi+ pi- pi- pi- pi0 (n pi+ pi+ pi+ pi- pi- pi- pi0)
597  {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
598  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
599  0.01, 0.02, 0.04, 0.07, 0.11, 0.13, 0.14, 0.13, 0.12, 0.11},
600 
601  // p pi+ pi+ pi- pi- pi0 pi0 pi0 (n pi+ pi+ pi- pi- pi0 pi0 pi0)
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.01, 0.01, 0.03, 0.05, 0.08, 0.10, 0.10, 0.09, 0.09, 0.08},
605 
606  // p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0 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.02, 0.04, 0.06, 0.07, 0.07, 0.06, 0.06, 0.05},
610 
611  // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 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.01, 0.02, 0.03, 0.04, 0.04, 0.04, 0.04, 0.04},
615 
616  // n pi+ pi+ pi+ pi+ pi- pi- pi- (p pi+ pi+ pi+ pi- pi- pi- pi-)
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.02, 0.04, 0.07, 0.12, 0.19, 0.26, 0.26, 0.24, 0.23, 0.21},
620 
621  // n pi+ pi+ pi+ pi- pi- pi0 pi0 (p pi+ pi+ pi- pi- pi- 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.02, 0.04, 0.07, 0.12, 0.19, 0.26, 0.25, 0.23, 0.23, 0.2},
625 
626  // n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p pi+ pi- pi- 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.01, 0.03, 0.05, 0.08, 0.13, 0.13, 0.12, 0.11, 0.09, 0.08},
630 
631  // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 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.01, 0.02, 0.04, 0.06, 0.10, 0.13, 0.12, 0.11, 0.11, 0.1},
635 
636  //
637  // multiplicity 9 (9 channels)
638  //
639  // p pi+ pi+ pi+ pi+ pi- pi- pi- pi- (n pi+ pi+ pi+ pi+ pi- pi- pi- pi-)
640  {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
641  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
642  0.0, 0.01, 0.02, 0.05, 0.11, 0.16, 0.22, 0.27, 0.27, 0.27},
643 
644  // p pi+ pi+ pi+ pi- pi- pi- pi0 pi0 (n pi+ pi+ pi+ pi- pi- pi- pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
647  0.0, 0.01, 0.02, 0.03, 0.07, 0.10, 0.13, 0.18, 0.18, 0.18},
648 
649  // p pi+ pi+ pi- pi- pi0 pi0 pi0 pi0 (n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
652  0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.09, 0.11, 0.11, 0.11},
653 
654  // p pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 (n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
657  0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.07, 0.07, 0.07},
658 
659  // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
662  0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.05, 0.07, 0.07, 0.07},
663 
664  // n pi+ pi+ pi+ pi+ pi- pi- pi- pi0 (p pi+ pi+ pi+ pi- pi- pi- pi- pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
667  0.0, 0.0, 0.01, 0.03, 0.06, 0.11, 0.14, 0.16, 0.16, 0.16},
668 
669  // n pi+ pi+ pi+ pi- pi- pi0 pi0 pi0 (p pi+ pi+ pi- pi- pi- pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
672  0.0, 0.0, 0.01, 0.02, 0.05, 0.07, 0.08, 0.1, 0.1, 0.1},
673 
674  // n pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0 (p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
677  0.0, 0.0, 0.0, 0.01, 0.02, 0.04, 0.06, 0.06, 0.06, 0.06},
678 
679  // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0)
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
682  0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.03, 0.04, 0.04, 0.04}};
683 }
684 
685 // Initialize both |T Tz> = |1/2 1/2> channels, using pizP cross-section table
686 
687 const G4CascadePiZeroPChannelData::data_t
688 G4CascadePiZeroPChannelData::data(pizP2bfs, pizP3bfs, pizP4bfs,
689  pizP5bfs, pizP6bfs, pizP7bfs,
690  pizP8bfs, pizP9bfs, pizPCrossSections,
691  pizPtotXSec, pi0*pro, "PiZeroP");
692 
693 const G4CascadePiZeroNChannelData::data_t
694 G4CascadePiZeroNChannelData::data(pizN2bfs, pizN3bfs, pizN4bfs,
695  pizN5bfs, pizN6bfs, pizN7bfs,
696  pizN8bfs, pizN9bfs, pizPCrossSections,
697  pizPtotXSec, pi0*neu, "PiZeroN");
G4InuclParticleNames.hh
G4int
int G4int
Definition: G4Types.hh:78
G4CascadePiZeroPChannel.hh
G4CascadePiZeroNChannel.hh
G4double
double G4double
Definition: G4Types.hh:76