source: trunk/source/processes/hadronic/models/im_r_matrix/src/G4XDeltaNstarTable.cc @ 819

Last change on this file since 819 was 819, checked in by garnier, 16 years ago

import all except CVS
File size: 21.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//      GEANT4 Class file
29//
30//
31//      File name:     G4XDeltaNstarTable
32//
33//      Author:      Maria Grazia Pia (MariaGrazia.Pia@genova.infn.it) 
34//
35//      Creation date: 4 June 1999
36//
37//      Modifications:
38//
39// Hadron Kinetic Model
40// p p -> Delta N* cross section tables
41//
42// -------------------------------------------------------------------
43
44#include "globals.hh"
45#include "G4ios.hh"
46#include "G4XDeltaNstarTable.hh"
47#include "G4PhysicsFreeVector.hh"
48
49
50const G4int G4XDeltaNstarTable::sizeDeltaNstar = 121;
51
52
53// Energies (GeV) corresponding to the cross section table
54// Units are assigned while filling the PhysicsVector
55
56const G4double G4XDeltaNstarTable::energyTable[121] =
57{
58  0.0, 
59  2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.031, 2.037, 2.044, 2.052, 
60  2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.135, 2.151, 2.168, 2.185, 
61  2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.335, 2.360, 2.386, 2.413, 
62  2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.629, 2.664, 2.699, 2.736,
63  2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.018, 3.062, 3.107, 3.153, 
64  3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.502, 3.555, 3.610, 3.666, 
65  3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.081, 4.144, 4.208, 4.273, 
66  4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.754, 4.827, 4.900, 4.975,
67  5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 11.807, 12.941, 14.076, 15.210, 
68 16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555, 
69 27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899, 
70 39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
71};
72
73// Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998
74// Units are assigned while filling the PhysicsVector
75
76const G4double G4XDeltaNstarTable::sigmaDN1440[121] = 
77{
78  0.0, 
79  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,
81  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
82  0.000, 0.000, 0.000, 0.001, 0.003, 0.006, 0.011, 0.021,
83  0.041, 0.000, 0.000, 0.002, 0.011, 0.047, 0.131, 0.257,
84  0.408, 0.568, 0.729, 0.886,  1.036,  1.178,  1.309,  1.431,
85  1.542,  1.644,  1.736,  1.819,  1.894,  1.960,  2.017,  2.068,
86  2.111,  2.148,  2.178,  2.203,  2.222,  2.237,  2.247,  2.253,
87  2.254,  2.252,  2.247,  2.239,  2.228,  2.215,  2.199,  2.181,
88  2.161,  2.140,  2.117,  2.093,  2.068,  2.042,  2.014,  1.986,
89  1.977,  1.558,  1.224, 0.976, 0.797, 0.654, 0.548, 0.465,
90  0.399, 0.346, 0.303, 0.268, 0.238, 0.213,  0.192,  0.173,
91  0.158,  0.144,  0.132,  0.121,  0.112,  0.103,  0.096,  0.089,
92  0.083,  0.078,  0.073,  0.069,  0.065,  0.061,  0.057,  0.054,
93  0.051,  0.049,  0.046,  0.044,  0.042,  0.040,  0.038,  0.036 
94};
95
96const G4double G4XDeltaNstarTable::sigmaDN1520[121] = 
97{
98  0.0,
99  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
100  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
101  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
102  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,
103  0.003,  0.006,  0.013,  0.027,  0.000,  0.001,  0.010,  0.048,
104  0.152,  0.318,  0.513,  0.713,  0.908,  1.091,  1.262,  1.417,
105  1.559,  1.686,  1.801,  1.903,  1.992,  2.071,  2.140,  2.198,
106  2.248,  2.289,  2.323,  2.349,  2.369,  2.383,  2.392,  2.395,
107  2.394,  2.389,  2.380,  2.368,  2.352,  2.334,  2.313,  2.290,
108  2.266,  2.239,  2.211,  2.182,  2.151,  2.119,  2.087,  2.054,
109  2.043,  1.571,  1.211,  0.953,  0.769,  0.627,  0.521,  0.440,
110  0.376,  0.325,  0.284,  0.250,  0.222,  0.198,  0.178,  0.160,
111  0.146,  0.133,  0.121,  0.112,  0.103,  0.095,  0.088,  0.082,
112  0.076,  0.071,  0.067,  0.063,  0.059,  0.056,  0.052,  0.050,
113  0.047,  0.045,  0.042,  0.040,  0.038,  0.036,  0.035,  0.033 
114};
115
116const G4double G4XDeltaNstarTable::sigmaDN1535[121] = 
117{
118  0.0,
119  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
120  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
121  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
122  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.006,  0.010,
123  0.014,  0.021,  0.030,  0.044,  0.000,  0.000,  0.002,  0.010,
124  0.039,  0.097,  0.174,  0.257,  0.339,  0.417,  0.490,  0.556,
125  0.615,  0.667,  0.714,  0.754,  0.790,  0.820,  0.846,  0.867,
126  0.885,  0.899,  0.910,  0.918,  0.924,  0.927,  0.928,  0.927,
127  0.924,  0.920,  0.914,  0.907,  0.899,  0.890,  0.880,  0.870,
128  0.858,  0.847,  0.834,  0.822,  0.809,  0.795,  0.782,  0.768,
129  0.764,  0.576,  0.439,  0.343,  0.275,  0.223,  0.185,  0.156,
130  0.133,  0.115,  0.100,  0.088,  0.078,  0.070,  0.062,  0.056,
131  0.051,  0.047,  0.043,  0.039,  0.036,  0.033,  0.031,  0.029,
132  0.027,  0.025,  0.023,  0.022,  0.021,  0.019,  0.018,  0.017,
133  0.016,  0.016,  0.015,  0.014,  0.013,  0.013,  0.012,  0.012 
134};
135
136const G4double G4XDeltaNstarTable::sigmaDN1650[121] = 
137{
138  0.0,
139  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
140  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
141  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
142  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.003,
143  0.004,  0.006,  0.008,  0.011,  0.014,  0.019,  0.026,  0.000,
144  0.000,  0.001,  0.004,  0.018,  0.047,  0.083,  0.122,  0.159,
145  0.194,  0.227,  0.256,  0.282,  0.305,  0.326,  0.343,  0.359,
146  0.372,  0.383,  0.393,  0.400,  0.406,  0.411,  0.415,  0.417,
147  0.418,  0.418,  0.418,  0.417,  0.415,  0.412,  0.409,  0.405,
148  0.401,  0.397,  0.392,  0.387,  0.382,  0.376,  0.371,  0.365,
149  0.363,  0.278,  0.214,  0.168,  0.135,  0.110,  0.091,  0.077,
150  0.066,  0.057,  0.049,  0.043,  0.038,  0.034,  0.031,  0.028,
151  0.025,  0.023,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,
152  0.013,  0.012,  0.012,  0.011,  0.010,  0.010,  0.009,  0.009,
153  0.008,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006 
154};
155
156const G4double G4XDeltaNstarTable::sigmaDN1675[121] = 
157{
158  0.0,
159  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
160  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
161  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
162  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
163  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.007,  0.013,
164  0.000,  0.000,  0.003,  0.021,  0.072,  0.155,  0.251,  0.350,
165  0.446,  0.538,  0.622,  0.699,  0.771,  0.835,  0.893,  0.945,
166  0.992,  1.033,  1.068,  1.099,  1.126,  1.148,  1.167,  1.182,
167  1.194,  1.203,  1.209,  1.212,  1.214,  1.213,  1.210,  1.205,
168  1.199,  1.191,  1.182,  1.172,  1.161,  1.148,  1.135,  1.122,
169  1.117,  0.891,  0.700,  0.556,  0.454,  0.371,  0.310,  0.262,
170  0.225,  0.195,  0.170,  0.150,  0.133,  0.119,  0.107,  0.097,
171  0.088,  0.080,  0.073,  0.067,  0.062,  0.057,  0.053,  0.050,
172  0.046,  0.043,  0.040,  0.038,  0.036,  0.034,  0.032,  0.030,
173  0.028,  0.027,  0.026,  0.024,  0.023,  0.022,  0.021,  0.020 
174};
175
176const G4double G4XDeltaNstarTable::sigmaDN1680[121] = 
177{
178  0.0,
179  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
180  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
181  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
182  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
183  0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.009,
184  0.000,  0.000,  0.003,  0.019,  0.070,  0.157,  0.258,  0.361,
185  0.461,  0.555,  0.641,  0.720,  0.792,  0.857,  0.915,  0.966,
186  1.012,  1.051,  1.086,  1.116,  1.141,  1.162,  1.179,  1.192,
187  1.202,  1.209,  1.214,  1.216,  1.215,  1.212,  1.208,  1.202,
188  1.194,  1.185,  1.175,  1.163,  1.151,  1.137,  1.123,  1.108,
189  1.103,  0.870,  0.678,  0.536,  0.437,  0.355,  0.296,  0.250,
190  0.214,  0.185,  0.162,  0.142,  0.126,  0.113,  0.101,  0.092,
191  0.083,  0.076,  0.069,  0.064,  0.059,  0.054,  0.050,  0.047,
192  0.044,  0.041,  0.038,  0.036,  0.034,  0.032,  0.030,  0.028,
193  0.027,  0.025,  0.024,  0.023,  0.022,  0.021,  0.020,  0.019
194};
195
196const G4double G4XDeltaNstarTable::sigmaDN1700[121] = 
197{
198  0.0,
199  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
200  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
201  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
202  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
203  0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.007,
204  0.013,  0.000,  0.001,  0.006,  0.028,  0.078,  0.142,  0.209,
205  0.273,  0.332,  0.385,  0.432,  0.475,  0.512,  0.545,  0.573,
206  0.598,  0.620,  0.638,  0.653,  0.665,  0.675,  0.683,  0.688,
207  0.692,  0.694,  0.695,  0.694,  0.692,  0.689,  0.685,  0.680,
208  0.674,  0.668,  0.660,  0.653,  0.645,  0.636,  0.627,  0.618,
209  0.615,  0.477,  0.368,  0.290,  0.234,  0.190,  0.158,  0.134,
210  0.114,  0.099,  0.086,  0.076,  0.067,  0.060,  0.054,  0.049,
211  0.044,  0.040,  0.037,  0.034,  0.031,  0.029,  0.027,  0.025,
212  0.023,  0.022,  0.020,  0.019,  0.018,  0.017,  0.016,  0.015,
213  0.014,  0.013,  0.013,  0.012,  0.012,  0.011,  0.010,  0.010 
214};
215
216const G4double G4XDeltaNstarTable::sigmaDN1710[121] = 
217{
218  0.0,
219  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
220  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
221  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
222  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
223  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,  0.003,
224  0.005,  0.000,  0.000,  0.002,  0.009,  0.028,  0.056,  0.086,
225  0.116,  0.144,  0.170,  0.194,  0.215,  0.234,  0.251,  0.266,
226  0.279,  0.290,  0.300,  0.308,  0.315,  0.321,  0.326,  0.329,
227  0.332,  0.334,  0.335,  0.336,  0.335,  0.334,  0.333,  0.331,
228  0.329,  0.327,  0.324,  0.320,  0.317,  0.313,  0.309,  0.305,
229  0.304,  0.239,  0.187,  0.161,  0.121,  0.098,  0.082,  0.069,
230  0.059,  0.051,  0.045,  0.040,  0.035,  0.031,  0.028,  0.025,
231  0.023,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,  0.013,
232  0.012,  0.011,  0.011,  0.010,  0.009,  0.009,  0.008,  0.008,
233  0.007,  0.007,  0.007,  0.006,  0.006,  0.006,  0.006,  0.005 
234};
235
236const G4double G4XDeltaNstarTable::sigmaDN1720[121] = 
237{
238  0.0,
239  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
240  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
241  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
242  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
243  0.000,  0.000,  0.001,  0.001,  0.002,  0.003,  0.005,  0.008,
244  0.014,  0.000,  0.000,  0.001,  0.009,  0.033,  0.075,  0.124,
245  0.175,  0.225,  0.271,  0.315,  0.355,  0.391,  0.423,  0.452,
246  0.478,  0.501,  0.521,  0.538,  0.553,  0.566,  0.576,  0.585,
247  0.592,  0.597,  0.601,  0.603,  0.604,  0.604,  0.603,  0.601,
248  0.598,  0.595,  0.590,  0.585,  0.580,  0.574,  0.568,  0.561,
249  0.559,  0.446,  0.351,  0.279,  0.228,  0.202,  0.156,  0.132,
250  0.113,  0.098,  0.086,  0.076,  0.067,  0.060,  0.054,  0.049,
251  0.044,  0.040,  0.037,  0.034,  0.031,  0.029,  0.027,  0.025,
252  0.023,  0.022,  0.020,  0.019,  0.018,  0.017,  0.016,  0.015,
253  0.014,  0.014,  0.013,  0.012,  0.012,  0.011,  0.011,  0.010 
254};
255
256const G4double G4XDeltaNstarTable::sigmaDN1900[121] = 
257{
258  0.0,
259  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
260  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
261  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
262  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
263  0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.002,  0.003,
264  0.005,  0.007,  0.009,  0.014,  0.021,  0.000,  0.000,  0.001,
265  0.003,  0.009,  0.017,  0.028,  0.038,  0.049,  0.060,  0.070,
266  0.080,  0.090,  0.099,  0.107,  0.115,  0.122,  0.128,  0.134,
267  0.140,  0.145,  0.149,  0.153,  0.156,  0.159,  0.162,  0.164,
268  0.166,  0.168,  0.169,  0.170,  0.170,  0.171,  0.171,  0.171,
269  0.171,  0.155,  0.131,  0.110,  0.095,  0.079,  0.067,  0.058,
270  0.051,  0.044,  0.039,  0.035,  0.031,  0.028,  0.026,  0.023,
271  0.021,  0.019,  0.018,  0.016,  0.016,  0.015,  0.013,  0.012,
272  0.011,  0.011,  0.010,  0.009,  0.009,  0.008,  0.008,  0.008,
273  0.007,  0.007,  0.006,  0.006,  0.006,  0.006,  0.005,  0.005 
274};
275
276const G4double G4XDeltaNstarTable::sigmaDN1990[121] = 
277{
278  0.0,
279  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
280  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
281  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
282  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
283  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
284  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.008,  0.016,
285  0.000,  0.000,  0.001,  0.005,  0.013,  0.024,  0.035,  0.046,
286  0.057,  0.068,  0.078,  0.088,  0.097,  0.105,  0.113,  0.120,
287  0.126,  0.132,  0.138,  0.143,  0.147,  0.151,  0.155,  0.158,
288  0.161,  0.164,  0.166,  0.168,  0.169,  0.171,  0.172,  0.172,
289  0.173,  0.164,  0.143,  0.122,  0.105,  0.089,  0.076,  0.066,
290  0.058,  0.053,  0.045,  0.040,  0.036,  0.032,  0.029,  0.027,
291  0.024,  0.022,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,
292  0.013,  0.012,  0.012,  0.011,  0.010,  0.010,  0.009,  0.009,
293  0.008,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006 
294};
295
296const G4double G4XDeltaNstarTable::sigmaDN2090[121] = 
297{
298  0.0,
299  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
300  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
301  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
302  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
303  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
304  0.000,  0.000,  0.001,  0.001,  0.002,  0.003,  0.004,  0.006,
305  0.010,  0.000,  0.000,  0.001,  0.005,  0.013,  0.022,  0.033,
306  0.043,  0.053,  0.063,  0.071,  0.079,  0.086,  0.093,  0.099,
307  0.104,  0.108,  0.112,  0.116,  0.119,  0.121,  0.123,  0.125,
308  0.126,  0.127,  0.128,  0.129,  0.129,  0.129,  0.129,  0.128,
309  0.128,  0.110,  0.090,  0.073,  0.061,  0.049,  0.041,  0.035,
310  0.030,  0.026,  0.023,  0.020,  0.018,  0.016,  0.015,  0.013,
311  0.012,  0.011,  0.010,  0.009,  0.008,  0.008,  0.007,  0.007,
312  0.006,  0.006,  0.006,  0.005,  0.005,  0.005,  0.004,  0.004,
313  0.004,  0.004,  0.004,  0.003,  0.003,  0.003,  0.003,  0.003 
314};
315
316const G4double G4XDeltaNstarTable::sigmaDN2190[121] = 
317{
318  0.0,
319  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
320  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
321  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
322  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
323  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
324  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,
325  0.002,  0.003,  0.006,  0.010,  0.000,  0.000,  0.000,  0.003,
326  0.007,  0.013,  0.019,  0.026,  0.032,  0.038,  0.044,  0.050,
327  0.055,  0.060,  0.065,  0.069,  0.073,  0.077,  0.080,  0.083,
328  0.086,  0.088,  0.090,  0.092,  0.094,  0.095,  0.096,  0.097,
329  0.098,  0.097,  0.086,  0.073,  0.062,  0.053,  0.046,  0.039,
330  0.034,  0.030,  0.027,  0.024,  0.021,  0.019,  0.017,  0.016,
331  0.015,  0.013,  0.012,  0.011,  0.010,  0.010,  0.009,  0.008,
332  0.008,  0.007,  0.007,  0.006,  0.006,  0.006,  0.005,  0.005,
333  0.005,  0.005,  0.004,  0.004,  0.004,  0.004,  0.004,  0.003 
334};
335
336const G4double G4XDeltaNstarTable::sigmaDN2220[121] = 
337{
338  0.0,
339  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
340  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
341  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
342  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
343  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
344  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
345  0.001,  0.001,  0.003,  0.006,  0.011,  0.000,  0.000,  0.001,
346  0.005,  0.010,  0.017,  0.024,  0.030,  0.037,  0.044,  0.050,
347  0.055,  0.061,  0.066,  0.070,  0.074,  0.078,  0.082,  0.085,
348  0.088,  0.091,  0.093,  0.095,  0.097,  0.099,  0.100,  0.101,
349  0.102,  0.103,  0.092,  0.079,  0.067,  0.058,  0.050,  0.043,
350  0.038,  0.033,  0.029,  0.027,  0.023,  0.021,  0.019,  0.017,
351  0.016,  0.014,  0.013,  0.012,  0.011,  0.011,  0.010,  0.009,
352  0.009,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006,
353  0.005,  0.005,  0.005,  0.005,  0.004,  0.004,  0.004,  0.004, 
354};
355
356const G4double G4XDeltaNstarTable::sigmaDN2250[121] = 
357{
358  0.0,
359  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
360  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
361  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
362  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
363  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
364  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
365  0.001,  0.001,  0.002,  0.003,  0.006,  0.000,  0.000,  0.000,
366  0.003,  0.007,  0.014,  0.021,  0.028,  0.035,  0.042,  0.049,
367  0.055,  0.060,  0.066,  0.071,  0.075,  0.079,  0.083,  0.087,
368  0.090,  0.093,  0.095,  0.098,  0.100,  0.101,  0.103,  0.104,
369  0.105,  0.105,  0.093,  0.079,  0.067,  0.057,  0.049,  0.042,
370  0.037,  0.032,  0.029,  0.025,  0.023,  0.020,  0.018,  0.017,
371  0.015,  0.014,  0.013,  0.012,  0.011,  0.010,  0.009,  0.009,
372  0.008,  0.008,  0.007,  0.007,  0.006,  0.006,  0.006,  0.005,
373  0.005,  0.005,  0.005,  0.004,  0.004,  0.004,  0.004,  0.004 
374};
375
376
377G4XDeltaNstarTable::G4XDeltaNstarTable() 
378{
379  xMap["N(1440)0"] = (G4double*) sigmaDN1440;
380  xMap["N(1440)+"] = (G4double*) sigmaDN1440;
381 
382  xMap["N(1520)0"] = (G4double*) sigmaDN1520;
383  xMap["N(1520)+"] = (G4double*) sigmaDN1520;
384 
385  xMap["N(1535)0"] = (G4double*) sigmaDN1535;
386  xMap["N(1535)+"] = (G4double*) sigmaDN1535;
387 
388  xMap["N(1650)0"] = (G4double*) sigmaDN1650;
389  xMap["N(1650)+"] = (G4double*) sigmaDN1650;
390 
391  xMap["N(1675)0"] = (G4double*) sigmaDN1675;
392  xMap["N(1675)+"] = (G4double*) sigmaDN1675;
393 
394  xMap["N(1680)0"] = (G4double*) sigmaDN1680;
395  xMap["N(1680)+"] = (G4double*) sigmaDN1680;
396 
397  xMap["N(1700)0"] = (G4double*) sigmaDN1700;
398  xMap["N(1700)+"] = (G4double*) sigmaDN1700;
399 
400  xMap["N(1710)0"] = (G4double*) sigmaDN1710;
401  xMap["N(1710)+"] = (G4double*) sigmaDN1710;
402 
403  xMap["N(1720)0"] = (G4double*) sigmaDN1720;
404  xMap["N(1720)+"] = (G4double*) sigmaDN1720;
405 
406  xMap["N(1900)0"] = (G4double*) sigmaDN1900;
407  xMap["N(1900)+"] = (G4double*) sigmaDN1900;
408   
409  xMap["N(1990)0"] = (G4double*) sigmaDN1990;
410  xMap["N(1990)+"] = (G4double*) sigmaDN1990;
411 
412  xMap["N(2090)0"] = (G4double*) sigmaDN2090;
413  xMap["N(2090)+"] = (G4double*) sigmaDN2090;
414 
415  xMap["N(2190)0"] = (G4double*) sigmaDN2190;
416  xMap["N(2190)+"] = (G4double*) sigmaDN2190;
417 
418  xMap["N(2220)0"] = (G4double*) sigmaDN2220;
419  xMap["N(2220)+"] = (G4double*) sigmaDN2220;
420 
421  xMap["N(2250)0"] = (G4double*) sigmaDN2250;
422  xMap["N(2250)+"] = (G4double*) sigmaDN2250; 
423}
424
425
426G4XDeltaNstarTable::~G4XDeltaNstarTable()
427{ }
428
429
430const G4PhysicsVector* G4XDeltaNstarTable::CrossSectionTable(const G4String& particleName) const
431{
432   // NOTE: the returned pointer is owned by the client
433
434  if (xMap.find(particleName) != xMap.end())
435    {
436      // Cross section table for the requested particle available in the Map
437      G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeDeltaNstar);
438      std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
439      G4double* sigmaPointer = 0;
440      for (iter = xMap.begin(); iter != xMap.end(); ++iter)
441        {
442          G4String str = (*iter).first;
443          if (str == particleName)
444            {
445              sigmaPointer = (*iter).second; 
446            }
447        }
448
449      G4int i;
450      for (i=0; i<sizeDeltaNstar; i++)
451        {
452          G4double value = *(sigmaPointer + i) * 0.5* millibarn;
453          G4double energy = energyTable[i] * GeV;
454          sigmaVector->PutValue(i,energy,value);
455        }           
456      return sigmaVector;
457    }
458  else
459    // No cross section table for the requested particle is available in the Map
460    return 0;
461}
462
463
464
465
466
467
468
469
470
471
472
473
Note: See TracBrowser for help on using the repository browser.