| 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 // 28 //
29 // GEANT4 Class file 29 // GEANT4 Class file
30 // 30 //
31 // Description: Data on stopping power 31 // Description: Data on stopping power
32 // 32 //
33 // Author: Anton Ivantchenko 18.04.2006 33 // Author: Anton Ivantchenko 18.04.2006
34 // 34 //
35 // Organisation: QinetiQ Ltd, UK 35 // Organisation: QinetiQ Ltd, UK
36 // Customer: ESA/ESTEC, NOORDWIJK 36 // Customer: ESA/ESTEC, NOORDWIJK
37 // Contract: CSMAN-5288 37 // Contract: CSMAN-5288
38 // 38 //
39 // Modifications: 39 // Modifications:
40 // 25.04.2010 V.Ivanchenko introduce G4LPhysic 40 // 25.04.2010 V.Ivanchenko introduce G4LPhysicsFreeVector
41 // 26.04.2010 A.Ivanchenko fixed data for G4_P 41 // 26.04.2010 A.Ivanchenko fixed data for G4_PROPANE and G4_Pt,
42 // checked all data 42 // checked all data
43 // 43 //
44 //-------------------------------------------- 44 //----------------------------------------------------------------------------
45 // 45 //
46 46
47 //....oooOO0OOooo........oooOO0OOooo........oo 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
48 48
49 #include "G4PSTARStopping.hh" 49 #include "G4PSTARStopping.hh"
50 #include "G4NISTStoppingData.hh" 50 #include "G4NISTStoppingData.hh"
51 #include "G4EmParameters.hh" <<
52 51
53 //....oooOO0OOooo........oooOO0OOooo........oo 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54 53
55 54
56 G4PSTARStopping::G4PSTARStopping() : nvectors( 55 G4PSTARStopping::G4PSTARStopping() : nvectors(0), emin(CLHEP::keV)
57 { 56 {
58 // 1st initialisation 57 // 1st initialisation
59 Initialise(); 58 Initialise();
60 } 59 }
61 60
62 //....oooOO0OOooo........oooOO0OOooo........oo 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63 62
64 G4PSTARStopping::~G4PSTARStopping() 63 G4PSTARStopping::~G4PSTARStopping()
65 { 64 {
66 if(0 < nvectors) { 65 if(0 < nvectors) {
67 for(G4int i=0; i<nvectors; ++i) { delete s 66 for(G4int i=0; i<nvectors; ++i) { delete sdata[i]; }
68 } 67 }
69 } 68 }
70 69
71 //....oooOO0OOooo........oooOO0OOooo........oo 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72 71
73 void G4PSTARStopping::PrintWarning(G4int i) co 72 void G4PSTARStopping::PrintWarning(G4int i) const
74 { 73 {
75 G4ExceptionDescription ed; 74 G4ExceptionDescription ed;
76 ed << "index of data " << i << " is <0 or >= 75 ed << "index of data " << i << " is <0 or >= " << nvectors
77 << " request ignored!"; 76 << " request ignored!";
78 G4Exception("G4PSTARStopping::PrintWarning() 77 G4Exception("G4PSTARStopping::PrintWarning()","em0033",JustWarning,ed);
79 } 78 }
80 79
81 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
82 81
83 void G4PSTARStopping::Initialise() 82 void G4PSTARStopping::Initialise()
84 { 83 {
85 // this method may be called several times d 84 // this method may be called several times during initialisation
86 G4int nmat = (G4int)G4Material::GetNumberOfM << 85 G4int nmat = G4Material::GetNumberOfMaterials();
87 if(nmat == nvectors) { return; } 86 if(nmat == nvectors) { return; }
88 87
89 // loop via material list to add extra data 88 // loop via material list to add extra data
90 G4int j; 89 G4int j;
91 for(G4int i=0; i<nmat; ++i) { 90 for(G4int i=0; i<nmat; ++i) {
92 const G4Material* mat = (*(G4Material::Get 91 const G4Material* mat = (*(G4Material::GetMaterialTable()))[i];
93 92
94 G4bool isThere = false; 93 G4bool isThere = false;
95 for(j=0; j<nvectors; ++j) { 94 for(j=0; j<nvectors; ++j) {
96 if(mat == materials[j]) { 95 if(mat == materials[j]) {
97 isThere = true; 96 isThere = true;
98 break; 97 break;
99 } 98 }
100 } 99 }
101 if(!isThere) { 100 if(!isThere) {
102 // check list of NIST materials 101 // check list of NIST materials
103 G4String mname = mat->GetName(); 102 G4String mname = mat->GetName();
104 for(j=0; j<74; ++j) { 103 for(j=0; j<74; ++j) {
105 if(mname == nameNIST[j]) { 104 if(mname == nameNIST[j]) {
106 FindData(j, mat); 105 FindData(j, mat);
107 isThere = true; 106 isThere = true;
108 break; 107 break;
109 } 108 }
110 } 109 }
111 if(!isThere) { 110 if(!isThere) {
112 G4String form = mat->GetChemicalFormula(); 111 G4String form = mat->GetChemicalFormula();
113 for(j=0; j<numberOfMolecula; ++j) { 112 for(j=0; j<numberOfMolecula; ++j) {
114 if(form == molecularName[j]) { 113 if(form == molecularName[j]) {
115 FindData(molecularIndex[j], mat); 114 FindData(molecularIndex[j], mat);
116 break; 115 break;
117 } 116 }
118 } 117 }
119 } 118 }
120 } 119 }
121 } 120 }
122 } 121 }
123 122
124 //....oooOO0OOooo........oooOO0OOooo........oo 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125 124
126 void G4PSTARStopping::FindData(G4int j, const 125 void G4PSTARStopping::FindData(G4int j, const G4Material* mat)
127 { 126 {
128 static const G4float e0[60] = { 176.8f, 216. 127 static const G4float e0[60] = { 176.8f, 216.6f, 250.1f, 279.6f, 306.3f, 351.1f, 390.9f, 426.8f, 459.8f, 489.9f, 517.8f, 544.f, 596.4f, 641.3f, 680.6f, 715.7f, 747.f, 775.f, 800.f, 822.6f, 861.7f, 894.2f, 920.7f, 941.5f, 957.3f, 968.6f, 976.2f, 980.5f, 982.f, 981.f, 978.f, 973.2f, 967.f, 959.5f, 910.4f, 853.6f, 797.9f, 746.9f, 701.5f, 661.5f, 626.4f, 595.6f, 544.8f, 504.4f, 471.1f, 442.8f, 416.4f, 392.f, 369.9f, 350.f, 332.2f, 316.4f, 302.2f, 289.5f, 278.1f, 267.9f, 228.9f, 200.8f, 179.5f, 162.6f };
129 128
130 static const G4float e1[60] = { 185.5f, 227. 129 static const G4float e1[60] = { 185.5f, 227.2f, 262.3f, 293.2f, 321.2f, 364.8f, 404.f, 439.3f, 471.6f, 500.3f, 526.4f, 550.6f, 600.2f, 641.6f, 677.4f, 709.2f, 737.3f, 762.f, 783.7f, 803.1f, 837.2f, 866.8f, 891.7f, 911.6f, 926.6f, 937.4f, 944.7f, 949.f, 950.8f, 950.3f, 947.9f, 943.9f, 938.4f, 931.8f, 887.1f, 833.7f, 779.9f, 729.5f, 686.8f, 650.4f, 618.5f, 590.2f, 541.4f, 500.4f, 465.3f, 434.8f, 408.f, 384.6f, 364.4f, 346.5f, 330.5f, 316.1f, 303.2f, 291.3f, 280.4f, 270.4f, 230.7f, 202.2f, 180.5f, 163.5f };
131 130
132 static const G4float e2[60] = { 195.9f, 239. 131 static const G4float e2[60] = { 195.9f, 239.9f, 277.f, 309.7f, 339.3f, 387.5f, 430.6f, 469.4f, 505.f, 537.2f, 566.9f, 594.6f, 649.9f, 696.7f, 737.4f, 773.3f, 805.2f, 833.4f, 858.3f, 880.5f, 918.9f, 951.f, 977.2f, 997.6f, 1013.f, 1023.f, 1030.f, 1034.f, 1034.f, 1033.f, 1029.f, 1023.f, 1016.f, 1008.f, 955.6f, 895.f, 835.1f, 779.6f, 729.7f, 685.6f, 646.6f, 612.3f, 555.f, 509.1f, 471.6f, 440.3f, 413.7f, 390.4f, 370.1f, 352.1f, 336.f, 321.4f, 308.3f, 296.3f, 285.3f, 275.2f, 234.9f, 206.f, 184.f, 166.7f };
133 132
134 static const G4float e3[60] = { 29.73f, 36.4 133 static const G4float e3[60] = { 29.73f, 36.41f, 42.05f, 47.01f, 51.5f, 59.46f, 66.48f, 72.83f, 78.66f, 84.09f, 89.2f, 94.02f, 103.7f, 112.1f, 119.7f, 126.6f, 132.8f, 138.5f, 143.7f, 148.5f, 157.1f, 164.3f, 170.5f, 175.7f, 180.1f, 183.7f, 186.7f, 189.1f, 191.1f, 192.5f, 193.6f, 194.4f, 194.8f, 195.f, 193.4f, 189.3f, 184.4f, 179.2f, 174.1f, 169.3f, 164.8f, 160.5f, 152.7f, 145.7f, 139.5f, 133.8f, 128.6f, 123.9f, 119.5f, 115.4f, 111.7f, 108.1f, 104.8f, 101.7f, 98.74f, 95.98f, 84.49f, 75.92f, 69.3f, 64.05f };
135 134
136 static const G4float e4[60] = { 119.7f, 146. 135 static const G4float e4[60] = { 119.7f, 146.7f, 169.3f, 189.3f, 207.4f, 239.5f, 267.8f, 293.3f, 316.8f, 338.7f, 359.2f, 378.7f, 417.f, 450.4f, 480.1f, 506.7f, 530.7f, 552.6f, 572.4f, 590.5f, 622.1f, 648.3f, 670.f, 687.7f, 702.f, 713.2f, 721.7f, 727.8f, 731.9f, 734.1f, 734.8f, 734.f, 732.f, 729.f, 702.9f, 667.2f, 629.1f, 592.2f, 558.3f, 527.8f, 500.6f, 476.3f, 434.9f, 401.2f, 373.3f, 349.8f, 329.7f, 312.1f, 296.4f, 282.4f, 269.9f, 258.7f, 248.5f, 239.1f, 230.6f, 222.7f, 191.1f, 168.2f, 150.8f, 137.f };
137 136
138 static const G4float e5[60] = { 92.38f, 113. 137 static const G4float e5[60] = { 92.38f, 113.1f, 130.6f, 146.1f, 160.f, 184.8f, 206.6f, 226.3f, 244.4f, 261.3f, 277.1f, 292.1f, 320.6f, 344.8f, 365.7f, 383.8f, 399.6f, 413.2f, 425.f, 435.1f, 451.f, 462.f, 469.2f, 473.4f, 475.2f, 475.1f, 473.7f, 471.2f, 468.f, 464.2f, 460.1f, 455.8f, 451.3f, 446.8f, 424.5f, 404.5f, 386.7f, 371.f, 356.8f, 344.f, 332.3f, 321.5f, 301.7f, 284.2f, 268.6f, 254.8f, 242.5f, 231.4f, 221.5f, 212.4f, 204.2f, 196.6f, 189.7f, 183.3f, 177.4f, 171.9f, 149.4f, 132.7f, 119.8f, 109.4f };
139 138
140 static const G4float e6[60] = { 73.51f, 90.0 139 static const G4float e6[60] = { 73.51f, 90.03f, 104.f, 116.2f, 127.3f, 147.f, 164.4f, 180.1f, 194.5f, 207.9f, 220.5f, 232.5f, 256.1f, 276.8f, 295.3f, 311.9f, 327.f, 340.8f, 353.4f, 364.9f, 385.3f, 402.5f, 417.f, 429.2f, 439.3f, 447.6f, 454.3f, 459.6f, 463.7f, 466.6f, 468.6f, 469.8f, 470.2f, 470.f, 461.9f, 447.2f, 430.1f, 412.7f, 396.5f, 381.4f, 367.5f, 354.5f, 331.1f, 310.6f, 292.6f, 276.5f, 262.2f, 249.4f, 238.1f, 227.9f, 218.7f, 210.2f, 202.6f, 195.7f, 189.2f, 183.2f, 158.9f, 140.9f, 127.f, 116.f };
141 140
142 static const G4float e7[60] = { 86.08f, 105. 141 static const G4float e7[60] = { 86.08f, 105.4f, 121.7f, 136.1f, 149.1f, 172.2f, 192.5f, 210.9f, 227.7f, 243.5f, 258.2f, 272.2f, 299.7f, 323.5f, 344.5f, 363.3f, 380.2f, 395.3f, 409.f, 421.4f, 442.5f, 459.4f, 472.8f, 483.1f, 490.7f, 496.f, 499.2f, 500.7f, 500.8f, 499.5f, 497.2f, 494.f, 490.f, 485.5f, 457.4f, 426.7f, 397.7f, 371.9f, 349.5f, 330.1f, 313.2f, 298.5f, 274.2f, 254.9f, 239.f, 225.6f, 214.4f, 204.7f, 196.1f, 188.4f, 181.3f, 174.9f, 168.9f, 163.4f, 158.2f, 153.3f, 133.f, 118.2f, 106.8f, 97.72f };
143 142
144 static const G4float e8[60] = { 14.76f, 18.0 143 static const G4float e8[60] = { 14.76f, 18.07f, 20.87f, 23.33f, 25.56f, 29.51f, 33.f, 36.15f, 39.04f, 41.74f, 44.27f, 46.66f, 51.51f, 55.8f, 59.66f, 63.18f, 66.43f, 69.43f, 72.23f, 74.84f, 79.58f, 83.77f, 87.5f, 90.81f, 93.77f, 96.4f, 98.74f, 100.8f, 102.6f, 104.3f, 105.7f, 106.9f, 108.f, 108.9f, 111.5f, 112.f, 111.2f, 109.6f, 107.6f, 105.4f, 103.2f, 101.f, 96.33f, 91.8f, 87.7f, 84.05f, 80.83f, 78.f, 75.47f, 73.18f, 71.11f, 69.21f, 67.46f, 65.83f, 64.32f, 62.91f, 56.98f, 52.39f, 48.68f, 45.58f };
145 144
146 static const G4float e9[60] = { 142.5f, 174. 145 static const G4float e9[60] = { 142.5f, 174.5f, 201.5f, 225.2f, 246.7f, 283.1f, 315.3f, 344.4f, 371.1f, 395.5f, 418.2f, 439.6f, 482.3f, 519.f, 551.3f, 580.1f, 606.f, 629.2f, 650.1f, 669.f, 701.9f, 729.5f, 752.2f, 770.4f, 784.5f, 795.1f, 802.6f, 807.4f, 810.f, 810.5f, 809.3f, 806.6f, 802.7f, 797.7f, 762.5f, 719.5f, 676.2f, 635.8f, 599.3f, 566.7f, 537.9f, 512.3f, 468.1f, 431.4f, 400.8f, 375.f, 353.f, 333.7f, 316.9f, 301.9f, 288.5f, 276.3f, 265.3f, 255.1f, 245.9f, 237.4f, 203.4f, 178.9f, 160.2f, 145.4f };
147 146
148 static const G4float e10[60] = { 149.7f, 183 147 static const G4float e10[60] = { 149.7f, 183.3f, 211.7f, 236.7f, 259.2f, 299.3f, 334.7f, 366.6f, 396.f, 423.3f, 449.f, 473.3f, 517.2f, 553.4f, 583.6f, 608.7f, 629.6f, 646.8f, 660.7f, 671.9f, 687.3f, 695.5f, 698.4f, 697.6f, 694.2f, 689.1f, 682.8f, 675.7f, 668.2f, 660.5f, 652.7f, 644.9f, 637.2f, 629.6f, 594.2f, 563.f, 535.4f, 510.7f, 488.4f, 468.2f, 449.8f, 432.8f, 404.4f, 380.1f, 358.1f, 337.6f, 318.8f, 302.2f, 287.4f, 274.2f, 262.3f, 251.5f, 241.7f, 232.7f, 224.4f, 216.8f, 186.f, 163.6f, 146.3f, 132.6f };
149 148
150 static const G4float e11[60] = { 141.6f, 173 149 static const G4float e11[60] = { 141.6f, 173.4f, 200.2f, 223.8f, 245.2f, 281.3f, 313.3f, 342.2f, 368.8f, 393.1f, 415.6f, 436.8f, 479.f, 515.2f, 547.f, 575.3f, 600.7f, 623.4f, 643.8f, 662.1f, 694.2f, 721.f, 743.2f, 760.9f, 774.8f, 785.2f, 792.6f, 797.6f, 800.3f, 801.1f, 800.3f, 797.9f, 794.4f, 789.9f, 756.4f, 714.5f, 671.4f, 630.5f, 593.2f, 559.7f, 530.f, 503.5f, 459.f, 423.f, 393.3f, 368.3f, 346.9f, 328.1f, 311.6f, 296.9f, 283.7f, 271.8f, 261.f, 251.1f, 242.f, 233.7f, 200.3f, 176.3f, 158.f, 143.5f };
151 150
152 static const G4float e12[60] = { 138.4f, 169 151 static const G4float e12[60] = { 138.4f, 169.5f, 195.7f, 218.8f, 239.7f, 270.2f, 297.8f, 322.7f, 345.4f, 365.1f, 382.6f, 398.7f, 433.5f, 462.4f, 487.4f, 509.9f, 530.f, 547.7f, 563.4f, 577.6f, 603.9f, 628.3f, 650.1f, 668.6f, 683.5f, 695.3f, 704.5f, 711.4f, 716.3f, 719.4f, 720.8f, 720.9f, 719.7f, 717.4f, 694.5f, 661.1f, 624.7f, 589.1f, 556.f, 526.f, 499.1f, 475.f, 435.8f, 404.4f, 377.9f, 354.6f, 334.f, 315.9f, 300.2f, 286.2f, 273.7f, 262.4f, 252.2f, 242.7f, 234.1f, 226.1f, 194.2f, 171.1f, 153.4f, 139.4f };
153 152
154 static const G4float e13[60] = { 142.6f, 175 153 static const G4float e13[60] = { 142.6f, 175.3f, 203.1f, 228.f, 250.5f, 283.4f, 313.3f, 340.4f, 365.1f, 386.7f, 405.9f, 423.5f, 462.f, 494.f, 521.9f, 546.9f, 569.3f, 589.2f, 606.9f, 622.9f, 652.4f, 679.9f, 704.5f, 725.3f, 742.2f, 755.6f, 766.1f, 774.f, 779.6f, 783.1f, 784.7f, 784.8f, 783.4f, 780.7f, 753.8f, 714.5f, 671.4f, 629.1f, 590.2f, 555.6f, 524.9f, 497.7f, 453.9f, 419.3f, 390.3f, 365.1f, 343.f, 323.7f, 307.1f, 292.3f, 279.1f, 267.3f, 256.6f, 246.8f, 237.8f, 229.5f, 196.7f, 173.f, 155.f, 140.8f };
155 154
156 static const G4float e14[60] = { 225.8f, 276 155 static const G4float e14[60] = { 225.8f, 276.5f, 319.3f, 357.f, 391.1f, 445.9f, 494.9f, 539.1f, 579.6f, 616.f, 649.5f, 680.6f, 742.9f, 795.1f, 840.1f, 879.9f, 914.8f, 945.5f, 972.4f, 996.1f, 1037.f, 1071.f, 1099.f, 1120.f, 1135.f, 1145.f, 1150.f, 1152.f, 1152.f, 1148.f, 1142.f, 1135.f, 1126.f, 1115.f, 1052.f, 981.5f, 912.8f, 849.7f, 795.6f, 749.f, 708.3f, 672.4f, 611.7f, 562.1f, 520.6f, 485.3f, 454.9f, 428.4f, 405.5f, 385.2f, 367.1f, 350.8f, 336.2f, 322.8f, 310.5f, 299.3f, 254.6f, 222.7f, 198.6f, 179.6f };
157 156
158 static const G4float e15[60] = { 119.2f, 146 157 static const G4float e15[60] = { 119.2f, 146.f, 168.6f, 188.5f, 206.5f, 235.1f, 260.7f, 283.8f, 305.f, 323.9f, 341.3f, 357.4f, 390.7f, 418.9f, 443.7f, 465.8f, 485.7f, 503.6f, 519.5f, 534.1f, 560.1f, 583.2f, 603.2f, 620.f, 633.7f, 644.7f, 653.4f, 660.1f, 665.f, 668.4f, 670.3f, 671.1f, 670.7f, 669.5f, 652.9f, 626.4f, 596.2f, 565.7f, 537.f, 510.7f, 486.8f, 465.1f, 427.5f, 396.1f, 369.6f, 347.f, 327.4f, 310.1f, 294.9f, 281.3f, 269.f, 258.f, 247.9f, 238.6f, 230.2f, 222.4f, 191.f, 168.2f, 150.8f, 137.f };
159 158
160 static const G4float e16[60] = { 106.4f, 130 159 static const G4float e16[60] = { 106.4f, 130.3f, 150.4f, 168.2f, 184.2f, 212.7f, 237.9f, 260.6f, 281.4f, 300.9f, 319.1f, 336.4f, 370.9f, 401.2f, 428.3f, 452.8f, 475.1f, 495.6f, 514.4f, 531.8f, 562.6f, 588.8f, 611.2f, 630.2f, 646.1f, 659.3f, 670.1f, 678.6f, 685.2f, 690.f, 693.2f, 694.9f, 695.4f, 694.8f, 679.f, 650.5f, 617.3f, 583.6f, 552.1f, 523.6f, 497.8f, 474.6f, 434.7f, 401.8f, 374.3f, 351.f, 330.9f, 313.2f, 297.6f, 283.7f, 271.2f, 259.9f, 249.6f, 240.3f, 231.7f, 223.7f, 191.9f, 168.9f, 151.4f, 137.5f };
161 160
162 static const G4float e17[60] = { 74.48f, 91. 161 static const G4float e17[60] = { 74.48f, 91.21f, 105.3f, 117.8f, 129.f, 149.f, 166.5f, 182.4f, 197.f, 210.6f, 223.4f, 235.5f, 259.6f, 280.7f, 299.6f, 316.6f, 332.f, 346.2f, 359.1f, 371.1f, 392.2f, 410.1f, 425.2f, 438.f, 448.7f, 457.5f, 464.6f, 470.3f, 474.7f, 477.9f, 480.1f, 481.4f, 481.9f, 481.7f, 473.f, 456.7f, 437.6f, 417.8f, 398.6f, 380.6f, 364.f, 348.7f, 322.9f, 301.7f, 283.5f, 267.5f, 253.4f, 240.9f, 229.8f, 219.9f, 210.8f, 202.5f, 194.8f, 187.7f, 181.3f, 175.5f, 151.6f, 134.4f, 121.1f, 110.6f };
163 162
164 static const G4float e18[60] = { 135.4f, 165 163 static const G4float e18[60] = { 135.4f, 165.8f, 191.5f, 214.1f, 234.5f, 270.8f, 302.8f, 331.7f, 358.2f, 383.f, 406.2f, 428.2f, 471.7f, 509.8f, 543.6f, 574.f, 601.5f, 626.5f, 649.3f, 670.1f, 706.4f, 736.5f, 761.2f, 781.4f, 797.4f, 809.8f, 819.f, 825.3f, 829.2f, 830.8f, 830.6f, 828.7f, 825.4f, 820.9f, 785.9f, 741.4f, 696.f, 653.6f, 618.4f, 588.8f, 562.8f, 539.6f, 499.3f, 464.8f, 434.6f, 407.9f, 384.2f, 363.4f, 344.9f, 328.4f, 313.6f, 300.3f, 288.2f, 277.2f, 267.f, 257.7f, 220.4f, 193.5f, 173.f, 156.9f };
165 164
166 static const G4float e19[60] = { 132.2f, 161 165 static const G4float e19[60] = { 132.2f, 161.9f, 187.f, 209.f, 229.f, 262.6f, 292.5f, 319.4f, 344.2f, 366.8f, 387.8f, 407.5f, 447.f, 481.f, 510.8f, 537.6f, 561.6f, 583.2f, 602.7f, 620.4f, 651.5f, 677.9f, 700.f, 718.2f, 732.7f, 744.1f, 752.7f, 759.f, 763.2f, 765.5f, 766.2f, 765.5f, 763.5f, 760.6f, 734.4f, 698.2f, 659.3f, 621.4f, 586.f, 553.8f, 524.9f, 498.9f, 452.4f, 413.9f, 382.9f, 358.f, 337.7f, 319.7f, 303.9f, 289.7f, 277.1f, 265.6f, 255.2f, 245.6f, 236.9f, 228.8f, 196.4f, 173.f, 155.f, 140.8f };
167 166
168 static const G4float e20[60] = { 123.5f, 151 167 static const G4float e20[60] = { 123.5f, 151.3f, 174.7f, 195.3f, 214.f, 246.1f, 274.5f, 300.2f, 323.8f, 345.6f, 365.9f, 385.f, 423.f, 455.7f, 484.7f, 510.5f, 533.8f, 554.8f, 573.7f, 590.9f, 620.8f, 645.9f, 666.5f, 683.3f, 696.6f, 706.8f, 714.4f, 719.7f, 723.f, 724.5f, 724.5f, 723.1f, 720.6f, 717.1f, 689.6f, 653.5f, 615.8f, 579.8f, 546.6f, 516.7f, 490.f, 466.2f, 426.1f, 393.6f, 366.6f, 343.7f, 324.f, 306.7f, 291.5f, 278.f, 265.8f, 254.7f, 244.7f, 235.5f, 227.1f, 219.4f, 188.4f, 166.1f, 149.f, 135.4f };
169 168
170 static const G4float e21[60] = { 34.67f, 42. 169 static const G4float e21[60] = { 34.67f, 42.46f, 49.03f, 54.82f, 60.05f, 69.34f, 77.52f, 84.92f, 91.73f, 98.06f, 104.f, 109.6f, 121.f, 130.9f, 139.9f, 148.f, 155.5f, 162.3f, 168.6f, 174.5f, 184.9f, 194.f, 201.8f, 208.5f, 214.2f, 219.f, 223.f, 226.3f, 228.9f, 231.f, 232.5f, 233.4f, 234.f, 234.1f, 230.4f, 221.9f, 211.4f, 200.6f, 190.1f, 180.6f, 171.9f, 164.1f, 151.1f, 140.6f, 132.1f, 125.1f, 119.1f, 114.f, 109.6f, 105.6f, 102.1f, 98.92f, 95.83f, 92.83f, 89.95f, 87.22f, 75.76f, 67.41f, 61.26f, 56.63f };
171 170
172 static const G4float e22[60] = { 34.9f, 42.7 171 static const G4float e22[60] = { 34.9f, 42.74f, 49.35f, 55.18f, 60.45f, 69.8f, 78.04f, 85.48f, 92.33f, 98.71f, 104.7f, 110.4f, 121.4f, 130.9f, 139.2f, 146.6f, 153.2f, 159.f, 164.3f, 168.9f, 176.9f, 183.2f, 188.4f, 192.5f, 195.9f, 198.7f, 201.f, 202.8f, 204.4f, 205.6f, 206.7f, 207.5f, 208.2f, 208.7f, 209.9f, 209.2f, 207.5f, 205.1f, 202.3f, 199.1f, 195.8f, 192.3f, 185.4f, 178.6f, 172.f, 165.8f, 159.6f, 153.6f, 148.f, 142.7f, 137.8f, 133.3f, 129.1f, 125.2f, 121.7f, 118.3f, 104.7f, 94.43f, 86.45f, 79.87f };
173 172
174 static const G4float e23[60] = { 37.81f, 46. 173 static const G4float e23[60] = { 37.81f, 46.31f, 53.47f, 59.78f, 65.49f, 75.62f, 84.54f, 92.61f, 100.f, 106.9f, 113.4f, 119.6f, 132.f, 143.f, 153.f, 162.f, 170.4f, 178.2f, 185.4f, 192.2f, 204.6f, 215.6f, 225.4f, 234.2f, 242.1f, 249.2f, 255.5f, 261.2f, 266.2f, 270.7f, 274.6f, 278.1f, 281.f, 283.6f, 290.9f, 291.5f, 287.6f, 281.f, 272.8f, 263.9f, 254.9f, 246.f, 229.5f, 214.8f, 202.1f, 191.1f, 181.6f, 173.4f, 166.2f, 159.7f, 153.8f, 148.5f, 143.6f, 139.1f, 135.f, 131.1f, 115.3f, 103.5f, 94.02f, 86.52f };
175 174
176 static const G4float e24[60] = { 135.1f, 165 175 static const G4float e24[60] = { 135.1f, 165.5f, 191.1f, 213.6f, 234.f, 270.2f, 302.1f, 330.9f, 357.4f, 382.1f, 405.3f, 427.2f, 470.7f, 508.7f, 542.4f, 572.8f, 600.3f, 625.3f, 648.f, 668.8f, 705.f, 735.1f, 759.9f, 780.1f, 796.1f, 808.5f, 817.7f, 824.1f, 827.9f, 829.6f, 829.4f, 827.5f, 824.2f, 819.7f, 784.9f, 740.5f, 695.1f, 652.8f, 617.8f, 588.4f, 562.8f, 539.9f, 499.9f, 465.6f, 435.5f, 408.8f, 385.f, 364.1f, 345.6f, 329.1f, 314.3f, 300.9f, 288.8f, 277.7f, 267.6f, 258.3f, 220.8f, 193.9f, 173.4f, 157.2f };
177 176
178 static const G4float e25[60] = { 25.71f, 31. 177 static const G4float e25[60] = { 25.71f, 31.49f, 36.37f, 40.66f, 44.54f, 51.43f, 57.5f, 62.99f, 68.03f, 72.73f, 77.14f, 81.32f, 89.78f, 97.24f, 104.f, 110.1f, 115.8f, 121.f, 125.8f, 130.4f, 138.5f, 145.7f, 152.1f, 157.7f, 162.6f, 166.9f, 170.7f, 173.9f, 176.7f, 179.f, 181.f, 182.5f, 183.8f, 184.7f, 185.6f, 182.1f, 176.3f, 169.4f, 162.2f, 155.2f, 148.7f, 142.6f, 132.f, 123.3f, 116.2f, 110.2f, 105.f, 100.6f, 96.67f, 93.2f, 90.08f, 87.27f, 84.72f, 82.38f, 80.24f, 78.26f, 70.14f, 63.87f, 58.87f, 54.76f };
179 178
180 static const G4float e26[60] = { 45.91f, 56. 179 static const G4float e26[60] = { 45.91f, 56.23f, 64.93f, 72.59f, 79.52f, 91.82f, 102.7f, 112.5f, 121.5f, 129.9f, 137.7f, 145.2f, 159.7f, 172.3f, 183.3f, 193.1f, 201.7f, 209.4f, 216.3f, 222.3f, 232.3f, 240.f, 245.6f, 249.6f, 252.3f, 253.9f, 254.5f, 254.5f, 253.9f, 252.9f, 251.5f, 250.f, 248.2f, 246.4f, 236.5f, 227.2f, 218.9f, 211.6f, 205.1f, 199.3f, 193.9f, 189.f, 180.1f, 172.3f, 165.2f, 158.8f, 152.6f, 146.6f, 140.9f, 135.7f, 130.8f, 126.3f, 122.1f, 118.3f, 114.7f, 111.3f, 97.72f, 87.77f, 80.18f, 74.18f };
181 180
182 static const G4float e27[60] = { 88.74f, 108 181 static const G4float e27[60] = { 88.74f, 108.7f, 125.5f, 140.3f, 153.7f, 177.5f, 198.4f, 217.4f, 234.8f, 251.f, 266.2f, 280.6f, 308.9f, 333.6f, 355.6f, 375.3f, 393.1f, 409.2f, 423.9f, 437.3f, 460.8f, 480.4f, 496.7f, 510.2f, 521.2f, 530.1f, 537.1f, 542.5f, 546.4f, 549.1f, 550.6f, 551.2f, 551.f, 550.1f, 537.6f, 518.1f, 496.1f, 473.9f, 451.2f, 429.2f, 408.8f, 390.1f, 357.9f, 331.5f, 309.8f, 291.7f, 276.2f, 262.4f, 250.1f, 239.1f, 229.2f, 220.2f, 212.f, 204.6f, 197.7f, 191.4f, 165.6f, 146.7f, 132.1f, 120.4f };
183 182
184 static const G4float e28[60] = { 746.5f, 914 183 static const G4float e28[60] = { 746.5f, 914.3f, 1056.f, 1180.f, 1293.f, 1493.f, 1669.f, 1829.f, 1975.f, 2112.f, 2240.f, 2361.f, 2586.f, 2777.f, 2942.f, 3084.f, 3208.f, 3316.f, 3409.f, 3490.f, 3618.f, 3710.f, 3771.f, 3808.f, 3823.f, 3822.f, 3807.f, 3781.f, 3745.f, 3702.f, 3654.f, 3601.f, 3545.f, 3487.f, 3183.f, 2890.f, 2629.f, 2402.f, 2204.f, 2031.f, 1882.f, 1753.f, 1544.f, 1384.f, 1258.f, 1159.f, 1077.f, 1007.f, 947.f, 894.2f, 847.5f, 805.9f, 768.6f, 734.9f, 704.3f, 676.4f, 566.7f, 489.7f, 432.5f, 388.1f };
185 184
186 static const G4float e29[60] = { 184.2f, 225 185 static const G4float e29[60] = { 184.2f, 225.6f, 260.6f, 291.3f, 319.1f, 368.5f, 412.f, 451.3f, 487.4f, 521.1f, 552.7f, 582.6f, 641.3f, 691.9f, 736.7f, 776.6f, 812.5f, 844.9f, 874.2f, 900.8f, 946.6f, 984.1f, 1014.f, 1039.f, 1058.f, 1072.f, 1083.f, 1089.f, 1093.f, 1094.f, 1093.f, 1090.f, 1085.f, 1079.f, 1032.f, 973.f, 911.2f, 852.1f, 797.5f, 748.3f, 704.5f, 665.8f, 600.4f, 547.7f, 504.5f, 468.3f, 437.7f, 411.3f, 388.2f, 368.f, 349.9f, 333.8f, 319.3f, 306.2f, 294.2f, 283.2f, 239.6f, 208.6f, 185.2f, 167.f };
187 186
188 static const G4float e30[60] = { 145.4f, 178 187 static const G4float e30[60] = { 145.4f, 178.1f, 205.6f, 229.9f, 251.8f, 286.2f, 317.1f, 344.9f, 370.4f, 393.1f, 413.7f, 432.9f, 472.6f, 506.1f, 535.2f, 561.3f, 584.6f, 605.3f, 623.8f, 640.4f, 670.2f, 696.5f, 719.1f, 737.7f, 752.5f, 763.9f, 772.5f, 778.6f, 782.6f, 784.6f, 784.9f, 783.8f, 781.4f, 778.f, 749.2f, 710.6f, 669.6f, 630.f, 593.5f, 560.4f, 530.8f, 504.4f, 459.6f, 423.3f, 393.3f, 368.2f, 346.8f, 328.f, 311.6f, 296.9f, 283.8f, 272.f, 261.3f, 251.5f, 242.4f, 234.1f, 200.8f, 176.8f, 158.4f, 143.8f };
189 188
190 static const G4float e31[60] = { 45.92f, 56. 189 static const G4float e31[60] = { 45.92f, 56.24f, 64.94f, 72.6f, 79.53f, 91.84f, 102.7f, 112.5f, 121.5f, 129.9f, 137.8f, 145.2f, 160.1f, 173.2f, 184.9f, 195.5f, 205.1f, 213.8f, 221.8f, 229.1f, 242.f, 252.7f, 261.6f, 268.8f, 274.6f, 279.f, 282.3f, 284.5f, 285.8f, 286.2f, 286.f, 285.1f, 283.8f, 281.9f, 268.6f, 252.2f, 236.2f, 221.7f, 209.f, 198.2f, 188.9f, 180.9f, 168.1f, 158.1f, 150.2f, 143.7f, 138.1f, 133.2f, 128.8f, 124.8f, 121.2f, 117.8f, 114.6f, 111.6f, 108.8f, 106.1f, 94.63f, 85.43f, 77.91f, 71.84f };
191 190
192 static const G4float e32[60] = { 111.f, 135. 191 static const G4float e32[60] = { 111.f, 135.9f, 157.f, 175.5f, 192.2f, 222.f, 248.2f, 271.9f, 293.6f, 313.9f, 333.f, 351.f, 386.5f, 417.4f, 444.8f, 469.4f, 491.6f, 511.7f, 530.1f, 546.7f, 575.8f, 599.9f, 619.8f, 636.f, 649.1f, 659.3f, 667.2f, 672.9f, 676.7f, 678.9f, 679.6f, 679.1f, 677.6f, 675.1f, 653.f, 622.4f, 589.5f, 557.3f, 525.5f, 495.6f, 468.4f, 444.1f, 403.1f, 370.6f, 344.6f, 323.5f, 305.9f, 290.3f, 276.4f, 263.9f, 252.7f, 242.4f, 233.2f, 224.7f, 216.9f, 209.6f, 180.4f, 159.2f, 142.8f, 129.9f };
193 192
194 static const G4float e33[60] = { 80.87f, 99. 193 static const G4float e33[60] = { 80.87f, 99.04f, 114.4f, 127.9f, 140.1f, 161.7f, 180.8f, 198.1f, 214.f, 228.7f, 242.6f, 255.7f, 281.9f, 304.9f, 325.4f, 344.f, 360.9f, 376.5f, 390.8f, 404.f, 427.6f, 447.8f, 465.3f, 480.3f, 493.1f, 504.1f, 513.3f, 521.1f, 527.4f, 532.5f, 536.5f, 539.6f, 541.7f, 543.f, 540.3f, 527.4f, 509.2f, 488.8f, 467.4f, 446.5f, 426.9f, 408.7f, 376.9f, 350.4f, 328.f, 309.1f, 292.7f, 278.1f, 265.1f, 253.3f, 242.5f, 232.8f, 223.9f, 215.6f, 208.2f, 201.2f, 173.1f, 152.5f, 136.8f, 124.4f };
195 194
196 static const G4float e34[60] = { 198.1f, 242 195 static const G4float e34[60] = { 198.1f, 242.6f, 280.2f, 313.3f, 343.1f, 391.9f, 435.4f, 474.6f, 510.6f, 543.1f, 573.f, 601.f, 656.7f, 703.7f, 744.4f, 780.4f, 812.2f, 840.3f, 865.f, 887.f, 925.f, 956.5f, 982.1f, 1002.f, 1016.f, 1026.f, 1032.f, 1035.f, 1036.f, 1033.f, 1029.f, 1023.f, 1016.f, 1007.f, 953.3f, 891.9f, 831.3f, 775.4f, 725.1f, 680.6f, 641.3f, 606.6f, 546.6f, 497.8f, 458.8f, 427.4f, 401.7f, 379.2f, 359.6f, 342.2f, 326.7f, 312.7f, 300.1f, 288.5f, 277.8f, 268.1f, 229.1f, 201.f, 179.7f, 162.8f };
197 196
198 static const G4float e35[60] = { 174.4f, 213 197 static const G4float e35[60] = { 174.4f, 213.6f, 246.7f, 275.8f, 302.1f, 345.f, 383.2f, 417.7f, 449.4f, 477.9f, 504.3f, 528.8f, 578.1f, 619.8f, 656.f, 688.2f, 716.7f, 742.f, 764.4f, 784.4f, 819.3f, 848.7f, 872.8f, 892.f, 906.4f, 916.8f, 923.7f, 927.8f, 929.4f, 928.8f, 926.4f, 922.3f, 916.9f, 910.4f, 866.6f, 814.5f, 762.f, 712.8f, 668.3f, 628.6f, 593.3f, 562.1f, 507.7f, 463.3f, 427.6f, 398.9f, 375.3f, 354.7f, 336.7f, 320.6f, 306.2f, 293.3f, 281.6f, 270.9f, 261.f, 252.f, 215.7f, 189.5f, 169.6f, 153.9f };
199 198
200 static const G4float e36[60] = { 272.8f, 334 199 static const G4float e36[60] = { 272.8f, 334.1f, 385.8f, 431.4f, 472.5f, 545.6f, 610.f, 668.2f, 721.8f, 771.6f, 818.4f, 862.7f, 949.2f, 1024.f, 1090.f, 1148.f, 1200.f, 1246.f, 1288.f, 1325.f, 1387.f, 1436.f, 1473.f, 1500.f, 1518.f, 1528.f, 1532.f, 1530.f, 1524.f, 1514.f, 1501.f, 1485.f, 1467.f, 1447.f, 1338.f, 1227.f, 1127.f, 1041.f, 967.3f, 904.6f, 850.4f, 803.2f, 724.6f, 661.8f, 610.5f, 567.7f, 531.4f, 499.8f, 472.5f, 448.3f, 426.7f, 407.4f, 390.f, 374.1f, 359.6f, 346.4f, 293.7f, 256.2f, 227.9f, 205.7f };
201 200
202 static const G4float e37[60] = { 40.18f, 49. 201 static const G4float e37[60] = { 40.18f, 49.21f, 56.82f, 63.53f, 69.59f, 80.35f, 89.84f, 98.41f, 106.3f, 113.6f, 120.5f, 127.f, 140.2f, 151.7f, 162.1f, 171.5f, 180.1f, 188.f, 195.4f, 202.1f, 214.3f, 224.8f, 233.8f, 241.6f, 248.3f, 253.9f, 258.7f, 262.6f, 265.7f, 268.2f, 270.f, 271.3f, 272.1f, 272.4f, 268.7f, 259.5f, 248.f, 235.8f, 224.f, 213.1f, 203.1f, 194.1f, 178.9f, 166.6f, 156.5f, 148.1f, 141.f, 134.9f, 129.5f, 124.7f, 120.4f, 116.4f, 112.6f, 109.1f, 105.8f, 102.6f, 89.8f, 80.34f, 73.15f, 67.51f };
203 202
204 static const G4float e38[60] = { 149.4f, 183 203 static const G4float e38[60] = { 149.4f, 183.f, 211.3f, 236.2f, 258.7f, 298.3f, 333.3f, 364.8f, 393.8f, 420.8f, 446.f, 469.8f, 516.6f, 557.2f, 593.1f, 625.2f, 654.f, 680.f, 703.5f, 724.8f, 761.7f, 792.f, 816.7f, 836.5f, 851.9f, 863.6f, 871.9f, 877.4f, 880.4f, 881.2f, 880.1f, 877.3f, 873.2f, 867.9f, 829.5f, 781.9f, 733.4f, 688.f, 647.9f, 612.8f, 581.8f, 554.2f, 507.1f, 468.4f, 435.9f, 408.3f, 384.4f, 363.5f, 344.9f, 328.4f, 313.7f, 300.3f, 288.2f, 277.2f, 267.f, 257.7f, 220.4f, 193.5f, 173.f, 156.9f };
205 204
206 static const G4float e39[60] = { 146.6f, 179 205 static const G4float e39[60] = { 146.6f, 179.5f, 207.3f, 231.8f, 253.9f, 292.9f, 327.3f, 358.4f, 386.9f, 413.5f, 438.4f, 461.9f, 508.2f, 548.4f, 584.f, 615.8f, 644.5f, 670.4f, 693.8f, 715.1f, 752.f, 782.3f, 807.1f, 827.f, 842.5f, 854.3f, 862.9f, 868.5f, 871.6f, 872.6f, 871.6f, 869.f, 865.f, 859.8f, 822.f, 774.9f, 726.9f, 682.f, 643.f, 609.2f, 579.4f, 552.9f, 507.4f, 469.6f, 437.7f, 410.1f, 386.1f, 365.1f, 346.5f, 329.9f, 315.f, 301.6f, 289.4f, 278.3f, 268.2f, 258.8f, 221.3f, 194.3f, 173.7f, 157.5f };
207 206
208 static const G4float e40[60] = { 141.3f, 173 207 static const G4float e40[60] = { 141.3f, 173.1f, 199.9f, 223.5f, 244.8f, 282.5f, 315.8f, 345.8f, 373.5f, 399.2f, 423.3f, 446.1f, 491.1f, 530.2f, 565.f, 596.1f, 624.2f, 649.7f, 672.8f, 693.8f, 730.4f, 760.6f, 785.3f, 805.3f, 821.1f, 833.1f, 841.9f, 847.9f, 851.4f, 852.6f, 852.f, 849.7f, 846.1f, 841.2f, 804.8f, 758.9f, 712.2f, 668.5f, 631.2f, 599.f, 570.8f, 545.6f, 502.2f, 465.8f, 434.6f, 407.5f, 383.7f, 362.8f, 344.4f, 327.9f, 313.1f, 299.8f, 287.7f, 276.7f, 266.6f, 257.3f, 220.f, 193.2f, 172.8f, 156.7f };
209 208
210 static const G4float e41[60] = { 141.7f, 173 209 static const G4float e41[60] = { 141.7f, 173.5f, 200.4f, 224.f, 245.4f, 283.2f, 316.5f, 346.6f, 374.3f, 400.f, 424.2f, 447.f, 492.1f, 531.3f, 566.1f, 597.3f, 625.4f, 650.9f, 674.f, 695.1f, 731.6f, 761.9f, 786.7f, 806.6f, 822.4f, 834.4f, 843.3f, 849.2f, 852.6f, 853.9f, 853.2f, 851.f, 847.3f, 842.4f, 805.8f, 759.8f, 713.f, 669.2f, 632.f, 600.1f, 572.1f, 547.1f, 504.f, 467.7f, 436.5f, 409.3f, 385.3f, 364.3f, 345.7f, 329.2f, 314.4f, 300.9f, 288.8f, 277.8f, 267.6f, 258.3f, 220.8f, 193.9f, 173.4f, 157.2f };
211 210
212 static const G4float e42[60] = { 126.5f, 155 211 static const G4float e42[60] = { 126.5f, 155.f, 179.f, 200.1f, 219.2f, 253.1f, 283.f, 310.f, 334.8f, 357.9f, 379.6f, 400.2f, 440.6f, 475.8f, 507.f, 535.f, 560.2f, 583.1f, 603.8f, 622.7f, 655.5f, 682.7f, 704.9f, 722.9f, 737.3f, 748.4f, 756.6f, 762.3f, 765.9f, 767.5f, 767.5f, 766.f, 763.3f, 759.4f, 729.3f, 689.8f, 648.6f, 609.1f, 572.9f, 540.4f, 511.6f, 486.2f, 443.4f, 408.8f, 380.1f, 356.f, 335.4f, 317.2f, 301.2f, 286.8f, 274.1f, 262.6f, 252.2f, 242.7f, 234.f, 225.9f, 193.8f, 170.5f, 152.8f, 138.8f };
213 212
214 static const G4float e43[60] = { 41.1f, 50.3 213 static const G4float e43[60] = { 41.1f, 50.34f, 58.12f, 64.98f, 71.19f, 82.2f, 91.9f, 100.7f, 108.7f, 116.2f, 123.3f, 130.f, 143.4f, 155.2f, 165.8f, 175.4f, 184.2f, 192.3f, 199.7f, 206.6f, 219.f, 229.6f, 238.9f, 246.8f, 253.5f, 259.2f, 264.f, 268.f, 271.1f, 273.6f, 275.4f, 276.7f, 277.4f, 277.7f, 273.7f, 264.2f, 252.3f, 239.7f, 227.6f, 216.2f, 205.9f, 196.6f, 180.7f, 167.8f, 157.3f, 148.5f, 141.1f, 134.7f, 129.1f, 124.1f, 119.7f, 115.7f, 111.9f, 108.4f, 105.1f, 101.9f, 88.79f, 78.99f, 71.48f, 65.54f };
215 214
216 static const G4float e44[60] = { 58.01f, 71. 215 static const G4float e44[60] = { 58.01f, 71.05f, 82.04f, 91.73f, 100.5f, 116.f, 129.7f, 142.1f, 153.5f, 164.1f, 174.f, 183.5f, 202.5f, 219.3f, 234.5f, 248.3f, 261.f, 272.8f, 283.8f, 294.f, 312.7f, 329.2f, 343.8f, 356.9f, 368.6f, 379.f, 388.3f, 396.5f, 403.8f, 410.2f, 415.8f, 420.7f, 424.8f, 428.4f, 438.1f, 437.9f, 431.3f, 420.8f, 408.2f, 394.6f, 380.7f, 367.1f, 339.f, 313.8f, 293.3f, 277.2f, 264.3f, 252.6f, 241.9f, 232.2f, 223.2f, 215.f, 207.4f, 200.4f, 193.8f, 187.7f, 163.f, 144.6f, 130.3f, 118.9f };
217 216
218 static const G4float e45[60] = { 191.f, 233. 217 static const G4float e45[60] = { 191.f, 233.9f, 270.1f, 302.f, 330.8f, 377.8f, 419.7f, 457.5f, 492.1f, 523.4f, 552.3f, 579.2f, 633.f, 678.5f, 717.9f, 752.9f, 783.8f, 811.1f, 835.3f, 856.8f, 894.2f, 925.4f, 951.f, 971.f, 985.8f, 996.3f, 1003.f, 1007.f, 1008.f, 1006.f, 1003.f, 997.7f, 991.2f, 983.5f, 933.5f, 875.4f, 817.5f, 763.6f, 715.f, 671.7f, 633.5f, 599.6f, 542.9f, 497.4f, 460.4f, 429.7f, 403.7f, 381.f, 361.2f, 343.6f, 327.9f, 313.8f, 301.f, 289.3f, 278.6f, 268.8f, 229.5f, 201.3f, 179.8f, 162.9f };
219 218
220 static const G4float e46[60] = { 99.46f, 121 219 static const G4float e46[60] = { 99.46f, 121.8f, 140.7f, 157.3f, 172.3f, 198.9f, 222.4f, 243.6f, 263.1f, 281.3f, 298.4f, 314.5f, 346.6f, 374.7f, 399.8f, 422.4f, 443.f, 461.9f, 479.1f, 495.f, 523.1f, 547.f, 567.2f, 584.4f, 598.7f, 610.6f, 620.2f, 628.f, 633.9f, 638.3f, 641.3f, 643.f, 643.6f, 643.3f, 630.5f, 606.7f, 578.3f, 549.1f, 520.9f, 494.6f, 470.9f, 449.7f, 413.2f, 383.1f, 357.8f, 336.3f, 317.6f, 301.1f, 286.3f, 273.2f, 261.3f, 250.5f, 240.8f, 231.9f, 223.7f, 216.1f, 185.7f, 163.7f, 146.7f, 133.4f };
221 220
222 static const G4float e47[60] = { 228.8f, 280 221 static const G4float e47[60] = { 228.8f, 280.2f, 323.5f, 361.7f, 396.2f, 451.9f, 501.6f, 546.5f, 587.6f, 624.6f, 658.6f, 690.2f, 753.4f, 806.4f, 852.1f, 892.4f, 927.9f, 959.1f, 986.3f, 1010.f, 1052.f, 1086.f, 1114.f, 1135.f, 1150.f, 1160.f, 1166.f, 1168.f, 1166.f, 1163.f, 1157.f, 1149.f, 1140.f, 1129.f, 1064.f, 992.4f, 922.6f, 858.6f, 801.5f, 751.1f, 706.8f, 667.8f, 601.8f, 548.9f, 506.2f, 471.3f, 442.2f, 416.9f, 394.9f, 375.4f, 358.f, 342.3f, 328.2f, 315.3f, 303.4f, 292.6f, 249.3f, 218.3f, 194.8f, 176.3f };
223 222
224 static const G4float e48[60] = { 15.4f, 18.8 223 static const G4float e48[60] = { 15.4f, 18.86f, 21.78f, 24.35f, 26.68f, 30.81f, 34.44f, 37.73f, 40.75f, 43.57f, 46.21f, 48.71f, 53.82f, 58.34f, 62.45f, 66.21f, 69.7f, 72.94f, 75.98f, 78.85f, 84.11f, 88.84f, 93.13f, 97.02f, 100.6f, 103.8f, 106.8f, 109.5f, 111.9f, 114.1f, 116.2f, 118.f, 119.6f, 121.1f, 126.1f, 128.1f, 127.9f, 126.3f, 123.6f, 120.4f, 116.9f, 113.3f, 106.4f, 99.98f, 94.37f, 89.5f, 85.18f, 81.37f, 78.01f, 75.05f, 72.42f, 70.08f, 67.99f, 66.12f, 64.44f, 62.92f, 56.88f, 52.22f, 48.46f, 45.34f };
225 224
226 static const G4float e49[60] = { 61.33f, 75. 225 static const G4float e49[60] = { 61.33f, 75.12f, 86.74f, 96.97f, 106.2f, 122.2f, 136.3f, 149.1f, 160.8f, 171.6f, 181.6f, 191.1f, 210.f, 226.3f, 240.6f, 253.5f, 265.f, 275.4f, 284.8f, 293.3f, 308.f, 320.3f, 330.3f, 338.4f, 344.7f, 349.4f, 352.8f, 355.1f, 356.4f, 356.9f, 356.6f, 355.7f, 354.3f, 352.4f, 338.9f, 322.3f, 305.5f, 289.9f, 275.7f, 263.f, 251.7f, 241.5f, 224.4f, 210.3f, 198.3f, 188.f, 179.f, 171.f, 163.9f, 157.5f, 151.6f, 146.3f, 141.3f, 136.6f, 132.1f, 128.f, 110.6f, 97.75f, 87.92f, 80.2f };
227 226
228 static const G4float e50[60] = { 190.1f, 232 227 static const G4float e50[60] = { 190.1f, 232.8f, 268.8f, 300.5f, 329.2f, 374.1f, 414.4f, 450.7f, 483.9f, 513.5f, 540.4f, 565.4f, 616.5f, 659.2f, 696.f, 728.7f, 757.6f, 783.f, 805.3f, 825.2f, 860.1f, 890.2f, 915.4f, 935.4f, 950.4f, 961.1f, 968.2f, 972.3f, 973.8f, 973.f, 970.2f, 965.7f, 959.9f, 952.8f, 906.f, 850.6f, 795.1f, 743.3f, 696.4f, 654.7f, 617.9f, 585.3f, 530.7f, 487.1f, 451.3f, 421.4f, 395.9f, 373.7f, 354.4f, 337.2f, 321.9f, 308.1f, 295.7f, 284.2f, 273.8f, 264.1f, 225.7f, 198.1f, 177.1f, 160.5f };
229 228
230 static const G4float e51[60] = { 164.8f, 201 229 static const G4float e51[60] = { 164.8f, 201.8f, 233.f, 260.5f, 285.4f, 324.6f, 359.7f, 391.4f, 420.3f, 446.2f, 469.7f, 491.7f, 536.5f, 574.3f, 607.1f, 636.2f, 662.2f, 685.1f, 705.4f, 723.6f, 755.9f, 783.9f, 807.7f, 826.9f, 841.7f, 852.8f, 860.8f, 866.f, 868.8f, 869.6f, 868.6f, 866.f, 862.1f, 857.1f, 820.5f, 774.5f, 727.2f, 682.2f, 641.1f, 604.1f, 571.2f, 542.f, 492.3f, 452.1f, 419.1f, 391.8f, 368.6f, 348.4f, 330.7f, 314.9f, 300.8f, 288.1f, 276.7f, 266.1f, 256.5f, 247.6f, 212.f, 186.3f, 166.8f, 151.3f };
231 230
232 static const G4float e52[60] = { 184.9f, 226 231 static const G4float e52[60] = { 184.9f, 226.5f, 261.5f, 292.4f, 320.3f, 369.9f, 413.5f, 453.f, 489.3f, 523.f, 554.8f, 584.8f, 643.9f, 695.3f, 740.8f, 781.4f, 817.9f, 850.9f, 880.6f, 907.5f, 953.7f, 990.9f, 1021.f, 1044.f, 1061.f, 1073.f, 1081.f, 1084.f, 1085.f, 1083.f, 1078.f, 1072.f, 1064.f, 1054.f, 994.1f, 927.8f, 865.f, 808.9f, 760.1f, 717.9f, 681.4f, 649.7f, 591.2f, 541.f, 499.2f, 464.9f, 436.3f, 411.4f, 389.8f, 370.6f, 353.4f, 338.1f, 324.2f, 311.4f, 299.8f, 289.1f, 246.4f, 215.9f, 192.7f, 174.4f };
233 232
234 static const G4float e53[60] = { 150.9f, 184 233 static const G4float e53[60] = { 150.9f, 184.9f, 213.5f, 238.7f, 261.4f, 297.8f, 330.2f, 359.6f, 386.4f, 410.5f, 432.5f, 453.f, 494.9f, 530.3f, 561.2f, 588.8f, 613.4f, 635.3f, 654.7f, 672.3f, 703.4f, 730.5f, 753.5f, 772.3f, 787.f, 798.3f, 806.6f, 812.4f, 816.f, 817.6f, 817.4f, 815.8f, 812.9f, 808.9f, 777.6f, 736.5f, 693.4f, 651.9f, 613.6f, 579.1f, 548.2f, 520.6f, 472.9f, 434.f, 402.2f, 376.1f, 354.3f, 335.f, 318.2f, 303.2f, 289.8f, 277.7f, 266.7f, 256.6f, 247.4f, 238.9f, 204.8f, 180.2f, 161.4f, 146.5f };
235 234
236 static const G4float e54[60] = { 174.9f, 214 235 static const G4float e54[60] = { 174.9f, 214.2f, 247.4f, 276.6f, 303.f, 345.9f, 384.1f, 418.7f, 450.4f, 479.f, 505.2f, 529.8f, 579.1f, 620.8f, 657.1f, 689.3f, 717.9f, 743.2f, 765.7f, 785.8f, 820.8f, 850.5f, 875.f, 894.5f, 909.3f, 920.1f, 927.4f, 931.9f, 933.9f, 933.7f, 931.6f, 928.f, 922.9f, 916.8f, 874.4f, 823.2f, 771.3f, 722.4f, 676.2f, 634.f, 596.2f, 562.8f, 507.2f, 463.4f, 428.5f, 400.2f, 376.7f, 356.1f, 338.f, 321.9f, 307.5f, 294.5f, 282.8f, 272.f, 262.1f, 253.f, 216.6f, 190.4f, 170.3f, 154.5f };
237 236
238 static const G4float e55[60] = { 197.6f, 242 237 static const G4float e55[60] = { 197.6f, 242.f, 279.4f, 312.4f, 342.2f, 395.1f, 441.8f, 483.9f, 522.7f, 558.8f, 592.7f, 624.8f, 687.3f, 741.3f, 789.f, 831.4f, 869.4f, 903.6f, 934.3f, 962.f, 1009.f, 1048.f, 1078.f, 1101.f, 1119.f, 1132.f, 1140.f, 1145.f, 1146.f, 1145.f, 1141.f, 1135.f, 1128.f, 1119.f, 1060.f, 991.7f, 924.3f, 862.f, 806.1f, 756.5f, 712.8f, 674.2f, 605.4f, 548.6f, 503.3f, 467.5f, 438.7f, 413.6f, 391.8f, 372.4f, 355.2f, 339.7f, 325.7f, 312.9f, 301.1f, 290.4f, 247.5f, 216.8f, 193.5f, 175.1f };
239 238
240 static const G4float e56[60] = { 186.f, 227. 239 static const G4float e56[60] = { 186.f, 227.8f, 263.1f, 294.1f, 322.2f, 372.1f, 416.f, 455.7f, 492.2f, 526.2f, 558.1f, 588.3f, 645.3f, 693.8f, 735.8f, 772.4f, 804.5f, 832.6f, 857.2f, 878.7f, 913.5f, 939.2f, 957.2f, 969.f, 975.6f, 978.f, 976.8f, 972.8f, 966.5f, 958.4f, 948.9f, 938.3f, 926.9f, 914.9f, 851.3f, 789.6f, 734.2f, 685.4f, 642.9f, 605.7f, 573.f, 543.9f, 499.5f, 465.2f, 435.8f, 409.1f, 384.6f, 363.2f, 344.6f, 328.1f, 313.3f, 300.f, 288.f, 277.f, 266.8f, 257.5f, 220.4f, 193.6f, 173.2f, 157.f };
241 240
242 static const G4float e57[60] = { 83.27f, 102 241 static const G4float e57[60] = { 83.27f, 102.f, 117.8f, 131.7f, 144.2f, 165.f, 183.4f, 200.1f, 215.3f, 229.2f, 242.f, 254.f, 278.7f, 300.f, 319.f, 336.2f, 351.9f, 366.1f, 379.2f, 391.2f, 413.1f, 432.7f, 450.f, 465.1f, 478.f, 489.1f, 498.5f, 506.4f, 513.f, 518.4f, 522.7f, 526.f, 528.4f, 530.1f, 528.6f, 516.7f, 499.4f, 479.6f, 460.4f, 442.3f, 425.4f, 409.6f, 381.2f, 356.4f, 334.6f, 315.4f, 298.3f, 283.1f, 269.6f, 257.4f, 246.5f, 236.5f, 227.4f, 219.1f, 211.5f, 204.4f, 175.9f, 155.2f, 139.3f, 126.7f };
243 242
244 static const G4float e58[60] = { 140.2f, 171 243 static const G4float e58[60] = { 140.2f, 171.7f, 198.3f, 221.7f, 242.8f, 277.9f, 309.f, 337.2f, 363.f, 386.4f, 408.1f, 428.4f, 468.9f, 503.4f, 533.6f, 560.3f, 584.1f, 605.3f, 624.1f, 640.9f, 670.1f, 694.5f, 714.3f, 729.8f, 741.4f, 749.6f, 755.f, 757.9f, 758.6f, 757.6f, 755.f, 751.f, 746.f, 740.f, 701.3f, 656.7f, 612.9f, 572.6f, 538.7f, 509.8f, 484.7f, 462.4f, 424.2f, 392.4f, 365.3f, 341.8f, 321.1f, 302.9f, 287.3f, 273.6f, 261.3f, 250.1f, 240.1f, 231.1f, 222.7f, 215.f, 184.5f, 162.6f, 145.9f, 132.6f };
245 244
246 static const G4float e59[60] = { 2.183E+02f, 245 static const G4float e59[60] = { 2.183E+02f, 2.673E+02f, 3.087E+02f, 3.451E+02f, 3.780E+02f, 4.365E+02f, 4.880E+02f, 5.346E+02f, 5.775E+02f, 6.173E+02f, 6.548E+02f, 6.902E+02f, 7.595E+02f, 8.190E+02f, 8.714E+02f, 9.180E+02f, 9.597E+02f, 9.970E+02f, 1.031E+03f, 1.061E+03f, 1.112E+03f, 1.153E+03f, 1.186E+03f, 1.211E+03f, 1.229E+03f, 1.242E+03f, 1.250E+03f, 1.254E+03f, 1.254E+03f, 1.252E+03f, 1.247E+03f, 1.239E+03f, 1.230E+03f, 1.220E+03f, 1.153E+03f, 1.078E+03f, 1.005E+03f, 9.379E+02f, 8.769E+02f, 8.222E+02f, 7.736E+02f, 7.305E+02f, 6.583E+02f, 6.006E+02f, 5.537E+02f, 5.150E+02f, 4.825E+02f, 4.542E+02f, 4.297E+02f, 4.079E+02f, 3.886E+02f, 3.712E+02f, 3.556E+02f, 3.413E+02f, 3.282E+02f, 3.163E+02f, 2.687E+02f, 2.348E+02f, 2.092E+02f, 1.890E+02f };
247 246
248 static const G4float e60[60] = { 1.377E+01f 247 static const G4float e60[60] = { 1.377E+01f, 1.687E+01f, 1.947E+01f, 2.177E+01f, 2.385E+01f, 2.754E+01f, 3.079E+01f, 3.373E+01f, 3.643E+01f, 3.895E+01f, 4.131E+01f, 4.355E+01f, 4.811E+01f, 5.216E+01f, 5.583E+01f, 5.920E+01f, 6.232E+01f, 6.523E+01f, 6.795E+01f, 7.051E+01f, 7.522E+01f, 7.946E+01f, 8.330E+01f, 8.679E+01f, 8.997E+01f, 9.287E+01f, 9.551E+01f, 9.792E+01f, 1.001E+02f, 1.021E+02f, 1.039E+02f, 1.055E+02f, 1.070E+02f, 1.083E+02f, 1.127E+02f, 1.145E+02f, 1.144E+02f, 1.130E+02f, 1.108E+02f, 1.081E+02f, 1.053E+02f, 1.024E+02f, 9.692E+01f, 9.212E+01f, 8.791E+01f, 8.422E+01f, 8.097E+01f, 7.811E+01f, 7.556E+01f, 7.326E+01f, 7.117E+01f, 6.926E+01f, 6.750E+01f, 6.587E+01f, 6.436E+01f, 6.294E+01f, 5.700E+01f, 5.239E+01f, 4.867E+01f, 4.557E+01f };
249 248
250 static const G4float e61[60] = { 105.f, 128. 249 static const G4float e61[60] = { 105.f, 128.6f, 148.5f, 166.f, 181.8f, 210.f, 234.8f, 257.2f, 277.8f, 296.9f, 315.f, 332.f, 364.6f, 392.5f, 416.7f, 437.8f, 456.3f, 472.5f, 486.6f, 498.9f, 518.4f, 532.3f, 541.6f, 547.f, 549.2f, 549.f, 546.8f, 543.f, 538.f, 532.2f, 525.7f, 518.8f, 511.7f, 504.4f, 468.7f, 437.3f, 410.7f, 388.4f, 369.4f, 353.f, 338.7f, 325.9f, 304.1f, 285.9f, 270.3f, 256.7f, 245.2f, 235.1f, 225.9f, 217.2f, 209.f, 201.1f, 193.6f, 187.f, 180.9f, 175.3f, 152.3f, 135.4f, 122.2f, 111.7f };
251 250
252 static const G4float e62[60] = { 85.64f, 104 251 static const G4float e62[60] = { 85.64f, 104.9f, 121.1f, 135.4f, 148.3f, 171.3f, 191.5f, 209.8f, 226.6f, 242.2f, 256.9f, 270.8f, 298.1f, 322.f, 343.1f, 362.1f, 379.3f, 394.9f, 409.1f, 422.1f, 444.8f, 463.8f, 479.7f, 492.8f, 503.6f, 512.3f, 519.2f, 524.5f, 528.5f, 531.2f, 532.9f, 533.7f, 533.6f, 532.9f, 521.9f, 503.9f, 483.5f, 462.8f, 442.9f, 424.2f, 406.7f, 390.6f, 358.9f, 330.8f, 307.6f, 288.9f, 273.6f, 259.9f, 247.8f, 236.9f, 227.1f, 218.2f, 210.1f, 202.7f, 195.9f, 189.7f, 164.2f, 145.5f, 131.1f, 119.6f };
253 252
254 static const G4float e63[60] = { 179.2f, 219 253 static const G4float e63[60] = { 179.2f, 219.5f, 253.4f, 283.3f, 310.4f, 352.2f, 389.8f, 423.8f, 454.7f, 482.3f, 507.2f, 530.3f, 577.9f, 617.7f, 652.1f, 682.6f, 709.7f, 733.5f, 754.4f, 773.f, 806.1f, 835.1f, 859.5f, 879.2f, 894.2f, 905.1f, 912.7f, 917.4f, 919.6f, 919.6f, 917.7f, 914.2f, 909.3f, 903.3f, 861.4f, 810.7f, 759.2f, 710.8f, 667.1f, 628.3f, 593.9f, 563.4f, 511.8f, 470.2f, 435.9f, 407.2f, 382.8f, 361.4f, 342.9f, 326.4f, 311.6f, 298.4f, 286.4f, 275.4f, 265.3f, 256.1f, 219.f, 192.4f, 172.1f, 156.1f };
255 254
256 static const G4float e64[60] = { 60.88f, 74. 255 static const G4float e64[60] = { 60.88f, 74.56f, 86.1f, 96.26f, 105.4f, 121.8f, 136.1f, 149.1f, 161.1f, 172.2f, 182.6f, 192.5f, 212.3f, 229.6f, 245.1f, 259.1f, 271.8f, 283.5f, 294.2f, 304.1f, 321.6f, 336.5f, 349.1f, 359.8f, 368.7f, 376.f, 382.f, 386.7f, 390.2f, 392.7f, 394.4f, 395.2f, 395.3f, 394.8f, 385.4f, 369.3f, 351.f, 332.6f, 315.4f, 299.6f, 285.5f, 272.8f, 251.3f, 233.9f, 219.6f, 207.6f, 197.4f, 188.7f, 181.f, 174.1f, 167.9f, 162.2f, 157.f, 152.2f, 147.7f, 143.5f, 126.f, 112.5f, 101.7f, 93.07f };
257 256
258 static const G4float e65[60] = { 32.4f, 39.6 257 static const G4float e65[60] = { 32.4f, 39.68f, 45.82f, 51.23f, 56.12f, 64.8f, 72.45f, 79.36f, 85.72f, 91.64f, 97.2f, 102.5f, 113.1f, 122.4f, 130.8f, 138.5f, 145.5f, 151.9f, 157.9f, 163.4f, 173.4f, 182.f, 189.5f, 196.f, 201.6f, 206.4f, 210.4f, 213.8f, 216.5f, 218.7f, 220.4f, 221.6f, 222.4f, 222.9f, 220.7f, 213.8f, 204.9f, 195.4f, 186.f, 177.3f, 169.4f, 162.3f, 150.1f, 140.3f, 132.2f, 125.5f, 119.7f, 114.7f, 110.1f, 106.1f, 102.5f, 99.14f, 96.12f, 93.34f, 90.77f, 88.4f, 78.69f, 71.47f, 65.7f, 60.79f };
259 258
260 static const G4float e66[60] = { 174.5f, 213 259 static const G4float e66[60] = { 174.5f, 213.8f, 246.8f, 276.f, 302.3f, 349.1f, 390.3f, 427.6f, 461.8f, 493.7f, 523.6f, 552.f, 607.8f, 656.2f, 699.1f, 737.4f, 771.9f, 803.f, 831.1f, 856.5f, 900.2f, 935.6f, 963.8f, 985.9f, 1003.f, 1015.f, 1022.f, 1027.f, 1028.f, 1027.f, 1023.f, 1018.f, 1011.f, 1002.f, 947.8f, 885.5f, 824.9f, 769.7f, 720.9f, 678.2f, 641.f, 608.6f, 553.2f, 507.9f, 470.5f, 439.1f, 412.4f, 389.f, 368.6f, 350.5f, 334.4f, 319.8f, 306.6f, 294.6f, 283.6f, 273.5f, 233.2f, 204.4f, 182.5f, 165.2f };
261 260
262 static const G4float e67[60] = { 169.9f, 208 261 static const G4float e67[60] = { 169.9f, 208.1f, 240.3f, 268.6f, 294.3f, 339.8f, 379.9f, 416.2f, 449.5f, 480.6f, 509.7f, 537.3f, 591.5f, 638.4f, 679.9f, 716.9f, 750.2f, 780.2f, 807.4f, 832.f, 874.4f, 909.1f, 937.1f, 959.3f, 976.6f, 989.5f, 998.7f, 1005.f, 1008.f, 1008.f, 1006.f, 1003.f, 998.f, 991.7f, 946.9f, 891.9f, 836.f, 783.4f, 735.7f, 693.2f, 655.4f, 622.f, 563.5f, 515.3f, 475.8f, 443.4f, 416.2f, 392.6f, 371.9f, 353.6f, 337.3f, 322.5f, 309.3f, 297.1f, 286.f, 275.8f, 235.1f, 206.f, 183.9f, 166.5f };
263 262
264 static const G4float e68[60] = { 191.6f, 234 263 static const G4float e68[60] = { 191.6f, 234.7f, 271.f, 302.9f, 331.9f, 377.2f, 417.8f, 454.5f, 488.f, 517.9f, 545.1f, 570.4f, 621.9f, 665.f, 702.2f, 735.2f, 764.3f, 790.f, 812.4f, 832.5f, 867.7f, 897.9f, 923.2f, 943.3f, 958.3f, 969.f, 976.f, 980.f, 981.4f, 980.4f, 977.5f, 972.9f, 966.9f, 959.8f, 912.2f, 856.2f, 800.1f, 747.8f, 700.6f, 658.5f, 621.4f, 588.5f, 535.1f, 492.8f, 457.8f, 427.9f, 401.8f, 379.1f, 359.4f, 341.9f, 326.3f, 312.3f, 299.6f, 287.9f, 277.3f, 267.5f, 228.5f, 200.4f, 179.1f, 162.3f };
265 264
266 static const G4float e69[60] = { 18.38f, 22. 265 static const G4float e69[60] = { 18.38f, 22.51f, 25.99f, 29.06f, 31.83f, 36.75f, 41.09f, 45.01f, 48.62f, 51.98f, 55.13f, 58.11f, 64.19f, 69.57f, 74.44f, 78.9f, 83.02f, 86.84f, 90.42f, 93.78f, 99.92f, 105.4f, 110.3f, 114.8f, 118.8f, 122.4f, 125.7f, 128.6f, 131.2f, 133.5f, 135.6f, 137.4f, 139.f, 140.4f, 144.4f, 144.8f, 142.7f, 139.1f, 134.5f, 129.6f, 124.2f, 118.8f, 108.7f, 100.f, 92.78f, 86.72f, 81.64f, 77.35f, 73.71f, 70.6f, 67.92f, 65.59f, 63.56f, 61.77f, 60.19f, 58.78f, 53.33f, 49.08f, 45.62f, 42.73f };
267 266
268 static const G4float e70[60] = { 14.97f, 18. 267 static const G4float e70[60] = { 14.97f, 18.34f, 21.17f, 23.67f, 25.93f, 29.94f, 33.47f, 36.67f, 39.61f, 42.34f, 44.91f, 47.34f, 52.3f, 56.7f, 60.68f, 64.33f, 67.7f, 70.84f, 73.78f, 76.55f, 81.61f, 86.15f, 90.24f, 93.95f, 97.3f, 100.3f, 103.1f, 105.6f, 107.8f, 109.9f, 111.7f, 113.3f, 114.7f, 116.f, 120.f, 121.2f, 120.5f, 118.5f, 115.9f, 112.9f, 109.8f, 106.7f, 100.2f, 94.14f, 89.08f, 85.f, 81.68f, 78.78f, 76.2f, 73.88f, 71.78f, 69.86f, 68.09f, 66.45f, 64.93f, 63.5f, 57.52f, 52.9f, 49.15f, 46.04f };
269 268
270 static const G4float e71[60] = { 133.7f, 163 269 static const G4float e71[60] = { 133.7f, 163.8f, 189.1f, 211.4f, 231.6f, 267.5f, 299.f, 327.6f, 353.8f, 378.2f, 401.2f, 422.9f, 466.f, 503.6f, 537.2f, 567.3f, 594.6f, 619.5f, 642.1f, 662.8f, 698.9f, 729.f, 753.8f, 774.f, 790.1f, 802.6f, 811.9f, 818.3f, 822.3f, 824.1f, 823.9f, 822.2f, 819.f, 814.5f, 780.1f, 736.f, 695.9f, 660.4f, 628.6f, 599.9f, 573.7f, 549.7f, 507.5f, 471.4f, 440.1f, 412.8f, 388.8f, 367.6f, 348.9f, 332.2f, 317.2f, 303.7f, 291.4f, 280.3f, 270.f, 260.6f, 222.8f, 195.5f, 174.8f, 158.5f };
271 270
272 static const G4float e72[60] = { 152.2f, 186 271 static const G4float e72[60] = { 152.2f, 186.5f, 215.3f, 240.7f, 263.7f, 304.5f, 340.4f, 372.9f, 402.8f, 430.6f, 456.7f, 481.4f, 530.5f, 573.2f, 611.3f, 645.4f, 676.4f, 704.4f, 730.f, 753.3f, 793.8f, 827.4f, 855.f, 877.2f, 894.8f, 908.2f, 918.1f, 924.7f, 928.5f, 929.9f, 929.1f, 926.5f, 922.3f, 916.9f, 876.3f, 825.9f, 775.3f, 728.5f, 686.8f, 650.1f, 618.f, 589.9f, 536.1f, 489.2f, 450.8f, 420.2f, 395.5f, 373.8f, 354.6f, 337.5f, 322.2f, 308.3f, 295.8f, 284.4f, 273.9f, 264.3f, 225.7f, 198.f, 176.9f, 160.4f };
273 272
274 static const G4float e73[60] = { 38.1f, 46.6 273 static const G4float e73[60] = { 38.1f, 46.67f, 53.88f, 60.24f, 65.99f, 76.2f, 85.2f, 93.33f, 100.8f, 107.8f, 114.3f, 120.5f, 132.8f, 143.6f, 153.2f, 161.8f, 169.6f, 176.7f, 183.2f, 189.1f, 199.4f, 208.f, 215.f, 220.7f, 225.2f, 228.7f, 231.2f, 232.9f, 233.9f, 234.3f, 234.1f, 233.5f, 232.5f, 231.1f, 221.3f, 209.3f, 197.4f, 186.5f, 177.f, 168.6f, 161.3f, 155.f, 144.4f, 135.9f, 128.8f, 122.8f, 117.5f, 112.7f, 108.4f, 104.5f, 101.f, 97.75f, 94.8f, 92.09f, 89.6f, 87.29f, 77.8f, 70.52f, 64.84f, 60.13f };
275 274
276 if(0 == j) { AddData(e0, mat); } 275 if(0 == j) { AddData(e0, mat); }
277 else if(1 == j) { AddData(e1, mat); } 276 else if(1 == j) { AddData(e1, mat); }
278 else if(2 == j) { AddData(e2, mat); } 277 else if(2 == j) { AddData(e2, mat); }
279 else if(3 == j) { AddData(e3, mat); } 278 else if(3 == j) { AddData(e3, mat); }
280 else if(4 == j) { AddData(e4, mat); } 279 else if(4 == j) { AddData(e4, mat); }
281 else if(5 == j) { AddData(e5, mat); } 280 else if(5 == j) { AddData(e5, mat); }
282 else if(6 == j) { AddData(e6, mat); } 281 else if(6 == j) { AddData(e6, mat); }
283 else if(7 == j) { AddData(e7, mat); } 282 else if(7 == j) { AddData(e7, mat); }
284 else if(8 == j) { AddData(e8, mat); } 283 else if(8 == j) { AddData(e8, mat); }
285 else if(9 == j) { AddData(e9, mat); } 284 else if(9 == j) { AddData(e9, mat); }
286 else if(10 == j) { AddData(e10, mat); } 285 else if(10 == j) { AddData(e10, mat); }
287 else if(11 == j) { AddData(e11, mat); } 286 else if(11 == j) { AddData(e11, mat); }
288 else if(12 == j) { AddData(e12, mat); } 287 else if(12 == j) { AddData(e12, mat); }
289 else if(13 == j) { AddData(e13, mat); } 288 else if(13 == j) { AddData(e13, mat); }
290 else if(14 == j) { AddData(e14, mat); } 289 else if(14 == j) { AddData(e14, mat); }
291 else if(15 == j) { AddData(e15, mat); } 290 else if(15 == j) { AddData(e15, mat); }
292 else if(16 == j) { AddData(e16, mat); } 291 else if(16 == j) { AddData(e16, mat); }
293 else if(17 == j) { AddData(e17, mat); } 292 else if(17 == j) { AddData(e17, mat); }
294 else if(18 == j) { AddData(e18, mat); } 293 else if(18 == j) { AddData(e18, mat); }
295 else if(19 == j) { AddData(e19, mat); } 294 else if(19 == j) { AddData(e19, mat); }
296 else if(20 == j) { AddData(e20, mat); } 295 else if(20 == j) { AddData(e20, mat); }
297 else if(21 == j) { AddData(e21, mat); } 296 else if(21 == j) { AddData(e21, mat); }
298 else if(22 == j) { AddData(e22, mat); } 297 else if(22 == j) { AddData(e22, mat); }
299 else if(23 == j) { AddData(e23, mat); } 298 else if(23 == j) { AddData(e23, mat); }
300 else if(24 == j) { AddData(e24, mat); } 299 else if(24 == j) { AddData(e24, mat); }
301 else if(25 == j) { AddData(e25, mat); } 300 else if(25 == j) { AddData(e25, mat); }
302 else if(26 == j) { AddData(e26, mat); } 301 else if(26 == j) { AddData(e26, mat); }
303 else if(27 == j) { AddData(e27, mat); } 302 else if(27 == j) { AddData(e27, mat); }
304 else if(28 == j) { AddData(e28, mat); } 303 else if(28 == j) { AddData(e28, mat); }
305 else if(29 == j) { AddData(e29, mat); } 304 else if(29 == j) { AddData(e29, mat); }
306 else if(30 == j) { AddData(e30, mat); } 305 else if(30 == j) { AddData(e30, mat); }
307 else if(31 == j) { AddData(e31, mat); } 306 else if(31 == j) { AddData(e31, mat); }
308 else if(32 == j) { AddData(e32, mat); } 307 else if(32 == j) { AddData(e32, mat); }
309 else if(33 == j) { AddData(e33, mat); } 308 else if(33 == j) { AddData(e33, mat); }
310 else if(34 == j) { AddData(e34, mat); } 309 else if(34 == j) { AddData(e34, mat); }
311 else if(35 == j) { AddData(e35, mat); } 310 else if(35 == j) { AddData(e35, mat); }
312 else if(36 == j) { AddData(e36, mat); } 311 else if(36 == j) { AddData(e36, mat); }
313 else if(37 == j) { AddData(e37, mat); } 312 else if(37 == j) { AddData(e37, mat); }
314 else if(38 == j) { AddData(e38, mat); } 313 else if(38 == j) { AddData(e38, mat); }
315 else if(39 == j) { AddData(e39, mat); } 314 else if(39 == j) { AddData(e39, mat); }
316 else if(40 == j) { AddData(e40, mat); } 315 else if(40 == j) { AddData(e40, mat); }
317 else if(41 == j) { AddData(e41, mat); } 316 else if(41 == j) { AddData(e41, mat); }
318 else if(42 == j) { AddData(e42, mat); } 317 else if(42 == j) { AddData(e42, mat); }
319 else if(43 == j) { AddData(e43, mat); } 318 else if(43 == j) { AddData(e43, mat); }
320 else if(44 == j) { AddData(e44, mat); } 319 else if(44 == j) { AddData(e44, mat); }
321 else if(45 == j) { AddData(e45, mat); } 320 else if(45 == j) { AddData(e45, mat); }
322 else if(46 == j) { AddData(e46, mat); } 321 else if(46 == j) { AddData(e46, mat); }
323 else if(47 == j) { AddData(e47, mat); } 322 else if(47 == j) { AddData(e47, mat); }
324 else if(48 == j) { AddData(e48, mat); } 323 else if(48 == j) { AddData(e48, mat); }
325 else if(49 == j) { AddData(e49, mat); } 324 else if(49 == j) { AddData(e49, mat); }
326 else if(50 == j) { AddData(e50, mat); } 325 else if(50 == j) { AddData(e50, mat); }
327 else if(51 == j) { AddData(e51, mat); } 326 else if(51 == j) { AddData(e51, mat); }
328 else if(52 == j) { AddData(e52, mat); } 327 else if(52 == j) { AddData(e52, mat); }
329 else if(53 == j) { AddData(e53, mat); } 328 else if(53 == j) { AddData(e53, mat); }
330 else if(54 == j) { AddData(e54, mat); } 329 else if(54 == j) { AddData(e54, mat); }
331 else if(55 == j) { AddData(e55, mat); } 330 else if(55 == j) { AddData(e55, mat); }
332 else if(56 == j) { AddData(e56, mat); } 331 else if(56 == j) { AddData(e56, mat); }
333 else if(57 == j) { AddData(e57, mat); } 332 else if(57 == j) { AddData(e57, mat); }
334 else if(58 == j) { AddData(e58, mat); } 333 else if(58 == j) { AddData(e58, mat); }
335 else if(59 == j) { AddData(e59, mat); } 334 else if(59 == j) { AddData(e59, mat); }
336 else if(60 == j) { AddData(e60, mat); } 335 else if(60 == j) { AddData(e60, mat); }
337 else if(61 == j) { AddData(e61, mat); } 336 else if(61 == j) { AddData(e61, mat); }
338 else if(62 == j) { AddData(e62, mat); } 337 else if(62 == j) { AddData(e62, mat); }
339 else if(63 == j) { AddData(e63, mat); } 338 else if(63 == j) { AddData(e63, mat); }
340 else if(64 == j) { AddData(e64, mat); } 339 else if(64 == j) { AddData(e64, mat); }
341 else if(65 == j) { AddData(e65, mat); } 340 else if(65 == j) { AddData(e65, mat); }
342 else if(66 == j) { AddData(e66, mat); } 341 else if(66 == j) { AddData(e66, mat); }
343 else if(67 == j) { AddData(e67, mat); } 342 else if(67 == j) { AddData(e67, mat); }
344 else if(68 == j) { AddData(e68, mat); } 343 else if(68 == j) { AddData(e68, mat); }
345 else if(69 == j) { AddData(e69, mat); } 344 else if(69 == j) { AddData(e69, mat); }
346 else if(70 == j) { AddData(e70, mat); } 345 else if(70 == j) { AddData(e70, mat); }
347 else if(71 == j) { AddData(e71, mat); } 346 else if(71 == j) { AddData(e71, mat); }
348 else if(72 == j) { AddData(e72, mat); } 347 else if(72 == j) { AddData(e72, mat); }
349 else if(73 == j) { AddData(e73, mat); } 348 else if(73 == j) { AddData(e73, mat); }
350 } 349 }
351 350
352 //....oooOO0OOooo........oooOO0OOooo........oo 351 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
353 352
354 void G4PSTARStopping::AddData(const G4float* s 353 void G4PSTARStopping::AddData(const G4float* stop, const G4Material* mat)
355 { 354 {
356 auto v = new G4PhysicsFreeVector(60, T0[0], << 355 G4LPhysicsFreeVector* v =
357 for(size_t i=0; i<60; ++i) { << 356 new G4LPhysicsFreeVector(60, T0[0], T0[59]);
358 v->PutValues(i, T0[i], ((G4double)stop[i]) << 357 for(size_t i=0; i<60; ++i) { v->PutValues(i, T0[i], ((G4double)stop[i])*fac); }
359 } << 358 v->SetSpline(true);
360 v->FillSecondDerivatives(); 359 v->FillSecondDerivatives();
361 v->EnableLogBinSearch(G4EmParameters::Instan <<
362 materials.push_back(mat); 360 materials.push_back(mat);
363 sdata.push_back(v); 361 sdata.push_back(v);
364 ++nvectors; 362 ++nvectors;
365 } 363 }
366 364
367 //....oooOO0OOooo........oooOO0OOooo........oo 365 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
368 366