00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #include "globals.hh"
00029 #include "G4ios.hh"
00030 #include "G4SystemOfUnits.hh"
00031 #include "G4XNNstarTable.hh"
00032 #include "G4PhysicsFreeVector.hh"
00033
00034 const G4int G4XNNstarTable::sizeNNstar = 121;
00035
00036
00037
00038
00039 const G4double G4XNNstarTable::energyTable[121] =
00040 {
00041 0.0,
00042 2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.031, 2.037, 2.044, 2.052,
00043 2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.135, 2.151, 2.168, 2.185,
00044 2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.335, 2.360, 2.386, 2.413,
00045 2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.629, 2.664, 2.699, 2.736,
00046 2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.018, 3.062, 3.107, 3.153,
00047 3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.502, 3.555, 3.610, 3.666,
00048 3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.081, 4.144, 4.208, 4.273,
00049 4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.754, 4.827, 4.900, 4.975,
00050 5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 11.807, 12.941, 14.076, 15.210,
00051 16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555,
00052 27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899,
00053 39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
00054 };
00055
00056
00057
00058
00059 const G4double G4XNNstarTable::sigmaNN1440[121] =
00060 {
00061 0.0,
00062 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00063 0.000, 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.004, 0.006, 0.009,
00064 0.014, 0.020, 0.031, 0.047, 0.072, 0.113, 0.178, 0.266, 0.370, 0.476,
00065 0.575, 0.665, 0.744, 0.814, 0.874, 0.926, 0.971, 1.009, 1.040, 1.066,
00066 1.087, 1.103, 1.115, 1.124, 1.129, 1.132, 1.131, 1.129, 1.124, 1.117,
00067 1.109, 1.099, 1.088, 1.075, 1.062, 1.048, 1.033, 1.017, 1.001, 0.984,
00068 0.967, 0.950, 0.932, 0.914, 0.896, 0.879, 0.861, 0.843, 0.825, 0.807,
00069 0.790, 0.773, 0.755, 0.738, 0.722, 0.705, 0.689, 0.673, 0.657, 0.641,
00070 0.636, 0.453, 0.336, 0.258, 0.204, 0.166, 0.137, 0.115, 0.098, 0.084,
00071 0.073, 0.064, 0.057, 0.051, 0.046, 0.041, 0.037, 0.034, 0.031, 0.029,
00072 0.026, 0.024, 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
00073 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008
00074 };
00075
00076 const G4double G4XNNstarTable::sigmaNN1520[121] =
00077 {
00078 0.0,
00079 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00080 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001, 0.001,
00081 0.001, 0.002, 0.003, 0.005, 0.008, 0.014, 0.024, 0.043, 0.079, 0.150,
00082 0.280, 0.465, 0.662, 0.841, 0.995, 1.125, 1.235, 1.327, 1.410, 1.474,
00083 1.525, 1.567, 1.600, 1.625, 1.644, 1.651, 1.657, 1.659, 1.657, 1.651,
00084 1.642, 1.629, 1.614, 1.597, 1.578, 1.558, 1.536, 1.512, 1.488, 1.463,
00085 1.437, 1.411, 1.384, 1.357, 1.329, 1.302, 1.274, 1.247, 1.219, 1.192,
00086 1.165, 1.138, 1.112, 1.086, 1.060, 1.035, 1.010, 0.985, 0.961, 0.938,
00087 0.930, 0.652, 0.479, 0.365, 0.287, 0.232, 0.191, 0.160, 0.135, 0.116,
00088 0.101, 0.089, 0.078, 0.070, 0.062, 0.056, 0.051, 0.046, 0.042, 0.039,
00089 0.036, 0.033, 0.031, 0.028, 0.026, 0.025, 0.023, 0.022, 0.020, 0.019,
00090 0.018, 0.017, 0.016, 0.015, 0.015, 0.014, 0.013, 0.013, 0.012, 0.011
00091 };
00092
00093 const G4double G4XNNstarTable::sigmaNN1535[121] =
00094 {
00095 0.0,
00096 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003,
00097 0.004, 0.005, 0.006, 0.008, 0.010, 0.012, 0.015, 0.019,
00098 0.024, 0.031, 0.039, 0.052, 0.069, 0.097, 0.145, 0.216,
00099 0.298, 0.378, 0.451, 0.513, 0.566, 0.610, 0.646, 0.675,
00100 0.699, 0.718, 0.732, 0.742, 0.749, 0.753, 0.754, 0.753,
00101 0.751, 0.746, 0.740, 0.733, 0.724, 0.715, 0.705, 0.694,
00102 0.683, 0.671, 0.659, 0.647, 0.634, 0.621, 0.608, 0.595,
00103 0.582, 0.569, 0.556, 0.543, 0.531, 0.518, 0.506, 0.493,
00104 0.481, 0.469, 0.458, 0.446, 0.435, 0.424, 0.413, 0.402,
00105 0.399, 0.276, 0.201, 0.153, 0.120, 0.096, 0.079, 0.066,
00106 0.056, 0.048, 0.042, 0.037, 0.032, 0.029, 0.026, 0.023,
00107 0.021, 0.019, 0.017, 0.016, 0.015, 0.014, 0.013, 0.012,
00108 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007, 0.007,
00109 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005
00110 };
00111
00112 const G4double G4XNNstarTable::sigmaNN1650[121] =
00113 {
00114 0.0,
00115 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00116 0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.001,
00117 0.002, 0.002, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007,
00118 0.008, 0.010, 0.012, 0.015, 0.019, 0.023, 0.029, 0.038,
00119 0.051, 0.071, 0.104, 0.150, 0.201, 0.249, 0.290, 0.326,
00120 0.354, 0.378, 0.397, 0.412, 0.424, 0.434, 0.440, 0.445,
00121 0.448, 0.449, 0.449, 0.448, 0.445, 0.442, 0.438, 0.433,
00122 0.428, 0.422, 0.416, 0.409, 0.403, 0.395, 0.388, 0.381,
00123 0.373, 0.366, 0.358, 0.350, 0.343, 0.335, 0.327, 0.320,
00124 0.312, 0.305, 0.298, 0.291, 0.284, 0.277, 0.270, 0.263,
00125 0.261, 0.182, 0.133, 0.102, 0.080, 0.064, 0.053, 0.044,
00126 0.037, 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015,
00127 0.014, 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008,
00128 0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
00129 0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
00130 };
00131
00132 const G4double G4XNNstarTable::sigmaNN1675[121] =
00133 {
00134 0.0,
00135 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00136 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00137 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
00138 0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.009, 0.015,
00139 0.026, 0.048, 0.095, 0.189, 0.324, 0.463, 0.586, 0.691,
00140 0.780, 0.855, 0.919, 0.972, 1.016, 1.052, 1.081, 1.105,
00141 1.123, 1.136, 1.145, 1.151, 1.153, 1.153, 1.149, 1.144,
00142 1.136, 1.127, 1.116, 1.104, 1.090, 1.076, 1.061, 1.045,
00143 1.028, 1.011, 0.993, 0.975, 0.957, 0.939, 0.921, 0.902,
00144 0.884, 0.865, 0.847, 0.828, 0.810, 0.792, 0.775, 0.757,
00145 0.751, 0.538, 0.399, 0.307, 0.242, 0.196, 0.162, 0.136,
00146 0.115, 0.099, 0.086, 0.076, 0.067, 0.060, 0.053, 0.048,
00147 0.044, 0.040, 0.036, 0.033, 0.031, 0.028, 0.026, 0.024,
00148 0.023, 0.021, 0.020, 0.019, 0.018, 0.016, 0.016, 0.015,
00149 0.014, 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010
00150 };
00151
00152 const G4double G4XNNstarTable::sigmaNN1680[121] =
00153 {
00154 0.0,
00155 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00156 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00157 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00158 0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.006, 0.010,
00159 0.018, 0.035, 0.073, 0.156, 0.294, 0.446, 0.580, 0.693,
00160 0.788, 0.867, 0.933, 0.988, 1.033, 1.070, 1.100, 1.124,
00161 1.142, 1.155, 1.163, 1.168, 1.170, 1.168, 1.164, 1.158,
00162 1.149, 1.139, 1.127, 1.114, 1.100, 1.085, 1.068, 1.052,
00163 1.034, 1.016, 0.998, 0.979, 0.960, 0.941, 0.922, 0.903,
00164 0.884, 0.865, 0.846, 0.827, 0.809, 0.790, 0.772, 0.754,
00165 0.748, 0.533, 0.394, 0.301, 0.238, 0.192, 0.158, 0.133,
00166 0.113, 0.097, 0.084, 0.074, 0.065, 0.058, 0.052, 0.047,
00167 0.042, 0.039, 0.035, 0.032, 0.030, 0.028, 0.026, 0.024,
00168 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
00169 0.013, 0.013, 0.012, 0.012, 0.011, 0.010, 0.010, 0.010
00170 };
00171
00172 const G4double G4XNNstarTable::sigmaNN1700[121] =
00173 {
00174 0.0,
00175 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00176 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00177 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00178 0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.005, 0.008,
00179 0.013, 0.022, 0.038, 0.070, 0.140, 0.245, 0.347, 0.431,
00180 0.500, 0.556, 0.601, 0.637, 0.666, 0.689, 0.708, 0.721,
00181 0.731, 0.738, 0.742, 0.743, 0.743, 0.740, 0.736, 0.730,
00182 0.723, 0.716, 0.707, 0.697, 0.687, 0.676, 0.665, 0.654,
00183 0.642, 0.630, 0.617, 0.605, 0.593, 0.580, 0.568, 0.555,
00184 0.543, 0.531, 0.519, 0.507, 0.495, 0.483, 0.471, 0.460,
00185 0.456, 0.322, 0.236, 0.180, 0.142, 0.114, 0.094, 0.079,
00186 0.067, 0.057, 0.050, 0.044, 0.039, 0.034, 0.031, 0.028,
00187 0.025, 0.023, 0.021, 0.019, 0.018, 0.016, 0.015, 0.014,
00188 0.013, 0.012, 0.011, 0.011, 0.010, 0.009, 0.009, 0.008,
00189 0.008, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006
00190 };
00191
00192 const G4double G4XNNstarTable::sigmaNN1710[121] =
00193 {
00194 0.0,
00195 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00196 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00197 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00198 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.003,
00199 0.005, 0.008, 0.013, 0.025, 0.052, 0.096, 0.144, 0.185,
00200 0.219, 0.246, 0.269, 0.288, 0.304, 0.316, 0.326, 0.334,
00201 0.340, 0.344, 0.347, 0.349, 0.349, 0.349, 0.348, 0.346,
00202 0.343, 0.340, 0.336, 0.332, 0.328, 0.323, 0.319, 0.313,
00203 0.308, 0.303, 0.297, 0.292, 0.286, 0.280, 0.275, 0.269,
00204 0.263, 0.257, 0.252, 0.246, 0.241, 0.235, 0.230, 0.224,
00205 0.223, 0.158, 0.117, 0.090, 0.071, 0.057, 0.047, 0.040,
00206 0.034, 0.029, 0.025, 0.022, 0.019, 0.017, 0.016, 0.014,
00207 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
00208 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
00209 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
00210 };
00211
00212 const G4double G4XNNstarTable::sigmaNN1720[121] =
00213 {
00214 0.0,
00215 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00216 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00217 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00218 0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.007, 0.010,
00219 0.015, 0.023, 0.036, 0.061, 0.106, 0.174, 0.250, 0.321,
00220 0.382, 0.434, 0.477, 0.513, 0.543, 0.568, 0.589, 0.605,
00221 0.617, 0.627, 0.634, 0.639, 0.641, 0.642, 0.641, 0.639,
00222 0.635, 0.630, 0.625, 0.618, 0.611, 0.603, 0.595, 0.586,
00223 0.576, 0.567, 0.557, 0.547, 0.537, 0.527, 0.516, 0.506,
00224 0.496, 0.485, 0.475, 0.465, 0.455, 0.444, 0.435, 0.425,
00225 0.421, 0.302, 0.224, 0.172, 0.136, 0.117, 0.091, 0.076,
00226 0.065, 0.056, 0.049, 0.043, 0.038, 0.034, 0.030, 0.027,
00227 0.025, 0.022, 0.020, 0.019, 0.017, 0.016, 0.015, 0.014,
00228 0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008,
00229 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006, 0.006
00230 };
00231
00232 const G4double G4XNNstarTable::sigmaNN1900[121] =
00233 {
00234 0.0,
00235 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00236 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00237 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
00238 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.003, 0.005,
00239 0.006, 0.008, 0.010, 0.014, 0.019, 0.026, 0.037, 0.054,
00240 0.074, 0.094, 0.114, 0.131, 0.147, 0.161, 0.173, 0.184,
00241 0.193, 0.201, 0.208, 0.213, 0.218, 0.221, 0.224, 0.226,
00242 0.228, 0.229, 0.229, 0.229, 0.229, 0.228, 0.227, 0.225,
00243 0.223, 0.221, 0.219, 0.217, 0.214, 0.212, 0.209, 0.206,
00244 0.203, 0.200, 0.197, 0.194, 0.190, 0.187, 0.184, 0.181,
00245 0.180, 0.137, 0.106, 0.083, 0.067, 0.056, 0.047, 0.039,
00246 0.034, 0.029, 0.026, 0.023, 0.020, 0.018, 0.016, 0.015,
00247 0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.008, 0.008,
00248 0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
00249 0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
00250 };
00251
00252 const G4double G4XNNstarTable::sigmaNN1990[121] =
00253 {
00254 0.0,
00255 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00256 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00257 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00258 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00259 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009,
00260 0.017, 0.030, 0.051, 0.076, 0.102, 0.127, 0.149, 0.168,
00261 0.185, 0.199, 0.211, 0.221, 0.230, 0.237, 0.243, 0.247,
00262 0.251, 0.254, 0.256, 0.258, 0.258, 0.259, 0.258, 0.258,
00263 0.257, 0.255, 0.254, 0.252, 0.249, 0.247, 0.245, 0.242,
00264 0.239, 0.236, 0.233, 0.230, 0.226, 0.223, 0.220, 0.216,
00265 0.215, 0.167, 0.131, 0.104, 0.085, 0.070, 0.059, 0.050,
00266 0.043, 0.038, 0.033, 0.029, 0.026, 0.023, 0.021, 0.019,
00267 0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.010, 0.010,
00268 0.009, 0.009, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006,
00269 0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004
00270 };
00271
00272 const G4double G4XNNstarTable::sigmaNN2090[121] =
00273 {
00274 0.0,
00275 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00276 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00277 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00278 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00279 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.005,
00280 0.007, 0.011, 0.016, 0.024, 0.036, 0.053, 0.071, 0.089,
00281 0.106, 0.120, 0.133, 0.142, 0.151, 0.158, 0.164, 0.169,
00282 0.172, 0.175, 0.178, 0.179, 0.180, 0.180, 0.180, 0.180,
00283 0.179, 0.178, 0.176, 0.175, 0.173, 0.171, 0.169, 0.166,
00284 0.164, 0.162, 0.159, 0.156, 0.154, 0.151, 0.148, 0.145,
00285 0.144, 0.107, 0.081, 0.063, 0.050, 0.041, 0.034, 0.028,
00286 0.024, 0.021, 0.018, 0.016, 0.014, 0.013, 0.011, 0.010,
00287 0.009, 0.008, 0.008, 0.007, 0.006, 0.006, 0.006, 0.005,
00288 0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003,
00289 0.003, 0.003, 0.003, 0.003, 0.002, 0.002, 0.002, 0.002
00290 };
00291
00292 const G4double G4XNNstarTable::sigmaNN2190[121] =
00293 {
00294 0.0,
00295 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00296 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00297 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00298 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
00299 0.001, 0.002, 0.003, 0.005, 0.009, 0.015, 0.024, 0.036,
00300 0.050, 0.064, 0.078, 0.090, 0.100, 0.110, 0.118, 0.125,
00301 0.131, 0.136, 0.140, 0.143, 0.146, 0.149, 0.150, 0.152,
00302 0.153, 0.153, 0.153, 0.153, 0.153, 0.152, 0.152, 0.151,
00303 0.150, 0.148, 0.147, 0.145, 0.144, 0.142, 0.140, 0.139,
00304 0.138, 0.110, 0.087, 0.069, 0.056, 0.047, 0.039, 0.033,
00305 0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.014, 0.013,
00306 0.012, 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006,
00307 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004,
00308 0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003
00309 };
00310
00311 const G4double G4XNNstarTable::sigmaNN2220[121] =
00312 {
00313 0.0,
00314 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00315 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00316 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00317 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00318 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00319 0.001, 0.001, 0.001, 0.003, 0.005, 0.008, 0.015, 0.026,
00320 0.041, 0.058, 0.074, 0.089, 0.101, 0.113, 0.122, 0.131,
00321 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161, 0.163,
00322 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164, 0.163,
00323 0.162, 0.161, 0.160, 0.158, 0.157, 0.155, 0.153, 0.151,
00324 0.150, 0.121, 0.096, 0.077, 0.062, 0.052, 0.043, 0.037,
00325 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
00326 0.013, 0.011, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
00327 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
00328 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003
00329 };
00330
00331 const G4double G4XNNstarTable::sigmaNN2250[121] =
00332 {
00333 0.0,
00334 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00335 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00336 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00337 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00338 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
00339 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009, 0.016,
00340 0.027, 0.043, 0.060, 0.076, 0.090, 0.103, 0.114, 0.123,
00341 0.131, 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161,
00342 0.163, 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164,
00343 0.163, 0.162, 0.161, 0.159, 0.158, 0.156, 0.154, 0.152,
00344 0.152, 0.122, 0.096, 0.077, 0.062, 0.051, 0.043, 0.037,
00345 0.031, 0.027, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
00346 0.012, 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007,
00347 0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004,
00348 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
00349 };
00350
00351
00352 G4XNNstarTable::G4XNNstarTable()
00353 {
00354 xMap["N(1440)0"] = (G4double*) sigmaNN1440;
00355 xMap["N(1440)+"] = (G4double*) sigmaNN1440;
00356
00357 xMap["N(1520)0"] = (G4double*) sigmaNN1520;
00358 xMap["N(1520)+"] = (G4double*) sigmaNN1520;
00359
00360 xMap["N(1535)0"] = (G4double*) sigmaNN1535;
00361 xMap["N(1535)+"] = (G4double*) sigmaNN1535;
00362
00363 xMap["N(1650)0"] = (G4double*) sigmaNN1650;
00364 xMap["N(1650)+"] = (G4double*) sigmaNN1650;
00365
00366 xMap["N(1675)0"] = (G4double*) sigmaNN1675;
00367 xMap["N(1675)+"] = (G4double*) sigmaNN1675;
00368
00369 xMap["N(1680)0"] = (G4double*) sigmaNN1680;
00370 xMap["N(1680)+"] = (G4double*) sigmaNN1680;
00371
00372 xMap["N(1700)0"] = (G4double*) sigmaNN1700;
00373 xMap["N(1700)+"] = (G4double*) sigmaNN1700;
00374
00375 xMap["N(1710)0"] = (G4double*) sigmaNN1710;
00376 xMap["N(1710)+"] = (G4double*) sigmaNN1710;
00377
00378 xMap["N(1720)0"] = (G4double*) sigmaNN1720;
00379 xMap["N(1720)+"] = (G4double*) sigmaNN1720;
00380
00381 xMap["N(1900)0"] = (G4double*) sigmaNN1900;
00382 xMap["N(1900)+"] = (G4double*) sigmaNN1900;
00383
00384 xMap["N(1990)0"] = (G4double*) sigmaNN1990;
00385 xMap["N(1990)+"] = (G4double*) sigmaNN1990;
00386
00387 xMap["N(2090)0"] = (G4double*) sigmaNN2090;
00388 xMap["N(2090)+"] = (G4double*) sigmaNN2090;
00389
00390 xMap["N(2190)0"] = (G4double*) sigmaNN2190;
00391 xMap["N(2190)+"] = (G4double*) sigmaNN2190;
00392
00393 xMap["N(2220)0"] = (G4double*) sigmaNN2220;
00394 xMap["N(2220)+"] = (G4double*) sigmaNN2220;
00395
00396 xMap["N(2250)0"] = (G4double*) sigmaNN2250;
00397 xMap["N(2250)+"] = (G4double*) sigmaNN2250;
00398 }
00399
00400
00401 G4XNNstarTable::~G4XNNstarTable()
00402 { }
00403
00404
00405 const G4PhysicsVector* G4XNNstarTable::CrossSectionTable(const G4String& particleName) const
00406 {
00407
00408
00409 if (xMap.find(particleName) != xMap.end())
00410 {
00411
00412 G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeNNstar);
00413 std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
00414 G4double* sigmaPointer = 0;
00415 for (iter = xMap.begin(); iter != xMap.end(); ++iter)
00416 {
00417 G4String str = (*iter).first;
00418 if (str == particleName)
00419 {
00420 sigmaPointer = (*iter).second;
00421 }
00422 }
00423 G4int i;
00424 for (i=0; i<sizeNNstar; i++)
00425 {
00426 G4double value = *(sigmaPointer + i) * millibarn;
00427 G4double energy = energyTable[i] * GeV;
00428 sigmaVector->PutValue(i,energy,value);
00429 }
00430 return sigmaVector;
00431 }
00432 else
00433
00434 return 0;
00435 }
00436
00437
00438