ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4CascadeXiZeroPChannel.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4CascadeXiZeroPChannel.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 // 20100804 M. Kelsey -- Add name string to ctor
28 // 20110719 M. Kelsey -- Add initial state code to ctor
29 // 20110725 M. Kelsey -- Instantiate cross-section object for self-registration
30 // 20110916 M. Kelsey -- Drop self-registration due to platform inconsistencies
31 
32 #include "G4CascadeXiZeroPChannel.hh"
33 #include "G4InuclParticleNames.hh"
34 using namespace G4InuclParticleNames;
35 
36 namespace {
37  // Outgoing particle types of a given multiplicity
38 
39  static const G4int x0p2bfs[3][2] =
40  {{1, 29}, {21,23}, {23,25}};
41 
42  static const G4int x0p3bfs[18][3] =
43  {{1,17,21}, {3,21,21}, {1,13,23}, {2,17,23}, {7,21,23}, {5,23,23},
44  {7,23,25}, {3,23,27}, {1,17,25}, {3,21,25}, {3,25,25}, {1,7,29},
45  {2,3,29}, {11,21,29},{11,25,29},{15,23,29},{1,3,31}, {11,23,31}};
46 
47  static const G4int x0p4bfs[53][4] =
48  {{1,1,13,17}, {1,2,17,17}, {1,7,17,21}, {1,3,13,21}, {1,7,17,25},
49  {1,3,13,25}, {1,7,13,23}, {1,5,17,23}, {1,3,17,27}, {1,7,7,29},
50  {1,3,5,29}, {1,11,13,29}, {1,15,17,29}, {1,3,7,31}, {1,11,17,31},
51  {2,3,17,21}, {2,3,17,25}, {2,7,17,23}, {2,3,13,23}, {2,3,7,29},
52  {2,11,17,29}, {2,3,3,31}, {3,7,21,21}, {11,17,21,21},{3,7,21,25},
53  {11,17,21,25},{7,7,21,23}, {3,5,21,23}, {11,13,21,23},{15,17,21,23},
54  {3,3,21,27}, {7,11,21,29}, {3,15,21,29}, {3,11,21,31}, {7,7,23,25},
55  {3,5,23,25}, {11,13,23,25},{15,17,23,25},{3,7,25,25}, {11,17,25,25},
56  {3,3,25,27}, {3,15,25,29}, {7,11,25,29}, {3,11,25,31}, {3,7,23,27},
57  {11,17,23,27},{5,7,23,23}, {13,15,23,23},{7,15,23,29}, {5,11,23,29},
58  {3,15,23,31}, {7,11,23,31}, {3,11,27,29}};
59 
60  static const G4int x0p5bfs[2][5] =
61  {{1,7,7,7,29}, {1,3,5,7,29}};
62 
63  static const G4int x0p6bfs[2][6] =
64  {{1,7,7,7,7,29}, {1,3,5,7,7,29}};
65 
66  static const G4int x0p7bfs[2][7] =
67  {{1,7,7,7,7,7,29}, {1,3,5,7,7,7,29}};
68 
69  // Cross sections for X0 p -> 2-7 body final states
70  //
71  // first index: 0-2: channels for mult = 2
72  // 3-20: channels for mult = 3
73  // 21-73: channels for mult = 4
74  // 74-75: channels for mult = 5
75  // 76-77: channels for mult = 6
76  // 78-79: channels for mult = 7
77  //
78  // second index: kinetic energy
79  //
80  static const G4double x0pCrossSections[80][31] = {
81  //
82  // multiplicity 2 (3 channels)
83  //
84  // X0 p
85  {16.00, 9.60, 7.00, 6.80, 6.60, 6.40, 6.20, 6.00, 5.82, 5.64,
86  5.45, 4.80, 4.58, 4.20, 3.80, 3.50, 3.20, 3.10, 2.95, 2.80,
87  2.65, 2.55, 2.45, 2.35, 2.25, 2.15, 2.10, 2.05, 2.04, 2.03, 1.95},
88 
89  // L S+
90  { 0.0, 0.0, 0.0, 1.80, 3.60, 5.40, 7.20, 9.00, 9.82,10.64,
91  11.45,14.40,12.97,11.31, 8.38, 5.45, 2.50, 2.10, 1.90, 1.70,
92  1.50, 1.40, 1.30, 1.20, 1.10, 1.05, 1.00, 0.95, 0.94, 0.93, 0.85},
93 
94  // S0 S+
95  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.55, 1.09,
96  1.64, 3.60, 4.32, 3.77, 2.80, 2.00, 1.50, 1.40, 1.30, 1.20,
97  1.10, 1.00, 0.95, 0.90, 0.85, 0.80, 0.75, 0.70, 0.69, 0.68, 0.60},
98 
99  //
100  // multiplicity 3 (18 channels)
101  //
102  // L p K0bar
103  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09, 0.18,
104  0.27, 0.60, 0.71, 0.90, 1.10, 1.00, 0.90, 0.80, 0.68, 0.56,
105  0.40, 0.30, 0.20, 0.10, 0.07, 0.05, 0.03, 0.01, 0.01, 0.01, 0.0},
106 
107  // L L pi+
108  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.46, 0.91,
109  1.36, 3.00, 1.93, 1.80, 1.70, 1.50, 1.30, 1.10, 1.00, 0.90,
110  0.75, 0.60, 0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.14, 0.13, 0.10},
111 
112  // S+ p K-
113  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
114  0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
115  0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
116 
117  // S+ n K0bar
118  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
119  0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
120  0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
121 
122  // S+ L pi0
123  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.55, 1.09,
124  1.64, 3.60, 3.36, 3.10, 3.00, 2.80, 2.60, 2.30, 2.00, 1.80,
125  1.60, 1.50, 1.40, 1.30, 1.20, 1.15, 1.10, 1.05, 1.04, 1.03, 0.90},
126 
127  // S+ S+ pi-
128  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
129  0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
130  0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
131 
132  // S+ S0 pi0
133  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
134  0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
135  0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
136 
137  // S+ S- pi+
138  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
139  0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
140  0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
141 
142  // S0 p K0bar
143  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
144  0.0, 0.0, 0.01, 0.05, 0.10, 0.15, 0.17, 0.17, 0.15, 0.13,
145  0.10, 0.08, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01, 0.0},
146 
147  // S0 L pi+
148  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
149  0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
150  0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
151 
152  // S0 S0 pi+
153  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
154  0.0, 0.0, 0.10, 0.15, 0.20, 0.22, 0.23, 0.20, 0.17, 0.15,
155  0.12, 0.10, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
156 
157  // X0 p pi0
158  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.14, 0.27,
159  0.41, 0.90, 1.02, 1.50, 1.40, 1.10, 1.00, 0.90, 0.80, 0.65,
160  0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.12, 0.10, 0.09, 0.09, 0.05},
161 
162  // X0 n pi+
163  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.14, 0.27,
164  0.41, 0.90, 1.02, 1.50, 1.40, 1.10, 1.00, 0.90, 0.80, 0.65,
165  0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.12, 0.10, 0.09, 0.09, 0.05},
166 
167  // X0 L K+
168  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
169  0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
170  0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
171 
172  // X0 S0 K+
173  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
174  0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
175  0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
176 
177  // X0 S+ K0
178  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
179  0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
180  0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
181 
182  // X- p pi+
183  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.46, 0.91,
184  1.36, 3.00, 3.56, 3.32, 3.10, 3.10, 3.00, 2.80, 2.60, 2.40,
185  2.20, 2.00, 1.90, 1.80, 1.70, 1.65, 1.60, 1.55, 1.54, 1.53, 1.40},
186 
187  // X- S+ K+
188  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
189  0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.02, 0.02, 0.03,
190  0.03, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01, 0.01, 0.01, 0.01},
191 
192  //
193  // multiplicity 4 (53 channels)
194  //
195  // p p K0bar K-
196  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
197  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
198  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
199 
200  // p n K0bar K0bar
201  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
202  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
203  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
204 
205  // p L K0bar pi0
206  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
207  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
208  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
209 
210  // p L K- pi+
211  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
212  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
213  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
214 
215  // p S0 K0bar pi0
216  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
217  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
218  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
219 
220  // p S0 K- pi+
221  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
222  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
223  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
224 
225  // p S+ K- pi0
226  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
227  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
228  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
229 
230  // p S+ K0bar pi-
231  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
232  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
233  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
234 
235  // p S- K0bar pi+
236  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
237  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.13, 0.12, 0.11,
238  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
239 
240  // p X0 2pi0
241  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
242  0.0, 0.0, 0.0, 0.04, 0.09, 0.14, 0.20, 0.25, 0.20, 0.17,
243  0.15, 0.12, 0.10, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.02, 0.01},
244 
245  // p X0 pi+ pi-
246  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
247  0.0, 0.0, 0.0, 0.06, 0.23, 0.40, 0.57, 0.75, 0.70, 0.60,
248  0.50, 0.40, 0.30, 0.25, 0.20, 0.15, 0.10, 0.08, 0.07, 0.07, 0.04},
249 
250  // p X0 K+ K-
251  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
252  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
253  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
254 
255  // p X0 K0bar K0
256  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
257  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
258  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
259 
260  // p X- pi+ pi0
261  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
262  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
263  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
264 
265  // p X- K+ K0bar
266  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
267  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
268  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
269 
270  // n L K0bar pi+
271  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
272  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
273  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
274 
275  // n S0 K0bar pi+
276  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
277  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
278  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
279 
280  // n S+ K0bar pi0
281  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
282  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
283  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
284 
285  // n S+ K- pi+
286  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
287  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.12, 0.11, 0.10,
288  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
289 
290  // n X0 pi+ pi0
291  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
292  0.0, 0.0, 0.0, 0.0, 0.05, 0.11, 0.14, 0.15, 0.14, 0.13,
293  0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.05, 0.04, 0.03},
294 
295  // n X0 K+ K0bar
296  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
297  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
298  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
299 
300  // n X- 2pi+
301  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
302  0.0, 0.0, 0.0, 0.0, 0.04, 0.10, 0.13, 0.13, 0.12, 0.11,
303  0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.03, 0.03, 0.01},
304 
305  // L L pi+ pi0
306  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
307  0.0, 0.0, 0.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
308  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
309 
310  // L L K+ K0bar
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.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
313  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
314 
315  // L S0 pi+ pi0
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.0, 0.0, 0.05, 0.10, 0.13, 0.12, 0.11, 0.10,
318  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
319 
320  // L S0 K+ K0bar
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.05, 0.10, 0.13, 0.12, 0.11, 0.10,
323  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
324 
325  // L S+ 2pi0
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.07, 0.34, 0.60, 0.87, 1.14, 1.50, 1.40,
328  1.30, 1.20, 1.10, 1.00, 0.85, 0.75, 0.60, 0.50, 0.48, 0.46, 0.30},
329 
330  // L S+ pi+ 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.12, 0.56, 1.00, 1.45, 1.90, 2.50, 2.40,
333  2.30, 2.10, 1.95, 1.80, 1.65, 1.50, 1.35, 1.20, 1.18, 1.16, 0.90},
334 
335  // L S+ K+ K-
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.05, 0.10, 0.13, 0.12, 0.11, 0.10,
338  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
339 
340  // L S+ K0bar K0
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.05, 0.10, 0.13, 0.12, 0.11, 0.10,
343  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01},
344 
345  // L S- 2pi+
346  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
347  0.0, 0.0, 0.0, 0.12, 0.56, 1.00, 1.45, 1.90, 2.50, 2.20,
348  2.05, 1.90, 1.75, 1.60, 1.45, 1.30, 1.10, 1.00, 0.88, 0.86, 0.70},
349 
350  // L X0 K+ pi0
351  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
352  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
353  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
354 
355  // L X0 K0 pi+
356  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
357  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
358  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
359 
360  // L X- K+ pi+
361  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
362  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
363  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
364 
365  // S0 S+ 2pi0
366  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
367  0.0, 0.0, 0.0, 0.04, 0.12, 0.20, 0.28, 0.36, 0.24, 0.20,
368  0.17, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
369 
370  // S0 S+ pi+ pi-
371  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
372  0.0, 0.0, 0.0, 0.06, 0.19, 0.33, 0.46, 0.60, 0.40, 0.25,
373  0.15, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
374 
375  // S0 S+ K+ K-
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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
378  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
379 
380  // S0 S+ K0 K0bar
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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
383  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
384 
385  // S0 S0 pi+ 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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
388  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
389 
390  // S0 S0 K+ K0bar
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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
393  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
394 
395  // S0 S- 2pi+
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.06, 0.19, 0.33, 0.46, 0.60, 0.40, 0.25,
398  0.15, 0.10, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
399 
400  // S0 X0 K0 pi+
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
404 
405  // S0 X0 K+ 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
409 
410  // S0 X- K+ 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
414 
415  // S+ S- pi+ 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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
418  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
419 
420  // S+ S- K+ K0bar
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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
423  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
424 
425  // S+ S+ pi- 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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
428  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
429 
430  // S+ S+ K- K0
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.04, 0.10, 0.13, 0.12, 0.11, 0.10,
433  0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.01, 0.01},
434 
435  // S+ X0 K0 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
439 
440  // S+ X0 K+ pi-
441  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
442  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
443  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
444 
445  // S+ X- K0 pi+
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
449 
450  // S+ X- K+ pi0
451  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
452  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
453  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
454 
455  // S- X0 K+ pi+
456  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
457  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
458  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
459 
460  //
461  // multiplicity 5 (2 channels)
462  //
463  // p X0 3pi0
464  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
465  0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
466  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
467 
468  // p X0 pi+ pi- pi0
469  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
470  0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01,
471  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
472 
473  //
474  // multiplicity 6 (2 channels)
475  //
476  // p X0 4pi0
477  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
478  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
479  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
480 
481  // p X0 pi+ pi- 2pi0
482  { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
483  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01, 0.01, 0.01, 0.01,
484  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
485 
486  //
487  // multiplicity 7 (2 channels)
488  //
489  // p X0 5pi0
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.01, 0.01,
492  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01},
493 
494  // p X0 pi+ pi- 3pi0
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.01, 0.01,
497  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01} };
498 }
499 
500 const G4CascadeXiZeroPChannelData::data_t
501 G4CascadeXiZeroPChannelData::data(x0p2bfs, x0p3bfs, x0p4bfs,
502  x0p5bfs, x0p6bfs, x0p7bfs,
503  x0pCrossSections, xi0*pro, "XiZeroP");