Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // 28 // 29 // History: 30 // ----------- 31 // 21 Apr 2009 ALF 1st implementation 32 // 29 Apr 2009 ALF Updated Desing for Integ 33 // 15 Mar 2011 ALF introduced the usage of 34 // 20 Oct 2011 ALF updated to take into acc 35 // 09 Mar 2012 LP update methods 36 // 09 Mar 2012 ALF update for M-shells Sim 37 // 10 Nov 2021 S. Guatelli & S. Bakr Added 38 // in Bakr et al, NIM B, vol. 4 39 // called here ECPSSR_ANSTO 40 41 #include "globals.hh" 42 #include "G4teoCrossSection.hh" 43 #include "G4Proton.hh" 44 #include "G4ecpssrBaseKxsModel.hh" 45 #include "G4ecpssrBaseLixsModel.hh" 46 #include "G4ecpssrFormFactorKxsModel.hh" 47 #include "G4ecpssrFormFactorLixsModel.hh" 48 #include "G4ecpssrFormFactorMixsModel.hh" 49 #include "G4ANSTOecpssrKxsModel.hh" 50 #include "G4ANSTOecpssrLixsModel.hh" 51 #include "G4ANSTOecpssrMixsModel.hh" 52 53 //....oooOO0OOooo........oooOO0OOooo........oo 54 G4teoCrossSection::G4teoCrossSection(const G4S 55 :G4VhShellCrossSection(nam),totalCS(0.0) 56 { 57 ecpssrShellMi = nullptr; 58 if (nam == "ECPSSR_Analytical") 59 { 60 ecpssrShellK = new G4ecpssrBaseKxsModel 61 ecpssrShellLi = new G4ecpssrBaseLixsMode 62 } 63 else if (nam == "ECPSSR_FormFactor") 64 { 65 ecpssrShellK = new G4ecpssrFormFactorKx 66 ecpssrShellLi = new G4ecpssrFormFactorLi 67 ecpssrShellMi = new G4ecpssrFormFactorMi 68 } 69 else if (nam == "ECPSSR_ANSTO") 70 { 71 ecpssrShellK = new G4ANSTOecpssrKxsMode 72 ecpssrShellLi = new G4ANSTOecpssrLixsMod 73 ecpssrShellMi = new G4ANSTOecpssrMixsMod 74 } 75 else 76 { 77 G4cout << "G4teoCrossSection::G4teoCross 78 << " in cross section name ECPSSR_Analy 79 << G4endl; 80 ecpssrShellK = new G4ecpssrBaseKxsModel 81 ecpssrShellLi = new G4ecpssrBaseLixsMode 82 } 83 } 84 85 //....oooOO0OOooo........oooOO0OOooo........oo 86 87 G4teoCrossSection::~G4teoCrossSection() 88 { 89 delete ecpssrShellK; 90 delete ecpssrShellLi; 91 delete ecpssrShellMi; 92 } 93 94 //....oooOO0OOooo........oooOO0OOooo........oo 95 96 std::vector<G4double> G4teoCrossSection::GetCr 97 G4double incidentEnergy, 98 G4double mass, 99 G4double, 100 const G4Material*) 101 { 102 std::vector<G4double> crossSections; 103 104 crossSections.push_back( ecpssrShellK->Calcu 105 106 crossSections.push_back( ecpssrShellLi->Calc 107 crossSections.push_back( ecpssrShellLi->Calc 108 crossSections.push_back( ecpssrShellLi->Calc 109 110 if (ecpssrShellMi) { 111 crossSections.push_back( ecpssrShellMi->Ca 112 crossSections.push_back( ecpssrShellMi->Ca 113 crossSections.push_back( ecpssrShellMi->Ca 114 crossSections.push_back( ecpssrShellMi->Ca 115 crossSections.push_back( ecpssrShellMi->Ca 116 } 117 return crossSections; 118 } 119 120 //....oooOO0OOooo........oooOO0OOooo........oo 121 122 G4double G4teoCrossSection::CrossSection(G4int 123 G4double incidentEnergy, 124 G4double mass, 125 const G4Material*) 126 { 127 G4double res = 0.0; 128 if(shell > 3 && !ecpssrShellMi) { 129 return res; 130 } 131 else if(shell > 8) { 132 return res; 133 } 134 else if(fKShell == shell) 135 { 136 res = ecpssrShellK->CalculateCrossSectio 137 } 138 else if(fL1Shell == shell) 139 { 140 res = ecpssrShellLi->CalculateL1CrossSec 141 } 142 else if(fL2Shell == shell) 143 { 144 res = ecpssrShellLi->CalculateL2CrossSec 145 } 146 else if(fL3Shell == shell) 147 { 148 res = ecpssrShellLi->CalculateL3CrossSec 149 } 150 else if(fM1Shell == shell) 151 { 152 res = ecpssrShellMi->CalculateM1CrossSec 153 } 154 else if(fM2Shell == shell) 155 { 156 res = ecpssrShellMi->CalculateM2CrossSec 157 } 158 else if(fM3Shell == shell) 159 { 160 res = ecpssrShellMi->CalculateM3CrossSec 161 } 162 else if(fM4Shell == shell) 163 { 164 res = ecpssrShellMi->CalculateM4CrossSec 165 } 166 else if(fM5Shell == shell) 167 { 168 res = ecpssrShellMi->CalculateM5CrossSec 169 } 170 return res; 171 } 172 173 //....oooOO0OOooo........oooOO0OOooo........oo 174 175 std::vector<G4double> G4teoCrossSection::Proba 176 G4double incidentEnergy, 177 G4double mass, 178 G4double deltaEnergy, 179 const G4Material*) 180 { 181 std::vector<G4double> crossSections = 182 GetCrossSection(Z, incidentEnergy, mass, d 183 184 for (size_t i=0; i<crossSections.size(); ++i 185 if (totalCS) { 186 crossSections[i] = crossSections[i]/tota 187 } 188 } 189 return crossSections; 190 } 191 192 //....oooOO0OOooo........oooOO0OOooo........oo 193 194 void G4teoCrossSection::SetTotalCS(G4double va 195 totalCS = val; 196 } 197 198 199 200