Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/lowenergy/src/G4ecpssrFormFactorMixsModel.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/lowenergy/src/G4ecpssrFormFactorMixsModel.cc (Version 11.3.0) and /processes/electromagnetic/lowenergy/src/G4ecpssrFormFactorMixsModel.cc (Version 11.1.3)


  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 // History:                                        26 // History:
 27 // -----------                                     27 // -----------
 28 //  01 Oct 2011   A.M., S.I. - 1st implementat     28 //  01 Oct 2011   A.M., S.I. - 1st implementation
 29 //                                                 29 // 
 30 // Class description                               30 // Class description
 31 // ----------------                                31 // ----------------
 32 //  Computation of K, L & M shell ECPSSR ionis     32 //  Computation of K, L & M shell ECPSSR ionisation cross sections for protons and alphas
 33 //  Based on the work of A. Taborda et al.         33 //  Based on the work of A. Taborda et al. 
 34 //  EXRS2012 proceedings                           34 //  EXRS2012 proceedings
 35 // -------------------------------------------     35 // ---------------------------------------------------------------------------------------
 36                                                    36 
 37 #include <fstream>                                 37 #include <fstream>
 38 #include <iomanip>                                 38 #include <iomanip>
 39 #include "globals.hh"                              39 #include "globals.hh"
 40 #include "G4ios.hh"                                40 #include "G4ios.hh"
 41 #include "G4SystemOfUnits.hh"                      41 #include "G4SystemOfUnits.hh"
 42 #include "G4EMDataSet.hh"                          42 #include "G4EMDataSet.hh"
 43 #include "G4LinInterpolation.hh"                   43 #include "G4LinInterpolation.hh"
 44 #include "G4Proton.hh"                             44 #include "G4Proton.hh"
 45 #include "G4Alpha.hh"                              45 #include "G4Alpha.hh"
 46 #include "G4ecpssrFormFactorMixsModel.hh"          46 #include "G4ecpssrFormFactorMixsModel.hh"
 47                                                    47 
 48 //....oooOO0OOooo........oooOO0OOooo........oo     48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 49                                                    49 
 50 G4ecpssrFormFactorMixsModel::G4ecpssrFormFacto     50 G4ecpssrFormFactorMixsModel::G4ecpssrFormFactorMixsModel()
 51 {                                                  51 { 
 52   interpolation = new G4LinInterpolation();        52   interpolation = new G4LinInterpolation();
 53                                                    53 
 54   for (G4int i=29; i<93; i++)                      54   for (G4int i=29; i<93; i++) 
 55   {                                                55   {
 56       protonM1DataSetMap[i] = new G4EMDataSet(     56       protonM1DataSetMap[i] = new G4EMDataSet(i,interpolation);
 57       protonM1DataSetMap[i]->LoadData("pixe/ec     57       protonM1DataSetMap[i]->LoadData("pixe/ecpssr/proton/m1-i01m001c01-");
 58                                                    58       
 59       protonM2DataSetMap[i] = new G4EMDataSet(     59       protonM2DataSetMap[i] = new G4EMDataSet(i,interpolation);
 60       protonM2DataSetMap[i]->LoadData("pixe/ec     60       protonM2DataSetMap[i]->LoadData("pixe/ecpssr/proton/m2-i01m001c01-");
 61                                                    61       
 62       protonM3DataSetMap[i] = new G4EMDataSet(     62       protonM3DataSetMap[i] = new G4EMDataSet(i,interpolation);
 63       protonM3DataSetMap[i]->LoadData("pixe/ec     63       protonM3DataSetMap[i]->LoadData("pixe/ecpssr/proton/m3-i01m001c01-");
 64                                                    64       
 65       protonM4DataSetMap[i] = new G4EMDataSet(     65       protonM4DataSetMap[i] = new G4EMDataSet(i,interpolation);
 66       protonM4DataSetMap[i]->LoadData("pixe/ec     66       protonM4DataSetMap[i]->LoadData("pixe/ecpssr/proton/m4-i01m001c01-");
 67                                                    67       
 68       protonM5DataSetMap[i] = new G4EMDataSet(     68       protonM5DataSetMap[i] = new G4EMDataSet(i,interpolation);
 69       protonM5DataSetMap[i]->LoadData("pixe/ec     69       protonM5DataSetMap[i]->LoadData("pixe/ecpssr/proton/m5-i01m001c01-");
 70   }                                                70   }
 71                                                    71 
 72   protonMiXsVector.push_back(protonM1DataSetMa     72   protonMiXsVector.push_back(protonM1DataSetMap);
 73   protonMiXsVector.push_back(protonM2DataSetMa     73   protonMiXsVector.push_back(protonM2DataSetMap);
 74   protonMiXsVector.push_back(protonM3DataSetMa     74   protonMiXsVector.push_back(protonM3DataSetMap);
 75   protonMiXsVector.push_back(protonM4DataSetMa     75   protonMiXsVector.push_back(protonM4DataSetMap);
 76   protonMiXsVector.push_back(protonM5DataSetMa     76   protonMiXsVector.push_back(protonM5DataSetMap);
 77                                                    77 
 78   for (G4int i=29; i<93; i++)                      78   for (G4int i=29; i<93; i++) 
 79   {                                                79   {
 80       alphaM1DataSetMap[i] = new G4EMDataSet(i     80       alphaM1DataSetMap[i] = new G4EMDataSet(i,interpolation);
 81       alphaM1DataSetMap[i]->LoadData("pixe/ecp     81       alphaM1DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m1-i02m004c02-");
 82                                                    82        
 83       alphaM2DataSetMap[i] = new G4EMDataSet(i     83       alphaM2DataSetMap[i] = new G4EMDataSet(i,interpolation);
 84       alphaM2DataSetMap[i]->LoadData("pixe/ecp     84       alphaM2DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m2-i02m004c02-");
 85                                                    85        
 86       alphaM3DataSetMap[i] = new G4EMDataSet(i     86       alphaM3DataSetMap[i] = new G4EMDataSet(i,interpolation);
 87       alphaM3DataSetMap[i]->LoadData("pixe/ecp     87       alphaM3DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m3-i02m004c02-");
 88                                                    88        
 89       alphaM4DataSetMap[i] = new G4EMDataSet(i     89       alphaM4DataSetMap[i] = new G4EMDataSet(i,interpolation);
 90       alphaM4DataSetMap[i]->LoadData("pixe/ecp     90       alphaM4DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m4-i02m004c02-");
 91                                                    91        
 92       alphaM5DataSetMap[i] = new G4EMDataSet(i     92       alphaM5DataSetMap[i] = new G4EMDataSet(i,interpolation);
 93       alphaM5DataSetMap[i]->LoadData("pixe/ecp     93       alphaM5DataSetMap[i]->LoadData("pixe/ecpssr/alpha/m5-i02m004c02-");
 94                                                    94 
 95   }                                                95   }
 96                                                    96 
 97   alphaMiXsVector.push_back(alphaM1DataSetMap)     97   alphaMiXsVector.push_back(alphaM1DataSetMap);
 98   alphaMiXsVector.push_back(alphaM2DataSetMap)     98   alphaMiXsVector.push_back(alphaM2DataSetMap);
 99   alphaMiXsVector.push_back(alphaM3DataSetMap)     99   alphaMiXsVector.push_back(alphaM3DataSetMap);
100   alphaMiXsVector.push_back(alphaM4DataSetMap)    100   alphaMiXsVector.push_back(alphaM4DataSetMap);
101   alphaMiXsVector.push_back(alphaM5DataSetMap)    101   alphaMiXsVector.push_back(alphaM5DataSetMap);
102 }                                                 102 }
103                                                   103 
104 //....oooOO0OOooo........oooOO0OOooo........oo    104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
105                                                   105 
106 G4ecpssrFormFactorMixsModel::~G4ecpssrFormFact    106 G4ecpssrFormFactorMixsModel::~G4ecpssrFormFactorMixsModel()
107 {                                                 107 { 
108   protonM1DataSetMap.clear();                     108   protonM1DataSetMap.clear();
109   alphaM1DataSetMap.clear();                      109   alphaM1DataSetMap.clear();
110                                                   110   
111   protonM2DataSetMap.clear();                     111   protonM2DataSetMap.clear();
112   alphaM2DataSetMap.clear();                      112   alphaM2DataSetMap.clear();
113                                                   113   
114   protonM3DataSetMap.clear();                     114   protonM3DataSetMap.clear();
115   alphaM3DataSetMap.clear();                      115   alphaM3DataSetMap.clear();
116                                                   116   
117   protonM4DataSetMap.clear();                     117   protonM4DataSetMap.clear();
118   alphaM4DataSetMap.clear();                      118   alphaM4DataSetMap.clear();
119                                                   119   
120   protonM5DataSetMap.clear();                     120   protonM5DataSetMap.clear();
121   alphaM5DataSetMap.clear();                      121   alphaM5DataSetMap.clear();
122                                                   122   
123   delete interpolation;                           123   delete interpolation;
124 }                                                 124 }
125                                                   125 
126 //....oooOO0OOooo........oooOO0OOooo........oo    126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127                                                   127 
128 G4double G4ecpssrFormFactorMixsModel::Calculat    128 G4double G4ecpssrFormFactorMixsModel::CalculateMiCrossSection(G4int zTarget,G4double massIncident, G4double energyIncident, G4int mShellId)
129 {                                                 129 {
130   G4Proton* aProton = G4Proton::Proton();         130   G4Proton* aProton = G4Proton::Proton();
131   G4Alpha* aAlpha = G4Alpha::Alpha();             131   G4Alpha* aAlpha = G4Alpha::Alpha();  
132   G4double sigma = 0;                             132   G4double sigma = 0;
133   G4int mShellIndex = mShellId -1;                133   G4int mShellIndex = mShellId -1;
134                                                   134 
135   if (energyIncident > 0.1*MeV && energyIncide    135   if (energyIncident > 0.1*MeV && energyIncident < 100*MeV && zTarget < 93 && zTarget > 28) {
136                                                   136 
137     if (massIncident == aProton->GetPDGMass())    137     if (massIncident == aProton->GetPDGMass())
138       {                                           138       {      
139   sigma = protonMiXsVector[mShellIndex][zTarge    139   sigma = protonMiXsVector[mShellIndex][zTarget]->FindValue(energyIncident/MeV);  
140         if (sigma !=0 && energyIncident > prot    140         if (sigma !=0 && energyIncident > protonMiXsVector[mShellIndex][zTarget]->GetEnergies(0).back()*MeV) return 0.;
141       }                                           141       }
142     else if (massIncident == aAlpha->GetPDGMas    142     else if (massIncident == aAlpha->GetPDGMass())
143       {                                           143       {
144         sigma = alphaMiXsVector[mShellIndex][z    144         sigma = alphaMiXsVector[mShellIndex][zTarget]->FindValue(energyIncident/MeV); 
145         if (sigma !=0 && energyIncident > alph    145         if (sigma !=0 && energyIncident > alphaMiXsVector[mShellIndex][zTarget]->GetEnergies(0).back()*MeV) return 0.;
146       }                                           146       }
147     else                                          147     else
148       {                                           148       { 
149   sigma = 0.;                                     149   sigma = 0.;
150       }                                           150       }
151   }                                               151   }
152                                                   152   
153   // sigma is in internal units: it has been c    153   // sigma is in internal units: it has been converted from 
154   // the input file in barns bt the EmDataset     154   // the input file in barns bt the EmDataset
155   return sigma;                                   155   return sigma;
156 }                                                 156 }
157                                                   157 
158 //....oooOO0OOooo........oooOO0OOooo........oo    158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159                                                   159 
160 G4double G4ecpssrFormFactorMixsModel::Calculat    160 G4double G4ecpssrFormFactorMixsModel::CalculateM1CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
161 {                                                 161 {
162   // mShellId                                     162   // mShellId
163   return  CalculateMiCrossSection (zTarget, ma    163   return  CalculateMiCrossSection (zTarget, massIncident, energyIncident, 1); 
164 }                                                 164 }
165                                                   165 
166 //....oooOO0OOooo........oooOO0OOooo........oo    166 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
167                                                   167 
168 G4double G4ecpssrFormFactorMixsModel::Calculat    168 G4double G4ecpssrFormFactorMixsModel::CalculateM2CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
169 {                                                 169 {
170   // mShellId                                     170   // mShellId
171   return  CalculateMiCrossSection (zTarget, ma    171   return  CalculateMiCrossSection (zTarget, massIncident, energyIncident, 2); 
172 }                                                 172 }
173                                                   173 
174 //....oooOO0OOooo........oooOO0OOooo........oo    174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
175                                                   175 
176 G4double G4ecpssrFormFactorMixsModel::Calculat    176 G4double G4ecpssrFormFactorMixsModel::CalculateM3CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
177 {                                                 177 {
178   return  CalculateMiCrossSection (zTarget, ma    178   return  CalculateMiCrossSection (zTarget, massIncident, energyIncident, 3); 
179 }                                                 179 }
180                                                   180 
181 //....oooOO0OOooo........oooOO0OOooo........oo    181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
182                                                   182 
183 G4double G4ecpssrFormFactorMixsModel::Calculat    183 G4double G4ecpssrFormFactorMixsModel::CalculateM4CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
184 {                                                 184 {
185   return  CalculateMiCrossSection (zTarget, ma    185   return  CalculateMiCrossSection (zTarget, massIncident, energyIncident, 4); 
186 }                                                 186 }
187                                                   187 
188 //....oooOO0OOooo........oooOO0OOooo........oo    188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
189                                                   189 
190 G4double G4ecpssrFormFactorMixsModel::Calculat    190 G4double G4ecpssrFormFactorMixsModel::CalculateM5CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident)
191 {                                                 191 {
192   return  CalculateMiCrossSection (zTarget, ma    192   return  CalculateMiCrossSection (zTarget, massIncident, energyIncident, 5); 
193 }                                                 193 }
194                                                   194