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