00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 #include <fstream>
00038 #include <iomanip>
00039
00040 #include "globals.hh"
00041 #include "G4ios.hh"
00042 #include "G4SystemOfUnits.hh"
00043
00044 #include "G4EMDataSet.hh"
00045 #include "G4LinInterpolation.hh"
00046 #include "G4Proton.hh"
00047 #include "G4Alpha.hh"
00048
00049 #include "G4ecpssrFormFactorMixsModel.hh"
00050
00051
00052
00053 G4ecpssrFormFactorMixsModel::G4ecpssrFormFactorMixsModel()
00054 {
00055 interpolation = new G4LinInterpolation();
00056
00057 for (G4int i=62; i<93; i++)
00058 {
00059 protonM1DataSetMap[i] = new G4EMDataSet(i,interpolation);
00060 protonM1DataSetMap[i]->LoadData("pixe/ecpssr/proton/m1-");
00061
00062 protonM2DataSetMap[i] = new G4EMDataSet(i,interpolation);
00063 protonM2DataSetMap[i]->LoadData("pixe/ecpssr/proton/m2-");
00064
00065 protonM3DataSetMap[i] = new G4EMDataSet(i,interpolation);
00066 protonM3DataSetMap[i]->LoadData("pixe/ecpssr/proton/m3-");
00067
00068 protonM4DataSetMap[i] = new G4EMDataSet(i,interpolation);
00069 protonM4DataSetMap[i]->LoadData("pixe/ecpssr/proton/m4-");
00070
00071 protonM5DataSetMap[i] = new G4EMDataSet(i,interpolation);
00072 protonM5DataSetMap[i]->LoadData("pixe/ecpssr/proton/m5-");
00073 }
00074
00075 protonMiXsVector.push_back(protonM1DataSetMap);
00076 protonMiXsVector.push_back(protonM2DataSetMap);
00077 protonMiXsVector.push_back(protonM3DataSetMap);
00078 protonMiXsVector.push_back(protonM4DataSetMap);
00079 protonMiXsVector.push_back(protonM5DataSetMap);
00080
00081
00082 for (G4int i=62; i<93; i++)
00083 {
00084 alphaM1DataSetMap[i] = new G4EMDataSet(i,interpolation);
00085 alphaM1DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m1-");
00086
00087 alphaM2DataSetMap[i] = new G4EMDataSet(i,interpolation);
00088 alphaM2DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m2-");
00089
00090 alphaM3DataSetMap[i] = new G4EMDataSet(i,interpolation);
00091 alphaM3DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m3-");
00092
00093 alphaM4DataSetMap[i] = new G4EMDataSet(i,interpolation);
00094 alphaM4DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m4-");
00095
00096 alphaM5DataSetMap[i] = new G4EMDataSet(i,interpolation);
00097 alphaM5DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m5-");
00098
00099 }
00100
00101 alphaMiXsVector.push_back(alphaM1DataSetMap);
00102 alphaMiXsVector.push_back(alphaM2DataSetMap);
00103 alphaMiXsVector.push_back(alphaM3DataSetMap);
00104 alphaMiXsVector.push_back(alphaM4DataSetMap);
00105 alphaMiXsVector.push_back(alphaM5DataSetMap);
00106
00107
00108
00109 }
00110
00111
00112
00113 G4ecpssrFormFactorMixsModel::~G4ecpssrFormFactorMixsModel()
00114 {
00115 protonM1DataSetMap.clear();
00116 alphaM1DataSetMap.clear();
00117
00118 protonM2DataSetMap.clear();
00119 alphaM2DataSetMap.clear();
00120
00121 protonM3DataSetMap.clear();
00122 alphaM3DataSetMap.clear();
00123
00124 protonM4DataSetMap.clear();
00125 alphaM4DataSetMap.clear();
00126
00127 protonM5DataSetMap.clear();
00128 alphaM5DataSetMap.clear();
00129
00130 delete interpolation;
00131 }
00132
00133
00134
00135 G4double G4ecpssrFormFactorMixsModel::CalculateMiCrossSection(G4int zTarget,G4double massIncident, G4double energyIncident, G4int mShellId)
00136 {
00137 G4Proton* aProton = G4Proton::Proton();
00138 G4Alpha* aAlpha = G4Alpha::Alpha();
00139 G4double sigma = 0;
00140 G4int mShellIndex = mShellId -1;
00141
00142 if (energyIncident > 0.1*MeV && energyIncident < 10*MeV && zTarget < 93 && zTarget > 61) {
00143
00144 if (massIncident == aProton->GetPDGMass())
00145 {
00146 sigma = protonMiXsVector[mShellIndex][zTarget]->FindValue(energyIncident/MeV);
00147 if (sigma !=0 && energyIncident > protonMiXsVector[mShellIndex][zTarget]->GetEnergies(0).back()*MeV) return 0.;
00148 }
00149 else if (massIncident == aAlpha->GetPDGMass())
00150 {
00151 sigma = alphaMiXsVector[mShellIndex][zTarget]->FindValue(energyIncident/MeV);
00152 if (sigma !=0 && energyIncident > protonMiXsVector[mShellIndex][zTarget]->GetEnergies(0).back()*MeV) return 0.;
00153 }
00154 else
00155 {
00156 sigma = 0.;
00157 }
00158 }
00159
00160
00161
00162 return sigma;
00163 }
00164
00165
00166
00167 G4double G4ecpssrFormFactorMixsModel::CalculateM1CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
00168 {
00169
00170
00171 return CalculateMiCrossSection (zTarget, massIncident, energyIncident, 1);
00172
00173 }
00174
00175
00176
00177 G4double G4ecpssrFormFactorMixsModel::CalculateM2CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
00178 {
00179
00180
00181 return CalculateMiCrossSection (zTarget, massIncident, energyIncident, 2);
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211 }
00212
00213
00214
00215 G4double G4ecpssrFormFactorMixsModel::CalculateM3CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
00216 {
00217
00218 return CalculateMiCrossSection (zTarget, massIncident, energyIncident, 3);
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248 }
00249
00250
00251
00252 G4double G4ecpssrFormFactorMixsModel::CalculateM4CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
00253 {
00254
00255 return CalculateMiCrossSection (zTarget, massIncident, energyIncident, 4);
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275
00276
00277
00278
00279
00280
00281
00282
00283 }
00284
00285
00286
00287 G4double G4ecpssrFormFactorMixsModel::CalculateM5CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
00288 {
00289
00290 return CalculateMiCrossSection (zTarget, massIncident, energyIncident, 5);
00291
00292
00293
00294
00295
00296
00297
00298
00299
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315
00316
00317
00318 }