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