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

Last change on this file since 836 was 819, checked in by garnier, 16 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
34const G4int G4XNNstarTable::sizeNNstar = 121;
35
36
37// Energies (GeV) corresponding to the cross section table
38// Units are assigned while filling the PhysicsVector
39
40const 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
60const 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
77const 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
94const 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
113const 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
133const 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
153const 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
173const 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
193const 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
213const 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
233const 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
253const 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
273const 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
293const 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
312const 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
332const 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
353G4XNNstarTable::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
402G4XNNstarTable::~G4XNNstarTable()
403{ }
404
405
406const 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.