Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 // History: 28 // ----------- 29 // 21 Apr 2008 H. Abdelohauwed - 1st implementation 30 // 29 Apr 2009 ALF Major Design Revision 31 // 22 Oct 2011 ALF Warning removals and Z checks fix for alphas 32 // 29 Oct 2011 ALF Changed name to G4PAulKxsModel 33 // 34 // ------------------------------------------------------------------- 35 36 // Class description: 37 // Low Energy Electromagnetic Physics, Cross section, p ionisation, K shell 38 // Further documentation available from http://www.ge.infn.it/geant4/lowE 39 40 // ------------------------------------------------------------------- 41 42 #include <fstream> 43 #include <iomanip> 44 45 #include "G4PaulKxsModel.hh" 46 47 #include "globals.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "G4ios.hh" 50 #include "G4EMDataSet.hh" 51 #include "G4LogLogInterpolation.hh" 52 #include "G4Proton.hh" 53 #include "G4Alpha.hh" 54 55 // ------------------------------------------------------------------- 56 57 G4PaulKxsModel::G4PaulKxsModel() 58 { 59 interpolation = new G4LogLogInterpolation(); 60 61 for (G4int i=4; i<93; i++) { 62 protonDataSetMap[i] = new G4EMDataSet(i,interpolation); 63 protonDataSetMap[i]->LoadData("pixe/kpcsPaul/kcs-"); 64 } 65 for (G4int i=6; i<93; i++) { 66 alphaDataSetMap[i] = new G4EMDataSet(i,interpolation); 67 alphaDataSetMap[i]->LoadData("pixe/kacsPaul/kacs-"); 68 } 69 } 70 71 G4PaulKxsModel::~G4PaulKxsModel() 72 { 73 protonDataSetMap.clear(); 74 alphaDataSetMap.clear(); 75 delete interpolation; 76 } 77 78 // ------------------------------------------------------------------- 79 G4double G4PaulKxsModel::CalculateKCrossSection(G4int zTarget, 80 G4double massIncident, G4double energyIncident) 81 { 82 83 G4Proton* aProtone = G4Proton::Proton(); 84 G4Alpha* aAlpha = G4Alpha::Alpha(); 85 G4double sigma = 0; 86 87 if (massIncident == aProtone->GetPDGMass() && zTarget < 93 && zTarget > 3) 88 { 89 if (energyIncident > protonDataSetMap[zTarget]->GetEnergies(0).back() || 90 energyIncident < protonDataSetMap[zTarget]->GetEnergies(0).front() ) 91 {sigma = 0;} 92 else { 93 sigma = protonDataSetMap[zTarget]->FindValue(energyIncident/MeV); 94 } 95 } 96 else 97 { 98 if (massIncident == aAlpha->GetPDGMass() && zTarget < 93 && zTarget > 5) 99 { 100 if (energyIncident > alphaDataSetMap[zTarget]->GetEnergies(0).back() || 101 energyIncident < alphaDataSetMap[zTarget]->GetEnergies(0).front() ) 102 {sigma = 0;} 103 else { 104 sigma = alphaDataSetMap[zTarget]->FindValue(energyIncident/MeV); 105 } 106 } 107 else 108 { 109 sigma = 0.; 110 } 111 } 112 113 // sigma is in internal units (mm^2) 114 return sigma; 115 } 116 117 118 119 120 121 122 123 124 125 126 127 128