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 // 28 // GEANT4 Class file 29 // 30 // 31 // File name: G4EmLowEParameters 32 // 33 // Author: Vladimir Ivanchenko 34 // 35 // Creation date: 07.05.2019 36 // 37 // ------------------------------------------------------------------- 38 // 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 41 42 #include "G4EmLowEParameters.hh" 43 #include "G4PhysicalConstants.hh" 44 #include "G4UnitsTable.hh" 45 #include "G4SystemOfUnits.hh" 46 #include "G4VAtomDeexcitation.hh" 47 #include "G4EmLowEParametersMessenger.hh" 48 #include "G4RegionStore.hh" 49 #include "G4Region.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 52 53 G4EmLowEParameters::G4EmLowEParameters() 54 { 55 theMessenger = new G4EmLowEParametersMessenger(this); 56 Initialise(); 57 } 58 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 60 61 G4EmLowEParameters::~G4EmLowEParameters() 62 { 63 delete theMessenger; 64 } 65 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 67 68 void G4EmLowEParameters::Initialise() 69 { 70 fluo = false; 71 auger = false; 72 pixe = false; 73 deexIgnoreCut = false; 74 75 dnaFast = false; 76 dnaStationary = false; 77 dnaMsc = false; 78 dnaElectronSolvation = fMeesungnoen2002eSolvation; 79 fTimeStepModel = G4ChemTimeStepModel::Unknown; 80 81 fFluoDirectory = fluoDefault; 82 namePIXE = "Empirical"; 83 nameElectronPIXE = "Livermore"; 84 livDataDir = "epics_2017"; 85 } 86 87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 88 89 void G4EmLowEParameters::SetFluo(G4bool val) 90 { 91 fluo = val; 92 } 93 94 G4bool G4EmLowEParameters::Fluo() const 95 { 96 return fluo; 97 } 98 99 G4EmFluoDirectory G4EmLowEParameters::FluoDirectory() const 100 { 101 return fFluoDirectory; 102 } 103 104 void G4EmLowEParameters::SetFluoDirectory(G4EmFluoDirectory val) 105 { 106 fFluoDirectory = fluoDefault; 107 if(fluoBearden == val) { fFluoDirectory = fluoBearden; } 108 else if(fluoANSTO == val) { fFluoDirectory = fluoANSTO; } 109 else if(fluoXDB_EADL == val) { fFluoDirectory = fluoXDB_EADL; } 110 } 111 112 void G4EmLowEParameters::SetBeardenFluoDir(G4bool val) 113 { 114 fFluoDirectory = val ? fluoBearden : fluoDefault; 115 } 116 117 void G4EmLowEParameters::SetANSTOFluoDir(G4bool val) 118 { 119 fFluoDirectory = val ? fluoANSTO : fluoDefault; 120 } 121 122 void G4EmLowEParameters::SetXDB_EADLFluoDir(G4bool val) 123 { 124 fFluoDirectory = val ? fluoXDB_EADL : fluoDefault; 125 } 126 127 void G4EmLowEParameters::SetAuger(G4bool val) 128 { 129 auger = val; 130 if(val) { fluo = true; } 131 } 132 133 G4bool G4EmLowEParameters::Auger() const 134 { 135 return auger; 136 } 137 138 void G4EmLowEParameters::SetPixe(G4bool val) 139 { 140 pixe = val; 141 if(val) { fluo = true; } 142 } 143 144 G4bool G4EmLowEParameters::Pixe() const 145 { 146 return pixe; 147 } 148 149 void G4EmLowEParameters::SetDeexcitationIgnoreCut(G4bool val) 150 { 151 deexIgnoreCut = val; 152 } 153 154 G4bool G4EmLowEParameters::DeexcitationIgnoreCut() const 155 { 156 return deexIgnoreCut; 157 } 158 159 void G4EmLowEParameters::SetDNAFast(G4bool val) 160 { 161 dnaFast = val; 162 } 163 164 G4bool G4EmLowEParameters::DNAFast() const 165 { 166 return dnaFast; 167 } 168 169 void G4EmLowEParameters::SetDNAStationary(G4bool val) 170 { 171 dnaStationary = val; 172 } 173 174 G4bool G4EmLowEParameters::DNAStationary() const 175 { 176 return dnaStationary; 177 } 178 179 void G4EmLowEParameters::SetDNAElectronMsc(G4bool val) 180 { 181 dnaMsc = val; 182 } 183 184 G4bool G4EmLowEParameters::DNAElectronMsc() const 185 { 186 return dnaMsc; 187 } 188 189 void G4EmLowEParameters::SetDNAeSolvationSubType(G4DNAModelSubType val) 190 { 191 dnaElectronSolvation = val; 192 } 193 194 G4DNAModelSubType G4EmLowEParameters::DNAeSolvationSubType() const 195 { 196 return dnaElectronSolvation; 197 } 198 199 void G4EmLowEParameters::SetChemTimeStepModel(G4ChemTimeStepModel val) 200 { 201 fTimeStepModel = val; 202 } 203 204 G4ChemTimeStepModel G4EmLowEParameters::GetChemTimeStepModel() const 205 { 206 return fTimeStepModel; 207 } 208 209 void G4EmLowEParameters::SetPIXECrossSectionModel(const G4String& sss) 210 { 211 namePIXE = sss; 212 } 213 214 const G4String& G4EmLowEParameters::PIXECrossSectionModel() 215 { 216 return namePIXE; 217 } 218 219 void G4EmLowEParameters::SetPIXEElectronCrossSectionModel(const G4String& sss) 220 { 221 nameElectronPIXE = sss; 222 } 223 224 const G4String& G4EmLowEParameters::PIXEElectronCrossSectionModel() 225 { 226 return nameElectronPIXE; 227 } 228 229 void G4EmLowEParameters::SetLivermoreDataDir(const G4String& sss) 230 { 231 livDataDir = sss; 232 } 233 234 const G4String& G4EmLowEParameters::LivermoreDataDir() 235 { 236 return livDataDir; 237 } 238 239 void G4EmLowEParameters::PrintWarning(G4ExceptionDescription& ed) const 240 { 241 G4Exception("G4EmLowEParameters", "em0044", JustWarning, ed); 242 } 243 244 G4String G4EmLowEParameters::CheckRegion(const G4String& reg) const 245 { 246 G4String r = reg; 247 if(r == "" || r == "world" || r == "World") { 248 r = "DefaultRegionForTheWorld"; 249 } 250 return r; 251 } 252 253 void G4EmLowEParameters::AddMicroElec(const G4String& region) 254 { 255 G4String r = CheckRegion(region); 256 std::size_t nreg = m_regnamesME.size(); 257 for(std::size_t i=0; i<nreg; ++i) { 258 if(r == m_regnamesME[i]) { return; } 259 } 260 m_regnamesME.push_back(std::move(r)); 261 } 262 263 const std::vector<G4String>& G4EmLowEParameters::RegionsMicroElec() const 264 { 265 return m_regnamesME; 266 } 267 268 void G4EmLowEParameters::AddDNA(const G4String& region, const G4String& type) 269 { 270 G4String r = CheckRegion(region); 271 std::size_t nreg = m_regnamesDNA.size(); 272 for(std::size_t i=0; i<nreg; ++i) { 273 if(r == m_regnamesDNA[i]) { return; } 274 } 275 m_regnamesDNA.push_back(std::move(r)); 276 m_typesDNA.push_back(type); 277 } 278 279 const std::vector<G4String>& G4EmLowEParameters::RegionsDNA() const 280 { 281 return m_regnamesDNA; 282 } 283 284 const std::vector<G4String>& G4EmLowEParameters::TypesDNA() const 285 { 286 return m_typesDNA; 287 } 288 289 void 290 G4EmLowEParameters::SetDeexActiveRegion(const G4String& region, G4bool fdeex, 291 G4bool fauger, G4bool fpixe) 292 { 293 if(fdeex) { fluo = true; } 294 G4String r = CheckRegion(region); 295 std::size_t nreg = m_regnamesDeex.size(); 296 if(0 == nreg && r != "DefaultRegionForTheWorld") { 297 m_regnamesDeex.push_back("DefaultRegionForTheWorld"); 298 m_fluo.push_back(false); 299 m_auger.push_back(false); 300 m_pixe.push_back(false); 301 nreg = 1; 302 } 303 for(std::size_t i=0; i<nreg; ++i) { 304 if(r == m_regnamesDeex[i]) { 305 m_fluo[i] = fdeex; 306 m_auger[i]= fauger; 307 m_pixe[i] = fpixe; 308 return; 309 } 310 } 311 m_regnamesDeex.push_back(std::move(r)); 312 m_fluo.push_back(fdeex); 313 m_auger.push_back(fauger); 314 m_pixe.push_back(fpixe); 315 } 316 317 void G4EmLowEParameters::DefineRegParamForDeex(G4VAtomDeexcitation* ptr) const 318 { 319 std::size_t n = m_regnamesDeex.size(); 320 for(std::size_t i=0; i<n; ++i) { 321 ptr->SetDeexcitationActiveRegion(m_regnamesDeex[i], 322 m_fluo[i], m_auger[i], m_pixe[i]); 323 } 324 } 325 326 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 327