Loading...
Searching...
No Matches
G4eplusTo3GammaOKVIModel Class Reference
#include <G4eplusTo3GammaOKVIModel.hh>
Inheritance diagram for G4eplusTo3GammaOKVIModel:
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 *) |
| 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 |
| size_t | currentCoupleIndex = 0 |
| size_t | basedCoupleIndex = 0 |
| G4bool | lossFlucFlag = true |
Detailed Description
Definition at line 53 of file G4eplusTo3GammaOKVIModel.hh.
Constructor & Destructor Documentation
◆ G4eplusTo3GammaOKVIModel() [1/2]
|
explicit |
Definition at line 60 of file G4eplusTo3GammaOKVIModel.cc.
62 : G4VEmModel(nam), fDelta(0.001)
63{
64 theGamma = G4Gamma::Gamma();
65 fParticleChange = nullptr;
66}
Definition: G4VEmModel.hh:103
◆ ~G4eplusTo3GammaOKVIModel()
|
overridedefault |
◆ G4eplusTo3GammaOKVIModel() [2/2]
|
delete |
Member Function Documentation
◆ ComputeCrossSectionPerAtom()
|
finalvirtual |
Reimplemented from G4VEmModel.
Definition at line 202 of file G4eplusTo3GammaOKVIModel.cc.
206{
207 // Calculates the cross section per atom of annihilation into two photons
208
209
210
212 return cross;
213}
G4double ComputeCrossSectionPerElectron(G4double kinEnergy)
Definition: G4eplusTo3GammaOKVIModel.cc:181
◆ ComputeCrossSectionPerElectron()
Definition at line 181 of file G4eplusTo3GammaOKVIModel.cc.
182{
183 // Calculates the cross section per electron of annihilation into 3 photons
184 // from the Heilter formula.
185
186 G4double ekin = std::max(eV,kinEnergy);
187 G4double tau = ekin/electron_mass_c2;
190 G4double bg2 = tau * (tau+2.0);
191 G4double bg = sqrt(bg2);
192
194 - (gam+3.)/(sqrt(gam*gam - 1.));
195
197 return cross;
198}
Referenced by ComputeCrossSectionPerAtom(), CrossSectionPerVolume(), G4eplusTo2GammaOKVIModel::Initialise(), and G4eplusTo2GammaOKVIModel::SampleSecondaries().
◆ ComputeF()
| G4double G4eplusTo3GammaOKVIModel::ComputeF | ( | G4double | fr1, |
| G4double | fr2, | ||
| G4double | fr3, | ||
| G4double | kinEnergy | ||
| ) |
Definition at line 85 of file G4eplusTo3GammaOKVIModel.cc.
87{
88 G4double ekin = std::max(eV,kinEnergy);
89 G4double tau = ekin/electron_mass_c2;
92 G4double bg2 = tau * (tau+2.0);
93 G4double bg = sqrt(bg2);
94
96 - (gam+3.)/(sqrt(gam*gam - 1.)) + 1.;
97 G4double border;
98
99 if(ekin < 500*MeV) {
100 border = 1. - (electron_mass_c2)/(2*(ekin + electron_mass_c2));
101 } else {
102 border = 1. - (100*electron_mass_c2)/(2*(ekin + electron_mass_c2));
103 }
104
105 border = std::min(border, 0.9999);
106
107 if (fr1>border) { fr1 = border; }
108 if (fr2>border) { fr2 = border; }
109 if (fr3>border) { fr3 = border; }
110
112 G4double fr2s = fr2*fr2;
113 G4double fr3s = fr3*fr3;
114
115 G4double aa = (1.-fr1)*(1.-fr2);
116 G4double ab = fr3s + (fr1-fr2)*(fr1-fr2);
117 G4double add= ((1.-fr1)*(1.-fr1) + (1.-fr2)*(1.-fr2))/(fr3s*aa);
118
119 G4double fres = -rho*(1./fr1s + 1./fr2s)
120 + (ab/(2.*(fr1*fr2*aa)))*(G4Log(2.*gam*aa/(fr1*fr2)))
121 + (ab/(2.*fr1*fr2*(1-fr3)))*G4Log(2.*gam*(1.-fr3)/(fr1*fr2)) - add;
122
123 return fres;
124}
Referenced by ComputeFS().
◆ ComputeF0()
Definition at line 129 of file G4eplusTo3GammaOKVIModel.cc.
131{
132 G4double tau = 0.0;
135 G4double bg2 = tau * (tau+2.0);
136 G4double bg = sqrt(bg2);
137
139 - (gam+3.)/(sqrt(gam*gam - 1.)) + 1.;
140 G4double border = 0.5;
141
142 if (fr1>border) { fr1 = border; }
143 if (fr2>border) { fr2 = border; }
144 if (fr3>border) { fr3 = border; }
145
147 G4double fr2s = fr2*fr2;
148 G4double fr3s = fr3*fr3;
149
150 G4double aa = (1.-fr1)*(1.-fr2);
151 G4double ab = fr3s + (fr1-fr2)*(fr1-fr2);
152 G4double add= ((1.-fr1)*(1.-fr1) + (1.-fr2)*(1.-fr2))/(fr3s*aa);
153
154 G4double fres = -rho*(1./fr1s + 1./fr2s)
155 + (ab/(2.*(fr1*fr2*aa)))*(G4Log(2.*gam*aa/(fr1*fr2)))
156 + (ab/(2.*fr1*fr2*(1-fr3)))*G4Log(2.*gam*(1.-fr3)/(fr1*fr2)) - add;
157
158 return fres;
159}
◆ ComputeFS()
| G4double G4eplusTo3GammaOKVIModel::ComputeFS | ( | G4double | fr1, |
| G4double | fr2, | ||
| G4double | fr3, | ||
| G4double | kinEnergy | ||
| ) |
Definition at line 162 of file G4eplusTo3GammaOKVIModel.cc.
164{
165 G4double ekin = std::max(eV,kinEnergy);
166 G4double tau = ekin/electron_mass_c2;
168
170 ComputeF(fr3,fr1,fr2,ekin) +
171 ComputeF(fr2,fr3,fr1,ekin));
172
174
175 return dcross;
176}
G4double ComputeF(G4double fr1, G4double fr2, G4double fr3, G4double kinEnergy)
Definition: G4eplusTo3GammaOKVIModel.cc:85
Referenced by SampleSecondaries().
◆ CrossSectionPerVolume()
|
finalvirtual |
Reimplemented from G4VEmModel.
Definition at line 217 of file G4eplusTo3GammaOKVIModel.cc.
222{
223 // Calculates the cross section per volume of annihilation into two photons
224
227 return cross;
228}
◆ Initialise()
|
finalvirtual |
Implements G4VEmModel.
Definition at line 74 of file G4eplusTo3GammaOKVIModel.cc.
76{
77 // here particle change is set for the triplet model
78 if(fParticleChange) { return; }
79 fParticleChange = GetParticleChangeForGamma();
80}
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:124
Referenced by G4eplusTo2GammaOKVIModel::Initialise().
◆ operator=()
|
delete |
◆ SampleSecondaries()
|
finalvirtual |
Implements G4VEmModel.
Definition at line 236 of file G4eplusTo3GammaOKVIModel.cc.
240{
241
243 G4DynamicParticle *aGamma1, *aGamma2;
245 G4double border;
246
247 if(posiKinEnergy < 500*MeV) {
248 border = 1. - (electron_mass_c2)/(2*(posiKinEnergy + electron_mass_c2));
249 } else {
250 border = 1. - (100*electron_mass_c2)/(2*(posiKinEnergy + electron_mass_c2));
251 }
252 border = std::min(border, 0.9999);
253
254 CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
255
256 // Case at rest
257 if(posiKinEnergy == 0.0) {
259 G4double sint = sqrt((1. - cost)*(1. + cost));
261 G4ThreeVector dir(sint*cos(phi), sint*sin(phi), cost);
262 phi = twopi * rndmEngine->flat();
263 G4double cosphi = cos(phi);
264 G4double sinphi = sin(phi);
265 G4ThreeVector pol(cosphi, sinphi, 0.0);
266 pol.rotateUz(dir);
268 aGamma1->SetPolarization(pol.x(),pol.y(),pol.z());
270 pol.set(-sinphi, cosphi, 0.0);
271 pol.rotateUz(dir);
272 aGamma2->SetPolarization(pol.x(),pol.y(),pol.z());
273
274 } else {
275
277
278 // (A.A.) LIMITS FOR 1st GAMMA
279 G4double xmin = 0.01;
281
282 G4double d1, d0, x1, x2, dmax, x2min;
283
284 // (A.A.) sampling of x1 x2 x3 (whole cycle of rejection)
285 do {
286 x1 = 1/((1/xmin) - ((1/xmin)-(1/xmax))*rndmEngine->flat());
287 dmax = ComputeFS(posiKinEnergy, x1,1.-x1,border);
288 x2min = 1.-x1;
289 x2 = 1 - rndmEngine->flat()*(1-x2min);
290 d1 = dmax*rndmEngine->flat();
291 d0 = ComputeFS(posiKinEnergy,x1,x2,2-x1-x2);
292 }
293 while(d0 < d1);
294
295 G4double x3 = 2 - x1 - x2;
296
297 //
298 // angles between Gammas
299 //
300
301 G4double psi13 = 2*asin(sqrt(std::abs((x1+x3-1)/(x1*x3))));
302 G4double psi12 = 2*asin(sqrt(std::abs((x1+x2-1)/(x1*x2))));
303
304 // sin^t
305
306 //G4double phi = twopi * rndmEngine->flat();
307 //G4double psi = acos(x3); // Angle of the plane
308
309 //
310 // kinematic of the created pair
311 //
312
313 G4double TotalAvailableEnergy = posiKinEnergy + 2.0*electron_mass_c2;
314
315 G4double phot1Energy = 0.5*x1*TotalAvailableEnergy;
316 G4double phot2Energy = 0.5*x2*TotalAvailableEnergy;
317 G4double phot3Energy = 0.5*x3*TotalAvailableEnergy;
318
319
320 // DIRECTIONS
321
322 // The azimuthal angles of ql and q3 with respect to some plane
323 // through the beam axis are generated at random.
324
325 G4ThreeVector phot1Direction(0, 0, 1);
326 G4ThreeVector phot2Direction(0, sin(psi12), cos(psi12));
327 G4ThreeVector phot3Direction(0, sin(psi13), cos(psi13));
328
329 phot1Direction.rotateUz(posiDirection);
330 phot2Direction.rotateUz(posiDirection);
331 phot3Direction.rotateUz(posiDirection);
332
336
337
338 //POLARIZATION - ???
339 /*
340
341
342 phi = twopi * rndmEngine->flat();
343 G4double cosphi = cos(phi);
344 G4double sinphi = sin(phi);
345 G4ThreeVector pol(cosphi, sinphi, 0.0);
346 pol.rotateUz(phot1Direction);
347 aGamma1->SetPolarization(pol.x(),pol.y(),pol.z());
348
349 G4double phot2Energy =(1.-epsil)*TotalAvailableEnergy;
350 G4double posiP= sqrt(posiKinEnergy*(posiKinEnergy+2.*electron_mass_c2));
351 G4ThreeVector dir = posiDirection*posiP - phot1Direction*phot1Energy;
352 G4ThreeVector phot2Direction = dir.unit();
353
354 // create G4DynamicParticle object for the particle2
355 aGamma2 = new G4DynamicParticle (theGamma,phot2Direction, phot2Energy);
356
357 //!!! likely problematic direction to be checked
358 pol.set(-sinphi, cosphi, 0.0);
359 pol.rotateUz(phot1Direction);
360 cost = pol*phot2Direction;
361 pol -= cost*phot2Direction;
362 pol = pol.unit();
363 aGamma2->SetPolarization(pol.x(),pol.y(),pol.z());
364
365 */
366
367 }
368 /*
369 G4cout << "Annihilation in fly: e0= " << posiKinEnergy
370 << " m= " << electron_mass_c2
371 << " e1= " << phot1Energy
372 << " e2= " << phot2Energy << " dir= " << dir
373 << " -> " << phot1Direction << " "
374 << phot2Direction << G4endl;
375 */
376
377 vdp->push_back(aGamma1);
378 vdp->push_back(aGamma2);
379 if(aGamma3 != nullptr) { vdp->push_back(aGamma3); }
380
381 // kill primary positron
382 fParticleChange->SetProposedKineticEnergy(0.0);
384}
Definition: ThreeVector.h:36
Definition: RandomEngine.h:53
virtual double flat()=0
Definition: G4DynamicParticle.hh:65
void SetPolarization(const G4ThreeVector &)
const G4ThreeVector & GetMomentumDirection() const
G4double GetKineticEnergy() const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForGamma.hh:106
void ProposeTrackStatus(G4TrackStatus status)
G4double ComputeFS(G4double fr1, G4double fr2, G4double fr3, G4double kinEnergy)
Definition: G4eplusTo3GammaOKVIModel.cc:162
Referenced by G4eplusTo2GammaOKVIModel::SampleSecondaries().
◆ SetDelta()
|
inline |
Definition at line 93 of file G4eplusTo3GammaOKVIModel.hh.
93{ if(val > 0.0) { fDelta = val; } };
Referenced by G4eplusTo2GammaOKVIModel::Initialise(), and G4eplusTo2GammaOKVIModel::SampleSecondaries().
The documentation for this class was generated from the following files:
