ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4CascadeT33piNChannel.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4CascadeT33piNChannel.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 
28 #include "G4CascadePiPlusPChannel.hh"
29 #include "G4CascadePiMinusNChannel.hh"
30 #include "G4InuclParticleNames.hh"
31 using namespace G4InuclParticleNames;
32 
33 namespace {
34  // pi+ p : Outgoing particle types of a given multiplicity
35 
36  static const G4int pipP2bfs[2][2] =
37  {{pip,pro}, {kpl,sp}};
38 
39  static const G4int pipP3bfs[7][3] =
40  {{pip,pro,pi0}, {pip,neu,pip}, {pi0,sp,kpl}, {pip,sp,k0},
41  {pip,s0,kpl}, {pip,lam,kpl}, {kpl,pro,k0b}};
42 
43  static const G4int pipP4bfs[15][4] =
44  {{pip,pro,pip,pim}, {pip,pro,pi0,pi0}, {pip,neu,pip,pi0},
45  {pip,sp,kpl,pim}, {pi0,sp,kpl,pi0}, {pip,sp,k0,pi0},
46  {pip,s0,k0,pip}, {pip,s0,kpl,pi0}, {pip,lam,kpl,pi0},
47  {pip,lam,k0,pip}, {pip,sm,kpl,pip}, {pip,pro,kpl,kmi},
48  {pip,pro,k0,k0b}, {pi0,pro,kpl,k0b}, {pip,neu,kpl,k0b}};
49 
50  static const G4int pipP5bfs[24][5] =
51  {{pip,pro,pip,pim,pi0}, {pip,pro,pi0,pi0,pi0}, {pip,neu,pip,pip,pim},
52  {pip,neu,pip,pi0,pi0}, {pip,sp,kpl,pim,pi0}, {pi0,sp,kpl,pi0,pi0},
53  {pip,sp,k0,pip,pim}, {pip,sp,k0,pi0,pi0}, {pip,lam,k0,pip,pi0},
54  {pip,lam,kpl,pip,pim}, {pip,lam,kpl,pi0,pi0}, {pip,s0,kpl,pip,pim},
55  {pip,s0,kpl,pi0,pi0}, {pip,s0,k0,pip,pi0}, {pip,sm,kpl,pip,pi0},
56  {pip,sm,k0,pip,pip}, {pip,pro,pim,kpl,k0b}, {pip,pro,pip,k0,kmi},
57  {pip,pro,pi0,kpl,kmi}, {pip,pro,pi0,k0,k0b}, {pi0,pro,pi0,kpl,k0b},
58  {pip,neu,pip,kpl,kmi}, {pip,neu,pip,k0,k0b}, {pip,neu,pi0,kpl,k0b}};
59 
60  static const G4int pipP6bfs[33][6] =
61  {{pip,pro,pip,pip,pim,pim}, {pip,pro,pip,pim,pi0,pi0},
62  {pip,pro,pi0,pi0,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0},
63  {pip,neu,pip,pip,pim,pi0}, {pip,pip,pim,pi0,kpl,lam},
64  {pip,pi0,pi0,pi0,kpl,lam}, {pip,pip,pi0,pi0,k0,lam},
65  {pip,pip,pip,pim,k0,lam}, {pip,pip,pim,pim,kpl,sp},
66  {pip,pim,pi0,pi0,kpl,sp}, {pi0,pi0,pi0,pi0,kpl,sp},
67  {pip,pip,pim,pi0,k0,sp}, {pip,pi0,pi0,pi0,k0,sp},
68  {pip,pip,pim,pi0,kpl,s0}, {pip,pi0,pi0,pi0,kpl,s0},
69  {pip,pip,pip,pim,k0,s0}, {pip,pip,pi0,pi0,k0,s0},
70  {pip,pip,pi0,pi0,kpl,sm}, {pip,pip,pip,pim,kpl,sm},
71  {pip,pip,pip,pi0,k0,sm}, {pip,pro,pim,pi0,kpl,k0b},
72  {pi0,pro,pi0,pi0,kpl,k0b}, {pip,pro,pi0,pi0,kpl,kmi},
73  {pip,pro,pip,pim,kpl,kmi}, {pip,pro,pi0,pi0,k0,k0b},
74  {pip,pro,pip,pim,k0,k0b}, {pip,pro,pip,pi0,kmi,k0},
75  {pip,neu,pi0,pi0,kpl,k0b}, {pip,neu,pip,pim,kpl,k0b},
76  {pip,neu,pip,pi0,kpl,kmi}, {pip,neu,pip,pim,kpl,k0b},
77  {pip,neu,pip,pip,k0,kmi}};
78 
79  static const G4int pipP7bfs[41][7] =
80  {{pip,pro,pip,pip,pim,pim,pi0}, {pip,pro,pip,pim,pi0,pi0,pi0},
81  {pip,pro,pi0,pi0,pi0,pi0,pi0}, {pip,neu,pip,pip,pip,pim,pim},
82  {pip,neu,pip,pip,pim,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0,pi0},
83  {pip,pip,pim,pi0,pi0,kpl,lam}, {pip,pip,pip,pim,pim,kpl,lam},
84  {pip,pi0,pi0,pi0,pi0,kpl,lam}, {pip,pip,pip,pim,pi0,k0,lam},
85  {pip,pip,pi0,pi0,pi0,k0,lam}, {pip,pip,pim,pim,pi0,kpl,sp},
86  {pip,pim,pi0,pi0,pi0,kpl,sp}, {pip,pip,pim,pi0,pi0,k0,sp},
87  {pip,pip,pip,pim,pim,k0,sp}, {pip,pi0,pi0,pi0,pi0,k0,sp},
88  {pip,pip,pim,pi0,pi0,kpl,s0}, {pip,pip,pip,pim,pim,kpl,s0},
89  {pip,pi0,pi0,pi0,pi0,kpl,s0}, {pip,pip,pip,pim,pi0,k0,s0},
90  {pip,pip,pi0,pi0,pi0,k0,s0}, {pip,pip,pip,pim,pi0,kpl,sm},
91  {pip,pip,pi0,pi0,pi0,kpl,sm}, {pip,pip,pip,pip,pim,k0,sm},
92  {pip,pip,pip,pi0,pi0,k0,sm}, {pip,pro,pim,pi0,pi0,kpl,k0b},
93  {pip,pro,pip,pim,pim,kpl,k0b}, {pi0,pro,pi0,pi0,pi0,kpl,k0b},
94  {pip,pro,pip,pim,pi0,kpl,kmi}, {pip,pro,pi0,pi0,pi0,kpl,kmi},
95  {pip,pro,pip,pim,pi0,k0,k0b}, {pip,pro,pi0,pi0,pi0,k0,k0b},
96  {pip,pro,pip,pi0,pi0,kmi,k0}, {pip,pro,pip,pip,pim,kmi,k0},
97  {pip,neu,pip,pim,pi0,kpl,k0b}, {pip,neu,pi0,pi0,pi0,kpl,k0b},
98  {pip,neu,pip,pi0,pi0,kpl,kmi}, {pip,neu,pip,pip,pim,kpl,kmi},
99  {pip,neu,pip,pi0,pi0,k0,k0b}, {pip,neu,pip,pip,pim,k0,k0b},
100  {pip,neu,pip,pip,pi0,kmi,k0}};
101 
102  static const G4int pipP8bfs[47][8] =
103  {{pip,pro,pip,pip,pip,pim,pim,pim}, {pip,pro,pip,pip,pim,pim,pi0,pi0},
104  {pip,pro,pip,pim,pi0,pi0,pi0,pi0}, {pip,pro,pi0,pi0,pi0,pi0,pi0,pi0},
105  {pip,neu,pip,pip,pip,pim,pim,pi0}, {pip,neu,pip,pip,pim,pi0,pi0,pi0},
106  {pip,neu,pip,pi0,pi0,pi0,pi0,pi0}, {pip,pip,pip,pim,pim,pi0,kpl,lam},
107  {pip,pi0,pi0,pi0,pi0,pi0,kpl,lam}, {pip,pip,pip,pim,pi0,pi0,k0,lam},
108  {pip,pip,pi0,pi0,pi0,pi0,k0,lam}, {pip,pip,pip,pip,pim,pim,k0,lam},
109  {pip,pip,pim,pim,pi0,pi0,kpl,sp}, {pip,pim,pi0,pi0,pi0,pi0,kpl,sp},
110  {pip,pip,pip,pim,pim,pim,kpl,sp}, {pip,pip,pip,pim,pim,pi0,k0,sp},
111  {pip,pip,pim,pi0,pi0,pi0,k0,sp}, {pip,pip,pip,pim,pim,pi0,kpl,s0},
112  {pip,pi0,pi0,pi0,pi0,pi0,kpl,s0}, {pip,pip,pip,pim,pi0,pi0,k0,s0},
113  {pip,pip,pi0,pi0,pi0,pi0,k0,s0}, {pip,pip,pip,pip,pim,pim,k0,s0},
114  {pip,pip,pip,pim,pi0,pi0,kpl,sm}, {pip,pip,pi0,pi0,pi0,pi0,kpl,sm},
115  {pip,pip,pip,pip,pim,pim,kpl,sm}, {pip,pip,pip,pip,pim,pi0,k0,sm},
116  {pip,pip,pip,pi0,pi0,pi0,k0,sm}, {pip,pro,pip,pim,pim,pi0,kpl,k0b},
117  {pip,pro,pim,pi0,pi0,pi0,kpl,k0b}, {pi0,pro,pi0,pi0,pi0,pi0,kpl,k0b},
118  {pip,pro,pip,pim,pi0,pi0,kpl,kmi}, {pip,pro,pip,pip,pim,pim,kpl,kmi},
119  {pip,pro,pi0,pi0,pi0,pi0,kpl,kmi}, {pip,pro,pip,pim,pi0,pi0,k0,k0b},
120  {pip,pro,pip,pip,pim,pim,k0,k0b}, {pip,pro,pi0,pi0,pi0,pi0,k0,k0b},
121  {pip,pro,pip,pip,pim,pi0,k0,kmi}, {pip,pro,pip,pi0,pi0,pi0,k0,kmi},
122  {pip,neu,pip,pip,pim,pim,kpl,k0b}, {pip,neu,pip,pim,pi0,pi0,kpl,k0b},
123  {pip,neu,pi0,pi0,pi0,pi0,kpl,k0b}, {pip,neu,pip,pip,pim,pi0,kpl,kmi},
124  {pip,neu,pip,pi0,pi0,pi0,kpl,kmi}, {pip,neu,pip,pip,pim,pi0,k0,k0b},
125  {pip,neu,pip,pi0,pi0,pi0,k0,k0b}, {pip,neu,pip,pip,pi0,pi0,k0,kmi},
126  {pip,neu,pip,pip,pip,pim,k0,kmi}};
127 
128  static const G4int pipP9bfs[55][9] =
129  {{pip,pro,pip,pip,pip,pim,pim,pim,pi0}, {pip,pro,pip,pip,pim,pim,pi0,pi0,pi0},
130  {pip,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pip,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
131  {pip,neu,pip,pip,pip,pip,pim,pim,pim}, {pip,neu,pip,pip,pip,pim,pim,pi0,pi0},
132  {pip,neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {pip,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
133  {pip,pip,pip,pim,pim,pi0,pi0,kpl,lam}, {pip,pip,pim,pi0,pi0,pi0,pi0,kpl,lam},
134  {pip,pip,pip,pip,pim,pim,pim,kpl,lam}, {pip,pi0,pi0,pi0,pi0,pi0,pi0,kpl,lam},
135  {pip,pip,pip,pim,pi0,pi0,pi0,k0,lam}, {pip,pip,pip,pip,pim,pim,pi0,k0,lam},
136  {pip,pip,pi0,pi0,pi0,pi0,pi0,k0,lam}, {pip,pip,pim,pim,pi0,pi0,pi0,kpl,sp},
137  {pip,pip,pip,pim,pim,pim,pi0,kpl,sp}, {pip,pim,pi0,pi0,pi0,pi0,pi0,kpl,sp},
138  {pip,pip,pip,pim,pim,pi0,pi0,k0,sp}, {pip,pip,pim,pi0,pi0,pi0,pi0,k0,sp},
139  {pip,pip,pip,pip,pim,pim,pim,k0,sp}, {pip,pip,pip,pim,pim,pi0,pi0,kpl,s0},
140  {pip,pip,pim,pi0,pi0,pi0,pi0,kpl,s0}, {pip,pip,pip,pip,pim,pim,pim,kpl,s0},
141  {pip,pip,pip,pim,pi0,pi0,pi0,k0,s0}, {pip,pip,pip,pip,pim,pim,pi0,k0,s0},
142  {pip,pip,pi0,pi0,pi0,pi0,pi0,k0,s0}, {pip,pip,pip,pim,pi0,pi0,pi0,kpl,sm},
143  {pip,pip,pip,pip,pim,pim,pi0,kpl,sm}, {pip,pip,pi0,pi0,pi0,pi0,pi0,kpl,sm},
144  {pip,pip,pip,pip,pim,pi0,pi0,k0,sm}, {pip,pip,pip,pi0,pi0,pi0,pi0,k0,sm},
145  {pip,pro,pip,pim,pim,pi0,pi0,kpl,k0b}, {pip,pro,pim,pi0,pi0,pi0,pi0,kpl,k0b},
146  {pip,pro,pip,pip,pim,pim,pim,kpl,k0b}, {pip,pro,pip,pim,pi0,pi0,pi0,kpl,kmi},
147  {pip,pro,pip,pip,pim,pim,pi0,kpl,kmi}, {pip,pro,pi0,pi0,pi0,pi0,pi0,kpl,kmi},
148  {pip,pro,pip,pim,pi0,pi0,pi0,k0,k0b}, {pip,pro,pip,pip,pim,pim,pi0,k0,k0b},
149  {pip,pro,pi0,pi0,pi0,pi0,pi0,k0,k0b}, {pip,pro,pip,pip,pim,pi0,pi0,k0,kmi},
150  {pip,pro,pip,pi0,pi0,pi0,pi0,k0,kmi}, {pip,pro,pip,pip,pip,pim,pim,k0,kmi},
151  {pip,neu,pip,pim,pi0,pi0,pi0,kpl,k0b}, {pip,neu,pip,pip,pim,pim,pi0,kpl,k0b},
152  {pip,neu,pi0,pi0,pi0,pi0,pi0,kpl,k0b}, {pip,neu,pip,pip,pim,pi0,pi0,kpl,kmi},
153  {pip,neu,pip,pi0,pi0,pi0,pi0,kpl,kmi}, {pip,neu,pip,pip,pip,pim,pim,kpl,kmi},
154  {pip,neu,pip,pip,pim,pi0,pi0,k0,k0b}, {pip,neu,pip,pi0,pi0,pi0,pi0,k0,k0b},
155  {pip,neu,pip,pip,pip,pim,pim,k0,k0b}, {pip,neu,pip,pip,pip,pim,pi0,kmi,k0},
156  {pip,neu,pip,pip,pi0,pi0,pi0,kmi,k0}};
157 }
158 
159 namespace {
160  // pi- n : Outgoing particle types of a given multiplicity
161 
162  static const G4int pimN2bfs[2][2] =
163  {{pim,neu}, {k0,sm}};
164 
165  static const G4int pimN3bfs[7][3] =
166  {{pim,neu,pi0}, {pim,pro,pim}, {pi0,sm,k0}, {pim,sm,kpl},
167  {pim,s0,k0}, {pim,lam,k0}, {k0,neu,kmi}};
168 
169  static const G4int pimN4bfs[15][4] =
170  {{pim,neu,pip,pim},{pim,neu,pi0,pi0},{pim,pro,pim,pi0},
171  {pim,sm,k0,pip}, {pi0,sm,k0,pi0}, {pim,sm,kpl,pi0},
172  {pim,s0,kpl,pim}, {pim,s0,k0,pi0}, {pim,lam,k0,pi0},
173  {pim,lam,kpl,pim}, {pim,sp,k0,pim}, {pim,neu,k0,k0b},
174  {pim,neu,kpl,kmi}, {pi0,neu,k0,kmi}, {pim,pro,k0,kmi}};
175 
176  static const G4int pimN5bfs[24][5] =
177  {{pim,neu,pip,pim,pi0}, {pim,neu,pi0,pi0,pi0}, {pim,pro,pip,pim,pim},
178  {pim,pro,pim,pi0,pi0}, {pim,sm,k0,pip,pi0}, {pi0,sm,k0,pi0,pi0},
179  {pim,sm,kpl,pip,pim}, {pim,sm,kpl,pi0,pi0}, {pim,lam,kpl,pim,pi0},
180  {pim,lam,k0,pip,pim}, {pim,lam,k0,pi0,pi0}, {pim,s0,k0,pip,pim},
181  {pim,s0,k0,pi0,pi0}, {pim,s0,kpl,pim,pi0}, {pim,sp,k0,pim,pi0},
182  {pim,sp,kpl,pim,pim}, {pim,neu,pip,k0,kmi}, {pim,neu,pim,kpl,k0b},
183  {pim,neu,pi0,k0,k0b}, {pim,neu,pi0,kpl,kmi}, {pi0,neu,pi0,k0,kmi},
184  {pim,pro,pim,k0,k0b}, {pim,pro,pim,kpl,kmi}, {pim,pro,pi0,k0,kmi}};
185 
186  static const G4int pimN6bfs[33][6] =
187  {{pim,neu,pip,pip,pim,pim}, {pim,neu,pip,pim,pi0,pi0},
188  {pim,neu,pi0,pi0,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0},
189  {pim,pro,pip,pim,pim,pi0}, {pim,pip,pim,pi0,k0,lam},
190  {pim,pi0,pi0,pi0,k0,lam}, {pim,pim,pi0,pi0,kpl,lam},
191  {pim,pip,pim,pim,kpl,lam}, {pim,pip,pip,pim,k0,sm},
192  {pim,pip,pi0,pi0,k0,sm}, {pi0,pi0,pi0,pi0,k0,sm},
193  {pim,pip,pim,pi0,kpl,sm}, {pim,pi0,pi0,pi0,kpl,sm},
194  {pim,pip,pim,pi0,k0,s0}, {pim,pi0,pi0,pi0,k0,s0},
195  {pim,pip,pim,pim,kpl,s0}, {pim,pim,pi0,pi0,kpl,s0},
196  {pim,pim,pi0,pi0,k0,sp}, {pim,pip,pim,pim,k0,sp},
197  {pim,pim,pim,pi0,kpl,sp}, {pim,neu,pip,pi0,kmi,k0},
198  {pi0,neu,pi0,pi0,kmi,k0}, {pim,neu,pi0,pi0,k0,k0b},
199  {pim,neu,pip,pim,k0,k0b}, {pim,neu,pi0,pi0,kpl,kmi},
200  {pim,neu,pip,pim,kpl,kmi}, {pim,neu,pim,pi0,kpl,k0b},
201  {pim,pro,pi0,pi0,kmi,k0}, {pim,pro,pip,pim,kmi,k0},
202  {pim,pro,pim,pi0,k0,k0b}, {pim,pro,pip,pim,kmi,k0},
203  {pim,pro,pim,pim,kpl,k0b}};
204 
205  static const G4int pimN7bfs[41][7] =
206  {{pim,neu,pip,pip,pim,pim,pi0}, {pim,neu,pip,pim,pi0,pi0,pi0},
207  {pim,neu,pi0,pi0,pi0,pi0,pi0}, {pim,pro,pip,pip,pim,pim,pim},
208  {pim,pro,pip,pim,pim,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0,pi0},
209  {pim,pip,pim,pi0,pi0,k0,lam}, {pim,pip,pip,pim,pim,k0,lam},
210  {pim,pi0,pi0,pi0,pi0,k0,lam}, {pim,pip,pim,pim,pi0,kpl,lam},
211  {pim,pim,pi0,pi0,pi0,kpl,lam}, {pim,pip,pip,pim,pi0,k0,sm},
212  {pim,pip,pi0,pi0,pi0,k0,sm}, {pim,pip,pim,pi0,pi0,kpl,sm},
213  {pim,pip,pip,pim,pim,kpl,sm}, {pim,pi0,pi0,pi0,pi0,kpl,sm},
214  {pim,pip,pim,pi0,pi0,k0,s0}, {pim,pip,pip,pim,pim,k0,s0},
215  {pim,pi0,pi0,pi0,pi0,k0,s0}, {pim,pip,pim,pim,pi0,kpl,s0},
216  {pim,pim,pi0,pi0,pi0,kpl,s0}, {pim,pip,pim,pim,pi0,k0,sp},
217  {pim,pim,pi0,pi0,pi0,k0,sp}, {pim,pip,pim,pim,pim,kpl,sp},
218  {pim,pim,pim,pi0,pi0,kpl,sp}, {pim,neu,pip,pi0,pi0,kmi,k0},
219  {pim,neu,pip,pip,pim,kmi,k0}, {pi0,neu,pi0,pi0,pi0,kmi,k0},
220  {pim,neu,pip,pim,pi0,k0,k0b}, {pim,neu,pi0,pi0,pi0,k0,k0b},
221  {pim,neu,pip,pim,pi0,kpl,kmi}, {pim,neu,pi0,pi0,pi0,kpl,kmi},
222  {pim,neu,pim,pi0,pi0,kpl,k0b}, {pim,neu,pip,pim,pim,kpl,k0b},
223  {pim,pro,pip,pim,pi0,kmi,k0}, {pim,pro,pi0,pi0,pi0,kmi,k0},
224  {pim,pro,pim,pi0,pi0,k0,k0b}, {pim,pro,pip,pim,pim,k0,k0b},
225  {pim,pro,pim,pi0,pi0,kpl,kmi}, {pim,pro,pip,pim,pim,kpl,kmi},
226  {pim,pro,pim,pim,pi0,kpl,k0b}};
227 
228  static const G4int pimN8bfs[47][8] =
229  {{pim,neu,pip,pip,pip,pim,pim,pim}, {pim,neu,pip,pip,pim,pim,pi0,pi0},
230  {pim,neu,pip,pim,pi0,pi0,pi0,pi0}, {pim,neu,pi0,pi0,pi0,pi0,pi0,pi0},
231  {pim,pro,pip,pip,pim,pim,pim,pi0}, {pim,pro,pip,pim,pim,pi0,pi0,pi0},
232  {pim,pro,pim,pi0,pi0,pi0,pi0,pi0}, {pim,pip,pip,pim,pim,pi0,k0,lam},
233  {pim,pi0,pi0,pi0,pi0,pi0,k0,lam}, {pim,pip,pim,pim,pi0,pi0,kpl,lam},
234  {pim,pim,pi0,pi0,pi0,pi0,kpl,lam}, {pim,pip,pip,pim,pim,pim,kpl,lam},
235  {pim,pip,pip,pim,pi0,pi0,k0,sm}, {pim,pip,pi0,pi0,pi0,pi0,k0,sm},
236  {pim,pip,pip,pip,pim,pim,k0,sm}, {pim,pip,pip,pim,pim,pi0,kpl,sm},
237  {pim,pip,pim,pi0,pi0,pi0,kpl,sm}, {pim,pip,pip,pim,pim,pi0,k0,s0},
238  {pim,pi0,pi0,pi0,pi0,pi0,k0,s0}, {pim,pip,pim,pim,pi0,pi0,kpl,s0},
239  {pim,pim,pi0,pi0,pi0,pi0,kpl,s0}, {pim,pip,pip,pim,pim,pim,kpl,s0},
240  {pim,pip,pim,pim,pi0,pi0,k0,sp}, {pim,pim,pi0,pi0,pi0,pi0,k0,sp},
241  {pim,pip,pip,pim,pim,pim,k0,sp}, {pim,pip,pim,pim,pim,pi0,kpl,sp},
242  {pim,pim,pim,pi0,pi0,pi0,kpl,sp}, {pim,neu,pip,pip,pim,pi0,kmi,k0},
243  {pim,neu,pip,pi0,pi0,pi0,kmi,k0}, {pi0,neu,pi0,pi0,pi0,pi0,kmi,k0},
244  {pim,neu,pip,pim,pi0,pi0,k0,k0b}, {pim,neu,pip,pip,pim,pim,k0,k0b},
245  {pim,neu,pi0,pi0,pi0,pi0,k0,k0b}, {pim,neu,pip,pim,pi0,pi0,kpl,kmi},
246  {pim,neu,pip,pip,pim,pim,kpl,kmi}, {pim,neu,pi0,pi0,pi0,pi0,kpl,kmi},
247  {pim,neu,pip,pim,pim,pi0,kpl,k0b}, {pim,neu,pim,pi0,pi0,pi0,kpl,k0b},
248  {pim,pro,pip,pip,pim,pim,kmi,k0}, {pim,pro,pip,pim,pi0,pi0,kmi,k0},
249  {pim,pro,pi0,pi0,pi0,pi0,kmi,k0}, {pim,pro,pip,pim,pim,pi0,k0,k0b},
250  {pim,pro,pim,pi0,pi0,pi0,k0,k0b}, {pim,pro,pip,pim,pim,pi0,kpl,kmi},
251  {pim,pro,pim,pi0,pi0,pi0,kpl,kmi}, {pim,pro,pim,pim,pi0,pi0,kpl,k0b},
252  {pim,pro,pip,pim,pim,pim,kpl,k0b}};
253 
254  static const G4int pimN9bfs[55][9] =
255  {{pim,neu,pip,pip,pip,pim,pim,pim,pi0}, {pim,neu,pip,pip,pim,pim,pi0,pi0,pi0},
256  {pim,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pim,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
257  {pim,pro,pip,pip,pip,pim,pim,pim,pim}, {pim,pro,pip,pip,pim,pim,pim,pi0,pi0},
258  {pim,pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pim,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
259  {pim,pip,pip,pim,pim,pi0,pi0,k0,lam}, {pim,pip,pim,pi0,pi0,pi0,pi0,k0,lam},
260  {pim,pip,pip,pip,pim,pim,pim,k0,lam}, {pim,pi0,pi0,pi0,pi0,pi0,pi0,k0,lam},
261  {pim,pip,pim,pim,pi0,pi0,pi0,kpl,lam}, {pim,pip,pip,pim,pim,pim,pi0,kpl,lam},
262  {pim,pim,pi0,pi0,pi0,pi0,pi0,kpl,lam}, {pim,pip,pip,pim,pi0,pi0,pi0,k0,sm},
263  {pim,pip,pip,pip,pim,pim,pi0,k0,sm}, {pim,pip,pi0,pi0,pi0,pi0,pi0,k0,sm},
264  {pim,pip,pip,pim,pim,pi0,pi0,kpl,sm}, {pim,pip,pim,pi0,pi0,pi0,pi0,kpl,sm},
265  {pim,pip,pip,pip,pim,pim,pim,kpl,sm}, {pim,pip,pip,pim,pim,pi0,pi0,k0,s0},
266  {pim,pip,pim,pi0,pi0,pi0,pi0,k0,s0}, {pim,pip,pip,pip,pim,pim,pim,k0,s0},
267  {pim,pip,pim,pim,pi0,pi0,pi0,kpl,s0}, {pim,pip,pip,pim,pim,pim,pi0,kpl,s0},
268  {pim,pim,pi0,pi0,pi0,pi0,pi0,kpl,s0}, {pim,pip,pim,pim,pi0,pi0,pi0,k0,sp},
269  {pim,pip,pip,pim,pim,pim,pi0,k0,sp}, {pim,pim,pi0,pi0,pi0,pi0,pi0,k0,sp},
270  {pim,pip,pim,pim,pim,pi0,pi0,kpl,sp}, {pim,pim,pim,pi0,pi0,pi0,pi0,kpl,sp},
271  {pim,neu,pip,pip,pim,pi0,pi0,kmi,k0}, {pim,neu,pip,pi0,pi0,pi0,pi0,kmi,k0},
272  {pim,neu,pip,pip,pip,pim,pim,kmi,k0}, {pim,neu,pip,pim,pi0,pi0,pi0,k0,k0b},
273  {pim,neu,pip,pip,pim,pim,pi0,k0,k0b}, {pim,neu,pi0,pi0,pi0,pi0,pi0,k0,k0b},
274  {pim,neu,pip,pim,pi0,pi0,pi0,kpl,kmi}, {pim,neu,pip,pip,pim,pim,pi0,kpl,kmi},
275  {pim,neu,pi0,pi0,pi0,pi0,pi0,kpl,kmi}, {pim,neu,pip,pim,pim,pi0,pi0,kpl,k0b},
276  {pim,neu,pim,pi0,pi0,pi0,pi0,kpl,k0b}, {pim,neu,pip,pip,pim,pim,pim,kpl,k0b},
277  {pim,pro,pip,pim,pi0,pi0,pi0,kmi,k0}, {pim,pro,pip,pip,pim,pim,pi0,kmi,k0},
278  {pim,pro,pi0,pi0,pi0,pi0,pi0,kmi,k0}, {pim,pro,pip,pim,pim,pi0,pi0,k0,k0b},
279  {pim,pro,pim,pi0,pi0,pi0,pi0,k0,k0b}, {pim,pro,pip,pip,pim,pim,pim,k0,k0b},
280  {pim,pro,pip,pim,pim,pi0,pi0,kpl,kmi}, {pim,pro,pim,pi0,pi0,pi0,pi0,kpl,kmi},
281  {pim,pro,pip,pip,pim,pim,pim,kpl,kmi}, {pim,pro,pip,pim,pim,pim,pi0,kpl,k0b},
282  {pim,pro,pim,pim,pi0,pi0,pi0,kpl,k0b}};
283 }
284 
285 namespace {
286  // Total pi+ p cross sections as a function of kinetic energy
287  // New cs after 7, 8, 9-body strange channel additions (27 Oct 17)
288  static const G4double pipPtotXSec[30] =
289  { 0.0, 1.2, 2.5, 3.8, 5.0, 7.0, 9.0, 15.0, 30.0, 64.0,
290  130.0, 190.0, 130.01, 55.85, 28.0, 17.14, 19.27, 27.58, 40.07, 32.51,
291  30.46, 28.61, 27.26, 25.81, 25.0, 24.38, 24.0, 23.5, 23.0, 23.0};
292 
293  // pi+ p cross sections as functions of kinetic energy and multiplicity
294  static const G4double pipPCrossSections[224][30] = {
295  //
296  // multiplicity 2 (2 channels)
297  //
298  // p pi+ (n pi-)
299  { 0.0, 1.2, 2.5, 3.8, 5.0, 7.00, 9.00, 15.0, 30.0, 64.0,
300  130.0, 190.0, 130.0, 55.7, 27.2, 13.95, 8.38, 12.98, 18.53, 11.81,
301  9.4, 7.7, 6.8, 5.9, 5.2, 4.7, 4.429, 4.275, 4.159, 4.109},
302 
303  // S+ K+ (S- K0)
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.16, 0.62, 0.34,
306  0.25, 0.2, 0.1, 0.06, 0.036, 0.025, 0.018, 0.014, 0.012, 0.011},
307  //
308  // multiplicity 3 (7 channels)
309  //
310  // p pi+ pi0 (n pi- pi0)
311  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
312  0.0, 0.0, 0.007, 0.11, 0.6, 2.39, 8.67, 9.99, 12.36, 6.66,
313  4.4, 2.85, 1.85, 1.25, 0.83, 0.57, 0.432, 0.326, 0.251, 0.194},
314 
315  // n pi+ pi+ (p pi- pi-)
316  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
317  0.0, 0.0, 0.003, 0.04, 0.2, 0.6, 1.48, 2.5, 3.71, 3.22,
318  2.7, 2.1, 1.47, 1.05, 0.73, 0.52, 0.391, 0.288, 0.219, 0.171},
319 
320  // S+ K+ pi0 (S- K0 pi0)
321  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
322  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
323  0.15, 0.12, 0.09, 0.06, 0.042, 0.03, 0.023, 0.016, 0.013, 0.011},
324 
325  // S+ K0 pi+
326  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
327  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
328  0.15, 0.12, 0.09, 0.065, 0.043, 0.03, 0.023, 0.016, 0.013, 0.011},
329 
330  // S0 K+ pi+
331  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
332  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
333  0.15, 0.12, 0.09, 0.065, 0.043, 0.03, 0.023, 0.016, 0.013, 0.011},
334 
335  // L K+ pi+
336  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
337  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.13,
338  0.15, 0.12, 0.09, 0.065, 0.043, 0.03, 0.023, 0.016, 0.013, 0.011},
339 
340  // p K+ K0bar
341  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
342  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04,
343  0.07, 0.062, 0.045, 0.035, 0.026, 0.021, 0.018, 0.015, 0.013, 0.011},
344  //
345  // multiplicity 4 (15 channels)
346  //
347  // p pi+ pi+ pi-
348  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
349  0.0, 0.0, 0.0, 0.0, 0.0, 0.05, 0.206, 0.736, 2.26, 3.43,
350  3.5, 3.05, 2.6, 2.25, 1.88, 1.6, 1.442, 1.23, 1.073, 0.936},
351 
352  // p pi+ pi0 pi0
353  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
354  0.0, 0.0, 0.0, 0.0, 0.0, 0.075, 0.265, 0.575, 1.045, 1.835,
355  2.125, 1.7, 1.135, 0.675, 0.43, 0.27, 0.185, 0.125, 0.092, 0.072},
356 
357  // n pi+ pi+ pi0
358  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
359  0.0, 0.0, 0.0, 0.0, 0.0, 0.075, 0.265, 0.575, 1.045, 1.835,
360  2.125, 1.7, 1.135, 0.675, 0.43, 0.27, 0.185, 0.125, 0.092, 0.072},
361 
362  // S+ K+ pi+ pi-
363  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.03, 0.08, 0.09, 0.06, 0.037, 0.025, 0.016, 0.011, 0.008},
366 
367  // S+ K+ pi0 pi0
368  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.03, 0.08, 0.09, 0.055, 0.035, 0.025, 0.016, 0.011, 0.008},
371 
372  // S+ K0 pi+ pi0
373  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
374  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
375  0.02, 0.06, 0.07, 0.053, 0.04, 0.031, 0.025, 0.018, 0.014, 0.011},
376 
377  // S0 K0 pi+ pi+
378  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
379  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
380  0.0, 0.007, 0.018, 0.02, 0.018, 0.0155, 0.014, 0.012, 0.011, 0.01},
381 
382  // S0 K+ pi+ pi0
383  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
384  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
385  0.0, 0.007, 0.018, 0.02, 0.018, 0.0155, 0.014, 0.012, 0.011, 0.01},
386 
387  // L K+ pi+ pi0
388  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
389  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
390  0.0, 0.013, 0.04, 0.04, 0.031, 0.025, 0.021, 0.017, 0.014, 0.011},
391 
392  // L K0 pi+ pi+
393  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
394  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
395  0.0, 0.013, 0.04, 0.04, 0.031, 0.025, 0.021, 0.017, 0.014, 0.011},
396 
397  // S- K+ pi+ pi+
398  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
399  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
400  0.02, 0.16, 0.18, 0.14, 0.1, 0.065, 0.05, 0.04, 0.033, 0.026},
401 
402  // p pi+ K+ K-
403  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
404  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
405  0.02, 0.13, 0.16, 0.12, 0.09, 0.073, 0.063, 0.056, 0.049, 0.043},
406 
407  // p pi+ K0 K0bar
408  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
409  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
410  0.01, 0.06, 0.15, 0.14, 0.11, 0.085, 0.067, 0.053, 0.049, 0.046},
411 
412  // p pi0 K+ K0bar
413  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
414  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
415  0.01, 0.05, 0.09, 0.085, 0.08, 0.07, 0.064, 0.053, 0.049, 0.046},
416 
417  // n pi+ K+ K0bar
418  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
419  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
420  0.01, 0.05, 0.07, 0.065, 0.06, 0.05, 0.043, 0.037, 0.033, 0.029},
421  //
422  // multiplicity 5 (24 channels)
423  //
424  // p pi+ pi+ pi- pi0
425  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
426  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03, 0.31, 2.14,
427  3.45, 4.0, 3.5, 2.75, 2.15, 1.73, 1.442, 1.175, 0.963, 0.8},
428 
429  // p pi+ pi0 pi0 pi0
430  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
431  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.004, 0.03, 0.2,
432  0.391, 0.45, 0.4, 0.31, 0.234, 0.19, 0.154, 0.118, 0.098, 0.08},
433 
434  // n pi+ pi+ pi+ pi-
435  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
436  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.002, 0.015, 0.06, 0.2,
437  0.4, 0.715, 0.99, 0.93, 0.75, 0.58, 0.453, 0.342, 0.262, 0.205},
438 
439  // n pi+ pi+ pi0 pi0
440  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
441  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.002, 0.015, 0.06, 0.2,
442  0.4, 0.715, 0.99, 0.93, 0.75, 0.58, 0.453, 0.342, 0.262, 0.205},
443 
444  // S+ K+ pi+ pi- pi0
445  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.008, 0.05, 0.09, 0.09, 0.075, 0.062, 0.053, 0.046, 0.04},
448 
449  // S+ K+ pi0 pi0 pi0
450  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.006, 0.025, 0.04, 0.04, 0.035, 0.031, 0.027, 0.023, 0.019},
453 
454  // S+ K0 pi+ pi+ pi-
455  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.005, 0.025, 0.04, 0.04, 0.035, 0.031, 0.027, 0.023, 0.019},
458 
459  // S+ K0 pi+ pi0 pi0
460  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.005, 0.025, 0.04, 0.04, 0.035, 0.031, 0.027, 0.023, 0.019},
463 
464  // L K0 pi+ pi+ pi0
465  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
466  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
467  0, 0.006, 0.035, 0.06, 0.06, 0.055, 0.049, 0.045, 0.039, 0.034},
468 
469  // L K+ pi+ pi+ pi-
470  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
471  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
472  0, 0.004, 0.02, 0.045, 0.05, 0.045, 0.041, 0.037, 0.033, 0.029},
473 
474  // L K+ pi+ pi0 pi0
475  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.004, 0.02, 0.04, 0.05, 0.045, 0.041, 0.037, 0.033, 0.029},
478 
479  // S0 K+ pi+ pi+ pi-
480  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.002, 0.01, 0.02, 0.019, 0.017, 0.016, 0.013, 0.011},
483 
484  // S0 K+ pi+ pi0 pi0
485  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.002, 0.01, 0.02, 0.019, 0.017, 0.016, 0.013, 0.011},
488 
489  // S0 K0 pi+ pi+ pi0
490  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.002, 0.01, 0.02, 0.019, 0.017, 0.016, 0.013, 0.011},
493 
494  // S- K+ pi+ pi+ pi0
495  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.012, 0.03, 0.027, 0.023, 0.02, 0.016, 0.013, 0.011},
498 
499  // S- K0 pi+ pi+ pi+
500  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.012, 0.03, 0.027, 0.023, 0.02, 0.016, 0.013, 0.011},
503 
504  // p pi+ pi- K+ K0bar
505  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.015, 0.06, 0.08, 0.075, 0.067, 0.061, 0.055, 0.051},
508 
509  // p pi+ pi+ K0 K-
510  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.015, 0.04, 0.05, 0.05, 0.046, 0.043, 0.038, 0.034},
513 
514  // p pi+ pi0 K+ K-
515  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.02, 0.05, 0.06, 0.055, 0.05, 0.046, 0.041, 0.036, 0.034},
518 
519  // p pi+ pi0 K0 K0bar
520  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.039, 0.05, 0.055, 0.05, 0.046, 0.041, 0.036, 0.034},
523 
524  // p pi0 pi0 K+ K0bar
525  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.01, 0.03, 0.04, 0.04, 0.038, 0.032, 0.027, 0.023},
528 
529  // n pi+ pi+ K+ K-
530  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
531  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
532  0.0, 0.02, 0.045, 0.055, 0.05, 0.045, 0.041, 0.037, 0.033, 0.029},
533 
534  // n pi+ pi+ K0 K0bar
535  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
536  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
537  0.0, 0.01, 0.035, 0.035, 0.03, 0.026, 0.024, 0.02, 0.017, 0.015},
538 
539  // n pi+ pi0 K+ K0bar
540  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
541  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
542  0.0, 0.01, 0.035, 0.035, 0.03, 0.026, 0.024, 0.02, 0.017, 0.015},
543  //
544  // multiplicity 6 (33 channels)
545  //
546  // p pi+ pi+ pi+ pi- pi-
547  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
548  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
549  0.091, 0.275, 0.397, 0.48, 0.5, 0.44, 0.411, 0.363, 0.317, 0.274},
550 
551  // p pi+ pi+ pi- pi0 pi0
552  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
553  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
554  0.091, 0.275, 0.397, 0.48, 0.5, 0.44, 0.411, 0.363, 0.317, 0.274},
555 
556  // p pi+ pi0 pi0 pi0 pi0
557  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
558  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
559  0.01, 0.04, 0.064, 0.085, 0.09, 0.085, 0.082, 0.075, 0.066, 0.057},
560 
561  // n pi+ pi+ pi0 pi0 pi0
562  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
563  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
564  0.091, 0.275, 0.397, 0.48, 0.5, 0.44, 0.411, 0.363, 0.317, 0.274},
565 
566  // n pi+ pi+ pi+ pi- pi0
567  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
568  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02,
569  0.091, 0.275, 0.397, 0.48, 0.5, 0.44, 0.411, 0.363, 0.317, 0.274},
570 
571  // L K+ pi+ pi+ pi- pi0
572  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
573  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
574  0.001, 0.008, 0.03, 0.06, 0.08, 0.1, 0.103, 0.096, 0.08, 0.068},
575 
576  // L K+ pi+ pi0 pi0 pi0
577  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
578  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
579  0.0, 0.003, 0.01, 0.02, 0.027, 0.033, 0.034, 0.032, 0.026, 0.023},
580 
581  // L K0 pi+ pi+ pi0 pi0
582  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
583  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
584  0.0, 0.004, 0.015, 0.03, 0.04, 0.05, 0.051, 0.048, 0.039, 0.034},
585 
586  // L K0 pi+ pi+ pi+ pi-
587  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
588  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
589  0.0, 0.003, 0.01, 0.02, 0.027, 0.033, 0.034, 0.032, 0.026, 0.023},
590 
591  // S+ K+ pi+ pi+ pi- pi-
592  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
593  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
594  0.001, 0.006, 0.021, 0.042, 0.06, 0.07, 0.072, 0.067, 0.056, 0.048},
595 
596  // S+ K+ pi+ pi- pi0 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.002, 0.012, 0.042, 0.084, 0.12, 0.14, 0.144, 0.135, 0.112, 0.096},
600 
601  // S+ K+ pi0 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.0, 0.001, 0.003, 0.005, 0.006, 0.006, 0.006, 0.005, 0.004, 0.003},
605 
606  // S+ K0 pi+ pi+ pi- 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.003, 0.01, 0.02, 0.028, 0.033, 0.034, 0.032, 0.026, 0.023},
610 
611  // S+ K0 pi+ 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.001, 0.003, 0.007, 0.01, 0.011, 0.011, 0.011, 0.01, 0.008},
615 
616  // S0 K+ pi+ pi+ pi- 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.001, 0.003, 0.006, 0.009, 0.011, 0.011, 0.011, 0.01, 0.009},
620 
621  // S0 K+ pi+ 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.001, 0.002, 0.003, 0.004, 0.004, 0.004, 0.004, 0.003},
625 
626  // S0 K0 pi+ pi+ pi+ pi-
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.001, 0.002, 0.003, 0.005, 0.006, 0.006, 0.006, 0.005, 0.005},
630 
631  // S0 K0 pi+ pi+ 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.001, 0.002, 0.004, 0.007, 0.008, 0.008, 0.007, 0.007, 0.006},
635 
636  // S- K+ pi+ pi+ pi0 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.001, 0.004, 0.009, 0.014, 0.015, 0.015, 0.014, 0.012, 0.01},
640 
641  // S- K+ pi+ pi+ pi+ pi-
642  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
643  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
644  0.0, 0.001, 0.003, 0.006, 0.008, 0.01, 0.01, 0.01, 0.009, 0.008},
645 
646  // S- K0 pi+ pi+ pi+ pi0
647  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
648  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
649  0.0, 0.001, 0.003, 0.006, 0.008, 0.01, 0.01, 0.01, 0.008, 0.007},
650 
651  // p K+ K0b pi+ pi- pi0
652  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
653  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
654  0.0, 0.003, 0.01, 0.03, 0.075, 0.14, 0.165, 0.156, 0.146, 0.137},
655 
656  // p K+ K0b pi0 pi0 pi0
657  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
658  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
659  0.0, 0.0, 0.002, 0.006, 0.014, 0.023, 0.028, 0.027, 0.024, 0.023},
660 
661  // p K+ K- pi+ pi0 pi0
662  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
663  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
664  0.0, 0.001, 0.005, 0.015, 0.037, 0.07, 0.082, 0.078, 0.073, 0.068},
665 
666  // p K+ K- pi+ pi+ pi-
667  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
668  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
669  0.0, 0.001, 0.005, 0.015, 0.037, 0.07, 0.082, 0.078, 0.073, 0.068},
670 
671  // p K0 K0b pi+ pi0 pi0
672  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
673  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
674  0.0, 0.001, 0.003, 0.009, 0.02, 0.035, 0.041, 0.038, 0.036, 0.034},
675 
676  // p K0 K0b pi+ pi+ pi-
677  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
678  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
679  0.0, 0.001, 0.003, 0.009, 0.02, 0.035, 0.041, 0.038, 0.036, 0.034},
680 
681  // p K0 K- pi+ pi+ pi0
682  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
683  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
684  0.0, 0.001, 0.005, 0.015, 0.037, 0.07, 0.082, 0.078, 0.073, 0.068},
685 
686  // n K+ K0b pi+ pi0 pi0
687  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
688  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
689  0.0, 0.001, 0.005, 0.02, 0.045, 0.07, 0.082, 0.078, 0.073, 0.068},
690 
691  // n K+ K0b pi+ pi+ pi-
692  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
693  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
694  0.0, 0.001, 0.005, 0.02, 0.045, 0.07, 0.082, 0.078, 0.073, 0.068},
695 
696  // n K+ K- pi+ p1i+ pi0
697  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
698  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
699  0.0, 0.001, 0.005, 0.02, 0.045, 0.07, 0.082, 0.078, 0.073, 0.068},
700 
701  // n K0 K0b pi+ pi+ pi0
702  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
703  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
704  0.0, 0.001, 0.004, 0.015, 0.035, 0.052, 0.062, 0.059, 0.055, 0.051},
705 
706  // n K0 K- pi+ pi+ pi+
707  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
708  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
709  0.0, 0.0, 0.001, 0.004, 0.01, 0.018, 0.021, 0.019, 0.018, 0.017},
710  //
711  // multiplicity 7 (41 channels)
712  //
713  // p pi+ pi+ pi+ pi- pi- pi0
714  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
715  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
716  0.041, 0.22, 0.55, 0.8, 0.9, 0.9, 0.927, 0.961, 0.986, 1.028},
717 
718  // p pi+ pi+ pi- pi0 pi0 pi0
719  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
720  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
721  0.041, 0.22, 0.55, 0.8, 0.9, 0.9, 0.927, 0.961, 0.986, 1.028},
722 
723  // p pi+ pi0 pi0 pi0 pi0 pi0
724  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
725  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
726  0.004, 0.025, 0.06, 0.09, 0.1, 0.095, 0.093, 0.096, 0.098, 0.103},
727 
728  // n pi+ pi+ pi+ pi+ pi- pi-
729  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
730  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
731  0.004, 0.024, 0.075, 0.15, 0.25, 0.3, 0.309, 0.32, 0.328, 0.342},
732 
733  // n pi+ pi+ pi+ pi- pi0 pi0
734  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
735  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
736  0.02, 0.089, 0.24, 0.45, 0.75, 0.9, 0.927, 0.961, 0.986, 1.028},
737 
738  // n pi+ pi+ pi0 pi0 pi0 pi0
739  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
740  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
741  0.004, 0.023, 0.075, 0.15, 0.25, 0.3, 0.309, 0.32, 0.328, 0.342},
742 
743  // L K+ pi+ pi+ pi- pi0 pi0
744  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
745  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
746  0.0, 0.002, 0.006, 0.014, 0.027, 0.039, 0.04, 0.04, 0.038, 0.038},
747 
748  // L K+ pi+ pi+ pi+ pi- pi-
749  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
750  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
751  0.0, 0.001, 0.003, 0.007, 0.011, 0.013, 0.014, 0.014, 0.013, 0.013},
752 
753  // L K+ pi+ pi0 pi0 pi0 pi0
754  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
755  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
756  0.0, 0.0, 0.001, 0.003, 0.005, 0.007, 0.007, 0.007, 0.008, 0.007},
757 
758  // L K0 pi+ pi+ pi+ pi- pi0
759  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
760  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
761  0.0, 0.002, 0.007, 0.02, 0.04, 0.07, 0.082, 0.08, 0.073, 0.068},
762 
763  // L K0 pi+ pi+ pi0 pi0 pi0
764  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
765  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
766  0.0, 0.001, 0.003, 0.01, 0.023, 0.035, 0.041, 0.038, 0.036, 0.034},
767 
768  // S+ K+ pi+ pi+ pi- pi- pi0
769  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
770  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
771  0.0, 0.001, 0.003, 0.007, 0.012, 0.017, 0.017, 0.017, 0.016, 0.016},
772 
773  // S+ K+ pi+ pi- pi0 pi0 pi0
774  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
775  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
776  0.0, 0.001, 0.002, 0.004, 0.008, 0.011, 0.011, 0.012, 0.011, 0.01},
777 
778  // S+ K+ pi0 pi0 pi0 pi0 pi0 negligible
779 
780  // S+ K0 pi+ pi+ pi- pi0 pi0
781  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
782  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
783  0.0, 0.001, 0.002, 0.004, 0.007, 0.01, 0.01, 0.01, 0.009, 0.009},
784 
785  // S+ K0 pi+ pi+ pi+ pi- pi-
786  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
787  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
788  0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.004, 0.004, 0.004, 0.005},
789 
790  // S+ K0 pi+ pi0 pi0 pi0 pi0
791  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
792  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
793  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002, 0.002},
794 
795  // S0 K+ pi+ pi+ pi- pi0 pi0
796  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
797  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
798  0.0, 0.001, 0.003, 0.006, 0.009, 0.013, 0.013, 0.013, 0.013, 0.013},
799 
800  // S0 K+ pi+ pi+ pi+ pi- pi-
801  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
802  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
803  0.0, 0.0, 0.001, 0.002, 0.003, 0.004, 0.004, 0.004, 0.004, 0.005},
804 
805  // S0 K+ pi+ pi0 pi0 pi0 pi0
806  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
807  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
808  0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.003, 0.003, 0.003, 0.003},
809 
810  // S0 K0 pi+ pi+ pi+ pi- pi0
811  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
812  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
813  0.0, 0.001, 0.003, 0.007, 0.014, 0.023, 0.028, 0.027, 0.026, 0.024},
814 
815  // S0 K0 pi+ pi+ pi0 pi0 pi0
816  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
817  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
818  0.0, 0.0, 0.001, 0.003, 0.008, 0.012, 0.013, 0.013, 0.012, 0.011},
819 
820  // S- K+ pi+ pi+ pi+ pi- pi0
821  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
822  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
823  0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.010, 0.012, 0.012, 0.013},
824 
825  // S- K+ pi+ pi+ pi0 pi0 pi0
826  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
827  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
828  0.0, 0.0, 0.0, 0.001, 0.003, 0.005, 0.006, 0.007, 0.008, 0.008},
829 
830  // S- K0 pi+ pi+ pi+ pi+ pi-
831  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
832  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
833  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.005, 0.006, 0.007, 0.007},
834 
835  // S- K0 pi+ pi+ pi+ pi0 pi0
836  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
837  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
838  0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.007, 0.008, 0.009, 0.009},
839 
840  // p K+ K0b pi+ pi- pi0 pi0
841  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
842  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
843  0.0, 0.001, 0.004, 0.015, 0.037, 0.05, 0.05, 0.05, 0.05, 0.051},
844 
845  // p K+ K0b pi+ pi+ pi- pi-
846  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
847  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
848  0.0, 0.0, 0.002, 0.007, 0.018, 0.025, 0.026, 0.026, 0.025, 0.025},
849 
850  // p K+ K0b pi0 pi0 pi0 pi0
851  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
852  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
853  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.005, 0.005, 0.005, 0.006},
854 
855  // p K+ K- pi+ pi+ pi- pi0
856  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
857  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
858  0.0, 0.001, 0.004, 0.017, 0.043, 0.058, 0.059, 0.059, 0.058, 0.058},
859 
860  // p K+ K- pi+ pi0 pi0 pi0
861  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
862  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
863  0.0, 0.0, 0.002, 0.006, 0.014, 0.019, 0.02, 0.02, 0.02, 0.019},
864 
865  // p K0 K0b pi+ pi+ pi- pi0
866  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
867  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
868  0.0, 0.001, 0.006, 0.025, 0.065, 0.088, 0.088, 0.088, 0.087, 0.086},
869 
870  // p K0 K0b pi+ pi0 pi0 pi0
871  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
872  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
873  0.0, 0.0, 0.002, 0.008, 0.022, 0.029, 0.029, 0.03, 0.03, 0.03},
874 
875  // p K0 K- pi+ pi+ pi0 pi0
876  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
877  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
878  0.0, 0.002, 0.01, 0.033, 0.087, 0.117, 0.116, 0.118, 0.117, 0.114},
879 
880  // p K0 K- pi+ pi+ pi+ pi-
881  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
882  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
883  0.0, 0.003, 0.015, 0.05, 0.13, 0.175, 0.175, 0.176, 0.175, 0.172},
884 
885  // n K+ K0b pi+ pi+ pi- pi0
886  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
887  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
888  0.0, 0.001, 0.005, 0.018, 0.043, 0.058, 0.059, 0.059, 0.058, 0.058},
889 
890  // n K+ K0b pi+ pi0 pi0 pi0
891  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
892  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
893  0.0, 0.0, 0.002, 0.006, 0.014, 0.019, 0.02, 0.02, 0.02, 0.019},
894 
895  // n K+ K- pi+ pi+ pi0 pi0
896  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
897  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
898  0.0, 0.0, 0.002, 0.008, 0.022, 0.029, 0.029, 0.03, 0.03, 0.029},
899 
900  // n K+ K- pi+ pi+ pi+ pi-
901  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
902  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
903  0.0, 0.0, 0.002, 0.006, 0.014, 0.019, 0.02, 0.02, 0.02, 0.019},
904 
905  // n K0 K0b pi+ pi+ pi0 pi0
906  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
907  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
908  0.0, 0.0, 0.002, 0.008, 0.022, 0.029, 0.029, 0.03, 0.03, 0.029},
909 
910  // n K0 K0b pi+ pi+ pi+ pi-
911  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
912  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
913  0.0, 0.0, 0.002, 0.006, 0.014, 0.019, 0.02, 0.02, 0.02, 0.019},
914 
915  // n K0 K- pi+ pi+ pi+ pi0
916  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
917  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
918  0.0, 0.0, 0.002, 0.006, 0.014, 0.019, 0.02, 0.02, 0.02, 0.019},
919  //
920  // multiplicity 8 (47 channels)
921  //
922  // p pi+ pi+ pi+ pi+ pi- pi- pi-
923  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
924  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
925  0.01, 0.028, 0.045, 0.07, 0.085, 0.095, 0.103, 0.112, 0.12, 0.131},
926 
927  // p pi+ pi+ pi+ pi- pi- pi0 pi0
928  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
929  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
930  0.01, 0.05, 0.15, 0.3, 0.45, 0.57, 0.618, 0.673, 0.723, 0.787},
931 
932  // p pi+ pi+ pi- pi0 pi0 pi0 pi0
933  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
934  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
935  0.005, 0.025, 0.075, 0.15, 0.22, 0.28, 0.309, 0.337, 0.361, 0.394},
936 
937  // p pi+ pi0 pi0 pi0 pi0 pi0 pi0
938  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
939  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
940  0.001, 0.005, 0.01, 0.015, 0.018, 0.02, 0.021, 0.022, 0.024, 0.026},
941 
942  // n pi+ pi+ pi+ pi+ pi- pi- pi0
943  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
944  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
945  0.01, 0.06, 0.12, 0.2, 0.27, 0.33, 0.36, 0.395, 0.426, 0.468},
946 
947  // n pi+ pi+ pi+ pi- pi0 pi0 pi0
948  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
949  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
950  0.013, 0.07, 0.16, 0.254, 0.36, 0.44, 0.48, 0.527, 0.57, 0.624},
951 
952  // n pi+ pi+ pi0 pi0 pi0 pi0 pi0
953  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
954  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
955  0.002, 0.014, 0.024, 0.04, 0.052, 0.066, 0.072, 0.079, 0.085, 0.094},
956 
957  // L K+ pi+ pi+ pi+ pi- pi- pi0
958  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
959  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
960  0.0, 0.0, 0.002, 0.006, 0.014, 0.026, 0.041, 0.053, 0.059, 0.068},
961 
962  // L K+ pi+ pi0 pi0 pi0 pi0 pi0
963  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
964  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
965  0.0, 0.0, 0.0, 0.0, 0.001, 0.003, 0.004, 0.005, 0.006, 0.007},
966 
967  // L K0 pi+ pi+ pi+ pi- pi0 pi0
968  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
969  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
970  0.0, 0.0, 0.001, 0.003, 0.008, 0.017, 0.027, 0.036, 0.04, 0.046},
971 
972  // L K0 pi+ pi+ pi0 pi0 pi0 pi0
973  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
974  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
975  0.0, 0.0, 0.001, 0.003, 0.008, 0.015, 0.021, 0.027, 0.031, 0.034},
976 
977  // L K0 pi+ pi+ pi+ pi+ pi- pi-
978  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
979  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
980  0.0, 0.0, 0.001, 0.003, 0.008, 0.015, 0.021, 0.027, 0.031, 0.034},
981 
982  // S+ K+ pi+ pi+ pi- pi- pi0 pi0
983  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
984  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
985  0.0, 0.0, 0.001, 0.003, 0.005, 0.009, 0.012, 0.016, 0.017, 0.019},
986 
987  // S+ K+ pi+ pi- pi0 pi0 pi0 pi0
988  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
989  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
990  0.0, 0.0, 0.0, 0.0, 0.001, 0.003, 0.004, 0.005, 0.006, 0.007},
991 
992  // S+ K+ pi+ pi+ pi+ pi- pi- pi-
993  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
994  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
995  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.007, 0.008, 0.009},
996 
997  // S+ K+ pi0 pi0 pi0 pi0 pi0 pi0 (negligible)
998 
999  // S+ K0 pi+ pi+ pi+ pi- pi- pi0
1000  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1001  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1002  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.01, 0.011, 0.012},
1003 
1004  // S+ K0 pi+ pi+ pi- pi0 pi0 pi0
1005  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1006  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1007  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.01, 0.011, 0.012},
1008 
1009  // S+ K0 pi+ pi0 pi0 pi0 pi0 pi0 (negligible)
1010 
1011  // S0 K+ pi+ pi+ pi+ pi- pi- pi0
1012  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1013  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1014  0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.013, 0.018, 0.021, 0.023},
1015 
1016  // S0 K+ pi+ pi0 pi0 pi0 pi0 pi0
1017  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1018  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1019  0.0, 0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002},
1020 
1021  // S0 K0 pi+ pi+ pi+ pi- pi0 pi0
1022  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1023  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1024  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.012, 0.013, 0.015},
1025 
1026  // S0 K0 pi+ pi+ pi0 pi0 pi0 pi0
1027  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1028  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1029  0.0, 0.0, 0.0, 0.001, 0.003, 0.005, 0.007, 0.009, 0.01, 0.011},
1030 
1031  // S0 K0 pi+ pi+ pi+ pi+ pi- pi-
1032  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1033  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1034  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.009, 0.01, 0.011},
1035 
1036  // S- K+ pi+ pi+ pi+ pi- pi0 pi0
1037  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1038  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1039  0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.012, 0.016, 0.017, 0.019},
1040 
1041  // S- K+ pi+ pi+ pi0 pi0 pi0 pi0
1042  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1043  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1044  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.004, 0.004, 0.005},
1045 
1046  // S- K+ pi+ pi+ pi+ pi+ pi- pi-
1047  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1048  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1049  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.004, 0.004, 0.005},
1050 
1051  // S- K0 pi+ pi+ pi+ pi+ pi- pi0
1052  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1053  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1054  0.0, 0.0, 0.0, 0.001, 0.003, 0.005, 0.006, 0.007, 0.008, 0.009},
1055 
1056  // S- K0 pi+ pi+ pi+ pi0 pi0 pi0
1057  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1058  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1059  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.005, 0.006, 0.007, 0.008},
1060 
1061  // p K+ K0b pi+ pi+ pi- pi- pi0
1062  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1063  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1064  0.0, 0.0, 0.002, 0.007, 0.019, 0.036, 0.051, 0.064, 0.077, 0.091},
1065 
1066  // p K+ K0b pi+ pi- pi0 pi0 pi0
1067  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1068  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1069  0.0, 0.0, 0.001, 0.004, 0.011, 0.022, 0.034, 0.043, 0.051, 0.06},
1070 
1071  // p K+ K0b pi0 pi0 pi0 pi0 pi0
1072  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1073  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1074  0.0, 0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002},
1075 
1076  // p K+ K- pi+ pi+ pi- pi0 pi0
1077  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1078  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1079  0.0, 0.0, 0.002, 0.007, 0.019, 0.036, 0.051, 0.064, 0.077, 0.091},
1080 
1081  // p K+ K- pi+ pi+ pi+ pi- pi-
1082  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1083  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1084  0.0, 0.0, 0.001, 0.003, 0.007, 0.012, 0.017, 0.024, 0.027, 0.031},
1085 
1086  // p K+ K- pi+ pi0 pi0 pi0 pi0
1087  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1088  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1089  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.011, 0.013, 0.015},
1090 
1091  // p K0 K0b pi+ pi+ pi- pi0 pi0
1092  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1093  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1094  0.0, 0.0, 0.002, 0.007, 0.019, 0.036, 0.051, 0.064, 0.077, 0.091},
1095 
1096  // p K0 K0b pi+ pi+ pi+ pi- pi-
1097  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1098  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1099  0.0, 0.0, 0.001, 0.003, 0.007, 0.012, 0.017, 0.024, 0.027, 0.031},
1100 
1101  // p K0 K0b pi+ pi0 pi0 pi0 pi0
1102  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1103  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1104  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.011, 0.013, 0.015},
1105 
1106  // p K0 K- pi+ pi+ pi+ pi- pi0
1107  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1108  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1109  0.0, 0.0, 0.001, 0.004, 0.01, 0.023, 0.034, 0.043, 0.051, 0.06},
1110 
1111  // p K0 K- pi+ pi+ pi0 pi0 pi0
1112  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1113  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1114  0.0, 0.0, 0.001, 0.003, 0.007, 0.012, 0.017, 0.024, 0.027, 0.031},
1115 
1116  // n K+ K0b pi+ pi+ pi+ pi- pi-
1117  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1118  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1119  0.0, 0.0, 0.0, 0.001, 0.006, 0.016, 0.022, 0.03, 0.035, 0.037},
1120 
1121  // n K+ K0b pi+ pi+ pi- pi0 pi0
1122  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1123  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1124  0.0, 0.0, 0.0, 0.004, 0.02, 0.044, 0.065, 0.09, 0.1, 0.11},
1125 
1126  // n K+ K0b pi+ pi0 pi0 pi0 pi0
1127  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1128  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1129  0.0, 0.0, 0.0, 0.0, 0.003, 0.008, 0.01, 0.015, 0.017, 0.018},
1130 
1131  // n K+ K- pi+ pi+ pi+ pi- pi0
1132  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1133  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1134  0.0, 0.0, 0.001, 0.003, 0.007, 0.012, 0.017, 0.022, 0.026, 0.031},
1135 
1136  // n K+ K- pi+ pi+ pi0 pi0 pi0
1137  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1138  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1139  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.011, 0.013, 0.015},
1140 
1141  // n K0 K0b pi+ pi+ pi+ pi- pi0
1142  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1143  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1144  0.0, 0.0, 0.0, 0.002, 0.007, 0.016, 0.021, 0.026, 0.031, 0.037},
1145 
1146  // n K0 K0b pi+ pi+ pi0 pi0 pi0
1147  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1148  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1149  0.0, 0.0, 0.0, 0.0, 0.003, 0.008, 0.01, 0.013, 0.015, 0.018},
1150 
1151  // n K0 K- pi+ pi+ pi+ pi0 pi0
1152  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1153  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1154  0.0, 0.0, 0.0, 0.0, 0.003, 0.008, 0.01, 0.013, 0.015, 0.018},
1155 
1156  // n K0 K- pi+ pi+ pi+ pi+ pi-
1157  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1158  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1159  0.0, 0.0, 0.0, 0.0, 0.001, 0.004, 0.005, 0.006, 0.008, 0.009},
1160  //
1161  // multiplicity 9 (55 channels)
1162  //
1163  // p pi+ pi+ pi+ pi+ pi- pi- pi- pi0
1164  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1165  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1166  0.003, 0.013, 0.037, 0.08, 0.14, 0.2, 0.257, 0.31, 0.35, 0.399},
1167 
1168  // p pi+ pi+ pi+ pi- pi- pi0 pi0 pi0
1169  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1170  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1171  0.006, 0.032, 0.085, 0.16, 0.27, 0.4, 0.515, 0.62, 0.701, 0.799},
1172 
1173  // p pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0
1174  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1175  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1176  0.002, 0.011, 0.027, 0.05, 0.085, 0.12, 0.154, 0.186, 0.21, 0.24},
1177 
1178  // p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
1179  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1180  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1181  0.0, 0.0, 0.001, 0.003, 0.005, 0.007, 0.008, 0.009, 0.01, 0.011},
1182 
1183  // n pi+ pi+ pi+ pi+ pi+ pi- pi- pi-
1184  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1185  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1186  0.0, 0.003, 0.01, 0.02, 0.036, 0.048, 0.066, 0.082, 0.093, 0.103},
1187 
1188  // n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
1189  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1190  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1191  0.003, 0.022, 0.075, 0.15, 0.24, 0.36, 0.495, 0.617, 0.698, 0.77},
1192 
1193  // n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
1194  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1195  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1196  0.002, 0.015, 0.05, 0.1, 0.17, 0.25, 0.329, 0.412, 0.465, 0.513},
1197 
1198  // n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
1199  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1200  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1201  0.0, 0.002, 0.005, 0.01, 0.017, 0.025, 0.033, 0.041, 0.046, 0.051},
1202 
1203  // L K+ pi+ pi+ pi+ pi- pi- pi0 pi0
1204  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1205  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1206  0.0, 0.0, 0.002, 0.006, 0.013, 0.026, 0.041, 0.053, 0.059, 0.068},
1207 
1208  // L K+ pi+ pi+ pi- pi0 pi0 pi0 pi0
1209  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1210  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1211  0.0, 0.0, 0.001, 0.003, 0.007, 0.014, 0.021, 0.027, 0.03, 0.034},
1212 
1213  // L K+ pi+ pi+ pi+ pi+ pi- pi- pi-
1214  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1215  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1216  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.008, 0.01, 0.011, 0.012},
1217 
1218  // L K+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
1219  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1220  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1221  0.0, 0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002},
1222 
1223  // L K0 pi+ pi+ pi+ pi- pi0 pi0 pi0
1224  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1225  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1226  0.0, 0.0, 0.001, 0.004, 0.009, 0.018, 0.028, 0.035, 0.039, 0.046},
1227 
1228  // L K0 pi+ pi+ pi+ pi+ pi- pi- pi0
1229  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1230  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1231  0.0, 0.0, 0.001, 0.003, 0.007, 0.014, 0.021, 0.027, 0.03, 0.034},
1232 
1233  // L K0 pi+ pi+ pi0 pi0 pi0 pi0 pi0
1234  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1235  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1236  0.0, 0.0, 0.0, 0.0, 0.001, 0.003, 0.004, 0.005, 0.006, 0.007},
1237 
1238  // S+ K+ pi+ pi+ pi- pi- pi0 pi0 pi0
1239  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1240  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1241  0.0, 0.0, 0.0, 0.002, 0.005, 0.009, 0.013, 0.018, 0.021, 0.023},
1242 
1243  // S+ K+ pi+ pi+ pi+ pi- pi- pi- pi0
1244  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1245  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1246  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.013, 0.014, 0.015},
1247 
1248  // S+ K+ pi+ pi- pi0 pi0 pi0 pi0 pi0
1249  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1250  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1251  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.004, 0.004, 0.005},
1252 
1253  // S+ K+ pi0 pi0 pi0 pi0 pi0 pi0 pi0 (negligible)
1254 
1255  // S+ K0 pi+ pi+ pi+ pi- pi- pi0 pi0
1256  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1257  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1258  0.0, 0.0, 0.0, 0.002, 0.005, 0.009, 0.013, 0.018, 0.021, 0.023},
1259 
1260  // S+ K0 pi+ pi+ pi- pi0 pi0 pi0 pi0
1261  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1262  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1263  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.009, 0.01, 0.011},
1264 
1265  // S+ K0 pi+ pi+ pi+ pi+ pi- pi- pi-
1266  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1267  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1268  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.003, 0.003, 0.003},
1269 
1270  // S+ K0 pi+ pi0 pi0 pi0 pi0 pi0 pi0 (negligible)
1271 
1272  // S0 K+ pi+ pi+ pi+ pi- pi- pi0 pi0
1273  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1274  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1275  0.0, 0.0, 0.0, 0.002, 0.005, 0.009, 0.013, 0.018, 0.021, 0.023},
1276 
1277  // S0 K+ pi+ pi+ pi- pi0 pi0 pi0 pi0
1278  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1279  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1280  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.009, 0.01, 0.011},
1281 
1282  // S0 K+ pi+ pi+ pi+ pi+ pi- pi- pi-
1283  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1284  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1285  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.003, 0.003, 0.003, 0.003},
1286 
1287  // S0 K+ pi+ pi0 pi0 pi0 pi0 pi0 pi0 (negligible)
1288 
1289  // S0 K0 pi+ pi+ pi+ pi- pi0 pi0 pi0
1290  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1291  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1292  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.012, 0.013, 0.015},
1293 
1294  // S0 K0 pi+ pi+ pi+ pi+ pi- pi- pi0
1295  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1296  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1297  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.009, 0.01, 0.011},
1298 
1299  // S0 K0 pi+ pi+ pi0 pi0 pi0 pi0 pi0
1300  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1301  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1302  0.0, 0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002},
1303 
1304  // S- K+ pi+ pi+ pi+ pi- pi0 pi0 pi0
1305  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1306  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1307  0.0, 0.0, 0.0, 0.001, 0.003, 0.006, 0.009, 0.013, 0.014, 0.015},
1308 
1309  // S- K+ pi+ pi+ pi+ pi+ pi- pi- pi0
1310  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1311  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1312  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.01, 0.012, 0.013},
1313 
1314  // S- K+ pi+ pi+ pi0 pi0 pi0 pi0 pi0
1315  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1316  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1317  0.0, 0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.002, 0.002},
1318 
1319  // S- K0 pi+ pi+ pi+ pi+ pi- pi0 pi0
1320  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1321  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1322  0.0, 0.0, 0.0, 0.001, 0.002, 0.004, 0.006, 0.01, 0.012, 0.013},
1323 
1324  // S- K0 pi+ pi+ pi+ pi0 pi0 pi0 pi0
1325  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1326  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1327  0.0, 0.0, 0.0, 0.0, 0.001, 0.002, 0.002, 0.003, 0.003, 0.003},
1328 
1329  // S- K0 pi+ pi+ pi+ pi+ pi+ pi- pi- (negligible)
1330 
1331  // p K+ K0b pi+ pi+ pi- pi- pi0 pi0
1332  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1333  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1334  0.0, 0.0, 0.0, 0.005, 0.027, 0.056, 0.078, 0.098, 0.114, 0.132},
1335 
1336  // p K+ K0b pi+ pi- pi0 pi0 pi0 pi0
1337  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1338  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1339  0.0, 0.0, 0.0, 0.003, 0.009, 0.018, 0.026, 0.033, 0.038, 0.044},
1340 
1341  // p K+ K0b pi+ pi+ pi+ pi- pi- pi-
1342  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1343  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1344  0.0, 0.0, 0.0, 0.002, 0.006, 0.012, 0.018, 0.022, 0.026, 0.03},
1345 
1346  // p K+ K0b pi0 pi0 pi0 pi0 pi0 pi0 (negligible)
1347 
1348  // p K+ K- pi+ pi+ pi- pi0 pi0 pi0
1349  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1350  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1351  0.0, 0.0, 0.0, 0.004, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1352 
1353  // p K+ K- pi+ pi+ pi+ pi- pi- pi0
1354  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1355  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1356  0.0, 0.0, 0.0, 0.004, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1357 
1358  // p K+ K- pi+ pi0 pi0 pi0 pi0 pi0
1359  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1360  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1361  0.0, 0.0, 0.0, 0.0, 0.002, 0.004, 0.005, 0.006, 0.008, 0.009},
1362 
1363  // p K0 K0b pi+ pi+ pi- pi0 pi0 pi0
1364  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1365  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1366  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1367 
1368  // p K0 K0b pi+ pi+ pi+ pi- pi- pi0
1369  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1370  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1371  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1372 
1373  // p K0 K0b pi+ pi0 pi0 pi0 pi0 pi0
1374  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1375  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1376  0.0, 0.0, 0.0, 0.0, 0.002, 0.004, 0.005, 0.006, 0.008, 0.009},
1377 
1378  // p K0 K- pi+ pi+ pi+ pi- pi0 pi0
1379  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1380  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1381  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1382 
1383  // p K0 K- pi+ pi+ pi0 pi0 pi0 pi0
1384  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1385  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1386  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1387 
1388  // p K0 K- pi+ pi+ pi+ pi+ pi- pi-
1389  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1390  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1391  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1392 
1393  // n K+ K0b pi+ pi+ pi- pi0 pi0 pi0
1394  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1395  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1396  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1397 
1398  // n K+ K0b pi+ pi+ pi+ pi- pi- pi0
1399  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1400  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1401  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1402 
1403  // n K+ K0b pi+ pi0 pi0 pi0 pi0 pi0
1404  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1405  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1406  0.0, 0.0, 0.0, 0.0, 0.002, 0.004, 0.005, 0.006, 0.008, 0.009},
1407 
1408  // n K+ K- pi+ pi+ pi+ pi- pi0 pi0
1409  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1410  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1411  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1412 
1413  // n K+ K- pi+ pi+ pi0 pi0 pi0 pi0
1414  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1415  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1416  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1417 
1418  // n K+ K- pi+ pi+ pi+ pi+ pi- pi-
1419  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1420  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1421  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1422 
1423  // n K0 K0b pi+ pi+ pi+ pi- pi0 pi0
1424  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1425  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1426  0.0, 0.0, 0.0, 0.003, 0.015, 0.035, 0.052, 0.065, 0.075, 0.088},
1427 
1428  // n K0 K0b pi+ pi+ pi0 pi0 pi0 pi0
1429  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1430  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1431  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1432 
1433  // n K0 K0b pi+ pi+ pi+ pi+ pi- pi-
1434  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1435  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1436  0.0, 0.0, 0.0, 0.001, 0.004, 0.009, 0.013, 0.016, 0.019, 0.022},
1437 
1438  // n K0 K- pi+ pi+ pi+ pi+ pi- pi0
1439  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1440  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1441  0.0, 0.0, 0.0, 0.002, 0.008, 0.018, 0.026, 0.033, 0.038, 0.044},
1442 
1443  // n K0 K- pi+ pi+ pi+ pi0 pi0 pi0
1444  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1445  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
1446  0.0, 0.0, 0.0, 0.001, 0.006, 0.012, 0.019, 0.022, 0.026, 0.03}};
1447 }
1448 
1449 // Initialize both |T Tz> = |3/2 3/2> channels, using pipP cross-section table
1450 
1451 const G4CascadePiPlusPChannelData::data_t
1452 G4CascadePiPlusPChannelData::data(pipP2bfs, pipP3bfs, pipP4bfs,
1453  pipP5bfs, pipP6bfs, pipP7bfs,
1454  pipP8bfs, pipP9bfs, pipPCrossSections,
1455  pipPtotXSec, pip*pro, "PiPlusP");
1456 
1457 const G4CascadePiMinusNChannelData::data_t
1458 G4CascadePiMinusNChannelData::data(pimN2bfs, pimN3bfs, pimN4bfs,
1459  pimN5bfs, pimN6bfs, pimN7bfs,
1460  pimN8bfs, pimN9bfs, pipPCrossSections,
1461  pipPtotXSec, pim*neu, "PiMinusN");