| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // ------------------------------------------- 27 // -------------------------------------------------------------------
28 // GEANT4 Class file 28 // GEANT4 Class file
29 // 29 //
30 // Hadron Kinetic Model 30 // Hadron Kinetic Model
31 // p p -> Delta Delta* cross section tables 31 // p p -> Delta Delta* cross section tables
32 // 32 //
33 // ------------------------------------------- 33 // -------------------------------------------------------------------
34 34
35 #include "globals.hh" 35 #include "globals.hh"
36 #include "G4ios.hh" 36 #include "G4ios.hh"
37 #include "G4SystemOfUnits.hh" 37 #include "G4SystemOfUnits.hh"
38 #include "G4XDeltaDeltastarTable.hh" 38 #include "G4XDeltaDeltastarTable.hh"
39 #include "G4PhysicsFreeVector.hh" 39 #include "G4PhysicsFreeVector.hh"
40 40
41 const G4int G4XDeltaDeltastarTable::sizeDeltaD 41 const G4int G4XDeltaDeltastarTable::sizeDeltaDeltastar = 121;
42 42
43 // Energies (GeV) corresponding to the cross s 43 // Energies (GeV) corresponding to the cross section table
44 // Units are assigned when filling the Physics 44 // Units are assigned when filling the PhysicsVector
45 45
46 const G4double G4XDeltaDeltastarTable::energyT 46 const G4double G4XDeltaDeltastarTable::energyTable[121] =
47 { 47 {
48 0.0, 48 0.0,
49 2.014, 2.014, 2.016, 2.018, 2.022, 2.02 49 2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.031, 2.037, 2.044, 2.052,
50 2.061, 2.071, 2.082, 2.094, 2.107, 2.12 50 2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.135, 2.151, 2.168, 2.185,
51 2.204, 2.223, 2.244, 2.265, 2.287, 2.31 51 2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.335, 2.360, 2.386, 2.413,
52 2.441, 2.470, 2.500, 2.531, 2.562, 2.59 52 2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.629, 2.664, 2.699, 2.736,
53 2.773, 2.812, 2.851, 2.891, 2.933, 2.97 53 2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.018, 3.062, 3.107, 3.153,
54 3.200, 3.248, 3.297, 3.347, 3.397, 3.44 54 3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.502, 3.555, 3.610, 3.666,
55 3.722, 3.779, 3.838, 3.897, 3.957, 4.01 55 3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.081, 4.144, 4.208, 4.273,
56 4.339, 4.406, 4.473, 4.542, 4.612, 4.68 56 4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.754, 4.827, 4.900, 4.975,
57 5.000, 6.134, 7.269, 8.403, 9.538, 10.67 57 5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 11.807, 12.941, 14.076, 15.210,
58 16.345, 17.479, 18.613, 19.748, 20.882, 22.01 58 16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555,
59 27.689, 28.824, 29.958, 31.092, 32.227, 33.36 59 27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899,
60 39.034, 40.168, 41.303, 42.437, 43.571, 44.70 60 39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
61 }; 61 };
62 62
63 // Cross-sections in mb, from S.A. Bass et al. 63 // Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998
64 // Units are assigned when filling the Physics 64 // Units are assigned when filling the PhysicsVector
65 65
66 const G4double G4XDeltaDeltastarTable::sigmaDD 66 const G4double G4XDeltaDeltastarTable::sigmaDD1600[121]=
67 { 67 {
68 0.0, 68 0.0,
69 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 69 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
70 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 70 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
71 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 71 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
72 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 72 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
73 0.000, 0.000, 0.001, 0.001, 0.002, 0.00 73 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.010,
74 0.018, 0.000, 0.000, 0.003, 0.015, 0.04 74 0.018, 0.000, 0.000, 0.003, 0.015, 0.045, 0.087, 0.134,
75 0.183, 0.231, 0.277, 0.319, 0.359, 0.39 75 0.183, 0.231, 0.277, 0.319, 0.359, 0.395, 0.429, 0.459,
76 0.486, 0.510, 0.532, 0.551, 0.568, 0.58 76 0.486, 0.510, 0.532, 0.551, 0.568, 0.582, 0.595, 0.605,
77 0.614, 0.621, 0.626, 0.630, 0.633, 0.63 77 0.614, 0.621, 0.626, 0.630, 0.633, 0.635, 0.635, 0.634,
78 0.633, 0.630, 0.627, 0.623, 0.619, 0.61 78 0.633, 0.630, 0.627, 0.623, 0.619, 0.614, 0.608, 0.602,
79 0.600, 0.490, 0.391, 0.314, 0.259, 0.21 79 0.600, 0.490, 0.391, 0.314, 0.259, 0.213, 0.178, 0.152,
80 0.130, 0.113, 0.099, 0.088, 0.078, 0.07 80 0.130, 0.113, 0.099, 0.088, 0.078, 0.070, 0.063, 0.057,
81 0.052, 0.047, 0.043, 0.040, 0.037, 0.03 81 0.052, 0.047, 0.043, 0.040, 0.037, 0.034, 0.031, 0.029,
82 0.027, 0.026, 0.024, 0.023, 0.021, 0.02 82 0.027, 0.026, 0.024, 0.023, 0.021, 0.020, 0.019, 0.018,
83 0.017, 0.016, 0.015, 0.014, 0.014, 0.01 83 0.017, 0.016, 0.015, 0.014, 0.014, 0.013, 0.013, 0.012
84 }; 84 };
85 85
86 const G4double G4XDeltaDeltastarTable::sigmaDD 86 const G4double G4XDeltaDeltastarTable::sigmaDD1620[121] =
87 { 87 {
88 0.0, 88 0.0,
89 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 89 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
90 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 90 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
91 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 91 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
92 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 92 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001, 0.001,
93 0.001, 0.002, 0.003, 0.004, 0.005, 0.00 93 0.001, 0.002, 0.003, 0.004, 0.005, 0.007, 0.009, 0.014,
94 0.000, 0.000, 0.001, 0.006, 0.020, 0.04 94 0.000, 0.000, 0.001, 0.006, 0.020, 0.043, 0.070, 0.100,
95 0.128, 0.156, 0.182, 0.207, 0.229, 0.24 95 0.128, 0.156, 0.182, 0.207, 0.229, 0.249, 0.267, 0.284,
96 0.299, 0.312, 0.324, 0.334, 0.343, 0.35 96 0.299, 0.312, 0.324, 0.334, 0.343, 0.351, 0.357, 0.362,
97 0.367, 0.370, 0.372, 0.374, 0.375, 0.37 97 0.367, 0.370, 0.372, 0.374, 0.375, 0.375, 0.374, 0.373,
98 0.372, 0.370, 0.367, 0.364, 0.361, 0.35 98 0.372, 0.370, 0.367, 0.364, 0.361, 0.358, 0.354, 0.350,
99 0.348, 0.280, 0.220, 0.176, 0.144, 0.11 99 0.348, 0.280, 0.220, 0.176, 0.144, 0.117, 0.098, 0.091,
100 0.071, 0.062, 0.054, 0.048, 0.042, 0.03 100 0.071, 0.062, 0.054, 0.048, 0.042, 0.038, 0.034, 0.031,
101 0.028, 0.025, 0.023, 0.021, 0.020, 0.01 101 0.028, 0.025, 0.023, 0.021, 0.020, 0.018, 0.017, 0.016,
102 0.015, 0.014, 0.013, 0.012, 0.011, 0.01 102 0.015, 0.014, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010,
103 0.009, 0.009, 0.008, 0.008, 0.007, 0.00 103 0.009, 0.009, 0.008, 0.008, 0.007, 0.007, 0.007, 0.006
104 }; 104 };
105 105
106 const G4double G4XDeltaDeltastarTable::sigmaDD 106 const G4double G4XDeltaDeltastarTable::sigmaDD1700[121] =
107 { 107 {
108 0.0, 108 0.0,
109 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 109 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
110 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 110 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
111 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 111 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
112 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 112 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
113 0.000, 0.000, 0.001, 0.001, 0.002, 0.00 113 0.000, 0.000, 0.001, 0.001, 0.002, 0.004, 0.007, 0.012,
114 0.019, 0.031, 0.000, 0.000, 0.001, 0.00 114 0.019, 0.031, 0.000, 0.000, 0.001, 0.006, 0.021, 0.044,
115 0.073, 0.103, 0.134, 0.164, 0.193, 0.22 115 0.073, 0.103, 0.134, 0.164, 0.193, 0.220, 0.246, 0.269,
116 0.291, 0.310, 0.328, 0.344, 0.358, 0.37 116 0.291, 0.310, 0.328, 0.344, 0.358, 0.370, 0.381, 0.391,
117 0.399, 0.405, 0.411, 0.416, 0.419, 0.42 117 0.399, 0.405, 0.411, 0.416, 0.419, 0.421, 0.423, 0.424,
118 0.424, 0.423, 0.422, 0.420, 0.418, 0.41 118 0.424, 0.423, 0.422, 0.420, 0.418, 0.415, 0.412, 0.408,
119 0.407, 0.335, 0.268, 0.215, 0.177, 0.14 119 0.407, 0.335, 0.268, 0.215, 0.177, 0.145, 0.121, 0.103,
120 0.088, 0.077, 0.067, 0.059, 0.053, 0.04 120 0.088, 0.077, 0.067, 0.059, 0.053, 0.047, 0.042, 0.038,
121 0.035, 0.032, 0.029, 0.027, 0.025, 0.02 121 0.035, 0.032, 0.029, 0.027, 0.025, 0.023, 0.021, 0.020,
122 0.018, 0.017, 0.016, 0.015, 0.014, 0.01 122 0.018, 0.017, 0.016, 0.015, 0.014, 0.013, 0.013, 0.012,
123 0.011, 0.011, 0.010, 0.010, 0.009, 0.00 123 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008, 0.008
124 }; 124 };
125 125
126 126
127 const G4double G4XDeltaDeltastarTable::sigmaDD 127 const G4double G4XDeltaDeltastarTable::sigmaDD1900[121] =
128 { 128 {
129 0.0, 129 0.0,
130 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 130 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
131 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 131 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
132 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 132 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
133 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 133 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
134 0.001, 0.001, 0.001, 0.001, 0.002, 0.00 134 0.001, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.003,
135 0.004, 0.006, 0.007, 0.010, 0.000, 0.00 135 0.004, 0.006, 0.007, 0.010, 0.000, 0.000, 0.000, 0.001,
136 0.006, 0.014, 0.025, 0.036, 0.047, 0.05 136 0.006, 0.014, 0.025, 0.036, 0.047, 0.058, 0.068, 0.077,
137 0.086, 0.094, 0.101, 0.107, 0.113, 0.11 137 0.086, 0.094, 0.101, 0.107, 0.113, 0.118, 0.122, 0.126,
138 0.129, 0.132, 0.134, 0.136, 0.137, 0.13 138 0.129, 0.132, 0.134, 0.136, 0.137, 0.138, 0.139, 0.140,
139 0.140, 0.140, 0.139, 0.139, 0.138, 0.13 139 0.140, 0.140, 0.139, 0.139, 0.138, 0.137, 0.136, 0.135,
140 0.135, 0.111, 0.088, 0.070, 0.059, 0.04 140 0.135, 0.111, 0.088, 0.070, 0.059, 0.047, 0.039, 0.033,
141 0.029, 0.025, 0.022, 0.019, 0.017, 0.01 141 0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.014, 0.012,
142 0.011, 0.010, 0.009, 0.009, 0.008, 0.00 142 0.011, 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006,
143 0.006, 0.005, 0.005, 0.005, 0.005, 0.00 143 0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004,
144 0.004, 0.003, 0.003, 0.003, 0.003, 0.00 144 0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003
145 }; 145 };
146 146
147 const G4double G4XDeltaDeltastarTable::sigmaDD 147 const G4double G4XDeltaDeltastarTable::sigmaDD1905[121] =
148 { 148 {
149 0.0, 149 0.0,
150 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 150 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
151 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 151 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
152 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 152 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
153 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 153 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
154 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 154 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
155 0.001, 0.001, 0.003, 0.005, 0.009, 0.00 155 0.001, 0.001, 0.003, 0.005, 0.009, 0.000, 0.000, 0.001,
156 0.005, 0.018, 0.038, 0.061, 0.086, 0.11 156 0.005, 0.018, 0.038, 0.061, 0.086, 0.110, 0.134, 0.157,
157 0.178, 0.197, 0.216, 0.232, 0.248, 0.26 157 0.178, 0.197, 0.216, 0.232, 0.248, 0.262, 0.274, 0.285,
158 0.295, 0.304, 0.312, 0.318, 0.324, 0.32 158 0.295, 0.304, 0.312, 0.318, 0.324, 0.329, 0.333, 0.336,
159 0.338, 0.340, 0.341, 0.342, 0.342, 0.34 159 0.338, 0.340, 0.341, 0.342, 0.342, 0.341, 0.340, 0.339,
160 0.338, 0.293, 0.240, 0.196, 0.166, 0.13 160 0.338, 0.293, 0.240, 0.196, 0.166, 0.135, 0.114, 0.098,
161 0.084, 0.073, 0.064, 0.057, 0.055, 0.04 161 0.084, 0.073, 0.064, 0.057, 0.055, 0.045, 0.041, 0.037,
162 0.034, 0.031, 0.028, 0.026, 0.024, 0.02 162 0.034, 0.031, 0.028, 0.026, 0.024, 0.022, 0.021, 0.019,
163 0.018, 0.017, 0.016, 0.015, 0.014, 0.01 163 0.018, 0.017, 0.016, 0.015, 0.014, 0.013, 0.012, 0.012,
164 0.011, 0.011, 0.010, 0.010, 0.009, 0.00 164 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008, 0.008
165 }; 165 };
166 166
167 const G4double G4XDeltaDeltastarTable::sigmaDD 167 const G4double G4XDeltaDeltastarTable::sigmaDD1910[121] =
168 { 168 {
169 0.0, 169 0.0,
170 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 170 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
171 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 171 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
172 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 172 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
173 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 173 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
174 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 174 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
175 0.001, 0.001, 0.002, 0.003, 0.005, 0.00 175 0.001, 0.001, 0.002, 0.003, 0.005, 0.008, 0.000, 0.000,
176 0.001, 0.004, 0.010, 0.018, 0.027, 0.03 176 0.001, 0.004, 0.010, 0.018, 0.027, 0.035, 0.044, 0.052,
177 0.059, 0.066, 0.073, 0.078, 0.083, 0.08 177 0.059, 0.066, 0.073, 0.078, 0.083, 0.088, 0.092, 0.096,
178 0.099, 0.102, 0.104, 0.106, 0.108, 0.11 178 0.099, 0.102, 0.104, 0.106, 0.108, 0.110, 0.111, 0.112,
179 0.112, 0.113, 0.113, 0.113, 0.112, 0.11 179 0.112, 0.113, 0.113, 0.113, 0.112, 0.112, 0.112, 0.111,
180 0.111, 0.094, 0.076, 0.062, 0.052, 0.04 180 0.111, 0.094, 0.076, 0.062, 0.052, 0.042, 0.036, 0.030,
181 0.026, 0.023, 0.020, 0.018, 0.016, 0.01 181 0.026, 0.023, 0.020, 0.018, 0.016, 0.014, 0.013, 0.011,
182 0.010, 0.009, 0.009, 0.008, 0.007, 0.00 182 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006, 0.006,
183 0.006, 0.005, 0.005, 0.005, 0.004, 0.00 183 0.006, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004, 0.004,
184 0.003, 0.003, 0.003, 0.003, 0.003, 0.00 184 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.002
185 }; 185 };
186 186
187 const G4double G4XDeltaDeltastarTable::sigmaDD 187 const G4double G4XDeltaDeltastarTable::sigmaDD1920[121] =
188 { 188 {
189 0.0, 189 0.0,
190 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 190 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
191 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 191 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
192 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 192 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
193 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 193 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
194 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 194 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
195 0.001, 0.001, 0.001, 0.002, 0.003, 0.00 195 0.001, 0.001, 0.001, 0.002, 0.003, 0.006, 0.000, 0.000,
196 0.001, 0.007, 0.020, 0.040, 0.060, 0.08 196 0.001, 0.007, 0.020, 0.040, 0.060, 0.080, 0.099, 0.117,
197 0.133, 0.147, 0.160, 0.171, 0.181, 0.19 197 0.133, 0.147, 0.160, 0.171, 0.181, 0.190, 0.198, 0.204,
198 0.210, 0.215, 0.219, 0.222, 0.225, 0.22 198 0.210, 0.215, 0.219, 0.222, 0.225, 0.227, 0.228, 0.229,
199 0.229, 0.229, 0.229, 0.228, 0.227, 0.22 199 0.229, 0.229, 0.229, 0.228, 0.227, 0.226, 0.224, 0.222,
200 0.222, 0.183, 0.146, 0.117, 0.098, 0.07 200 0.222, 0.183, 0.146, 0.117, 0.098, 0.078, 0.066, 0.056,
201 0.048, 0.041, 0.036, 0.032, 0.028, 0.02 201 0.048, 0.041, 0.036, 0.032, 0.028, 0.025, 0.023, 0.021,
202 0.019, 0.017, 0.016, 0.014, 0.013, 0.01 202 0.019, 0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.011,
203 0.010, 0.009, 0.009, 0.008, 0.008, 0.00 203 0.010, 0.009, 0.009, 0.008, 0.008, 0.007, 0.007, 0.006,
204 0.006, 0.006, 0.005, 0.005, 0.005, 0.00 204 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.004
205 }; 205 };
206 206
207 const G4double G4XDeltaDeltastarTable::sigmaDD 207 const G4double G4XDeltaDeltastarTable::sigmaDD1930[121] =
208 { 208 {
209 0.0, 209 0.0,
210 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 210 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
211 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 211 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
212 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 212 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
213 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 213 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
214 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 214 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
215 0.000, 0.001, 0.001, 0.002, 0.004, 0.00 215 0.000, 0.001, 0.001, 0.002, 0.004, 0.007, 0.000, 0.000,
216 0.001, 0.004, 0.015, 0.033, 0.053, 0.07 216 0.001, 0.004, 0.015, 0.033, 0.053, 0.075, 0.096, 0.117,
217 0.136, 0.154, 0.171, 0.187, 0.201, 0.21 217 0.136, 0.154, 0.171, 0.187, 0.201, 0.214, 0.225, 0.236,
218 0.245, 0.253, 0.261, 0.267, 0.272, 0.27 218 0.245, 0.253, 0.261, 0.267, 0.272, 0.277, 0.281, 0.284,
219 0.287, 0.288, 0.290, 0.290, 0.291, 0.29 219 0.287, 0.288, 0.290, 0.290, 0.291, 0.291, 0.290, 0.289,
220 0.289, 0.251, 0.206, 0.168, 0.142, 0.11 220 0.289, 0.251, 0.206, 0.168, 0.142, 0.115, 0.097, 0.083,
221 0.071, 0.062, 0.059, 0.048, 0.043, 0.03 221 0.071, 0.062, 0.059, 0.048, 0.043, 0.038, 0.035, 0.031,
222 0.028, 0.026, 0.024, 0.022, 0.020, 0.01 222 0.028, 0.026, 0.024, 0.022, 0.020, 0.019, 0.017, 0.016,
223 0.015, 0.014, 0.013, 0.012, 0.012, 0.01 223 0.015, 0.014, 0.013, 0.012, 0.012, 0.011, 0.010, 0.010,
224 0.009, 0.009, 0.008, 0.008, 0.008, 0.00 224 0.009, 0.009, 0.008, 0.008, 0.008, 0.007, 0.007, 0.007
225 }; 225 };
226 226
227 const G4double G4XDeltaDeltastarTable::sigmaDD 227 const G4double G4XDeltaDeltastarTable::sigmaDD1950[121] =
228 { 228 {
229 0.0, 229 0.0,
230 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 230 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
231 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 231 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
232 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 232 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
233 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 233 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
234 0.000, 0.000, 0.000, 0.000, 0.000, 0.00 234 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
235 0.000, 0.001, 0.001, 0.002, 0.003, 0.00 235 0.000, 0.001, 0.001, 0.002, 0.003, 0.006, 0.011, 0.000,
236 0.000, 0.002, 0.011, 0.029, 0.053, 0.07 236 0.000, 0.002, 0.011, 0.029, 0.053, 0.078, 0.103, 0.128,
237 0.151, 0.173, 0.194, 0.213, 0.230, 0.24 237 0.151, 0.173, 0.194, 0.213, 0.230, 0.246, 0.261, 0.274,
238 0.285, 0.296, 0.305, 0.313, 0.320, 0.32 238 0.285, 0.296, 0.305, 0.313, 0.320, 0.326, 0.331, 0.335,
239 0.339, 0.342, 0.344, 0.345, 0.346, 0.34 239 0.339, 0.342, 0.344, 0.345, 0.346, 0.346, 0.345, 0.344,
240 0.344, 0.301, 0.247, 0.202, 0.170, 0.13 240 0.344, 0.301, 0.247, 0.202, 0.170, 0.139, 0.117, 0.099,
241 0.086, 0.075, 0.065, 0.058, 0.051, 0.04 241 0.086, 0.075, 0.065, 0.058, 0.051, 0.046, 0.042, 0.038,
242 0.034, 0.031, 0.029, 0.026, 0.024, 0.02 242 0.034, 0.031, 0.029, 0.026, 0.024, 0.022, 0.021, 0.019,
243 0.018, 0.017, 0.016, 0.015, 0.014, 0.01 243 0.018, 0.017, 0.016, 0.015, 0.014, 0.013, 0.013, 0.012,
244 0.011, 0.011, 0.010, 0.010, 0.009, 0.00 244 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008, 0.008
245 }; 245 };
246 246
247 247
248 248
249 G4XDeltaDeltastarTable::G4XDeltaDeltastarTable 249 G4XDeltaDeltastarTable::G4XDeltaDeltastarTable()
250 { 250 {
251 xMap["delta(1600)-"] = (G4double*) sigmaDD16 251 xMap["delta(1600)-"] = (G4double*) sigmaDD1600;
252 xMap["delta(1600)0"] = (G4double*) sigmaDD16 252 xMap["delta(1600)0"] = (G4double*) sigmaDD1600;
253 xMap["delta(1600)+"] = (G4double*) sigmaDD16 253 xMap["delta(1600)+"] = (G4double*) sigmaDD1600;
254 xMap["delta(1600)++"] = (G4double*) sigmaDD1 254 xMap["delta(1600)++"] = (G4double*) sigmaDD1600;
255 255
256 xMap["delta(1620)-"] = (G4double*) sigmaDD16 256 xMap["delta(1620)-"] = (G4double*) sigmaDD1620;
257 xMap["delta(1620)0"] = (G4double*) sigmaDD16 257 xMap["delta(1620)0"] = (G4double*) sigmaDD1620;
258 xMap["delta(1620)+"] = (G4double*) sigmaDD16 258 xMap["delta(1620)+"] = (G4double*) sigmaDD1620;
259 xMap["delta(1620)++"] = (G4double*) sigmaDD1 259 xMap["delta(1620)++"] = (G4double*) sigmaDD1620;
260 260
261 xMap["delta(1700)-"] = (G4double*) sigmaDD17 261 xMap["delta(1700)-"] = (G4double*) sigmaDD1700;
262 xMap["delta(1700)0"] = (G4double*) sigmaDD17 262 xMap["delta(1700)0"] = (G4double*) sigmaDD1700;
263 xMap["delta(1700)+"] = (G4double*) sigmaDD17 263 xMap["delta(1700)+"] = (G4double*) sigmaDD1700;
264 xMap["delta(1700)++"] = (G4double*) sigmaDD1 264 xMap["delta(1700)++"] = (G4double*) sigmaDD1700;
265 265
266 xMap["delta(1900)-"] = (G4double*) sigmaDD19 266 xMap["delta(1900)-"] = (G4double*) sigmaDD1900;
267 xMap["delta(1900)0"] = (G4double*) sigmaDD19 267 xMap["delta(1900)0"] = (G4double*) sigmaDD1900;
268 xMap["delta(1900)+"] = (G4double*) sigmaDD19 268 xMap["delta(1900)+"] = (G4double*) sigmaDD1900;
269 xMap["delta(1900)++"] = (G4double*) sigmaDD1 269 xMap["delta(1900)++"] = (G4double*) sigmaDD1900;
270 270
271 xMap["delta(1905)-"] = (G4double*) sigmaDD19 271 xMap["delta(1905)-"] = (G4double*) sigmaDD1905;
272 xMap["delta(1905)0"] = (G4double*) sigmaDD19 272 xMap["delta(1905)0"] = (G4double*) sigmaDD1905;
273 xMap["delta(1905)+"] = (G4double*) sigmaDD19 273 xMap["delta(1905)+"] = (G4double*) sigmaDD1905;
274 xMap["delta(1905)++"] = (G4double*) sigmaDD1 274 xMap["delta(1905)++"] = (G4double*) sigmaDD1905;
275 275
276 xMap["delta(1910)-"] = (G4double*) sigmaDD19 276 xMap["delta(1910)-"] = (G4double*) sigmaDD1910;
277 xMap["delta(1910)0"] = (G4double*) sigmaDD19 277 xMap["delta(1910)0"] = (G4double*) sigmaDD1910;
278 xMap["delta(1910)+"] = (G4double*) sigmaDD19 278 xMap["delta(1910)+"] = (G4double*) sigmaDD1910;
279 xMap["delta(1910)++"] = (G4double*) sigmaDD1 279 xMap["delta(1910)++"] = (G4double*) sigmaDD1910;
280 280
281 xMap["delta(1920)-"] = (G4double*) sigmaDD19 281 xMap["delta(1920)-"] = (G4double*) sigmaDD1920;
282 xMap["delta(1920)0"] = (G4double*) sigmaDD19 282 xMap["delta(1920)0"] = (G4double*) sigmaDD1920;
283 xMap["delta(1920)+"] = (G4double*) sigmaDD19 283 xMap["delta(1920)+"] = (G4double*) sigmaDD1920;
284 xMap["delta(1920)++"] = (G4double*) sigmaDD1 284 xMap["delta(1920)++"] = (G4double*) sigmaDD1920;
285 285
286 xMap["delta(1930)-"] = (G4double*) sigmaDD19 286 xMap["delta(1930)-"] = (G4double*) sigmaDD1930;
287 xMap["delta(1930)0"] = (G4double*) sigmaDD19 287 xMap["delta(1930)0"] = (G4double*) sigmaDD1930;
288 xMap["delta(1930)+"] = (G4double*) sigmaDD19 288 xMap["delta(1930)+"] = (G4double*) sigmaDD1930;
289 xMap["delta(1930)++"] = (G4double*) sigmaDD1 289 xMap["delta(1930)++"] = (G4double*) sigmaDD1930;
290 290
291 xMap["delta(1950)-"] = (G4double*) sigmaDD19 291 xMap["delta(1950)-"] = (G4double*) sigmaDD1950;
292 xMap["delta(1950)0"] = (G4double*) sigmaDD19 292 xMap["delta(1950)0"] = (G4double*) sigmaDD1950;
293 xMap["delta(1950)+"] = (G4double*) sigmaDD19 293 xMap["delta(1950)+"] = (G4double*) sigmaDD1950;
294 xMap["delta(1950)++"] = (G4double*) sigmaDD1 294 xMap["delta(1950)++"] = (G4double*) sigmaDD1950;
295 } 295 }
296 296
297 297
298 G4XDeltaDeltastarTable::~G4XDeltaDeltastarTabl 298 G4XDeltaDeltastarTable::~G4XDeltaDeltastarTable()
299 { } 299 { }
300 300
301 301
302 const G4PhysicsVector* G4XDeltaDeltastarTable: 302 const G4PhysicsVector* G4XDeltaDeltastarTable::CrossSectionTable(const G4String& particleName) const
303 { 303 {
304 // NOTE: the returned pointer is owned by th 304 // NOTE: the returned pointer is owned by the client
305 305
306 if (xMap.find(particleName) != xMap.end()) 306 if (xMap.find(particleName) != xMap.end())
307 { 307 {
308 // Cross section table for the requested 308 // Cross section table for the requested particle available in the Map
309 G4PhysicsFreeVector* sigmaVector = new G 309 G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeDeltaDeltastar);
310 std::map <G4String, G4double*, std::less 310 std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
311 G4double* sigmaPointer = 0; 311 G4double* sigmaPointer = 0;
312 for (iter = xMap.begin(); iter != xMap.e 312 for (iter = xMap.begin(); iter != xMap.end(); ++iter)
313 { 313 {
314 G4String str = (*iter).first; 314 G4String str = (*iter).first;
315 if (str == particleName) 315 if (str == particleName)
316 { 316 {
317 sigmaPointer = (*iter).second; 317 sigmaPointer = (*iter).second;
318 } 318 }
319 } 319 }
320 320
321 G4int i; 321 G4int i;
322 for (i=0; i<sizeDeltaDeltastar; i++) 322 for (i=0; i<sizeDeltaDeltastar; i++)
323 { 323 {
324 G4double value = *(sigmaPointer + i) * 0.5 324 G4double value = *(sigmaPointer + i) * 0.5* millibarn;
325 G4double energy = energyTable[i] * GeV; 325 G4double energy = energyTable[i] * GeV;
326 sigmaVector->PutValue(i,energy,value); 326 sigmaVector->PutValue(i,energy,value);
327 } 327 }
328 return sigmaVector; 328 return sigmaVector;
329 } 329 }
330 else 330 else
331 // No cross section table for the requeste 331 // No cross section table for the requested particle is available in the Map
332 return 0; 332 return 0;
333 } 333 }
334 334
335 335
336 336
337 337