Geant4-11
G4CascadeT1GamNChannel.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//
27
28#include "G4CascadeGamPChannel.hh"
29#include "G4CascadeGamNChannel.hh"
30#include "G4InuclParticleNames.hh"
31using namespace G4InuclParticleNames;
32
33namespace {
34 // gamma p : Outgoing particle types of a given multiplicity
35 static const G4int gamp2bfs[6][2] =
36 {{gam,pro}, {pi0,pro}, {pip,neu}, {kpl,lam}, {kpl,s0}, {k0,sp}};
37
38 static const G4int gamp3bfs[6][3] =
39 {{pip,pro,pim}, {pi0,pro,pi0}, {pip,neu,pi0},
40 {kpl,pro,kmi}, {k0,pro,k0b}, {kpl,neu,k0b}};
41
42 static const G4int gamp4bfs[4][4] =
43 {{pro,pip,pim,pi0}, {pro,pi0,pi0,pi0}, {neu,pip,pip,pim}, {neu,pip,pi0,pi0}};
44
45 static const G4int gamp5bfs[5][5] =
46 {{pro,pip,pip,pim,pim}, {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0},
47 {neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}};
48
49 static const G4int gamp6bfs[6][6] =
50 {{pro,pip,pip,pim,pim,pi0}, {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0},
51 {neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0}};
52
53 static const G4int gamp7bfs[7][7] =
54 {{pro,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pim,pim,pi0,pi0},
55 {pro,pip,pim,pi0,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0,pi0},
56 {neu,pip,pip,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pi0},
57 {neu,pip,pi0,pi0,pi0,pi0,pi0}};
58
59 static const G4int gamp8bfs[7][8] =
60 {{pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0},
61 {pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
62 {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim,pim},
63 {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}};
64
65 static const G4int gamp9bfs[7][9] =
66 {{pro,pip,pip,pip,pim,pim,pim,pi0,pi0}, {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0},
67 {pro,pip,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
68 {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim,pim,pi0},
69 {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}};
70}
71
72namespace {
73 // gamma n : Outgoing particle types of a given multiplicity
74 static const G4int gamn2bfs[6][2] =
75 {{gam,neu}, {neu,pi0}, {pro,pim}, {lam,k0}, {s0, k0}, {sm,kpl}};
76
77 static const G4int gamn3bfs[6][3] =
78 {{neu,pip,pim}, {neu,pi0,pi0}, {pro,pim,pi0},
79 {neu,k0,k0b}, {neu,kpl,kmi}, {pro,k0,kmi}};
80
81 static const G4int gamn4bfs[4][4] =
82 {{neu,pip,pim,pi0}, {neu,pi0,pi0,pi0}, {pro,pip,pim,pim}, {pro,pim,pi0,pi0}};
83
84 static const G4int gamn5bfs[5][5] =
85 {{neu,pip,pip,pim,pim}, {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0},
86 {pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}};
87
88 static const G4int gamn6bfs[6][6] =
89 {{neu,pip,pip,pim,pim,pi0}, {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0},
90 {pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0}};
91
92 static const G4int gamn7bfs[7][7] =
93 {{neu,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pim,pim,pi0,pi0},
94 {neu,pip,pim,pi0,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0,pi0},
95 {pro,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pi0},
96 {pro,pim,pi0,pi0,pi0,pi0,pi0}};
97
98 static const G4int gamn8bfs[7][8] =
99 {{neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0},
100 {neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
101 {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pim},
102 {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}};
103
104 static const G4int gamn9bfs[7][9] =
105 {{neu,pip,pip,pip,pim,pim,pim,pi0,pi0}, {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0},
106 {neu,pip,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
107 {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pim,pi0},
108 {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}};
109}
110
111
112namespace {
113 // Total gamma p cross section as a function of kinetic energy
114 static const G4double gampTotXSec[30] =
115 {0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002,
116 0.0002, 0.0002, 0.0723, 0.2327, 0.54, 0.24, 0.19, 0.28, 0.223,
117 0.18, 0.146, 0.138, 0.1296, 0.1275, 0.124, 0.122, 0.12, 0.1185,
118 0.117, 0.115, 0.115};
119
120 static const G4double gampCrossSections[48][30] = {
121 //
122 // multiplicity 2 (6 channels)
123 //
124 // gamma p elastic ( = gamma n elastic )
125 {0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002,
126 0.0002, 0.0002, 0.0003, 0.0007, 0.0027, 0.0014, 0.0010, 0.0012, 0.0005,
127 0.0003, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0001, 0.0001, 0.0001,
128 0.0001, 0.0001, 0.0001},
129
130 // p pi0 ( = n pi0 ) with new data above 500 MeV
131 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
132 0.0, 0.01, 0.08, 0.3, 0.1143, 0.03, 0.029, 0.0229, 0.0106, 0.005,
133 0.0018, 0.0009, 0.0004, 0.0002, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
134
135 // n pi+ ( = p pi- ) with new data above 800 MeV
136 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
137 0.0, 0.062, 0.152, 0.2373, 0.1143, 0.084, 0.0983, 0.0419, 0.015, 0.006,
138 0.0035, 0.0022, 0.0014, 0.0007, 0.0003, 0.0, 0.0, 0.0, 0.0, 0.0},
139
140 // L K+ ( = L K0 )
141 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
142 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0012, 0.0019, 0.0007,
143 0.0004, 0.0003, 0.0002, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
144
145 // S0 K+ ( = S0 K0 )
146 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
147 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0012, 0.0006,
148 0.0003, 0.0002, 0.0001, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
149
150 // S+ K0 ( = S- K+ )
151 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
152 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0014, 0.0006,
153 0.0003, 0.0002, 0.0001, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
154 //
155 // multiplicity 3 (6 channels)
156 //
157 // p pi+ pi- ( = n pi+ pi- )
158 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
159 0.0, 0.0, 0.0, 0.0, 0.004, 0.045, 0.072, 0.0728, 0.060, 0.041,
160 0.033, 0.026, 0.022, 0.018, 0.016, 0.015, 0.014, 0.0135, 0.013, 0.012},
161
162 // p pi0 pi0 ( = n pi0 pi0 )
163 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
164 0.0, 0.0, 0.0, 0.0, 0.002, 0.009, 0.0299, 0.0265, 0.025, 0.016,
165 0.011, 0.0089, 0.0076, 0.0071, 0.006, 0.0065, 0.0066, 0.007, 0.006, 0.006},
166
167 // n pi+ pi0 ( = p pi- pi0 )
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.004, 0.021, 0.0457, 0.0498, 0.0418, 0.029,
170 0.017, 0.013, 0.0081, 0.0043, 0.0017, 0.0006, 0.0002, 0.0, 0.0, 0.0},
171
172 // p K+ K- ( = n K0 K0b )
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.0, 0.0, 0.0, 0.0, 0.0,
175 0.0004, 0.001, 0.001, 0.0008, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},
176
177 // p K0 K0b (= n K+ K- )
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.0, 0.0, 0.0, 0.0, 0.0,
180 0.0004, 0.001, 0.0009, 0.0006, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},
181
182 // n K+ K0b (= p K0 K- )
183 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
184 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
185 0.0004, 0.001, 0.0009, 0.0006, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},
186 //
187 // multiplicity 4 (4 channels)
188 //
189 // p pi+ pi- pi0 ( = n pi+ pi- pi0 )
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.0021, 0.0041, 0.0132, 0.024,
192 0.025, 0.0225, 0.0196, 0.0162, 0.0125, 0.0088, 0.0059, 0.003, 0.0013, 0.0004},
193
194 // p pi0 pi0 pi0 ( = n pi0 pi0 pi0 )
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.0004, 0.0007, 0.0015, 0.0025,
197 0.004, 0.0038, 0.0033, 0.0027, 0.0021, 0.0015, 0.0009, 0.0005, 0.0002, 0.0001},
198
199 // n pi+ pi+ pi- ( = p pi+ pi- pi- )
200 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
201 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0007, 0.0013, 0.0038, 0.0088,
202 0.011, 0.010, 0.0081, 0.0064, 0.0046, 0.003, 0.0018, 0.0008, 0.0003, 0.0001},
203
204 // n pi+ pi0 pi0 (= p pi- pi0 pi0 )
205 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
206 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0007, 0.0013, 0.0029, 0.006,
207 0.011, 0.008, 0.006, 0.0053, 0.0037, 0.0027, 0.0017, 0.0008, 0.0003, 0.0001},
208 //
209 // multiplicity 5 (5 channels)
210 //
211 // p pi+ pi+ pi- pi- ( = n pi+ pi+ pi- pi- )
212 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
213 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0004, 0.0012,
214 0.0032, 0.005, 0.0055, 0.0052, 0.005, 0.0048, 0.0046, 0.0042, 0.0038, 0.0033},
215
216 // p pi+ pi- pi0 pi0 ( = n pi+ pi- pi0 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.0005, 0.0016,
219 0.0041, 0.0064, 0.0076, 0.0081, 0.0075, 0.008, 0.0086, 0.0084, 0.0076, 0.0066},
220
221 // p pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 )
222 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
223 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0001,
224 0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0008, 0.0007, 0.0006, 0.0005},
225
226 // n pi+ pi+ pi- pi0 ( = p pi+ pi- pi- 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.0005, 0.0016,
229 0.006, 0.0064, 0.0076, 0.0081, 0.009, 0.0093, 0.0086, 0.0084, 0.0076, 0.0066},
230
231 // n pi+ pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 )
232 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
233 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0006,
234 0.0013, 0.0022, 0.0026, 0.0027, 0.0025, 0.0027, 0.0028, 0.0028, 0.0026, 0.0022},
235 //
236 // multiplicity 6 (6 channels)
237 //
238 // p pi+ pi+ pi- pi- pi0 ( = n pi+ pi+ pi- pi- pi0 )
239 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
240 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0002,
241 0.0014, 0.0038, 0.0062, 0.0083, 0.0090, 0.0086, 0.0082, 0.0075, 0.0068, 0.0059},
242
243 // p pi+ pi- pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 )
244 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
245 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
246 0.0009, 0.0017, 0.0029, 0.005, 0.0045, 0.0049, 0.0051, 0.0050, 0.0046, 0.0040},
247
248 // p pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 )
249 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
250 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
251 0.0, 0.0001, 0.0002, 0.0002, 0.0002, 0.0003, 0.0003, 0.0002, 0.0002, 0.0002},
252
253 // n pi+ pi+ pi+ pi- pi- ( = p pi+ pi+ pi- pi- pi- )
254 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
255 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0001,
256 0.0003, 0.0009, 0.0026, 0.0032, 0.0035, 0.0033, 0.0031, 0.0028, 0.0025, 0.0022},
257
258 // n pi+ pi+ pi- pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 )
259 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
260 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0002,
261 0.0005, 0.0017, 0.0054, 0.0073, 0.0079, 0.0085, 0.0088, 0.0084, 0.0075, 0.0066},
262
263 // n pi+ pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 )
264 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
265 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
266 0.0001, 0.0003, 0.0009, 0.0012, 0.0013, 0.0013, 0.0014, 0.0014, 0.0012, 0.0011},
267 //
268 // multiplicity 7 (7 channels)
269 //
270 // p pi+ pi+ pi+ pi- pi- pi- ( = n pi+ pi+ pi+ pi- pi- 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
273 0.0, 0.0001, 0.0004, 0.0006, 0.0009, 0.0009, 0.0009, 0.0009, 0.0008, 0.0008},
274
275 // p pi+ pi+ pi- pi- pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 )
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
278 0.0, 0.0003, 0.0012, 0.0021, 0.004, 0.004, 0.0047, 0.005, 0.0045, 0.004},
279
280 // p pi+ pi- pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
283 0.0, 0.0001, 0.0004, 0.0007, 0.0013, 0.0014, 0.0016, 0.0017, 0.0015, 0.0013},
284
285 // p pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 )
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
288 0.0, 0.0, 0.0, 0.0, 0.0, 0.0001, 0.0001, 0.0001, 0.0, 0.0},
289
290 // n pi+ pi+ pi- pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
293 0.0, 0.0003, 0.0008, 0.0014, 0.002, 0.0023, 0.0025, 0.0027, 0.0025, 0.0025},
294
295 // n pi+ pi+ pi+ pi- pi- pi0 ( = p pi+ pi+ pi- pi- pi- pi0 )
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.0003, 0.0008, 0.0014, 0.002, 0.0023, 0.0025, 0.0027, 0.0025, 0.0025},
299
300 // n pi+ pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 )
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.0, 0.0, 0.0, 0.0, 0.0, 0.0,
303 0.0, 0.0, 0.0001, 0.0002, 0.0002, 0.0003, 0.0003, 0.0003, 0.0002, 0.0002},
304 //
305 // multiplicity 8 (7 channels)
306 //
307 // p pi+ pi+ pi- pi- pi0 pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 pi0 )
308 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.0001, 0.0003, 0.0006, 0.0015, 0.0026, 0.0033, 0.0034, 0.0034, 0.0033},
311
312 // p pi+ pi+ pi+ pi- pi- pi- pi0 ( = n pi+ pi+ pi+ pi- pi- pi- pi0 )
313 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.0001, 0.0002, 0.0005, 0.0011, 0.0017, 0.0023, 0.0022, 0.0022, 0.0022},
316
317 // p pi+ pi- pi0 pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 pi0 pi0 )
318 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
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.0001, 0.0002, 0.0004, 0.0005, 0.0007, 0.0007, 0.0007, 0.0007},
321
322 // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
323 // negligible
324
325 // n pi+ pi+ pi+ pi- pi- pi0 pi0 ( = p pi+ pi+ pi- pi- pi- pi0 pi0 )
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.0, 0.0,
328 0.0, 0.0001, 0.0003, 0.0006, 0.002, 0.0026, 0.0033, 0.0034, 0.0034, 0.0033},
329
330 // n pi+ pi+ pi- pi0 pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 pi0 pi0 )
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.0, 0.0,
333 0.0, 0.0, 0.0001, 0.0003, 0.001, 0.0013, 0.0017, 0.0017, 0.0017, 0.0016},
334
335 // n pi+ pi+ pi+ pi+ pi- pi- pi- ( = p pi+ pi+ pi+ pi- pi- pi- 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.0, 0.0,
338 0.0, 0.0, 0.0, 0.0001, 0.0002, 0.0004, 0.0006, 0.0005, 0.0005, 0.0005},
339
340 // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 pi0 )
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.0,
343 0.0, 0.0, 0.0, 0.0, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001},
344 //
345 // multiplicity 9 (7 channels)
346 //
347 // p pi+ pi+ pi+ pi- pi- pi- pi0 pi0 ( = n pi+ pi+ pi+ pi- pi- pi- pi0 pi0 )
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.0, 0.0, 0.0, 0.0, 0.0,
350 0.0, 0.0, 0.0002, 0.0005, 0.0013, 0.0019, 0.0023, 0.0026, 0.0026, 0.0024},
351
352 // p pi+ pi+ pi- pi- pi0 pi0 pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 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.0, 0.0, 0.0, 0.0, 0.0,
355 0.0, 0.0, 0.0002, 0.0003, 0.0011, 0.0014, 0.0018, 0.0020, 0.0020, 0.0018},
356
357 // p pi+ pi+ pi+ pi+ pi- pi- pi- pi- ( = n pi+ pi+ pi+ pi+ pi- pi- pi- pi- )
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.0, 0.0, 0.0, 0.0, 0.0,
360 0.0, 0.0, 0.0, 0.0001, 0.0002, 0.0003, 0.0003, 0.0003, 0.0003, 0.0003},
361
362 // p pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 )
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.0, 0.0, 0.0001, 0.0001, 0.0002, 0.0003, 0.0002, 0.0002, 0.0002},
366
367 // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
368 // negligible
369
370 // n pi+ pi+ pi+ pi- pi- pi0 pi0 pi0 ( = p pi+ pi+ pi- pi- pi- pi0 pi0 pi0 )
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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
373 0.0, 0.0, 0.0002, 0.0005, 0.0013, 0.0018, 0.0023, 0.0026, 0.0026, 0.0024},
374
375 // n pi+ pi+ pi+ pi+ pi- pi- pi- pi0 ( = p pi+ pi+ pi+ pi- pi- pi- pi- pi0 )
376 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
377 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
378 0.0, 0.0, 0.0001, 0.0002, 0.0007, 0.0009, 0.0011, 0.0013, 0.0013, 0.0012},
379
380 // n pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0 )
381 {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
382 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
383 0.0, 0.0, 0.0001, 0.0002, 0.0004, 0.0006, 0.0008, 0.0009, 0.0009, 0.0008}};
384
385 // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
386 // negligible
387}
388
389// Initialize both |T Tz> = |1/2 +-1/2> channels using gamma-p cross section
390// tables
391
392const G4CascadeGamPChannelData::data_t
393G4CascadeGamPChannelData::data(gamp2bfs, gamp3bfs, gamp4bfs, gamp5bfs,
394 gamp6bfs, gamp7bfs, gamp8bfs, gamp9bfs,
395 gampCrossSections, gampTotXSec, gam*pro,
396 "GammaProton");
397
398const G4CascadeGamNChannelData::data_t
399G4CascadeGamNChannelData::data(gamn2bfs, gamn3bfs, gamn4bfs, gamn5bfs,
400 gamn6bfs, gamn7bfs, gamn8bfs, gamn9bfs,
401 gampCrossSections, gampTotXSec, gam*neu,
402 "GammaNeutron");
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp5bfs
static const G4int gamp5bfs[5][5]
Definition: G4CascadeT1GamNChannel.cc:45
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn5bfs
static const G4int gamn5bfs[5][5]
Definition: G4CascadeT1GamNChannel.cc:84
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn2bfs
static const G4int gamn2bfs[6][2]
Definition: G4CascadeT1GamNChannel.cc:74
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn9bfs
static const G4int gamn9bfs[7][9]
Definition: G4CascadeT1GamNChannel.cc:104
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp3bfs
static const G4int gamp3bfs[6][3]
Definition: G4CascadeT1GamNChannel.cc:38
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn8bfs
static const G4int gamn8bfs[7][8]
Definition: G4CascadeT1GamNChannel.cc:98
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gampTotXSec
static const G4double gampTotXSec[30]
Definition: G4CascadeT1GamNChannel.cc:114
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn6bfs
static const G4int gamn6bfs[6][6]
Definition: G4CascadeT1GamNChannel.cc:88
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn4bfs
static const G4int gamn4bfs[4][4]
Definition: G4CascadeT1GamNChannel.cc:81
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn3bfs
static const G4int gamn3bfs[6][3]
Definition: G4CascadeT1GamNChannel.cc:77
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp7bfs
static const G4int gamp7bfs[7][7]
Definition: G4CascadeT1GamNChannel.cc:53
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gampCrossSections
static const G4double gampCrossSections[48][30]
Definition: G4CascadeT1GamNChannel.cc:120
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp2bfs
static const G4int gamp2bfs[6][2]
Definition: G4CascadeT1GamNChannel.cc:35
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp8bfs
static const G4int gamp8bfs[7][8]
Definition: G4CascadeT1GamNChannel.cc:59
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp6bfs
static const G4int gamp6bfs[6][6]
Definition: G4CascadeT1GamNChannel.cc:49
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamn7bfs
static const G4int gamn7bfs[7][7]
Definition: G4CascadeT1GamNChannel.cc:92
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp4bfs
static const G4int gamp4bfs[4][4]
Definition: G4CascadeT1GamNChannel.cc:42
anonymous_namespace{G4CascadeT1GamNChannel.cc}::gamp9bfs
static const G4int gamp9bfs[7][9]
Definition: G4CascadeT1GamNChannel.cc:65
Generated by doxygen 1.9.3