| 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 // author: Vladimir.Grichine@cern.ch 26 // author: Vladimir.Grichine@cern.ch
27 // 27 //
28 // Implements data from: Barashenkov V.S., Nuc 28 // Implements data from: Barashenkov V.S., Nucleon-Nucleus Cross Section,
29 // Preprint JINR P2-89-770, p. 12, Dubna 1989 29 // Preprint JINR P2-89-770, p. 12, Dubna 1989 (scanned version from KEK)
30 // Based on G. Folger version of G4PiNuclearCr 30 // Based on G. Folger version of G4PiNuclearCrossSection class
31 // 31 //
32 // Modified: 16.08.2018 V.Ivanchenko - major u << 32 // Modified:
>> 33 // 05.03.07 V.Ivanchenko - add IfZAApplicable, remove "debug"
33 // 34 //
34 35
35 #include "G4NucleonNuclearCrossSection.hh" 36 #include "G4NucleonNuclearCrossSection.hh"
36 #include "G4DynamicParticle.hh" << 37 //#include "G4HadronicException.hh"
37 #include "G4ParticleDefinition.hh" << 38 //#include "G4HadTmpUtil.hh"
38 #include "G4ComponentBarNucleonNucleusXsc.hh" << 39 #include "G4Neutron.hh"
>> 40 #include "G4Proton.hh"
>> 41 //#include "G4ping.hh"
39 42
40 ////////////////////////////////////////////// << 43 // Group 1: He, Be, C for 44 energies
>> 44
>> 45 const G4double G4NucleonNuclearCrossSection::e1[44] =
>> 46 {
>> 47 0.014, 0.015, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.035, .04, 0.045, 0.05, .06, 0.07,
>> 48 .08, 0.09, .1, .12, .14, .15, .16, .18, .20, .25, .30, .35, .4 , 0.5, 0.6, 0.7, 0.8,
>> 49 0.9, 1, 1.5, 2, 3, 5, 7, 10,
>> 50 20, 50, 100, 500, 1000
>> 51 };
>> 52
>> 53 const G4double G4NucleonNuclearCrossSection::he_m_t[44] =
>> 54 {
>> 55 1090, 1020, 915, 800, 710, 640, 600, 560, 500, 440, 390, 360, 295, 256, 220, 192,
>> 56 168, 136, 120, 116, 114, 110, 107, 104, 106, 108, 110, 120, 126, 135, 140, 144, 146,
>> 57 148, 152, 150, 146, 142, 138, 132, 129, 126, 127, 128
>> 58 };
>> 59 const G4double G4NucleonNuclearCrossSection::he_m_in[44] =
>> 60 {
>> 61 0, 5, 10, 20, 35, 55, 70, 80, 90, 105, 115, 115, 100, 90, 86, 84, 84, 82, 80, 80, 80, 80,
>> 62 79, 78, 80, 84, 88, 94, 100, 105, 108, 108, 108, 112, 114, 114, 112, 110, 108, 106, 104,
>> 63 101, 102, 102
>> 64 };
>> 65 const G4double G4NucleonNuclearCrossSection::he_p_in[44] =
>> 66 {
>> 67 0, 2, 3, 13, 30, 50, 65, 77, 90, 105, 115, 115, 100, 90, 86, 84, 84, 82, 80, 80, 80, 80,
>> 68 79, 78, 80, 84, 88, 94, 100, 105, 108, 108, 108, 112, 114, 114, 112, 110, 108, 106, 104,
>> 69 101, 102, 102
>> 70 };
>> 71
>> 72 const G4double G4NucleonNuclearCrossSection::be_m_t[44] =
>> 73 {
>> 74 1490, 1460, 1400, 1350, 1270, 1200, 1160, 1100, 1000, 910, 810, 740, 625, 575, 455, 406,
>> 75 365, 310, 275, 262, 255, 240, 235, 225, 225, 230, 238, 252, 270, 282, 288, 290, 294, 303,
>> 76 303, 300, 292, 284, 277, 267, 263, 264, 268, 268
>> 77 };
>> 78 const G4double G4NucleonNuclearCrossSection::be_m_in[44] =
>> 79 {
>> 80 650, 640, 617, 595, 555, 520, 495, 470, 430, 385, 350, 320, 270, 250, 210, 190, 185, 178,
>> 81 175, 175, 175, 175, 175, 170, 170, 172, 176, 184, 194, 200, 209, 213, 214, 216, 216, 212,
>> 82 210, 210, 210, 210, 210, 210, 210, 210
>> 83 };
>> 84 const G4double G4NucleonNuclearCrossSection::be_p_in[44] =
>> 85 {
>> 86 490, 540, 580, 545, 525, 495, 470, 450, 420, 370, 340, 310, 262, 242, 205, 185, 180, 175,
>> 87 172, 175, 175, 175, 175, 170, 170, 172, 176, 184, 194, 200, 209, 213, 214, 216, 216, 212,
>> 88 210, 210, 210, 210, 210, 210, 210, 210
>> 89 };
>> 90
>> 91 const G4double G4NucleonNuclearCrossSection::c_m_t[44] =
>> 92 {
>> 93 1240, 1370, 1450, 1455, 1445, 1385, 1345, 1290, 1210, 1110, 1020, 940, 800, 700, 604, 530,
>> 94 475, 396, 350, 336, 320, 303, 294, 280, 280, 286, 296, 314, 330, 344, 356, 360, 364, 384,
>> 95 388, 384, 364, 352, 344, 330, 324, 324, 332, 332
>> 96 };
>> 97 const G4double G4NucleonNuclearCrossSection::c_m_in[44] =
>> 98 {
>> 99 590, 570, 542, 510, 500, 460, 445, 430, 395, 380, 350, 330, 295, 270, 255, 240, 228, 222,
>> 100 216, 216, 210, 210, 210, 208, 210, 214, 216, 228, 240, 248, 254, 257, 260, 262, 260, 256,
>> 101 252, 252, 250, 250, 248, 248, 248, 248
>> 102 };
>> 103 const G4double G4NucleonNuclearCrossSection::c_p_in[44] =
>> 104 {
>> 105 310, 330, 400, 440, 450, 435, 430, 420, 385, 370, 340, 320, 288, 263, 249, 234, 222, 216,
>> 106 210, 211, 205, 208, 210, 208, 210, 214, 216, 228, 240, 248, 254, 257, 260, 262, 260, 256,
>> 107 252, 252, 250, 250, 248, 248, 248, 248
>> 108 };
>> 109
>> 110 // Group 2: N, O, Na for 44 energies (e1=e2)
>> 111
>> 112 const G4double G4NucleonNuclearCrossSection::e2[44] =
>> 113 {
>> 114 0.014, 0.015, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.035, .04, 0.045, 0.05, .06, 0.07,
>> 115 .08, 0.09, .1, .12, .14, .15, .16, .18, .20, .25, .30, .35, .4 , 0.5, 0.6, 0.7, 0.8,
>> 116 0.9, 1, 1.5, 2, 3, 5, 7, 10,
>> 117 20, 50, 100, 500, 1000
>> 118 };
>> 119
>> 120 const G4double G4NucleonNuclearCrossSection::n_m_t[44] =
>> 121 {
>> 122 1420,1480, 1537, 1550, 1525, 1500, 1480, 1425, 1340, 1260, 1175, 1090, 930, 805, 690, 612,
>> 123 552, 462, 402, 384, 372, 350, 345, 326, 324, 328, 336, 356, 372, 388, 400, 408, 415, 430,
>> 124 435, 432, 415, 402, 390, 375, 367, 370, 382, 385
>> 125 };
>> 126 const G4double G4NucleonNuclearCrossSection::n_m_in[44] =
>> 127 {
>> 128 680, 665, 625, 580, 562, 525, 510, 485, 450, 435, 410, 387, 340, 310, 290, 280, 276, 274,
>> 129 260, 258, 254, 247, 245, 240, 240, 244, 250, 260, 268, 275, 280, 285, 290, 295, 300, 294,
>> 130 292, 290, 285, 285, 282, 282, 282, 282
>> 131 };
>> 132 const G4double G4NucleonNuclearCrossSection::n_p_in[44] =
>> 133 {
>> 134 420, 440, 470, 490, 497, 500, 480, 462, 440, 425, 400, 377, 333, 303, 284, 274, 270, 268,
>> 135 254, 252, 247, 245, 245, 240, 240, 244, 250, 260, 268, 275, 280, 285, 290, 295, 300, 294,
>> 136 292, 290, 285, 285, 282, 282, 282, 282
>> 137 };
>> 138
>> 139 const G4double G4NucleonNuclearCrossSection::o_m_t[44] =
>> 140 {
>> 141 1520, 1570, 1630, 1660, 1647, 1623, 1595, 1555, 1475, 1395, 1290, 1207, 1035, 925, 816,
>> 142 720, 645, 540, 462, 438, 415, 392, 378, 362, 361, 381, 390, 403, 417, 440, 460, 470,
>> 143 479, 498, 504, 498, 477, 457, 443, 427, 420, 425, 429, 430
>> 144 };
>> 145 const G4double G4NucleonNuclearCrossSection::o_m_in[44] =
>> 146 {
>> 147 750, 740, 700, 650, 620, 575, 555, 530, 505, 462, 435, 420, 375, 345, 320, 310, 300, 293,
>> 148 288, 282, 282, 280, 276, 270, 271, 275, 280, 290, 295, 304, 310, 315, 318, 332, 335, 330,
>> 149 323, 320, 317, 315, 315, 315, 315, 315
>> 150 };
>> 151 const G4double G4NucleonNuclearCrossSection::o_p_in[44] =
>> 152 {
>> 153 460, 485, 510, 535, 537, 532, 520, 500, 460, 432, 405, 390, 350, 320, 310, 304, 293, 287,
>> 154 283, 279, 279, 278, 276, 270, 271, 275, 280, 290, 295, 304, 310, 315, 318, 332, 335, 330,
>> 155 323, 320, 317, 315, 315, 315, 315, 315
>> 156 };
>> 157
>> 158 const G4double G4NucleonNuclearCrossSection::na_m_t[44] =
>> 159 {
>> 160 1570, 1620, 1695, 1730, 1750, 1760, 1755, 1740, 1710, 1643, 1560, 1480, 1343, 1220, 1073,
>> 161 953, 860, 720, 618, 582, 546, 522, 504, 484, 492, 500, 512, 538, 560, 586, 608, 622, 632,
>> 162 660, 668, 664, 640, 616, 596, 568, 568, 568, 568, 568
>> 163 };
>> 164 const G4double G4NucleonNuclearCrossSection::na_m_in[44] =
>> 165 {
>> 166 960, 930, 890, 822, 790, 750, 725, 686, 620, 600, 575, 540, 497, 450, 414, 390, 380, 372,
>> 167 354, 360, 355, 354, 350, 350, 350, 356, 364, 384, 392, 400, 408, 410, 420, 408, 412, 420,
>> 168 411, 409, 407, 403, 400, 400, 400, 400
>> 169 };
>> 170 const G4double G4NucleonNuclearCrossSection::na_p_in[44] =
>> 171 {
>> 172 600, 617, 660, 675, 680, 680, 670, 650, 575, 550, 525, 490, 450, 420, 385, 367, 360, 350,
>> 173 350, 350, 345, 347, 350, 350, 350, 356, 364, 384, 392, 400, 408, 410, 420, 408, 412, 420,
>> 174 411, 409, 407, 403, 400, 400, 400, 400
>> 175 };
>> 176
>> 177 // Al, Si, Ca for 45 energies
>> 178
>> 179 const G4double G4NucleonNuclearCrossSection::e3[45] =
>> 180 {
>> 181 0.014, 0.015, 0.016, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.035, .04, 0.045, 0.05,
>> 182 .06, 0.07, .08, 0.09, .1, .12, .14, .15, .16, .18, .20, .25, .30, .35, .4 , 0.5, 0.6,
>> 183 0.7, 0.8, 0.9, 1, 1.5, 2, 3, 5, 7, 10, 20, 50, 100, 500, 1000
>> 184 };
>> 185
>> 186 const G4double G4NucleonNuclearCrossSection::al_m_t[45] =
>> 187 {
>> 188 1735, 1750, 1760, 1795, 1830, 1855, 1885, 1895, 1900, 1870, 1835, 1785, 1710, 1522,
>> 189 1350, 1212, 1080, 972, 816, 720, 678, 642, 600, 567, 558, 560, 578, 592, 616, 644,
>> 190 672, 688, 708, 720, 736, 754, 736, 706, 680, 672, 646, 632, 632, 632, 632
>> 191 };
>> 192 const G4double G4NucleonNuclearCrossSection::al_m_in[45] =
>> 193 {
>> 194 1000, 990, 975, 950, 905, 875, 825, 800, 762, 690, 652, 610, 570, 495, 480, 456, 444,
>> 195 432, 420, 420, 420, 420, 410, 410, 400, 402, 404, 408, 424, 438, 448, 450, 454, 456,
>> 196 472, 480, 466, 456, 452, 448, 444, 440, 440, 440, 440
>> 197 };
>> 198 const G4double G4NucleonNuclearCrossSection::al_p_in[45] =
>> 199 {
>> 200 650, 682, 690, 715, 750, 762, 750, 740, 720, 655, 617, 575, 540, 470, 455, 532, 420,
>> 201 408, 400, 403, 403, 408, 406, 404, 400, 402, 404, 408, 424, 438, 448, 450, 454, 456,
>> 202 472, 480, 466, 456, 452, 448, 444, 440, 440, 440, 440
>> 203 };
>> 204
>> 205 const G4double G4NucleonNuclearCrossSection::si_m_t[45] =
>> 206 {
>> 207 1810, 1833, 1850, 1872, 1920, 1950, 1995, 2020, 2035, 2000, 1930, 1850, 1760, 1570,
>> 208 1400, 1255, 1110, 1008, 846, 742, 696, 671, 623, 588, 584, 584, 602, 618, 645, 679,
>> 209 708, 727, 746, 757, 769, 782, 771, 734, 710, 698, 672, 654, 650, 650, 650
>> 210 };
>> 211 const G4double G4NucleonNuclearCrossSection::si_m_in[45] =
>> 212 {
>> 213 1060, 1035, 1015, 990, 935, 900, 860, 830, 790, 725, 665, 630, 600, 520, 504, 486,
>> 214 470, 456, 444, 432, 432, 432, 418, 418, 415, 412, 416, 422, 440, 460, 472, 476, 479,
>> 215 480, 492, 496, 488, 472, 472, 464, 460, 452, 448, 448, 448
>> 216 };
>> 217 const G4double G4NucleonNuclearCrossSection::si_p_in[45] =
>> 218 {
>> 219 670, 700, 725, 750, 780, 780, 770, 757, 735, 690, 635, 585, 570, 490, 475, 460, 446,
>> 220 431, 423, 425, 425, 425, 425, 422, 422, 412, 416, 422, 440, 460, 472, 476, 479,
>> 221 480, 492, 496, 488, 472, 472, 464, 460, 452, 448, 448, 448
>> 222 };
>> 223
>> 224 const G4double G4NucleonNuclearCrossSection::ca_m_t[45] =
>> 225 {
>> 226 2180, 2130, 2095, 2075, 2115, 2150, 2220, 2250, 2300, 2365, 2360, 2280, 2180, 2000,
>> 227 1805, 1650, 1500, 1340, 1140, 990, 940, 890, 825, 790, 770, 773, 787, 800, 830, 870,
>> 228 905, 930, 950, 965, 990, 1002, 990, 965, 945, 925, 892, 860, 860, 860, 860
>> 229 };
>> 230 const G4double G4NucleonNuclearCrossSection::ca_m_in[45] =
>> 231 {
>> 232 1240, 1225, 1200, 1180, 1125, 1090, 1045, 1020, 980, 925, 880, 825, 770, 680, 640,
>> 233 620, 615, 600, 580, 565, 560, 560, 560, 550, 535, 530, 540, 550, 570, 595, 610, 615,
>> 234 620, 622, 629, 630, 620, 612, 607, 592, 587, 580, 580, 580, 580
>> 235 };
>> 236 const G4double G4NucleonNuclearCrossSection::ca_p_in[45] =
>> 237 {
>> 238 770, 800, 823, 850, 900, 925, 935, 920, 895, 835, 800, 750, 715, 640, 605, 590, 588,
>> 239 573, 555, 543, 540, 540, 540, 535, 530, 530, 540, 550, 570, 595, 610, 615,
>> 240 620, 622, 629, 630, 620, 612, 607, 592, 587, 580, 580, 580, 580
>> 241 };
>> 242
>> 243 // Fe, Cu, Mo for 47 energies
>> 244
>> 245 const G4double G4NucleonNuclearCrossSection::e4[47] =
>> 246 {
>> 247 0.014, 0.015, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.033, 0.035, 0.037, .04, 0.045,
>> 248 0.05, 0.055, .06, 0.07, .08, 0.09, .1, .12, .14, .15, .16, .18, .20, .25, .30, .35,
>> 249 .4 , 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5, 2, 3, 5, 7, 10,
>> 250 20, 50, 100, 500, 1000
>> 251 };
>> 252
>> 253 const G4double G4NucleonNuclearCrossSection::fe_m_t[47] =
>> 254 {
>> 255 2580, 2490, 2370, 2282, 2275, 2285, 2320, 2370, 2432, 2445, 2460, 2485, 2530, 2540,
>> 256 2517, 2480, 2290, 2110, 1940, 1790, 1510, 1290, 1220, 1150, 1070, 1030, 1013, 1020,
>> 257 1030, 1043, 1075, 1110, 1133, 1163, 1185, 1225, 1252, 1260, 1260, 1233, 1207, 1185,
>> 258 1140, 1110, 1110, 1110, 1110
>> 259 };
>> 260 const G4double G4NucleonNuclearCrossSection::fe_m_in[47] =
>> 261 {
>> 262 1440, 1433, 1390, 1325, 1280, 1260, 1215, 1180, 1140, 1110, 1080, 1040, 990, 955, 920,
>> 263 885, 835, 800, 780, 765, 750, 725, 720, 720, 710, 700, 700, 700, 712, 705, 735, 750,
>> 264 765, 775, 780, 795, 810, 813, 810, 784, 757, 743, 735, 720, 720, 720, 720
>> 265 };
>> 266 const G4double G4NucleonNuclearCrossSection::fe_p_in[47] =
>> 267 {
>> 268 900, 960, 1070, 1090, 1115, 1120, 1115, 1080, 1045, 1025, 1000, 960, 900, 885, 865, 790,
>> 269 765, 740, 720, 700, 697, 697, 697, 697, 695, 690, 688, 690, 712, 705, 735, 750,
>> 270 765, 775, 780, 795, 810, 813, 810, 784, 757, 743, 735, 720, 720, 720, 720
>> 271 };
>> 272
>> 273 const G4double G4NucleonNuclearCrossSection::cu_m_t[47] =
>> 274 {
>> 275 2920, 2800, 2615, 2480, 2455, 2430, 2440, 2460, 2500, 2530, 2560, 2615, 2690, 2720,
>> 276 2700, 2645, 2500, 2320, 2140, 1970, 1670, 1460, 1380, 1285, 1200, 1160, 1140, 1147,
>> 277 1163, 1170, 1200, 1237, 1265, 1285, 1305, 1328, 1375, 1390, 1395, 1370, 1335, 1315,
>> 278 1270, 1230, 1230, 1230, 1230
>> 279 };
>> 280 const G4double G4NucleonNuclearCrossSection::cu_m_in[47] =
>> 281 {
>> 282 1540, 1535, 1500, 1445, 1407, 1380, 1330, 1300, 1285, 1270, 1240, 1190, 1090, 1010,
>> 283 940, 920, 860, 835, 820, 810, 800, 780, 775, 770, 760, 760, 758, 765, 765, 770, 795,
>> 284 810, 825, 830, 840, 848, 870, 870, 868, 840, 825, 810, 803, 795, 795, 795, 795
>> 285 };
>> 286 const G4double G4NucleonNuclearCrossSection::cu_p_in[47] =
>> 287 {
>> 288 935, 1000, 1060, 1190, 1220, 1250, 1240, 1210, 1150, 1130, 1115, 1050, 985, 950, 890,
>> 289 870, 820, 800, 785, 780, 770, 750, 745, 740, 735, 735, 745, 760, 762, 770, 795,
>> 290 810, 825, 830, 840, 848, 870, 870, 868, 840, 825, 810, 803, 795, 795, 795, 795
>> 291 };
>> 292
>> 293 const G4double G4NucleonNuclearCrossSection::mo_m_t[47] =
>> 294 {
>> 295 4150, 4040, 3800, 3490, 3300, 3060, 2960, 2845, 2785, 2820, 2850, 2980, 3170, 3230,
>> 296 3270, 3280, 3225, 3075, 2895, 2710, 2355, 2060, 1925, 1800, 1630, 1560, 1540, 1550,
>> 297 1570, 1590, 1650, 1685, 1715, 1740, 1760, 1780, 1850, 1880, 1858, 1815, 1790, 1782,
>> 298 1720, 1690, 1690, 1690, 1690
>> 299 };
>> 300 const G4double G4NucleonNuclearCrossSection::mo_m_in[47] =
>> 301 {
>> 302 1790, 1775, 1740, 1680, 1640, 1580, 1550, 1510, 1460, 1440, 1418, 1380, 1330, 1280,
>> 303 1240, 1200, 1155, 1140, 1110, 1110, 1080, 1065, 1050, 1050, 1025, 1020, 1015, 1020,
>> 304 1022, 1026, 1060, 1085, 1100, 1110, 1120, 1127, 1150, 1160, 1140, 1100, 1085, 1080,
>> 305 1070, 1070, 1070, 1070, 1070
>> 306 };
>> 307 const G4double G4NucleonNuclearCrossSection::mo_p_in[47] =
>> 308 {
>> 309 1025, 1080, 1190, 1380, 1440, 1495, 1475, 1420, 1350, 1310, 1300, 1290, 1250, 1200,
>> 310 1170, 1130, 1095, 1060, 1040, 1022, 1020, 1016, 1016, 1016, 1016, 1012, 1005, 1005,
>> 311 1005, 1010, 1060, 1085, 1100, 1110, 1120, 1127, 1150, 1160, 1140, 1100, 1085, 1080,
>> 312 1070, 1070, 1070, 1070, 1070
>> 313 };
>> 314
>> 315 // Cd, Sn, W for 48 energies
>> 316
>> 317 const G4double G4NucleonNuclearCrossSection::e5[48] =
>> 318 {
>> 319 0.014, 0.015, 0.017, 0.018, .02, 0.022, 0.025, 0.027, 0.03, 0.033, 0.035, .04,
>> 320 0.045, 0.05, 0.055, .06, .065, 0.07, .08, 0.09, .1, .12, .14, .15, .16, .18,
>> 321 .20, .25, .30, .35, .4 , 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5, 2, 3, 5, 7, 10, 20,
>> 322 50, 100, 500, 1000
>> 323 };
>> 324
>> 325 const G4double G4NucleonNuclearCrossSection::cd_m_t[48] =
>> 326 {
>> 327 4420, 4280, 4170, 4070, 3860, 3680, 3420, 3280, 3125, 3060, 3080, 3190, 3350, 3445,
>> 328 3510, 3540, 3560, 3550, 3460, 3300, 3030, 2640, 2340, 2190, 2070, 1950, 1770, 1732,
>> 329 1740, 1760, 1780, 1832, 1885, 1925, 1945, 1960, 1980, 2070, 2080, 2065, 2040, 2022,
>> 330 1980, 1940, 1870, 1870, 1870, 1870
>> 331 };
>> 332 const G4double G4NucleonNuclearCrossSection::cd_m_in[48]=
>> 333 {
>> 334 1920, 1910, 1880, 1860, 1840, 1800, 1760, 1720, 1675, 1630, 1600, 1520, 1465, 1420,
>> 335 1390, 1340, 1310, 1280, 1275, 1235, 1225, 1200, 1170, 1170, 1170, 1165, 1145, 1140,
>> 336 1140, 1135, 1160, 1180, 1220, 1240, 1250, 1260, 1265, 1270, 1275, 1250, 1222, 1222,
>> 337 1220, 1215, 1190, 1190, 1190, 1190
>> 338 };
>> 339 const G4double G4NucleonNuclearCrossSection::cd_p_in[48] =
>> 340 {
>> 341 1020, 1100, 1225, 1290, 1440, 1520, 1575, 1560, 1518, 1460, 1420, 1400, 1365, 1340,
>> 342 1300, 1280, 1260, 1200, 1190, 1160, 1125, 1125, 1125, 1125, 1125, 1125, 1120, 1120,
>> 343 1120, 1118, 1146, 1180, 1220, 1240, 1250, 1260, 1265, 1270, 1275, 1250, 1222, 1222,
>> 344 1220, 1215, 1190, 1190, 1190, 1190
>> 345 };
>> 346
>> 347 const G4double G4NucleonNuclearCrossSection::sn_m_t[48] =
>> 348 {
>> 349 4420, 4400, 4260, 4150, 3980, 3770, 3530, 3370, 3245, 3180, 3170, 3260, 3400, 3500,
>> 350 3560, 3610, 3650, 3680, 3580, 3390, 3190, 2760, 2430, 2295, 2175, 1990, 1880, 1810,
>> 351 1820, 1840, 1865, 1940, 1985, 2020, 2040, 2060, 2080, 2160, 2185, 2180, 2110, 2105,
>> 352 2080, 2050, 1980, 1980, 1980, 1980
>> 353 };
>> 354 const G4double G4NucleonNuclearCrossSection::sn_m_in[48] =
>> 355 {
>> 356 1945, 1940, 1905, 1890, 1860, 1830, 1780, 1755, 1717, 1680, 1645, 1570, 1500, 1455,
>> 357 1410, 1370, 1340, 1320, 1290, 1285, 1260, 1240, 1235, 1212, 1200, 1200, 1200, 1190,
>> 358 1190, 1200, 1210, 1240, 1270, 1285, 1300, 1300, 1310, 1320, 1320, 1290, 1240, 1240,
>> 359 1240, 1240, 1240, 1240, 1240, 1240
>> 360 };
>> 361 const G4double G4NucleonNuclearCrossSection::sn_p_in[48] =
>> 362 {
>> 363 1020, 1080, 1270, 1335, 1465, 1505, 1610, 1610, 1550, 1535, 1500, 1440, 1407, 1370,
>> 364 1340, 1300, 1285, 1260, 1230, 1215, 1200, 1180, 1170, 1170, 1165, 1165, 1170, 1165,
>> 365 1165, 1183, 1195, 1240, 1270, 1285, 1300, 1300, 1310, 1320, 1320, 1290, 1240, 1240,
>> 366 1240, 1240, 1240, 1240, 1240, 1240
>> 367 };
>> 368
>> 369 const G4double G4NucleonNuclearCrossSection::w_m_t[48] =
>> 370 {
>> 371 5320, 5430, 5480, 5450, 5330, 5190, 4960, 4790, 4550, 4340, 4200, 4070, 4000, 4030,
>> 372 4125, 4220, 4270, 4390, 4440, 4360, 4200, 3800, 3380, 3200, 3040, 2790, 2660, 2575,
>> 373 2575, 2600, 2640, 2690, 2755, 2790, 2812, 2837, 2850, 2950, 3000, 2970, 2940, 2910,
>> 374 2880, 2820, 2730, 2730, 2730, 2730
>> 375 };
>> 376 const G4double G4NucleonNuclearCrossSection::w_m_in[48] =
>> 377 {
>> 378 2440, 2400, 2370, 2350, 2310, 2270, 2220, 2195, 2150, 2100, 2070, 2010, 1945, 1900,
>> 379 1850, 1820, 1780, 1760, 1730, 1720, 1680, 1680, 1660, 1660, 1650, 1650, 1640, 1640,
>> 380 1612, 1615, 1625, 1640, 1700, 1720, 1730, 1740, 1750, 1780, 1780, 1750, 1740, 1735,
>> 381 1710, 1695, 1680, 1680, 1680, 1680
>> 382 };
>> 383 const G4double G4NucleonNuclearCrossSection::w_p_in[48] =
>> 384 {
>> 385 950, 1020, 1240, 1400, 1560, 1670, 1760, 1830, 1850, 1855, 1870, 1840, 1800, 1770,
>> 386 1740, 1715, 1680, 1670, 1650, 1620, 1610, 1600, 1600, 1600, 1600, 1600, 1600, 1595,
>> 387 1585, 1595, 1615, 1640, 1700, 1720, 1730, 1740, 1750, 1780, 1780, 1750, 1740, 1735,
>> 388 1710, 1695, 1680, 1680, 1680, 1680
>> 389 };
>> 390
>> 391 // Pb, U for 46 energies
>> 392
>> 393 const G4double G4NucleonNuclearCrossSection::e6[46] =
>> 394 {
>> 395 0.014, 0.015, 0.017, 0.019, .02, 0.022, 0.025, 0.027, 0.03, 0.035, .04, 0.045, 0.05,
>> 396 0.055, .06, 0.07, .08, 0.09, .1, .12, .14, .15, .16, .18, .20, .25, .30, .35, .4 ,
>> 397 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5, 2, 3, 5, 7, 10, 20, 50, 100, 500, 1000
>> 398 };
>> 399
>> 400 const G4double G4NucleonNuclearCrossSection::pb_m_t[46] =
>> 401 {
>> 402 5300, 5440, 5720, 5880, 5765, 5745, 5480, 5280, 4970, 4550, 4390, 4300, 4265, 4325,
>> 403 4450, 4540, 4740, 4710, 4600, 4100, 3660, 3480, 3300, 3000, 2890, 2865, 2855, 2850,
>> 404 2865, 2920, 2955, 3000, 3030, 3060, 3105, 3240, 3290, 3270, 3240, 3180, 3090, 3060,
>> 405 2970, 2970, 2970, 2970
>> 406
>> 407 };
>> 408 const G4double G4NucleonNuclearCrossSection::pb_m_in[46] =
>> 409 {
>> 410 2580, 2550, 2505, 2462, 2460, 2435, 2380, 2355, 2280, 2180, 2170, 2130, 2080, 2035,
>> 411 1980, 1940, 1900, 1870, 1840, 1800, 1800, 1800, 1780, 1760, 1760, 1740, 1730, 1725,
>> 412 1740, 1785, 1815, 1835, 1860, 1890, 1895, 1920, 1920, 1890, 1850, 1835, 1830, 1830,
>> 413 1830, 1830, 1830, 1830
>> 414 };
>> 415 const G4double G4NucleonNuclearCrossSection::pb_p_in[46] =
>> 416 {
>> 417 900, 1060, 1200, 1420, 1515, 1620, 1750, 1800, 1915, 2030, 1960, 1940, 1910, 1860,
>> 418 1840, 1780, 1770, 1760, 1740, 1720, 1725, 1740, 1740, 1730, 1720, 1700, 1710, 1720,
>> 419 1730, 1740, 1815, 1835, 1860, 1890, 1895, 1920, 1920, 1890, 1850, 1835, 1830, 1830,
>> 420 1830, 1830, 1830, 1830
>> 421 };
>> 422
>> 423 const G4double G4NucleonNuclearCrossSection::u_m_t[46] =
>> 424 {
>> 425 5800, 5940, 6160, 6345, 6360, 6350, 6170, 6020, 5760, 5350, 4990, 4800, 4710, 4690,
>> 426 4760, 5040, 5190, 5200, 5080, 4600, 4120, 3920, 3720, 3420, 3240, 3150, 3160, 3180,
>> 427 3210, 3240, 3280, 3350, 3390, 3435, 3480, 3560, 3585, 3580, 3540, 3500, 3470, 3410,
>> 428 3335, 3335, 3335, 3335
>> 429 };
>> 430 const G4double G4NucleonNuclearCrossSection::u_m_in[46] =
>> 431 {
>> 432 2820, 2770, 2700, 2660, 2645, 2620, 2580, 2550, 2515, 2450, 2390, 2320, 2260, 2225,
>> 433 2200, 2140, 2080, 2060, 2040, 2000, 1980, 1965, 1960, 1930, 1920, 1890, 1905, 1920,
>> 434 1945, 1970, 1985, 2010, 2040, 2070, 2080, 2090, 2095, 2080, 2063, 2060, 2050, 2040,
>> 435 2005, 2005, 2005, 2005
>> 436 };
>> 437 const G4double G4NucleonNuclearCrossSection::u_p_in[46] =
>> 438 {
>> 439 800, 900, 1100, 1300, 1410, 1510, 1680, 1800, 2000, 2200, 2080, 2060, 2035, 2100,
>> 440 2030, 2030, 2000, 1960, 1960, 1960, 1940, 1925, 1920, 1905, 1890, 1860, 1880, 1910,
>> 441 1930, 1945, 1985, 2010, 2040, 2070, 2080, 2090, 2095, 2080, 2063, 2060, 2050, 2040,
>> 442 2005, 2005, 2005, 2005
>> 443 };
>> 444
>> 445 using namespace std;
>> 446
>> 447 //////////////////////////////////////////////////////////////////////////////////
>> 448 //
>> 449 //
41 450
42 G4NucleonNuclearCrossSection::G4NucleonNuclear 451 G4NucleonNuclearCrossSection::G4NucleonNuclearCrossSection()
43 : G4VCrossSectionDataSet(Default_Name()) <<
44 { 452 {
45 fBarash = new G4ComponentBarNucleonNucleusXs << 453 theNeutron = G4Neutron::Neutron();
>> 454 theProton = G4Proton::Proton();
>> 455
>> 456 // He, Be, C
>> 457
>> 458 thePimData.push_back(new G4PiData(he_m_t, he_m_in, e1, 44));
>> 459 thePipData.push_back(new G4PiData(he_m_t, he_p_in, e1, 44));
>> 460
>> 461 thePimData.push_back(new G4PiData(be_m_t, be_m_in, e1, 44));
>> 462 thePipData.push_back(new G4PiData(be_m_t, be_p_in, e1, 44));
>> 463
>> 464 thePimData.push_back(new G4PiData(c_m_t, c_m_in, e1, 44));
>> 465 thePipData.push_back(new G4PiData(c_m_t, c_p_in, e1, 44));
>> 466
>> 467 // N, O, Na
>> 468
>> 469 thePimData.push_back(new G4PiData(n_m_t, n_m_in, e2, 44));
>> 470 thePipData.push_back(new G4PiData(n_m_t, n_p_in, e2, 44));
>> 471
>> 472 thePimData.push_back(new G4PiData(o_m_t, o_m_in, e2, 44));
>> 473 thePipData.push_back(new G4PiData(o_m_t, o_p_in, e2, 44));
>> 474
>> 475 thePimData.push_back(new G4PiData(na_m_t, na_m_in, e2, 44));
>> 476 thePipData.push_back(new G4PiData(na_m_t, na_p_in, e2, 44));
>> 477
>> 478 // Al, Si, Ca
>> 479
>> 480 thePimData.push_back(new G4PiData(al_m_t, al_m_in, e3, 45));
>> 481 thePipData.push_back(new G4PiData(al_m_t, al_p_in, e3, 45));
>> 482
>> 483 thePimData.push_back(new G4PiData(si_m_t, si_m_in, e3, 45));
>> 484 thePipData.push_back(new G4PiData(si_m_t, si_p_in, e3, 45));
>> 485
>> 486 thePimData.push_back(new G4PiData(ca_m_t, ca_m_in, e3, 45));
>> 487 thePipData.push_back(new G4PiData(ca_m_t, ca_p_in, e3, 45));
>> 488
>> 489 // Fe, Cu, Mo
>> 490
>> 491 thePimData.push_back(new G4PiData(fe_m_t, fe_m_in, e4, 47));
>> 492 thePipData.push_back(new G4PiData(fe_m_t, fe_p_in, e4, 47));
>> 493
>> 494 thePimData.push_back(new G4PiData(cu_m_t, cu_m_in, e4, 47));
>> 495 thePipData.push_back(new G4PiData(cu_m_t, cu_p_in, e4, 47));
>> 496
>> 497 thePimData.push_back(new G4PiData(mo_m_t, mo_m_in, e4, 47));
>> 498 thePipData.push_back(new G4PiData(mo_m_t, mo_p_in, e4, 47));
>> 499
>> 500 // Cd, Sn, W
>> 501
>> 502 thePimData.push_back(new G4PiData(cd_m_t, cd_m_in, e5, 48));
>> 503 thePipData.push_back(new G4PiData(cd_m_t, cd_p_in, e5, 48));
>> 504
>> 505 thePimData.push_back(new G4PiData(sn_m_t, sn_m_in, e5, 48));
>> 506 thePipData.push_back(new G4PiData(sn_m_t, sn_p_in, e5, 48));
>> 507
>> 508 thePimData.push_back(new G4PiData(w_m_t, w_m_in, e5, 48));
>> 509 thePipData.push_back(new G4PiData(w_m_t, w_p_in, e5, 48));
>> 510
>> 511 // Pb, U
>> 512
>> 513 thePimData.push_back(new G4PiData(pb_m_t, pb_m_in, e6, 46));
>> 514 thePipData.push_back(new G4PiData(pb_m_t, pb_p_in, e6, 46));
>> 515
>> 516 thePimData.push_back(new G4PiData(u_m_t, u_m_in, e6, 46));
>> 517 thePipData.push_back(new G4PiData(u_m_t, u_p_in, e6, 46));
>> 518
>> 519 theZ.push_back(2); // He
>> 520 theZ.push_back(4); // Be
>> 521 theZ.push_back(6); // C
>> 522 theZ.push_back(7); // N
>> 523 theZ.push_back(8); // O
>> 524 theZ.push_back(11); // Na
>> 525 theZ.push_back(13); // Al
>> 526 theZ.push_back(14); // Si
>> 527 theZ.push_back(20); // Ca
>> 528 theZ.push_back(26); // Fe
>> 529 theZ.push_back(29); // Cu
>> 530 theZ.push_back(42); // Mo
>> 531 theZ.push_back(48); // Cd
>> 532 theZ.push_back(50); // Sn
>> 533 theZ.push_back(74); // W
>> 534 theZ.push_back(82); // Pb
>> 535 theZ.push_back(92); // U
>> 536
>> 537 }
>> 538
>> 539 ///////////////////////////////////////////////////////////////////////////////////
>> 540 //
>> 541 //
>> 542
>> 543 G4NucleonNuclearCrossSection::~G4NucleonNuclearCrossSection()
>> 544 {
>> 545 std::for_each(thePimData.begin(), thePimData.end(), G4PiData::Delete());
>> 546 std::for_each(thePipData.begin(), thePipData.end(), G4PiData::Delete());
46 } 547 }
47 548
48 ////////////////////////////////////////////// 549 ////////////////////////////////////////////////////////////////////////////
>> 550 //
>> 551 //
>> 552
>> 553 G4double G4NucleonNuclearCrossSection::
>> 554 GetCrossSection( const G4DynamicParticle* aParticle,
>> 555 const G4Element* anElement,
>> 556 G4double )
49 557
50 G4bool G4NucleonNuclearCrossSection::IsElement <<
51 const G4DynamicParticle*, G4int Z, cons <<
52 { 558 {
53 return (Z > 1); << 559 return GetIsoZACrossSection(aParticle, anElement->GetZ(), anElement->GetN(), 0.);
54 } 560 }
55 561
56 ////////////////////////////////////////////// 562 ////////////////////////////////////////////////////////////////////////////
>> 563 //
>> 564 //
57 565
58 G4double G4NucleonNuclearCrossSection::GetElem << 566 G4bool G4NucleonNuclearCrossSection::IsApplicable(const G4DynamicParticle* aParticle,
59 const G4DynamicParticle* dp, G4int Z, << 567 const G4Element* anElement)
60 { 568 {
61 ComputeCrossSections(dp->GetDefinition(), dp << 569 return IsZAApplicable(aParticle, anElement->GetZ(), anElement->GetN());
62 return fInelasticXsc; <<
63 } 570 }
64 571
65 ////////////////////////////////////////////// 572 ////////////////////////////////////////////////////////////////////////////
>> 573 //
>> 574 //
66 575
67 void G4NucleonNuclearCrossSection::ComputeCros << 576 G4bool G4NucleonNuclearCrossSection::IsZAApplicable(const G4DynamicParticle* aParticle,
68 const G4ParticleDefinition* pd, << 577 G4double Z, G4double)
69 G4double kinEnergy, G4int Z) << 578 {
70 { << 579 G4bool result = false;
71 fBarash->ComputeCrossSections(pd, kinEnergy, << 580 if(aParticle->GetDefinition() == theNeutron ) result = true;
72 fTotalXsc = fBarash->GetTotalXsc(); << 581 if(aParticle->GetDefinition() == theProton) result = true;
73 fInelasticXsc = fBarash->GetInelasticXsc(); << 582 if(Z < 1.5) result = false;
74 fElasticXsc = fBarash->GetElasticXsc(); << 583 if(aParticle->GetKineticEnergy() > 999.9*GeV) result = false;
>> 584 return result;
75 } 585 }
76 586
77 ////////////////////////////////////////////// 587 ////////////////////////////////////////////////////////////////////////////
>> 588 //
>> 589 //
78 590
79 void << 591 G4double G4NucleonNuclearCrossSection::
80 G4NucleonNuclearCrossSection::CrossSectionDesc << 592 GetIsoZACrossSection( const G4DynamicParticle* aParticle,
>> 593 G4double zElement, G4double, G4double )
81 { 594 {
82 outFile << "G4NucleonNuclearCrossSection is << 595 G4double kineticEnergy = aParticle->GetKineticEnergy();
83 << "cross section parameterization t << 596
84 << "nucleons on targets heavier than << 597 G4double result = 0;
85 << "use as a cross section component << 598 G4int Z = G4int(zElement + 0.5);
86 << "G4BGGNucleonInelasticXS. It is << 599
87 << "to 1 TeV.\n"; << 600 // G4cout<<"Z = "<<Z<<G4endl;
>> 601
>> 602 size_t it = 0;
>> 603
>> 604 while( it < theZ.size() && Z > theZ[it] ) it++;
>> 605
>> 606 // if( Z == 29 ) tuning=.96;
>> 607
>> 608 if( Z > theZ[it] ) Z = theZ[it];
>> 609 G4int Z1, Z2;
>> 610 G4double x1, x2, xt1, xt2;
>> 611
>> 612 if(aParticle->GetDefinition() == theNeutron )
>> 613 {
>> 614 if( theZ[it] == Z )
>> 615 {
>> 616 result = thePimData[it]->ReactionXSection(kineticEnergy);
>> 617 fTotalXsc = thePimData[it]->TotalXSection(kineticEnergy);
>> 618 }
>> 619 else
>> 620 {
>> 621 x1 = thePimData[it-1]->ReactionXSection(kineticEnergy);
>> 622 xt1 = thePimData[it-1]->TotalXSection(kineticEnergy);
>> 623 Z1 = theZ[it-1];
>> 624 x2 = thePimData[it]->ReactionXSection(kineticEnergy);
>> 625 xt2 = thePimData[it]->TotalXSection(kineticEnergy);
>> 626 Z2 = theZ[it];
>> 627
>> 628 result = Interpolate(Z1, Z2, Z, x1, x2);
>> 629 fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2);
>> 630 }
>> 631 }
>> 632 else if(aParticle->GetDefinition() == theProton)
>> 633 {
>> 634 if( theZ[it] == Z )
>> 635 {
>> 636 // at high energies, when no data for proton, use neutron
>> 637
>> 638 std::vector<G4PiData *> * theData = &thePimData;
>> 639
>> 640 if( thePipData[it]->AppliesTo(kineticEnergy) )
>> 641 {
>> 642 theData = &thePipData;
>> 643 }
>> 644 result = theData->operator[](it)->ReactionXSection(kineticEnergy);
>> 645 fTotalXsc = theData->operator[](it)->TotalXSection(kineticEnergy);
>> 646
>> 647 }
>> 648 else
>> 649 {
>> 650 std::vector<G4PiData*>* theLData = &thePimData;
>> 651
>> 652 if(thePipData[it-1]->AppliesTo(kineticEnergy))
>> 653 {
>> 654 theLData = &thePipData;
>> 655 }
>> 656 std::vector<G4PiData *> * theHData = &thePimData;
>> 657
>> 658 if( thePipData[it]->AppliesTo(kineticEnergy) )
>> 659 {
>> 660 theHData = &thePipData;
>> 661 }
>> 662 x1 = theLData->operator[](it-1)->ReactionXSection(kineticEnergy);
>> 663 xt1 = theLData->operator[](it-1)->TotalXSection(kineticEnergy);
>> 664 Z1 = theZ[it-1];
>> 665 x2 = theHData->operator[](it)->ReactionXSection(kineticEnergy);
>> 666 xt2 = theHData->operator[](it)->TotalXSection(kineticEnergy);
>> 667 Z2 = theZ[it];
>> 668
>> 669 result = Interpolate(Z1, Z2, Z, x1, x2);
>> 670 fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2);
>> 671 }
>> 672 }
>> 673 fElasticXsc = fTotalXsc - result;
>> 674 if( fElasticXsc < 0.) fElasticXsc = 0.;
>> 675
>> 676 return result;
88 } 677 }
89 678
90 ////////////////////////////////////////////// << 679 /////////////////////////////////////////////////////////////////////////////
>> 680 //
>> 681 //
91 682
>> 683 G4double G4NucleonNuclearCrossSection::
>> 684 Interpolate(G4int Z1, G4int Z2, G4int Z, G4double x1, G4double x2)
>> 685 {
>> 686 // Nucleon numbers obtained from G4NistManager G4 8.0
>> 687
>> 688 static G4double alpha = 2./3.;
>> 689
>> 690 static const G4double A[92] =
>> 691 {
>> 692 1.0001, 4.0000, 6.9241, 9.0000, 10.801, 12.011, 14.004, 16.004, 19.000, 20.188,
>> 693 23.000, 24.320, 27.000, 28.109, 31.000, 32.094, 35.484, 39.985, 39.135, 40.116,
>> 694 45.000, 47.918, 50.998, 52.055, 55.000, 55.910, 59.000, 58.760, 63.617, 65.468,
>> 695 69.798, 72.691, 75.000, 79.042, 79.986, 83.887, 85.557, 87.710, 89.000, 91.318,
>> 696 93.000, 96.025, 98.000, 101.16, 103.00, 106.51, 107.96, 112.51, 114.91, 118.81,
>> 697 121.86, 127.70, 127.00, 131.39, 133.00, 137.42, 139.00, 140.21, 141.00, 144.32,
>> 698 145.00, 150.45, 152.04, 157.33, 159.00, 162.57, 165.00, 167.32, 169.00, 173.10,
>> 699 175.03, 178.54, 181.00, 183.89, 186.25, 190.27, 192.25, 195.11, 197.00, 200.63,
>> 700 204.41, 207.24, 209.00, 209.00, 210.00, 222.00, 223.00, 226.00, 227.00, 232.00,
>> 701 231.00, 237.98
>> 702 };
>> 703 static G4bool NeedInit = true;
>> 704
>> 705 static G4double A75[92];
>> 706
>> 707 if ( NeedInit )
>> 708 {
>> 709 for (G4int i=0; i<92; ++i)
>> 710 {
>> 711 A75[i] = std::pow(A[i], alpha); // interpolate by square ~ A^(2/3)
>> 712 }
>> 713 NeedInit=false;
>> 714 }
>> 715
>> 716 // for tabulated data, cross section scales with A^(2/3)
>> 717 G4double r1 = x1 / A75[Z1-1] * A75[Z-1];
>> 718 G4double r2 = x2 / A75[Z2-1] * A75[Z-1];
>> 719 G4double result = 0.5*(r1+r2);
>> 720
>> 721 // More precise average
>> 722 if(Z1 != Z2) {
>> 723 G4double alp1 = G4double(Z - Z1);
>> 724 G4double alp2 = G4double(Z2 - Z);
>> 725 result = (r1*alp2 + r2*alp1)/(alp1 + alp2);
>> 726 }
>> 727 // G4cout << "x1/2, z1/2 z" <<x1<<" "<<x2<<" "<<Z1<<" "<<Z2<<" "<<Z<<G4endl;
>> 728 // G4cout << "res1/2 " << r1 <<" " << r2 <<" " << result<< G4endl;
>> 729 return result;
>> 730 }
92 731