Loading...
Searching...
No Matches
G4PenelopeGammaConversionModel Class Reference
#include <G4PenelopeGammaConversionModel.hh>
Inheritance diagram for G4PenelopeGammaConversionModel:
Protected Attributes | |
| G4ParticleChangeForGamma * | fParticleChange |
| Protected Attributes inherited from G4VEmModel | |
| G4VParticleChange * | pParticleChange |
| G4PhysicsTable * | xSectionTable |
| const std::vector< G4double > * | theDensityFactor |
| const std::vector< G4int > * | theDensityIdx |
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 56 of file G4PenelopeGammaConversionModel.hh.
Constructor & Destructor Documentation
◆ G4PenelopeGammaConversionModel()
| G4PenelopeGammaConversionModel::G4PenelopeGammaConversionModel | ( | const G4ParticleDefinition * | p = 0, |
| const G4String & | processName = "PenConversion" |
||
| ) |
Definition at line 52 of file G4PenelopeGammaConversionModel.cc.
55 fEffectiveCharge(0),fMaterialInvScreeningRadius(0),
56 fScreeningFunction(0),isInitialised(false)
57{
58 fIntrinsicLowEnergyLimit = 2.0*electron_mass_c2;
59 fIntrinsicHighEnergyLimit = 100.0*GeV;
60 fSmallEnergy = 1.1*MeV;
61 InitializeScreeningRadii();
62
63 // SetLowEnergyLimit(fIntrinsicLowEnergyLimit);
64 SetHighEnergyLimit(fIntrinsicHighEnergyLimit);
65 //
66 verboseLevel= 0;
67 // Verbosity scale:
68 // 0 = nothing
69 // 1 = warning for energy non-conservation
70 // 2 = details of energy budget
71 // 3 = calculation of cross sections, file openings, sampling of atoms
72 // 4 = entering in methods
73}
G4ParticleChangeForGamma * fParticleChange
Definition: G4PenelopeGammaConversionModel.hh:86
Definition: G4VEmModel.hh:102
◆ ~G4PenelopeGammaConversionModel()
|
virtual |
Definition at line 77 of file G4PenelopeGammaConversionModel.cc.
78{
79 std::map <G4int,G4PhysicsFreeVector*>::iterator i;
80 if (logAtomicCrossSection)
81 {
82 for (i=logAtomicCrossSection->begin();i != logAtomicCrossSection->end();i++)
83 if (i->second) delete i->second;
84 delete logAtomicCrossSection;
85 }
86 if (fEffectiveCharge)
87 delete fEffectiveCharge;
88 if (fMaterialInvScreeningRadius)
89 delete fMaterialInvScreeningRadius;
90 if (fScreeningFunction)
91 delete fScreeningFunction;
92}
Member Function Documentation
◆ ComputeCrossSectionPerAtom()
|
virtual |
Reimplemented from G4VEmModel.
Definition at line 143 of file G4PenelopeGammaConversionModel.cc.
148{
149 //
150 // Penelope model v2008.
151 // Cross section (including triplet production) read from database and managed
152 // through the G4CrossSectionHandler utility. Cross section data are from
153 // M.J. Berger and J.H. Hubbel (XCOM), Report NBSIR 887-3598
154 //
155
156 if (energy < fIntrinsicLowEnergyLimit)
157 return 0;
158
160
161 //read data files
162 if (!logAtomicCrossSection->count(iZ))
163 ReadDataFile(iZ);
164 //now it should be ok
165 if (!logAtomicCrossSection->count(iZ))
166 {
167 G4ExceptionDescription ed;
171 }
172
173 G4double cs = 0;
174 G4double logene = std::log(energy);
175 G4PhysicsFreeVector* theVec = logAtomicCrossSection->find(iZ)->second;
176
178 cs = std::exp(logXS);
179
180 if (verboseLevel > 2)
183 return cs;
184}
G4DLLIMPORT std::ostream G4cout
Definition: G4PhysicsFreeVector.hh:63
G4double Value(G4double theEnergy)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
◆ GetVerbosityLevel()
|
inline |
Definition at line 83 of file G4PenelopeGammaConversionModel.hh.
83{return verboseLevel;};
◆ Initialise()
|
virtual |
Implements G4VEmModel.
Definition at line 97 of file G4PenelopeGammaConversionModel.cc.
99{
100 if (verboseLevel > 3)
102
103 // logAtomicCrossSection is created only once, since it is never cleared
104 if (!logAtomicCrossSection)
105 logAtomicCrossSection = new std::map<G4int,G4PhysicsFreeVector*>;
106
107 //delete old material data...
108 if (fEffectiveCharge)
109 {
110 delete fEffectiveCharge;
111 fEffectiveCharge = 0;
112 }
113 if (fMaterialInvScreeningRadius)
114 {
115 delete fMaterialInvScreeningRadius;
116 fMaterialInvScreeningRadius = 0;
117 }
118 if (fScreeningFunction)
119 {
120 delete fScreeningFunction;
121 fScreeningFunction = 0;
122 }
123 //and create new ones
124 fEffectiveCharge = new std::map<const G4Material*,G4double>;
125 fMaterialInvScreeningRadius = new std::map<const G4Material*,G4double>;
126 fScreeningFunction = new std::map<const G4Material*,std::pair<G4double,G4double> >;
127
128 if (verboseLevel > 0) {
130 << "Energy range: "
133 << G4endl;
134 }
135
136 if(isInitialised) return;
138 isInitialised = true;
139}
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:109
◆ SampleSecondaries()
|
virtual |
Implements G4VEmModel.
Definition at line 189 of file G4PenelopeGammaConversionModel.cc.
194{
195 //
196 // Penelope model v2008.
197 // Final state is sampled according to the Bethe-Heitler model with Coulomb
198 // corrections, according to the semi-empirical model of
199 // J. Baro' et al., Radiat. Phys. Chem. 44 (1994) 531.
200 //
201 // The model uses the high energy Coulomb correction from
202 // H. Davies et al., Phys. Rev. 93 (1954) 788
203 // and atomic screening radii tabulated from
204 // J.H. Hubbel et al., J. Phys. Chem. Ref. Data 9 (1980) 1023
205 // for Z= 1 to 92.
206 //
207 if (verboseLevel > 3)
209
211
212 // Always kill primary
215
216 if (photonEnergy <= fIntrinsicLowEnergyLimit)
217 {
219 return ;
220 }
221
224
225 //check if material data are available
226 if (!fEffectiveCharge->count(mat))
227 InitializeScreeningFunctions(mat);
228 if (!fEffectiveCharge->count(mat))
229 {
230 G4ExceptionDescription ed;
231 ed << "Unable to allocate the EffectiveCharge data for " <<
235 }
236
237 // eps is the fraction of the photon energy assigned to e- (including rest mass)
238 G4double eps = 0;
239 G4double eki = electron_mass_c2/photonEnergy;
240
241 //Do it fast for photon energy < 1.1 MeV (close to threshold)
242 if (photonEnergy < fSmallEnergy)
243 eps = eki + (1.0-2.0*eki)*G4UniformRand();
244 else
245 {
246 //Complete calculation
247 G4double effC = fEffectiveCharge->find(mat)->second;
248 G4double alz = effC*fine_structure_const;
249 G4double T = std::sqrt(2.0*eki);
251 +(8.523+7.326e1*alz-4.441e1*alz*alz)*T*T
252 -(1.352e1+1.211e2*alz-9.641e1*alz*alz)*T*T*T
253 +(8.946+6.205e1*alz-6.341e1*alz*alz)*T*T*T*T;
254
255 G4double F0b = fScreeningFunction->find(mat)->second.second;
257 G4double invRad = fMaterialInvScreeningRadius->find(mat)->second;
258 G4double bmin = 4.0*eki/invRad;
259 std::pair<G4double,G4double> scree = GetScreeningFunctions(bmin);
260 G4double g1 = scree.first;
261 G4double g2 = scree.second;
262 G4double g1min = g1+g0;
263 G4double g2min = g2+g0;
264 G4double xr = 0.5-eki;
265 G4double a1 = 2.*g1min*xr*xr/3.;
266 G4double p1 = a1/(a1+g2min);
267
269 //Random sampling of eps
270 do{
271 loopAgain = false;
273 {
275 if (ru2m1 < 0)
276 eps = 0.5-xr*std::pow(std::abs(ru2m1),1./3.);
277 else
278 eps = 0.5+xr*std::pow(ru2m1,1./3.);
279 G4double B = eki/(invRad*eps*(1.0-eps));
280 scree = GetScreeningFunctions(B);
281 g1 = scree.first;
282 g1 = std::max(g1+g0,0.);
284 loopAgain = true;
285 }
286 else
287 {
288 eps = eki+2.0*xr*G4UniformRand();
289 G4double B = eki/(invRad*eps*(1.0-eps));
290 scree = GetScreeningFunctions(B);
291 g2 = scree.second;
292 g2 = std::max(g2+g0,0.);
294 loopAgain = true;
295 }
296 }while(loopAgain);
297
298 }
299 if (verboseLevel > 4)
301
302 G4double electronTotEnergy = eps*photonEnergy;
303 G4double positronTotEnergy = (1.0-eps)*photonEnergy;
304
305 // Scattered electron (positron) angles. ( Z - axis along the parent photon)
306
307 //electron kinematics
308 G4double electronKineEnergy = std::max(0.,electronTotEnergy - electron_mass_c2) ;
310 G4double kk = std::sqrt(electronKineEnergy*(electronKineEnergy+2.*electron_mass_c2));
311 costheta_el = (costheta_el*electronTotEnergy+kk)/(electronTotEnergy+costheta_el*kk);
313 G4double dirX_el = std::sqrt(1.-costheta_el*costheta_el) * std::cos(phi_el);
314 G4double dirY_el = std::sqrt(1.-costheta_el*costheta_el) * std::sin(phi_el);
315 G4double dirZ_el = costheta_el;
316
317 //positron kinematics
318 G4double positronKineEnergy = std::max(0.,positronTotEnergy - electron_mass_c2) ;
320 kk = std::sqrt(positronKineEnergy*(positronKineEnergy+2.*electron_mass_c2));
321 costheta_po = (costheta_po*positronTotEnergy+kk)/(positronTotEnergy+costheta_po*kk);
323 G4double dirX_po = std::sqrt(1.-costheta_po*costheta_po) * std::cos(phi_po);
324 G4double dirY_po = std::sqrt(1.-costheta_po*costheta_po) * std::sin(phi_po);
325 G4double dirZ_po = costheta_po;
326
327 // Kinematics of the created pair:
328 // the electron and positron are assumed to have a symetric angular
329 // distribution with respect to the Z axis along the parent photon
330 G4double localEnergyDeposit = 0. ;
331
332 if (electronKineEnergy > 0.0)
333 {
334 G4ThreeVector electronDirection ( dirX_el, dirY_el, dirZ_el);
335 electronDirection.rotateUz(photonDirection);
337 electronDirection,
338 electronKineEnergy);
339 fvect->push_back(electron);
340 }
341 else
342 {
343 localEnergyDeposit += electronKineEnergy;
344 electronKineEnergy = 0;
345 }
346
347 //Generate the positron. Real particle in any case, because it will annihilate. If below
348 //threshold, produce it at rest
349 if (positronKineEnergy < 0.0)
350 {
351 localEnergyDeposit += positronKineEnergy;
352 positronKineEnergy = 0; //produce it at rest
353 }
354 G4ThreeVector positronDirection(dirX_po,dirY_po,dirZ_po);
355 positronDirection.rotateUz(photonDirection);
357 positronDirection, positronKineEnergy);
358 fvect->push_back(positron);
359
360 //Add rest of energy to the local energy deposit
362
363 if (verboseLevel > 1)
364 {
369 if (electronKineEnergy)
371 << G4endl;
372 if (positronKineEnergy)
375 if (localEnergyDeposit)
378 localEnergyDeposit+2.0*electron_mass_c2)/keV <<
381 }
382 if (verboseLevel > 0)
383 {
384 G4double energyDiff = std::fabs(electronKineEnergy+positronKineEnergy+
385 localEnergyDeposit+2.0*electron_mass_c2-photonEnergy);
386 if (energyDiff > 0.05*keV)
388 << (electronKineEnergy+positronKineEnergy+
389 localEnergyDeposit+2.0*electron_mass_c2)/keV
391 }
392}
Definition: ThreeVector.h:41
Definition: G4DynamicParticle.hh:74
const G4ThreeVector & GetMomentumDirection() const
G4double GetKineticEnergy() const
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150
Definition: G4Material.hh:119
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForGamma.hh:129
void ProposeTrackStatus(G4TrackStatus status)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
◆ SetVerbosityLevel()
|
inline |
Definition at line 82 of file G4PenelopeGammaConversionModel.hh.
82{verboseLevel = lev;};
Member Data Documentation
◆ fParticleChange
|
protected |
Definition at line 86 of file G4PenelopeGammaConversionModel.hh.
Referenced by Initialise(), and SampleSecondaries().
The documentation for this class was generated from the following files:
