Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

G4XNNstarTable.cc@ 1007

Last change on this file since 1007 was 819, checked in by garnier, 17 years ago
import all except CVS
File size: 19.1 KB

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

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: trunk/source/processes/hadronic/models/im_r_matrix/src/G4XNNstarTable.cc@ 1007

Download in other formats: