include/Geant4/G4BarashenkovData.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 #ifndef G4BarashenkovData_h
0027 #define G4BarashenkovData_h
0028
0029 //
0030 // author: Vladimir Ivanchenko 16 August 2018
0031 //
0032 //         on base of Vladimir Grichine compilation of cross section data:
0033 //         Barashenkov V.S., Nucleon-Nucleus Cross Section, Preprint
0034 //         JINR P2-89-770, p. 12, Dubna 1989 (scanned version from KEK)
0035 //
0036 //         Evaluation of proton and neutron x-sections
0037 //
0038
0039 // Group 1: He, Be, C for 44 energies
0040
0041 static const G4double e1[44] =
0042 {
0043   0.014, 0.015, 0.017, 0.02, 0.022, 0.025, 0.027, 0.03, 0.035, 0.04,
0044   0.045, 0.05,  0.06,  0.07, 0.08,  0.09,  0.1,   0.12, 0.14,  0.15,
0045   0.16,  0.18,  0.20,  0.25, 0.30,  0.35,  0.4,   0.5,  0.6,   0.7,
0046   0.8,   0.9,   1.0,   1.5,  2.0,   3.0,   5.0,   7.0, 10.0,  20.0,
0047  50.0, 100.0, 500.0, 1000.0
0048 };
0049
0050 static const G4double he_m_t[44] =
0051 {
0052   1090, 1020, 915, 800, 710, 640, 600, 560, 500, 440, 390, 360, 295, 256, 220, 192,
0053   168, 136, 120, 116, 114, 110, 107, 104, 106, 108, 110, 120, 126, 135, 140, 144, 146,
0054   148, 152, 150, 146, 142, 138, 132, 129, 126, 127, 128
0055 };
0056 static const G4double he_m_in[44] =
0057 {
0058   0, 5, 10, 20, 35, 55, 70, 80, 90, 105, 115, 115, 100, 90, 86, 84, 84, 82, 80, 80, 80, 80,
0059   79, 78, 80, 84, 88, 94, 100, 105, 108, 108, 108, 112, 114, 114, 112, 110, 108, 106, 104,
0060   101, 102, 102
0061 };
0062 static const G4double he_p_in[44] =
0063 {
0064   0, 2, 3, 13, 30, 50, 65, 77, 90, 105, 115, 115, 100, 90, 86, 84, 84, 82, 80, 80, 80, 80,
0065   79, 78, 80, 84, 88, 94, 100, 105, 108, 108, 108, 112, 114, 114, 112, 110, 108, 106, 104,
0066   101, 102, 102
0067 };
0068
0069 static const G4double be_m_t[44] =
0070 {
0071   1490, 1460, 1400, 1350, 1270, 1200, 1160, 1100, 1000, 910, 810, 740, 625, 575, 455, 406,
0072   365, 310, 275, 262, 255, 240, 235, 225, 225, 230, 238, 252, 270, 282, 288, 290, 294, 303,
0073   303, 300, 292, 284, 277, 267, 263, 264, 268, 268
0074 };
0075 static const G4double be_m_in[44] =
0076 {
0077   650, 640, 617, 595, 555, 520, 495, 470, 430, 385, 350, 320, 270, 250, 210, 190, 185, 178,
0078   175, 175, 175, 175, 175, 170, 170, 172, 176, 184, 194, 200, 209, 213, 214, 216, 216, 212,
0079   210, 210, 210, 210, 210, 210, 210, 210
0080 };
0081 static const G4double be_p_in[44] =
0082 {
0083   490, 540, 580, 545, 525, 495, 470, 450, 420, 370, 340, 310, 262, 242, 205, 185, 180, 175,
0084   172, 175, 175, 175, 175, 170, 170, 172, 176, 184, 194, 200, 209, 213, 214, 216, 216, 212,
0085   210, 210, 210, 210, 210, 210, 210, 210
0086 };
0087
0088 static const G4double c_m_t[44] =
0089 {
0090   1240, 1370, 1450, 1455, 1445, 1385, 1345, 1290, 1210, 1110, 1020, 940, 800, 700, 604, 530,
0091   475, 396, 350, 336, 320, 303, 294, 280, 280, 286, 296, 314, 330, 344, 356, 360, 364, 384,
0092   388, 384, 364, 352, 344, 330, 324, 324, 332, 332
0093 };
0094 static const G4double c_m_in[44] =
0095 {
0096   590, 570, 542, 510, 500, 460, 445, 430, 395, 380, 350, 330, 295, 270, 255, 240, 228, 222,
0097   216, 216, 210, 210, 210, 208, 210, 214, 216, 228, 240, 248, 254, 257, 260, 262, 260, 256,
0098   252, 252, 250, 250, 248, 248, 248, 248
0099 };
0100 static const G4double c_p_in[44] =
0101 {
0102   310, 330, 400, 440, 450, 435, 430, 420, 385, 370, 340, 320, 288, 263, 249, 234, 222, 216,
0103   210, 211, 205, 208, 210, 208, 210, 214, 216, 228, 240, 248, 254, 257, 260, 262, 260, 256,
0104   252, 252, 250, 250, 248, 248, 248, 248
0105 };
0106
0107 // Group 2: N, O, Na for 44 energies (e1=e2)
0108
0109 static const G4double e2[44] =
0110 {
0111  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,
0112   .08, 0.09,  .1, .12, .14, .15, .16, .18, .20, .25, .30, .35, .4 , 0.5, 0.6, 0.7,  0.8,
0113   0.9,   1, 1.5,   2,   3,   5,  7, 10,
0114   20,   50,  100,  500, 1000
0115 };
0116
0117 static const G4double n_m_t[44] =
0118 {
0119   1420,1480, 1537, 1550, 1525, 1500, 1480, 1425, 1340, 1260, 1175, 1090, 930, 805, 690, 612,
0120   552, 462, 402, 384, 372, 350, 345, 326, 324, 328, 336, 356, 372, 388, 400, 408, 415, 430,
0121   435, 432, 415, 402, 390, 375, 367, 370, 382, 385
0122 };
0123 static const G4double n_m_in[44] =
0124 {
0125   680, 665, 625, 580, 562, 525, 510, 485, 450, 435, 410, 387, 340, 310, 290, 280, 276, 274,
0126   260, 258, 254, 247, 245, 240, 240, 244, 250, 260, 268, 275, 280, 285, 290, 295, 300, 294,
0127   292, 290, 285, 285, 282, 282, 282, 282
0128 };
0129 static const G4double n_p_in[44] =
0130 {
0131   420, 440, 470, 490, 497, 500, 480, 462, 440, 425, 400, 377, 333, 303, 284, 274, 270, 268,
0132   254, 252, 247, 245, 245, 240, 240, 244, 250, 260, 268, 275, 280, 285, 290, 295, 300, 294,
0133   292, 290, 285, 285, 282, 282, 282, 282
0134 };
0135
0136 static const G4double o_m_t[44] =
0137 {
0138   1520, 1570, 1630, 1660, 1647, 1623, 1595, 1555, 1475, 1395, 1290, 1207, 1035, 925, 816,
0139   720, 645, 540, 462, 438, 415, 392, 378, 362, 361, 381, 390, 403, 417, 440, 460, 470,
0140   479, 498, 504, 498, 477, 457, 443, 427, 420, 425, 429, 430
0141 };
0142 static const G4double o_m_in[44] =
0143 {
0144   750, 740, 700, 650, 620, 575, 555, 530, 505, 462, 435, 420, 375, 345, 320, 310, 300, 293,
0145   288, 282, 282, 280, 276, 270, 271, 275, 280, 290, 295, 304, 310, 315, 318, 332, 335, 330,
0146   323, 320, 317, 315, 315, 315, 315, 315
0147 };
0148 static const G4double o_p_in[44] =
0149 {
0150   460, 485, 510, 535, 537, 532, 520, 500, 460, 432, 405, 390, 350, 320, 310, 304, 293, 287,
0151   283, 279, 279, 278, 276, 270, 271, 275, 280, 290, 295, 304, 310, 315, 318, 332, 335, 330,
0152   323, 320, 317, 315, 315, 315, 315, 315
0153 };
0154
0155 static const G4double na_m_t[44] =
0156 {
0157   1570, 1620, 1695, 1730, 1750, 1760, 1755, 1740, 1710, 1643, 1560, 1480, 1343, 1220, 1073,
0158   953, 860, 720, 618, 582, 546, 522, 504, 484, 492, 500, 512, 538, 560, 586, 608, 622, 632,
0159   660, 668, 664, 640, 616, 596, 568, 568, 568, 568, 568
0160 };
0161 static const G4double na_m_in[44] =
0162 {
0163   960, 930, 890, 822, 790, 750, 725, 686, 620, 600, 575, 540, 497, 450, 414, 390, 380, 372,
0164   354, 360, 355, 354, 350, 350, 350, 356, 364, 384, 392, 400, 408, 410, 420, 408, 412, 420,
0165   411, 409, 407, 403, 400, 400, 400, 400
0166 };
0167 static const G4double na_p_in[44] =
0168 {
0169   600, 617, 660, 675, 680, 680, 670, 650, 575, 550, 525, 490, 450, 420, 385, 367, 360, 350,
0170   350, 350, 345, 347, 350, 350, 350, 356, 364, 384, 392, 400, 408, 410, 420, 408, 412, 420,
0171   411, 409, 407, 403, 400, 400, 400, 400
0172 };
0173
0174 // Al, Si, Ca for 45 energies
0175
0176 static const G4double e3[45] =
0177 {
0178   0.014, 0.015, 0.016, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.035, .04, 0.045, 0.05, .06, 0.07,
0179     .08, 0.09,   .1,    .12,  .14,  .15,   .16,   .18,   .20,  .25,  .30,  .35,  0.4, 0.5,  0.6,
0180     0.7, 0.8,   0.9,    1,   1.5,   2,     3,     5,     7,   10,   20,   50,  100,   500, 1000
0181 };
0182
0183 static const G4double al_m_t[45] =
0184 {
0185   1735, 1750, 1760, 1795, 1830, 1855, 1885, 1895, 1900, 1870, 1835, 1785, 1710, 1522, 1350,
0186   1212, 1080,  972,  816,  720,  678,  642,  600,  567,  558,  560,  578,  592,  616,  644,
0187    672,  688,  708,  720,  736,  754,  736,  706,  680,  672,  646,  632,  632,  632,  632
0188 };
0189 static const G4double al_m_in[45] =
0190 {
0191   1000,  990,  975,  950,  905,  875,  825,  800,  762,  690,  652,  610,  570,  495,  480,
0192    456,  444,  432,  420,  420,  420,  420,  410,  410,  400,  402,  404,  408,  424,  438,
0193    448,  450,  454,  456,  472,  480,  466,  456,  452,  448,  444,  440,  440,  440,  440
0194 };
0195 static const G4double al_p_in[45] =
0196 {
0197    650,  682,  690,  715,  750,  762,  750,  740,  720,  655,  617,  575,  540,  470,  455,
0198    //   532,  420,  408,  400,  403,  403,  408,  406,  404,  400,  402,  404,  408,  424,  438,
0199    432,  420,  408,  400,  403,  403,  408,  406,  404,  400,  402,  404,  408,  424,  438,
0200    448,  450,  454,  456,  472,  480,  466,  456,  452,  448,  444,  440,  440,  440,  440
0201 };
0202
0203 static const G4double si_m_t[45] =
0204 {
0205   1810, 1833, 1850, 1872, 1920, 1950, 1995, 2020, 2035, 2000, 1930, 1850, 1760, 1570, 1400,
0206   1255, 1110, 1008,  846,  742,  696,  671,  623,  588,  584,  584,  602,  618,  645,  679,
0207    708,  727,  746,  757,  769,  782,  771,  734,  710,  698,  672,  654,  650,  650,  650
0208 };
0209 static const G4double si_m_in[45] =
0210 {
0211   1060, 1035, 1015,  990,  935,  900,  860,  830,  790,  725,  665,  630,  600,  520,  504,
0212    486,  470,  456,  444,  432,  432,  432,  418,  418,  415,  412,  416,  422,  440,  460,
0213    472,  476,  479,  480,  492,  496,  488,  472,  472,  464,  460,  452,  448,  448,  448
0214 };
0215 static const G4double si_p_in[45] =
0216 {
0217    670,  700,  725,  750,  780,  780,  770,  757,  735,  690,  635,  585,  570,  490,  475,
0218    460,  446,  431,  423,  425,  425,  425,  425,  422,  422,  412,  416,  422,  440,  460,
0219    472,  476,  479,  480,  492,  496,  488,  472,  472,  464,  460,  452,  448,  448,  448
0220 };
0221
0222 static const G4double ca_m_t[45] =
0223 {
0224   2180, 2130, 2095, 2075, 2115, 2150, 2220, 2250, 2300, 2365, 2360, 2280, 2180, 2000,
0225   1805, 1650, 1500, 1340, 1140, 990, 940, 890, 825, 790, 770, 773, 787, 800, 830, 870,
0226   905, 930, 950, 965, 990, 1002, 990, 965, 945, 925, 892, 860, 860, 860, 860
0227 };
0228 static const G4double ca_m_in[45] =
0229 {
0230   1240, 1225, 1200, 1180, 1125, 1090, 1045, 1020, 980, 925, 880, 825, 770, 680, 640,
0231   620, 615, 600, 580, 565, 560, 560, 560, 550, 535, 530, 540, 550, 570, 595, 610, 615,
0232   620, 622, 629, 630, 620, 612, 607, 592, 587, 580, 580, 580, 580
0233 };
0234 static const G4double ca_p_in[45] =
0235 {
0236   770, 800, 823, 850, 900, 925, 935, 920, 895, 835, 800, 750, 715, 640, 605, 590, 588,
0237   573, 555, 543, 540, 540, 540, 535, 530, 530, 540, 550, 570, 595, 610, 615,
0238   620, 622, 629, 630, 620, 612, 607, 592, 587, 580, 580, 580, 580
0239 };
0240
0241 // Fe, Cu, Mo for 47 energies
0242
0243 static const G4double e4[47] =
0244 {
0245   0.014, 0.015, 0.017, .02, 0.022, 0.025, 0.027, 0.03, 0.033, 0.035, 0.037, .04, 0.045,
0246   0.05, 0.055, .06, 0.07, .08, 0.09,  .1, .12, .14, .15, .16, .18, .20, .25, .30, .35,
0247   .4 , 0.5, 0.6, 0.7,  0.8,  0.9,   1, 1.5,   2,   3,   5,  7, 10,
0248   20,   50,  100,  500, 1000
0249 };
0250
0251 static const G4double fe_m_t[47] =
0252 {
0253   2580, 2490, 2370, 2282, 2275, 2285, 2320, 2370, 2432, 2445, 2460, 2485, 2530, 2540,
0254   2517, 2480, 2290, 2110, 1940, 1790, 1510, 1290, 1220, 1150, 1070, 1030, 1013, 1020,
0255   1030, 1043, 1075, 1110, 1133, 1163, 1185, 1225, 1252, 1260, 1260, 1233, 1207, 1185,
0256   1140, 1110, 1110, 1110, 1110
0257 };
0258 static const G4double fe_m_in[47] =
0259 {
0260   1440, 1433, 1390, 1325, 1280, 1260, 1215, 1180, 1140, 1110, 1080, 1040, 990, 955, 920,
0261   885, 835, 800, 780, 765, 750, 725, 720, 720, 710, 700, 700, 700, 712, 705, 735, 750,
0262   765, 775, 780, 795, 810, 813, 810, 784, 757, 743, 735, 720, 720, 720, 720
0263 };
0264 static const G4double fe_p_in[47] =
0265 {
0266   900, 960, 1070, 1090, 1115, 1120, 1115, 1080, 1045, 1025, 1000, 960, 900, 885, 865, 790,
0267   765, 740, 720, 700, 697, 697, 697, 697, 695, 690, 688, 690, 712, 705, 735, 750,
0268   765, 775, 780, 795, 810, 813, 810, 784, 757, 743, 735, 720, 720, 720, 720
0269 };
0270
0271 static const G4double cu_m_t[47] =
0272 {
0273   2920, 2800, 2615, 2480, 2455, 2430, 2440, 2460, 2500, 2530, 2560, 2615, 2690, 2720,
0274   2700, 2645, 2500, 2320, 2140, 1970, 1670, 1460, 1380, 1285, 1200, 1160, 1140, 1147,
0275   1163, 1170, 1200, 1237, 1265, 1285, 1305, 1328, 1375, 1390, 1395, 1370, 1335, 1315,
0276   1270, 1230, 1230, 1230, 1230
0277 };
0278 static const G4double cu_m_in[47] =
0279 {
0280   1540, 1535, 1500, 1445, 1407, 1380, 1330, 1300, 1285, 1270, 1240, 1190, 1090, 1010,
0281   940, 920, 860, 835, 820, 810, 800, 780, 775, 770, 760, 760, 758, 765, 765, 770, 795,
0282   810, 825, 830, 840, 848, 870, 870, 868, 840, 825, 810, 803, 795, 795, 795, 795
0283 };
0284 static const G4double cu_p_in[47] =
0285 {
0286   935, 1000, 1060, 1190, 1220, 1250, 1240, 1210, 1150, 1130, 1115, 1050, 985, 950, 890,
0287   870, 820, 800, 785, 780, 770, 750, 745, 740, 735, 735, 745, 760, 762, 770, 795,
0288   810, 825, 830, 840, 848, 870, 870, 868, 840, 825, 810, 803, 795, 795, 795, 795
0289 };
0290
0291 static const G4double mo_m_t[47] =
0292 {
0293   4150, 4040, 3800, 3490, 3300, 3060, 2960, 2845, 2785, 2820, 2850, 2980, 3170, 3230,
0294   3270, 3280, 3225, 3075, 2895, 2710, 2355, 2060, 1925, 1800, 1630, 1560, 1540, 1550,
0295   1570, 1590, 1650, 1685, 1715, 1740, 1760, 1780, 1850, 1880, 1858, 1815, 1790, 1782,
0296   1720, 1690, 1690, 1690, 1690
0297 };
0298 static const G4double mo_m_in[47] =
0299 {
0300   1790, 1775, 1740, 1680, 1640, 1580, 1550, 1510, 1460, 1440, 1418, 1380, 1330, 1280,
0301   1240, 1200, 1155, 1140, 1110, 1110, 1080, 1065, 1050, 1050, 1025, 1020, 1015, 1020,
0302   1022, 1026, 1060, 1085, 1100, 1110, 1120, 1127, 1150, 1160, 1140, 1100, 1085, 1080,
0303   1070, 1070, 1070, 1070, 1070
0304 };
0305 static const G4double mo_p_in[47] =
0306 {
0307   1025, 1080, 1190, 1380, 1440, 1495, 1475, 1420, 1350, 1310, 1300, 1290, 1250, 1200,
0308   1170, 1130, 1095, 1060, 1040, 1022, 1020, 1016, 1016, 1016, 1016, 1012, 1005, 1005,
0309   1005, 1010, 1060, 1085, 1100, 1110, 1120, 1127, 1150, 1160, 1140, 1100, 1085, 1080,
0310   1070, 1070, 1070, 1070, 1070
0311 };
0312
0313 // Cd, Sn, W for 48 energies
0314
0315 static const G4double e5[48] =
0316 {
0317    0.014, 0.015, 0.017, 0.018, .02, 0.022, 0.025, 0.027, 0.03, 0.033, 0.035, .04,
0318    0.045, 0.05, 0.055, .06, .065, 0.07, .08, 0.09,  .1, .12, .14, .15, .16, .18,
0319    .20, .25, .30, .35, .4 , 0.5, 0.6, 0.7,  0.8, 0.9, 1, 1.5, 2, 3, 5, 7, 10,  20,
0320    50, 100, 500, 1000
0321 };
0322
0323 static const G4double cd_m_t[48] =
0324 {
0325   4420, 4280, 4170, 4070, 3860, 3680, 3420, 3280, 3125, 3060, 3080, 3190, 3350, 3445,
0326   3510, 3540, 3560, 3550, 3460, 3300, 3030, 2640, 2340, 2190, 2070, 1950, 1770, 1732,
0327   1740, 1760, 1780, 1832, 1885, 1925, 1945, 1960, 1980, 2070, 2080, 2065, 2040, 2022,
0328   1980, 1940, 1870, 1870, 1870, 1870
0329 };
0330 static const G4double cd_m_in[48]=
0331 {
0332   1920, 1910, 1880, 1860, 1840, 1800, 1760, 1720, 1675, 1630, 1600, 1520, 1465, 1420,
0333   1390, 1340, 1310, 1280, 1275, 1235, 1225, 1200, 1170, 1170, 1170, 1165, 1145, 1140,
0334   1140, 1135, 1160, 1180, 1220, 1240, 1250, 1260, 1265, 1270, 1275, 1250, 1222, 1222,
0335   1220, 1215, 1190, 1190, 1190, 1190
0336 };
0337 static const G4double cd_p_in[48] =
0338 {
0339   1020, 1100, 1225, 1290, 1440, 1520, 1575, 1560, 1518, 1460, 1420, 1400, 1365, 1340,
0340   1300, 1280, 1260, 1200, 1190, 1160, 1125, 1125, 1125, 1125, 1125, 1125, 1120, 1120,
0341   1120, 1118, 1146, 1180, 1220, 1240, 1250, 1260, 1265, 1270, 1275, 1250, 1222, 1222,
0342   1220, 1215, 1190, 1190, 1190, 1190
0343 };
0344
0345 static const G4double sn_m_t[48] =
0346 {
0347   4420, 4400, 4260, 4150, 3980, 3770, 3530, 3370, 3245, 3180, 3170, 3260, 3400, 3500,
0348   3560, 3610, 3650, 3680, 3580, 3390, 3190, 2760, 2430, 2295, 2175, 1990, 1880, 1810,
0349   1820, 1840, 1865, 1940, 1985, 2020, 2040, 2060, 2080, 2160, 2185, 2180, 2110, 2105,
0350   2080, 2050, 1980, 1980, 1980, 1980
0351 };
0352 static const G4double sn_m_in[48] =
0353 {
0354   1945, 1940, 1905, 1890, 1860, 1830, 1780, 1755, 1717, 1680, 1645, 1570, 1500, 1455,
0355   1410, 1370, 1340, 1320, 1290, 1285, 1260, 1240, 1235, 1212, 1200, 1200, 1200, 1190,
0356   1190, 1200, 1210, 1240, 1270, 1285, 1300, 1300, 1310, 1320, 1320, 1290, 1240, 1240,
0357   1240, 1240, 1240, 1240, 1240, 1240
0358 };
0359 static const G4double sn_p_in[48] =
0360 {
0361   1020, 1080, 1270, 1335, 1465, 1505, 1610, 1610, 1550, 1535, 1500, 1440, 1407, 1370,
0362   1340, 1300, 1285, 1260, 1230, 1215, 1200, 1180, 1170, 1170, 1165, 1165, 1170, 1165,
0363   1165, 1183, 1195, 1240, 1270, 1285, 1300, 1300, 1310, 1320, 1320, 1290, 1240, 1240,
0364   1240, 1240, 1240, 1240, 1240, 1240
0365 };
0366
0367 static const G4double w_m_t[48] =
0368 {
0369   5320, 5430, 5480, 5450, 5330, 5190, 4960, 4790, 4550, 4340, 4200, 4070, 4000, 4030,
0370   4125, 4220, 4270, 4390, 4440, 4360, 4200, 3800, 3380, 3200, 3040, 2790, 2660, 2575,
0371   2575, 2600, 2640, 2690, 2755, 2790, 2812, 2837, 2850, 2950, 3000, 2970, 2940, 2910,
0372   2880, 2820, 2730, 2730, 2730, 2730
0373 };
0374 static const G4double w_m_in[48] =
0375 {
0376   2440, 2400, 2370, 2350, 2310, 2270, 2220, 2195, 2150, 2100, 2070, 2010, 1945, 1900,
0377   1850, 1820, 1780, 1760, 1730, 1720, 1680, 1680, 1660, 1660, 1650, 1650, 1640, 1640,
0378   1612, 1615, 1625, 1640, 1700, 1720, 1730, 1740, 1750, 1780, 1780, 1750, 1740, 1735,
0379   1710, 1695, 1680, 1680, 1680, 1680
0380 };
0381 static const G4double w_p_in[48] =
0382 {
0383   950,  1020, 1240, 1400, 1560, 1670, 1760, 1830, 1850, 1855, 1870, 1840, 1800, 1770,
0384   1740, 1715, 1680, 1670, 1650, 1620, 1610, 1600, 1600, 1600, 1600, 1600, 1600, 1595,
0385   1585, 1595, 1615, 1640, 1700, 1720, 1730, 1740, 1750, 1780, 1780, 1750, 1740, 1735,
0386   1710, 1695, 1680, 1680, 1680, 1680
0387 };
0388
0389 // Pb, U for 46 energies
0390
0391 static const G4double e6[46] =
0392 {
0393   0.014, 0.015, 0.017, 0.019, 0.02, 0.022, 0.025, 0.027, 0.03, 0.035,
0394   0.04,  0.045, 0.05,  0.055, 0.06, 0.07,  0.08,  0.09,  0.1,  0.12,
0395   0.14,  0.15,  0.16,  0.18,  0.20, 0.25,  0.30,  0.35,  0.4 , 0.5,
0396   0.6,   0.7,   0.8,   0.9,   1.0,  1.5,   2.0,   3.0,   5.0,  7.0,
0397  10.0,  20.0,  50.0, 100.0, 500.0, 1000.0
0398 };
0399
0400 static const G4double pb_m_t[46] =
0401 {
0402   5300, 5440, 5720, 5880, 5765, 5745, 5480, 5280, 4970, 4550, 4390, 4300, 4265, 4325,
0403   4450, 4540, 4740, 4710, 4600, 4100, 3660, 3480, 3300, 3000, 2890, 2865, 2855, 2850,
0404   2865, 2920, 2955, 3000, 3030, 3060, 3105, 3240, 3290, 3270, 3240, 3180, 3090, 3060,
0405   2970, 2970, 2970, 2970
0406
0407 };
0408 static const G4double pb_m_in[46] =
0409 {
0410   2580, 2550, 2505, 2462, 2460, 2435, 2380, 2355, 2280, 2180, 2170, 2130, 2080, 2035,
0411   1980, 1940, 1900, 1870, 1840, 1800, 1800, 1800, 1780, 1760, 1760, 1740, 1730, 1725,
0412   1740, 1785, 1815, 1835, 1860, 1890, 1895, 1920, 1920, 1890, 1850, 1835, 1830, 1830,
0413   1830, 1830, 1830, 1830
0414 };
0415 static const G4double pb_p_in[46] =
0416 {
0417   900,  1060, 1200, 1420, 1515, 1620, 1750, 1800, 1915, 2030, 1960, 1940, 1910, 1860,
0418   1840, 1780, 1770, 1760, 1740, 1720, 1725, 1740, 1740, 1730, 1720, 1700, 1710, 1720,
0419   1730, 1740, 1815, 1835, 1860, 1890, 1895, 1920, 1920, 1890, 1850, 1835, 1830, 1830,
0420   1830, 1830, 1830, 1830
0421 };
0422
0423 static const G4double u_m_t[46] =
0424 {
0425   5800, 5940, 6160, 6345, 6360, 6350, 6170, 6020, 5760, 5350, 4990, 4800, 4710, 4690,
0426   4760, 5040, 5190, 5200, 5080, 4600, 4120, 3920, 3720, 3420, 3240, 3150, 3160, 3180,
0427   3210, 3240, 3280, 3350, 3390, 3435, 3480, 3560, 3585, 3580, 3540, 3500, 3470, 3410,
0428   3335, 3335, 3335, 3335
0429 };
0430 static const G4double u_m_in[46] =
0431 {
0432   2820, 2770, 2700, 2660, 2645, 2620, 2580, 2550, 2515, 2450, 2390, 2320, 2260, 2225,
0433   2200, 2140, 2080, 2060, 2040, 2000, 1980, 1965, 1960, 1930, 1920, 1890, 1905, 1920,
0434   1945, 1970, 1985, 2010, 2040, 2070, 2080, 2090, 2095, 2080, 2063, 2060, 2050, 2040,
0435   2005, 2005, 2005, 2005
0436 };
0437 static const G4double u_p_in[46] =
0438 {
0439   800,  900,  1100, 1300, 1410, 1510, 1680, 1800, 2000, 2200, 2080, 2060, 2035, 2100,
0440   2030, 2030, 2000, 1960, 1960, 1960, 1940, 1925, 1920, 1905, 1890, 1860, 1880, 1910,
0441   1930, 1945, 1985, 2010, 2040, 2070, 2080, 2090, 2095, 2080, 2063, 2060, 2050, 2040,
0442   2005, 2005, 2005, 2005
0443 };
0444
0445 #endif