Loading...
Searching...
No Matches
G4MicroElecElasticModel_new Class Reference
#include <G4MicroElecElasticModel_new.hh>
Inheritance diagram for G4MicroElecElasticModel_new:
Protected Attributes | |
| G4ParticleChangeForGamma * | fParticleChangeForGamma |
| Protected Attributes inherited from G4VEmModel | |
| G4ElementData * | fElementData = nullptr |
| G4VParticleChange * | pParticleChange = nullptr |
| G4PhysicsTable * | xSectionTable = nullptr |
| const G4Material * | pBaseMaterial = nullptr |
| const std::vector< G4double > * | theDensityFactor = nullptr |
| const std::vector< G4int > * | theDensityIdx = nullptr |
| G4double | inveplus |
| G4double | pFactor = 1.0 |
| std::size_t | currentCoupleIndex = 0 |
| std::size_t | basedCoupleIndex = 0 |
| G4bool | lossFlucFlag = true |
Additional Inherited Members | |
| Protected Member Functions inherited from G4VEmModel | |
| G4ParticleChangeForLoss * | GetParticleChangeForLoss () |
| G4ParticleChangeForGamma * | GetParticleChangeForGamma () |
| virtual G4double | MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kineticEnergy) |
| const G4MaterialCutsCouple * | CurrentCouple () const |
| void | SetCurrentElement (const G4Element *) |
Detailed Description
Definition at line 90 of file G4MicroElecElasticModel_new.hh.
Constructor & Destructor Documentation
◆ G4MicroElecElasticModel_new()
| G4MicroElecElasticModel_new::G4MicroElecElasticModel_new | ( | const G4ParticleDefinition * | p = 0, |
| const G4String & | nam = "MicroElecElasticModel" ) |
Definition at line 89 of file G4MicroElecElasticModel_new.cc.
92{
93 killBelowEnergy = 0.1*eV; // Minimum e- energy for energy loss by excitation
94 lowEnergyLimit = 0.1 * eV;
95 lowEnergyLimitOfModel = 10 * eV; // The model lower energy is 10 eV
96 highEnergyLimit = 500. * keV;
97 SetLowEnergyLimit(lowEnergyLimit);
98 SetHighEnergyLimit(highEnergyLimit);
99
100 verboseLevel= 0;
101 // Verbosity scale:
102 // 0 = nothing
103 // 1 = warning for energy non-conservation
104 // 2 = details of energy budget
105 // 3 = calculation of cross sections, file openings, sampling of atoms
106 // 4 = entering in methods
107
108 if( verboseLevel>0 )
109 {
111 << "Energy range: "
112 << lowEnergyLimit / eV << " eV - "
113 << highEnergyLimit / MeV << " MeV"
114 << G4endl;
115 }
116 fParticleChangeForGamma = 0;
117
118 killElectron = false;
119 acousticModelEnabled = false;
120 currentMaterialName = "";
121}
G4GLOB_DLL std::ostream G4cout
G4ParticleChangeForGamma * fParticleChangeForGamma
Definition G4MicroElecElasticModel_new.hh:123
◆ ~G4MicroElecElasticModel_new()
|
override |
Definition at line 125 of file G4MicroElecElasticModel_new.cc.
126{
127 // For total cross section
128 TCSMap::iterator pos2;
129 for (pos2 = tableTCS.begin(); pos2 != tableTCS.end(); ++pos2) {
130 MapData* tableData = pos2->second;
131 std::map< G4String, G4MicroElecCrossSectionDataSet_new*, std::less<G4String> >::iterator pos;
132 for (pos = tableData->begin(); pos != tableData->end(); ++pos)
133 {
136 }
137 delete tableData;
138 }
139
140 //Clearing DCS maps
141
142 ThetaMap::iterator iterator_angle;
143 for (iterator_angle = thetaDataStorage.begin(); iterator_angle != thetaDataStorage.end(); ++iterator_angle) {
144 TriDimensionMap* eDiffCrossSectionData = iterator_angle->second;
145 eDiffCrossSectionData->clear();
146 delete eDiffCrossSectionData;
147 }
148
149 energyMap::iterator iterator_energy;
150 for (iterator_energy = eIncidentEnergyStorage.begin(); iterator_energy != eIncidentEnergyStorage.end(); ++iterator_energy) {
151 std::vector<G4double>* eTdummyVec = iterator_energy->second;
152 eTdummyVec->clear();
153 delete eTdummyVec;
154 }
155
156 ProbaMap::iterator iterator_proba;
157 for (iterator_proba = eProbaStorage.begin(); iterator_proba != eProbaStorage.end(); ++iterator_proba) {
158 VecMap* eVecm = iterator_proba->second;
159 eVecm->clear();
160 delete eVecm;
161 }
162
163}
Member Function Documentation
◆ AcousticCrossSectionPerVolume()
| G4double G4MicroElecElasticModel_new::AcousticCrossSectionPerVolume | ( | G4double | ekin, |
| G4double | kbz, | ||
| G4double | rho, | ||
| G4double | cs, | ||
| G4double | Aac, | ||
| G4double | Eac, | ||
| G4double | prefactor ) |
Definition at line 405 of file G4MicroElecElasticModel_new.cc.
412{
413
414 G4double e = 1.6e-19,
415 m0 = 9.10938356e-31,
416 h = 1.0546e-34,
417 kb = 1.38e-23;
418
419 G4double E = (ekin / eV) * e;
420 G4double D = (2 / (std::sqrt(2) * std::pow(pi, 2) * std::pow(h, 3))) * (1 + 2 * E) * std::pow(m0, 1.5) * std::sqrt(E);
421
422 // Parametres SiO2
423 G4double T = 300,
424 Ebz = (std::pow(h, 2) * std::pow(kbz, 2)) / (2 * m0),
425 hwbz = cs * kbz * h,
426 nbz = 1.0 / (exp(hwbz / (kb * T)) - 1),
427 Pac;
428
429 if (E < Ebz / 4.0)
430 {
432 }
433
434 else if (E > Ebz) //Screened relationship
435 {
436 Pac = ((2 * pi * m0 * (2 * nbz + 1)) / (h * rho * hwbz)) * std::pow(Eac, 2) * D * E * 2 * std::pow((Aac / E), 2) * (((-E / Aac) / (1 + (E / Aac))) + log(1 + (E / Aac)));
437 }
438 else //Linear interpolation
439 {
440 G4double fEbz = ((2 * pi * m0 * (2 * nbz + 1)) / (h * rho * hwbz)) * std::pow(Eac, 2) * D * Ebz * 2 * std::pow((Aac / Ebz), 2) * (((-Ebz / Aac) / (1 + (Ebz / Aac))) + log(1 + (Ebz / Aac)));
441 G4double fEbz4 = ((pi * kb * T) / (h * std::pow(cs, 2) * rho)) * (std::pow(Eac, 2) * D) / (1 + ((Ebz / 4) / Aac));
444 }
445
446 G4double MFP = (std::sqrt(2 * E / m0) / (prefactor * Pac)) * m;
447
448 return 1 / MFP;
449}
Referenced by CrossSectionPerVolume().
◆ CrossSectionPerVolume()
|
overridevirtual |
Reimplemented from G4VEmModel.
Definition at line 283 of file G4MicroElecElasticModel_new.cc.
288{
289 if (verboseLevel > 3)
291
293 const G4DataVector cuts;
294 // Calculate total cross section for model
295 MapEnergy::iterator lowEPos;
296 lowEPos = lowEnergyLimitTable.find(currentMaterialName);
297
298 MapEnergy::iterator highEPos;
299 highEPos = highEnergyLimitTable.find(currentMaterialName);
300
301 MapEnergy::iterator killEPos;
302 killEPos = workFunctionTable.find(currentMaterialName);
303
304 if (lowEPos == lowEnergyLimitTable.end() || highEPos == highEnergyLimitTable.end() || killEPos == workFunctionTable.end())
305 {
306 G4String str = "Material ";
307 str += currentMaterialName + " not found!";
309 return 0;
310 }
311 else {
312 // G4cout << "normal elastic " << G4endl;
313 lowEnergyLimit = lowEPos->second;
314 highEnergyLimit = highEPos->second;
315 killBelowEnergy = killEPos->second;
316 }
317
319
320 G4double sigma=0;
321
322 //Phonon for SiO2
323 if (currentMaterialName == "SILICON_DIOXIDE" && ekin < 100 * eV) {
324 acousticModelEnabled = true;
325
326 //Values for SiO2
327 G4double kbz = 11.54e9,
328 rho = 2.2 * 1000, // [g/cm3] * 1000
329 cs = 3560, //Sound speed
330 Ebz = 5.1 * 1.6e-19,
331 Aac = 17 * Ebz, //A screening parameter
332 Eac = 3.5 * 1.6e-19, //C deformation potential
333 prefactor = 2.2;// Facteur pour modifier les MFP
334
336 }
337
338 else if (currentMaterialName == "ALUMINUM_OXIDE" && ekin < 20 * eV) {
339 acousticModelEnabled = true;
340
341 //Values for Al2O3
342 G4double kbz = 8871930614.247564,
343 rho = 3.97 * 1000, // [g/cm3] * 1000
344 cs = 233329.07733059773, //Sound speed
345 Aac = 2.9912494342262614e-19, //A screening parameter
346 Eac = 2.1622471654789847e-18, //C deformation potential
347 prefactor = 1;
349 }
350 //Elastic
351 else {
352 acousticModelEnabled = false;
353
356
357 TCSMap::iterator tablepos;
358 tablepos = tableTCS.find(currentMaterialName);
359
360 if (tablepos != tableTCS.end())
361 {
362 MapData* tableData = tablepos->second;
363
364 if (ekin >= lowEnergyLimit && ekin < highEnergyLimit)
365 {
366 std::map< G4String, G4MicroElecCrossSectionDataSet_new*, std::less<G4String> >::iterator pos;
367 pos = tableData->find(particleName);
368
369 if (pos != tableData->end()){
371 if (table != 0)
372 {
373 sigma = table->FindValue(ekin);
374 }
375 }
376 else
377 {
378 G4Exception("G4MicroElecElasticModel_new::ComputeCrossSectionPerVolume", "em0002", FatalException, "Model not applicable to particle type.");
379 }
380 }
382 }
383 else {
384 G4String str = "Material ";
385 str += currentMaterialName + " TCS table not found!";
386 G4Exception("G4MicroElecElasticModel_new::ComputeCrossSectionPerVolume", "em0002", FatalException, str);
387 }
388
389 if (verboseLevel > 3) {
393 }
394
395 // Hsing-YinChangaAndrewAlvaradoaTreyWeberaJaimeMarianab Monte Carlo modeling of low-energy electron-induced secondary electron emission yields in micro-architected boron nitride surfaces - ScienceDirect, (n.d.). https://www.sciencedirect.com/science/article/pii/S0168583X19304069 (accessed April 1, 2022).
396 if (currentMaterialName == "BORON_NITRIDE") {
397 sigma = sigma * tanh(0.5 * pow(ekin / 5.2e-6, 2));
398 }
399 return sigma*density;
400 }
401}
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition G4Exception.cc:59
G4double GetTotNbOfAtomsPerVolume() const
Definition G4Material.hh:197
G4double AcousticCrossSectionPerVolume(G4double ekin, G4double kbz, G4double rho, G4double cs, G4double Aac, G4double Eac, G4double prefactor)
Definition G4MicroElecElasticModel_new.cc:405
const G4String & GetParticleName() const
Definition G4ParticleDefinition.hh:96
◆ DamageEnergy()
Definition at line 504 of file G4MicroElecElasticModel_new.cc.
505{
506 //.................. T in eV!!!!!!!!!!!!!
507 G4double Z2= Z;
509 G4double k_d;
510 G4double epsilon_d;
511 G4double g_epsilon_d;
512 G4double E_nu;
513
514 k_d=0.1334*std::pow(Z2,(2./3.))*std::pow(M2,(-1./2.));
515 epsilon_d=0.01014*std::pow(Z2,(-7./3.))*(T/eV);
516 g_epsilon_d= epsilon_d+0.40244*std::pow(epsilon_d,(3./4.))+3.4008*std::pow(epsilon_d,(1./6.));
517
518 E_nu=1./(1.+ k_d*g_epsilon_d);
519
520 return E_nu;
521}
◆ GetKillBelowThreshold()
|
inline |
Definition at line 118 of file G4MicroElecElasticModel_new.hh.
118{ return killBelowEnergy; }
◆ Initialise()
|
overridevirtual |
Implements G4VEmModel.
Definition at line 167 of file G4MicroElecElasticModel_new.cc.
169{
170 if (verboseLevel > -1)
172 // Energy limits
173 // Reading of data files
174
175 G4double scaleFactor = 1e-18 * cm * cm;
176
177 G4ProductionCutsTable* theCoupleTable = G4ProductionCutsTable::GetProductionCutsTable();
179
182
183 G4cout << "MicroElasticModel, Material " << i + 1 << " / " << numOfCouples << " : " << material->GetName() << G4endl;
185
188
189 currentMaterialStructure = new G4MicroElecMaterialStructure(matName);
190 lowEnergyLimitTable[matName]=currentMaterialStructure->GetElasticModelLowLimit();
191 highEnergyLimitTable[matName]=currentMaterialStructure->GetElasticModelHighLimit();
192 workFunctionTable[matName] = currentMaterialStructure->GetWorkFunction();
193
194 delete currentMaterialStructure;
195
197 G4String fileElectron = "Elastic/elsepa_elastic_cross_e_" + matName;
199
200 G4ParticleDefinition* electronDef = G4Electron::ElectronDefinition();
202
203 // For total cross section
204 MapData* tableData = new MapData();
205
206 G4MicroElecCrossSectionDataSet_new* tableE = new G4MicroElecCrossSectionDataSet_new(new G4LogLogInterpolation, eV, scaleFactor);
207 tableE->LoadData(fileElectron);
208 tableData->insert(make_pair(electron, tableE));
209 tableTCS[matName] = tableData; //Storage of TCS
210
211 // For final state
213 if (!path)
214 {
215 G4Exception("G4MicroElecElasticModel_new::Initialise","em0006",FatalException,"G4LEDATA environment variable not set.");
216 return;
217 }
218
219 //Reading DCS file
220 std::ostringstream eFullFileName;
221 eFullFileName << path << "/microelec/Elastic/elsepa_elastic_cumulated_diffcross_e_" + matName + ".dat";
223 std::ifstream eDiffCrossSection(eFullFileName.str().c_str());
224
225 if (!eDiffCrossSection)
226 G4Exception("G4MicroElecElasticModel_new::Initialise", "em0003", FatalException, "Missing data file: /microelec/sigmadiff_cumulated_elastic_e_Si.dat");
227
228 // October 21th, 2014 - Melanie Raine
229 // Added clear for MT
230 // Diff Cross Sections in cumulated mode
231 TriDimensionMap* eDiffCrossSectionData = new TriDimensionMap(); //Angles
232 std::vector<G4double>* eTdummyVec = new std::vector<G4double>; //Incident energy vector
233 VecMap* eProbVec = new VecMap; //Probabilities
234
235 eTdummyVec->push_back(0.);
236
237 while (!eDiffCrossSection.eof())
238 {
241 eDiffCrossSection >> tDummy >> eProb;
242
243 // SI : mandatory eVecm initialization
244 if (tDummy != eTdummyVec->back())
245 {
246 eTdummyVec->push_back(tDummy); //adding values for incident energy points
247 (*eProbVec)[tDummy].push_back(0.); //adding probability for the first angle, equal to 0
248 }
249
250 eDiffCrossSection >> (*eDiffCrossSectionData)[tDummy][eProb]; //adding Angle Value to map
251
252 if (eProb != (*eProbVec)[tDummy].back()) {
253 (*eProbVec)[tDummy].push_back(eProb); //Adding cumulated proba to map
254 }
255
256 }
257
258 //Filling maps for the material
259 thetaDataStorage[matName] = eDiffCrossSectionData;
260 eIncidentEnergyStorage[matName] = eTdummyVec;
261 eProbaStorage[matName] = eProbVec;
262 }
263 // End final state
264
265 if (verboseLevel > 2)
267
268 if (verboseLevel > 0) {
270 << "Energy range: "
274 }
275
276 if (isInitialised) { return; }
278 isInitialised = true;
279}
const G4Material * GetMaterial() const
Definition G4MaterialCutsCouple.hh:166
G4bool LoadData(const G4String &argFileName) override
Definition G4MicroElecCrossSectionDataSet_new.cc:59
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
Definition G4ProductionCutsTable.hh:256
std::size_t GetTableSize() const
Definition G4ProductionCutsTable.hh:249
static G4ProductionCutsTable * GetProductionCutsTable()
Definition G4ProductionCutsTable.cc:58
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition G4VEmModel.cc:121
◆ SampleSecondaries()
|
overridevirtual |
Implements G4VEmModel.
Definition at line 454 of file G4MicroElecElasticModel_new.cc.
459{
460 if (verboseLevel > 3)
462
464
465 if (electronEnergy0 < killBelowEnergy)
466 {
467 fParticleChangeForGamma->SetProposedKineticEnergy(0.);
469 fParticleChangeForGamma->ProposeLocalEnergyDeposit(electronEnergy0);
470 return;
471 }
472
473 if (electronEnergy0 < highEnergyLimit)
474 {
475 G4double cosTheta = 0;
476 if (acousticModelEnabled)
477 {
479 }
480 else if (electronEnergy0 >= lowEnergyLimit)
481 {
482 cosTheta = RandomizeCosTheta(electronEnergy0);
483 }
485
489
490 G4double xDir = std::sqrt(1. - cosTheta*cosTheta);
491 G4double yDir = xDir;
492 xDir *= std::cos(phi);
493 yDir *= std::sin(phi);
494
495 G4ThreeVector zPrimeVers((xDir*xVers + yDir*yVers + cosTheta*zVers));
496
497 fParticleChangeForGamma->ProposeMomentumDirection(zPrimeVers.unit());
498 fParticleChangeForGamma->SetProposedKineticEnergy(electronEnergy0);
499 }
500}
Hep3Vector orthogonal() const
Hep3Vector cross(const Hep3Vector &) const
const G4ThreeVector & GetMomentumDirection() const
G4double GetKineticEnergy() const
◆ SetKillBelowThreshold()
| void G4MicroElecElasticModel_new::SetKillBelowThreshold | ( | G4double | threshold | ) |
Definition at line 679 of file G4MicroElecElasticModel_new.cc.
680{
681 killBelowEnergy = threshold;
682
683 if (threshold < 5*CLHEP::eV)
684 {
685 G4Exception ("*** WARNING : the G4MicroElecElasticModel class is not validated below 5 eV !","",JustWarning,"") ;
686 threshold = 5*CLHEP::eV;
687 }
688}
Member Data Documentation
◆ fParticleChangeForGamma
|
protected |
Definition at line 123 of file G4MicroElecElasticModel_new.hh.
Referenced by G4MicroElecElasticModel_new(), Initialise(), and SampleSecondaries().
The documentation for this class was generated from the following files:
