Geant4-11
G4XNNstarTable.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 "globals.hh"
29#include "G4ios.hh"
30#include "G4SystemOfUnits.hh"
31#include "G4XNNstarTable.hh"
32#include "G4PhysicsFreeVector.hh"
33
34const G4int G4XNNstarTable::sizeNNstar = 121;
35
36// Energies (GeV) corresponding to the cross section table
37// Units are assigned when filling the PhysicsVector
38
39const G4double G4XNNstarTable::energyTable[121] =
40{
41 0.0,
42 2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.031, 2.037, 2.044, 2.052,
43 2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.135, 2.151, 2.168, 2.185,
44 2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.335, 2.360, 2.386, 2.413,
45 2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.629, 2.664, 2.699, 2.736,
46 2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.018, 3.062, 3.107, 3.153,
47 3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.502, 3.555, 3.610, 3.666,
48 3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.081, 4.144, 4.208, 4.273,
49 4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.754, 4.827, 4.900, 4.975,
50 5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 11.807, 12.941, 14.076, 15.210,
51 16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555,
52 27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899,
53 39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
54};
55
56// Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998
57// Units are assigned when filling the PhysicsVector
58
59const G4double G4XNNstarTable::sigmaNN1440[121] =
60{
61 0.0,
62 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
63 0.000, 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.004, 0.006, 0.009,
64 0.014, 0.020, 0.031, 0.047, 0.072, 0.113, 0.178, 0.266, 0.370, 0.476,
65 0.575, 0.665, 0.744, 0.814, 0.874, 0.926, 0.971, 1.009, 1.040, 1.066,
66 1.087, 1.103, 1.115, 1.124, 1.129, 1.132, 1.131, 1.129, 1.124, 1.117,
67 1.109, 1.099, 1.088, 1.075, 1.062, 1.048, 1.033, 1.017, 1.001, 0.984,
68 0.967, 0.950, 0.932, 0.914, 0.896, 0.879, 0.861, 0.843, 0.825, 0.807,
69 0.790, 0.773, 0.755, 0.738, 0.722, 0.705, 0.689, 0.673, 0.657, 0.641,
70 0.636, 0.453, 0.336, 0.258, 0.204, 0.166, 0.137, 0.115, 0.098, 0.084,
71 0.073, 0.064, 0.057, 0.051, 0.046, 0.041, 0.037, 0.034, 0.031, 0.029,
72 0.026, 0.024, 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
73 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008
74};
75
76const G4double G4XNNstarTable::sigmaNN1520[121] =
77{
78 0.0,
79 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
80 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001, 0.001,
81 0.001, 0.002, 0.003, 0.005, 0.008, 0.014, 0.024, 0.043, 0.079, 0.150,
82 0.280, 0.465, 0.662, 0.841, 0.995, 1.125, 1.235, 1.327, 1.410, 1.474,
83 1.525, 1.567, 1.600, 1.625, 1.644, 1.651, 1.657, 1.659, 1.657, 1.651,
84 1.642, 1.629, 1.614, 1.597, 1.578, 1.558, 1.536, 1.512, 1.488, 1.463,
85 1.437, 1.411, 1.384, 1.357, 1.329, 1.302, 1.274, 1.247, 1.219, 1.192,
86 1.165, 1.138, 1.112, 1.086, 1.060, 1.035, 1.010, 0.985, 0.961, 0.938,
87 0.930, 0.652, 0.479, 0.365, 0.287, 0.232, 0.191, 0.160, 0.135, 0.116,
88 0.101, 0.089, 0.078, 0.070, 0.062, 0.056, 0.051, 0.046, 0.042, 0.039,
89 0.036, 0.033, 0.031, 0.028, 0.026, 0.025, 0.023, 0.022, 0.020, 0.019,
90 0.018, 0.017, 0.016, 0.015, 0.015, 0.014, 0.013, 0.013, 0.012, 0.011
91};
92
93const G4double G4XNNstarTable::sigmaNN1535[121] =
94{
95 0.0,
96 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003,
97 0.004, 0.005, 0.006, 0.008, 0.010, 0.012, 0.015, 0.019,
98 0.024, 0.031, 0.039, 0.052, 0.069, 0.097, 0.145, 0.216,
99 0.298, 0.378, 0.451, 0.513, 0.566, 0.610, 0.646, 0.675,
100 0.699, 0.718, 0.732, 0.742, 0.749, 0.753, 0.754, 0.753,
101 0.751, 0.746, 0.740, 0.733, 0.724, 0.715, 0.705, 0.694,
102 0.683, 0.671, 0.659, 0.647, 0.634, 0.621, 0.608, 0.595,
103 0.582, 0.569, 0.556, 0.543, 0.531, 0.518, 0.506, 0.493,
104 0.481, 0.469, 0.458, 0.446, 0.435, 0.424, 0.413, 0.402,
105 0.399, 0.276, 0.201, 0.153, 0.120, 0.096, 0.079, 0.066,
106 0.056, 0.048, 0.042, 0.037, 0.032, 0.029, 0.026, 0.023,
107 0.021, 0.019, 0.017, 0.016, 0.015, 0.014, 0.013, 0.012,
108 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007, 0.007,
109 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005
110};
111
112const G4double G4XNNstarTable::sigmaNN1650[121] =
113{
114 0.0,
115 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
116 0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.001,
117 0.002, 0.002, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007,
118 0.008, 0.010, 0.012, 0.015, 0.019, 0.023, 0.029, 0.038,
119 0.051, 0.071, 0.104, 0.150, 0.201, 0.249, 0.290, 0.326,
120 0.354, 0.378, 0.397, 0.412, 0.424, 0.434, 0.440, 0.445,
121 0.448, 0.449, 0.449, 0.448, 0.445, 0.442, 0.438, 0.433,
122 0.428, 0.422, 0.416, 0.409, 0.403, 0.395, 0.388, 0.381,
123 0.373, 0.366, 0.358, 0.350, 0.343, 0.335, 0.327, 0.320,
124 0.312, 0.305, 0.298, 0.291, 0.284, 0.277, 0.270, 0.263,
125 0.261, 0.182, 0.133, 0.102, 0.080, 0.064, 0.053, 0.044,
126 0.037, 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015,
127 0.014, 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008,
128 0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
129 0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
130};
131
132const G4double G4XNNstarTable::sigmaNN1675[121] =
133{
134 0.0,
135 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
136 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
137 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
138 0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.009, 0.015,
139 0.026, 0.048, 0.095, 0.189, 0.324, 0.463, 0.586, 0.691,
140 0.780, 0.855, 0.919, 0.972, 1.016, 1.052, 1.081, 1.105,
141 1.123, 1.136, 1.145, 1.151, 1.153, 1.153, 1.149, 1.144,
142 1.136, 1.127, 1.116, 1.104, 1.090, 1.076, 1.061, 1.045,
143 1.028, 1.011, 0.993, 0.975, 0.957, 0.939, 0.921, 0.902,
144 0.884, 0.865, 0.847, 0.828, 0.810, 0.792, 0.775, 0.757,
145 0.751, 0.538, 0.399, 0.307, 0.242, 0.196, 0.162, 0.136,
146 0.115, 0.099, 0.086, 0.076, 0.067, 0.060, 0.053, 0.048,
147 0.044, 0.040, 0.036, 0.033, 0.031, 0.028, 0.026, 0.024,
148 0.023, 0.021, 0.020, 0.019, 0.018, 0.016, 0.016, 0.015,
149 0.014, 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010
150};
151
152const G4double G4XNNstarTable::sigmaNN1680[121] =
153{
154 0.0,
155 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
156 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
157 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
158 0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.006, 0.010,
159 0.018, 0.035, 0.073, 0.156, 0.294, 0.446, 0.580, 0.693,
160 0.788, 0.867, 0.933, 0.988, 1.033, 1.070, 1.100, 1.124,
161 1.142, 1.155, 1.163, 1.168, 1.170, 1.168, 1.164, 1.158,
162 1.149, 1.139, 1.127, 1.114, 1.100, 1.085, 1.068, 1.052,
163 1.034, 1.016, 0.998, 0.979, 0.960, 0.941, 0.922, 0.903,
164 0.884, 0.865, 0.846, 0.827, 0.809, 0.790, 0.772, 0.754,
165 0.748, 0.533, 0.394, 0.301, 0.238, 0.192, 0.158, 0.133,
166 0.113, 0.097, 0.084, 0.074, 0.065, 0.058, 0.052, 0.047,
167 0.042, 0.039, 0.035, 0.032, 0.030, 0.028, 0.026, 0.024,
168 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
169 0.013, 0.013, 0.012, 0.012, 0.011, 0.010, 0.010, 0.010
170};
171
172const G4double G4XNNstarTable::sigmaNN1700[121] =
173{
174 0.0,
175 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
176 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
177 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
178 0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.005, 0.008,
179 0.013, 0.022, 0.038, 0.070, 0.140, 0.245, 0.347, 0.431,
180 0.500, 0.556, 0.601, 0.637, 0.666, 0.689, 0.708, 0.721,
181 0.731, 0.738, 0.742, 0.743, 0.743, 0.740, 0.736, 0.730,
182 0.723, 0.716, 0.707, 0.697, 0.687, 0.676, 0.665, 0.654,
183 0.642, 0.630, 0.617, 0.605, 0.593, 0.580, 0.568, 0.555,
184 0.543, 0.531, 0.519, 0.507, 0.495, 0.483, 0.471, 0.460,
185 0.456, 0.322, 0.236, 0.180, 0.142, 0.114, 0.094, 0.079,
186 0.067, 0.057, 0.050, 0.044, 0.039, 0.034, 0.031, 0.028,
187 0.025, 0.023, 0.021, 0.019, 0.018, 0.016, 0.015, 0.014,
188 0.013, 0.012, 0.011, 0.011, 0.010, 0.009, 0.009, 0.008,
189 0.008, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006
190};
191
192const G4double G4XNNstarTable::sigmaNN1710[121] =
193{
194 0.0,
195 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
196 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
197 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
198 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.003,
199 0.005, 0.008, 0.013, 0.025, 0.052, 0.096, 0.144, 0.185,
200 0.219, 0.246, 0.269, 0.288, 0.304, 0.316, 0.326, 0.334,
201 0.340, 0.344, 0.347, 0.349, 0.349, 0.349, 0.348, 0.346,
202 0.343, 0.340, 0.336, 0.332, 0.328, 0.323, 0.319, 0.313,
203 0.308, 0.303, 0.297, 0.292, 0.286, 0.280, 0.275, 0.269,
204 0.263, 0.257, 0.252, 0.246, 0.241, 0.235, 0.230, 0.224,
205 0.223, 0.158, 0.117, 0.090, 0.071, 0.057, 0.047, 0.040,
206 0.034, 0.029, 0.025, 0.022, 0.019, 0.017, 0.016, 0.014,
207 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
208 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
209 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
210};
211
212const G4double G4XNNstarTable::sigmaNN1720[121] =
213{
214 0.0,
215 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
216 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
217 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
218 0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.007, 0.010,
219 0.015, 0.023, 0.036, 0.061, 0.106, 0.174, 0.250, 0.321,
220 0.382, 0.434, 0.477, 0.513, 0.543, 0.568, 0.589, 0.605,
221 0.617, 0.627, 0.634, 0.639, 0.641, 0.642, 0.641, 0.639,
222 0.635, 0.630, 0.625, 0.618, 0.611, 0.603, 0.595, 0.586,
223 0.576, 0.567, 0.557, 0.547, 0.537, 0.527, 0.516, 0.506,
224 0.496, 0.485, 0.475, 0.465, 0.455, 0.444, 0.435, 0.425,
225 0.421, 0.302, 0.224, 0.172, 0.136, 0.117, 0.091, 0.076,
226 0.065, 0.056, 0.049, 0.043, 0.038, 0.034, 0.030, 0.027,
227 0.025, 0.022, 0.020, 0.019, 0.017, 0.016, 0.015, 0.014,
228 0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008,
229 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006, 0.006
230};
231
232const G4double G4XNNstarTable::sigmaNN1900[121] =
233{
234 0.0,
235 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
236 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
237 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
238 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.003, 0.005,
239 0.006, 0.008, 0.010, 0.014, 0.019, 0.026, 0.037, 0.054,
240 0.074, 0.094, 0.114, 0.131, 0.147, 0.161, 0.173, 0.184,
241 0.193, 0.201, 0.208, 0.213, 0.218, 0.221, 0.224, 0.226,
242 0.228, 0.229, 0.229, 0.229, 0.229, 0.228, 0.227, 0.225,
243 0.223, 0.221, 0.219, 0.217, 0.214, 0.212, 0.209, 0.206,
244 0.203, 0.200, 0.197, 0.194, 0.190, 0.187, 0.184, 0.181,
245 0.180, 0.137, 0.106, 0.083, 0.067, 0.056, 0.047, 0.039,
246 0.034, 0.029, 0.026, 0.023, 0.020, 0.018, 0.016, 0.015,
247 0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.008, 0.008,
248 0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
249 0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
250};
251
252const G4double G4XNNstarTable::sigmaNN1990[121] =
253{
254 0.0,
255 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
256 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
257 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
258 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
259 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009,
260 0.017, 0.030, 0.051, 0.076, 0.102, 0.127, 0.149, 0.168,
261 0.185, 0.199, 0.211, 0.221, 0.230, 0.237, 0.243, 0.247,
262 0.251, 0.254, 0.256, 0.258, 0.258, 0.259, 0.258, 0.258,
263 0.257, 0.255, 0.254, 0.252, 0.249, 0.247, 0.245, 0.242,
264 0.239, 0.236, 0.233, 0.230, 0.226, 0.223, 0.220, 0.216,
265 0.215, 0.167, 0.131, 0.104, 0.085, 0.070, 0.059, 0.050,
266 0.043, 0.038, 0.033, 0.029, 0.026, 0.023, 0.021, 0.019,
267 0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.010, 0.010,
268 0.009, 0.009, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006,
269 0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004
270};
271
272const G4double G4XNNstarTable::sigmaNN2090[121] =
273{
274 0.0,
275 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
276 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
277 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
278 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
279 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.005,
280 0.007, 0.011, 0.016, 0.024, 0.036, 0.053, 0.071, 0.089,
281 0.106, 0.120, 0.133, 0.142, 0.151, 0.158, 0.164, 0.169,
282 0.172, 0.175, 0.178, 0.179, 0.180, 0.180, 0.180, 0.180,
283 0.179, 0.178, 0.176, 0.175, 0.173, 0.171, 0.169, 0.166,
284 0.164, 0.162, 0.159, 0.156, 0.154, 0.151, 0.148, 0.145,
285 0.144, 0.107, 0.081, 0.063, 0.050, 0.041, 0.034, 0.028,
286 0.024, 0.021, 0.018, 0.016, 0.014, 0.013, 0.011, 0.010,
287 0.009, 0.008, 0.008, 0.007, 0.006, 0.006, 0.006, 0.005,
288 0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003,
289 0.003, 0.003, 0.003, 0.003, 0.002, 0.002, 0.002, 0.002
290};
291
292const G4double G4XNNstarTable::sigmaNN2190[121] =
293{
294 0.0,
295 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
296 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
297 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
298 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
299 0.001, 0.002, 0.003, 0.005, 0.009, 0.015, 0.024, 0.036,
300 0.050, 0.064, 0.078, 0.090, 0.100, 0.110, 0.118, 0.125,
301 0.131, 0.136, 0.140, 0.143, 0.146, 0.149, 0.150, 0.152,
302 0.153, 0.153, 0.153, 0.153, 0.153, 0.152, 0.152, 0.151,
303 0.150, 0.148, 0.147, 0.145, 0.144, 0.142, 0.140, 0.139,
304 0.138, 0.110, 0.087, 0.069, 0.056, 0.047, 0.039, 0.033,
305 0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.014, 0.013,
306 0.012, 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006,
307 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004,
308 0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003
309};
310
311const G4double G4XNNstarTable::sigmaNN2220[121] =
312{
313 0.0,
314 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
315 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
316 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
317 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
318 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
319 0.001, 0.001, 0.001, 0.003, 0.005, 0.008, 0.015, 0.026,
320 0.041, 0.058, 0.074, 0.089, 0.101, 0.113, 0.122, 0.131,
321 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161, 0.163,
322 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164, 0.163,
323 0.162, 0.161, 0.160, 0.158, 0.157, 0.155, 0.153, 0.151,
324 0.150, 0.121, 0.096, 0.077, 0.062, 0.052, 0.043, 0.037,
325 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
326 0.013, 0.011, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
327 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
328 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003
329};
330
331const G4double G4XNNstarTable::sigmaNN2250[121] =
332{
333 0.0,
334 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
335 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
336 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
337 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
338 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
339 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009, 0.016,
340 0.027, 0.043, 0.060, 0.076, 0.090, 0.103, 0.114, 0.123,
341 0.131, 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161,
342 0.163, 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164,
343 0.163, 0.162, 0.161, 0.159, 0.158, 0.156, 0.154, 0.152,
344 0.152, 0.122, 0.096, 0.077, 0.062, 0.051, 0.043, 0.037,
345 0.031, 0.027, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
346 0.012, 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007,
347 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004,
348 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
349};
350
351
352G4XNNstarTable::G4XNNstarTable()
353{
354 xMap["N(1440)0"] = (G4double*) sigmaNN1440;
355 xMap["N(1440)+"] = (G4double*) sigmaNN1440;
356
357 xMap["N(1520)0"] = (G4double*) sigmaNN1520;
358 xMap["N(1520)+"] = (G4double*) sigmaNN1520;
359
360 xMap["N(1535)0"] = (G4double*) sigmaNN1535;
361 xMap["N(1535)+"] = (G4double*) sigmaNN1535;
362
363 xMap["N(1650)0"] = (G4double*) sigmaNN1650;
364 xMap["N(1650)+"] = (G4double*) sigmaNN1650;
365
366 xMap["N(1675)0"] = (G4double*) sigmaNN1675;
367 xMap["N(1675)+"] = (G4double*) sigmaNN1675;
368
369 xMap["N(1680)0"] = (G4double*) sigmaNN1680;
370 xMap["N(1680)+"] = (G4double*) sigmaNN1680;
371
372 xMap["N(1700)0"] = (G4double*) sigmaNN1700;
373 xMap["N(1700)+"] = (G4double*) sigmaNN1700;
374
375 xMap["N(1710)0"] = (G4double*) sigmaNN1710;
376 xMap["N(1710)+"] = (G4double*) sigmaNN1710;
377
378 xMap["N(1720)0"] = (G4double*) sigmaNN1720;
379 xMap["N(1720)+"] = (G4double*) sigmaNN1720;
380
381 xMap["N(1900)0"] = (G4double*) sigmaNN1900;
382 xMap["N(1900)+"] = (G4double*) sigmaNN1900;
383
384 xMap["N(1990)0"] = (G4double*) sigmaNN1990;
385 xMap["N(1990)+"] = (G4double*) sigmaNN1990;
386
387 xMap["N(2090)0"] = (G4double*) sigmaNN2090;
388 xMap["N(2090)+"] = (G4double*) sigmaNN2090;
389
390 xMap["N(2190)0"] = (G4double*) sigmaNN2190;
391 xMap["N(2190)+"] = (G4double*) sigmaNN2190;
392
393 xMap["N(2220)0"] = (G4double*) sigmaNN2220;
394 xMap["N(2220)+"] = (G4double*) sigmaNN2220;
395
396 xMap["N(2250)0"] = (G4double*) sigmaNN2250;
397 xMap["N(2250)+"] = (G4double*) sigmaNN2250;
398}
399
400
401G4XNNstarTable::~G4XNNstarTable()
402{ }
403
404
405const G4PhysicsVector* G4XNNstarTable::CrossSectionTable(const G4String& particleName) const
406{
407 // NOTE: the returned pointer is owned by the client
408
409 if (xMap.find(particleName) != xMap.end())
410 {
411 // Cross section table for the requested particle available in the Map
412 G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeNNstar);
413 std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
414 G4double* sigmaPointer = 0;
415 for (iter = xMap.begin(); iter != xMap.end(); ++iter)
416 {
417 G4String str = (*iter).first;
418 if (str == particleName)
419 {
420 sigmaPointer = (*iter).second;
421 }
422 }
423 G4int i;
424 for (i=0; i<sizeNNstar; i++)
425 {
426 G4double value = *(sigmaPointer + i) * millibarn;
427 G4double energy = energyTable[i] * GeV;
428 sigmaVector->PutValue(i,energy,value);
429 }
430 return sigmaVector;
431 }
432 else
433 // No cross section table for the requested particle is available in the Map
434 return 0;
435}
436
437
438
millibarn
static constexpr double millibarn
Definition: G4SIunits.hh:86
GeV
static constexpr double GeV
Definition: G4SIunits.hh:203
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4PhysicsFreeVector::PutValue
void PutValue(const std::size_t index, const G4double e, const G4double value)
Definition: G4PhysicsFreeVector.hh:95
G4XNNstarTable::sigmaNN1710
static const G4double sigmaNN1710[121]
Definition: G4XNNstarTable.hh:78
G4XNNstarTable::sigmaNN1650
static const G4double sigmaNN1650[121]
Definition: G4XNNstarTable.hh:74
G4XNNstarTable::sigmaNN1990
static const G4double sigmaNN1990[121]
Definition: G4XNNstarTable.hh:81
G4XNNstarTable::sigmaNN1700
static const G4double sigmaNN1700[121]
Definition: G4XNNstarTable.hh:77
G4XNNstarTable::sigmaNN2250
static const G4double sigmaNN2250[121]
Definition: G4XNNstarTable.hh:85
G4XNNstarTable::sigmaNN2220
static const G4double sigmaNN2220[121]
Definition: G4XNNstarTable.hh:84
G4XNNstarTable::sigmaNN1720
static const G4double sigmaNN1720[121]
Definition: G4XNNstarTable.hh:79
G4XNNstarTable::sigmaNN1680
static const G4double sigmaNN1680[121]
Definition: G4XNNstarTable.hh:76
G4XNNstarTable::sigmaNN1900
static const G4double sigmaNN1900[121]
Definition: G4XNNstarTable.hh:80
G4XNNstarTable::CrossSectionTable
virtual const G4PhysicsVector * CrossSectionTable(const G4String &particleName) const
Definition: G4XNNstarTable.cc:405
G4XNNstarTable::sizeNNstar
static const G4int sizeNNstar
Definition: G4XNNstarTable.hh:65
G4XNNstarTable::~G4XNNstarTable
virtual ~G4XNNstarTable()
Definition: G4XNNstarTable.cc:401
G4XNNstarTable::sigmaNN1675
static const G4double sigmaNN1675[121]
Definition: G4XNNstarTable.hh:75
G4XNNstarTable::sigmaNN1535
static const G4double sigmaNN1535[121]
Definition: G4XNNstarTable.hh:73
G4XNNstarTable::sigmaNN1520
static const G4double sigmaNN1520[121]
Definition: G4XNNstarTable.hh:72
G4XNNstarTable::xMap
std::map< G4String, G4double *, std::less< G4String > > xMap
Definition: G4XNNstarTable.hh:63
G4XNNstarTable::sigmaNN2190
static const G4double sigmaNN2190[121]
Definition: G4XNNstarTable.hh:83
G4XNNstarTable::sigmaNN1440
static const G4double sigmaNN1440[121]
Definition: G4XNNstarTable.hh:71
G4XNNstarTable::energyTable
static const G4double energyTable[121]
Definition: G4XNNstarTable.hh:68
G4XNNstarTable::sigmaNN2090
static const G4double sigmaNN2090[121]
Definition: G4XNNstarTable.hh:82
G4INCL::KinematicsUtils::energy
G4double energy(const ThreeVector &p, const G4double m)
Definition: G4INCLKinematicsUtils.cc:158
Generated by doxygen 1.9.3