PENELOPE 2003 - OECD Nuclear Energy Agency
PENELOPE 2003 - OECD Nuclear Energy Agency
PENELOPE 2003 - OECD Nuclear Energy Agency
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2.4. Electron-positron pair production 55<br />
electron shell i is sampled with relative probability f i and binding effects are accounted<br />
for by simply rejecting E ′ -values such that E − E ′ < U i .<br />
The azimuthal scattering angle φ of the photon is sampled uniformly in the interval<br />
(0, 2π). We assume that the Compton electron is emitted with energy E e = E − E ′ − U i<br />
in the direction of the momentum transfer vector q = ¯hk − ¯hk ′ , with polar angle θ e and<br />
azimuthal angle φ e = φ + π, relative to the direction of the incident photon. cos θ e is<br />
given by<br />
E − E ′ cos θ<br />
cos θ e = √<br />
E2 + E ′2 − 2EE ′ cos θ . (2.75)<br />
When E ′ = E C , this expression simplifies to<br />
cos θ e = E + m ec 2 (<br />
)<br />
E − E 1/2<br />
C<br />
, (2.76)<br />
E 2m e c 2 + E − E C<br />
which coincides with the result (A.20). Since the active electron shell is known, characteristic<br />
x rays and Auger electrons emitted in the de-excitation of the ionized atom<br />
can also be followed. This is important, for instance, to account for escape peaks in<br />
scintillation or solid state detectors<br />
Table 2.1: Average number n r of random numbers ξ needed to simulate a single incoherent<br />
scattering event for photons with energy E in aluminium, silver and gold.<br />
E (eV) Al Ag Au<br />
10 3 16.6 11.9 13.4<br />
10 4 11.0 11.4 11.5<br />
10 5 9.5 9.8 10.0<br />
10 6 8.2 8.2 8.3<br />
10 7 7.5 7.5 7.5<br />
As a measure of the efficiency of the sampling algorithm, we may consider the average<br />
number n r of random numbers ξ required to simulate an incoherent scattering event.<br />
n r is practically independent of the atomic number and decreases with photon energy<br />
(see table 2.1). The increase of n r at low energies stems from the loss of efficiency<br />
of the algorithm used to sample cos θ. Although the simulation of incoherent events<br />
becomes more laborious as the photon energy decreases, this has only a small influence<br />
on the speed of practical photon transport simulations since low-energy photons interact<br />
predominantly via photoelectric absorption (see fig. 2.10 below).<br />
2.4 Electron-positron pair production<br />
Electron-positron pairs can be created by absorption of a photon in the vicinity of a<br />
massive particle, a nucleus or an electron, which absorbs energy and momentum so that