| 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 |
|
|---|
| 50 | const G4int G4XDeltaNstarTable::sizeDeltaNstar = 121;
|
|---|
| 51 |
|
|---|
| 52 |
|
|---|
| 53 | // Energies (GeV) corresponding to the cross section table
|
|---|
| 54 | // Units are assigned while filling the PhysicsVector
|
|---|
| 55 |
|
|---|
| 56 | const 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 |
|
|---|
| 76 | const 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 |
|
|---|
| 96 | const 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 |
|
|---|
| 116 | const 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 |
|
|---|
| 136 | const 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 |
|
|---|
| 156 | const 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 |
|
|---|
| 176 | const 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 |
|
|---|
| 196 | const 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 |
|
|---|
| 216 | const 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 |
|
|---|
| 236 | const 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 |
|
|---|
| 256 | const 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 |
|
|---|
| 276 | const 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 |
|
|---|
| 296 | const 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 |
|
|---|
| 316 | const 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 |
|
|---|
| 336 | const 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 |
|
|---|
| 356 | const 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 |
|
|---|
| 377 | G4XDeltaNstarTable::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 |
|
|---|
| 426 | G4XDeltaNstarTable::~G4XDeltaNstarTable()
|
|---|
| 427 | { }
|
|---|
| 428 |
|
|---|
| 429 |
|
|---|
| 430 | const 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 |
|
|---|
