Livermore Model

G4LivermoreGammaConversionModel and G4LivermoreGammaConversion5DModel are available.

Total cross-section

The total cross-section of the Gamma Conversion process is based on either EPDL97 or EPICS2017 data as described in Low Energy Livermore Model, and determined from the data as described in Generic Calculation of Total Cross Sections.

Sampling of the final state

For low energy incident photons, the simulation of the Gamma Conversion final state is performed according to Final State.

The secondary \(e^\pm\) energies are sampled using the Bethe-Heitler cross-sections with Coulomb correction.

The Bethe-Heitler differential cross-section with the Coulomb correction for a photon of energy \(E\) to produce a pair with one of the particles having energy \(\epsilon E\) (\(\epsilon\) is the fraction of the photon energy carried by one particle of the pair) is given by [FN78]:

(13)\[\frac{d \sigma(Z,E,\epsilon)}{d \epsilon} = \frac{r_0^2\alpha Z (Z + \xi(Z))}{E^2}\left[ (\epsilon^2 + ( 1 -\epsilon)^2) \left( \Phi_1(\delta) - \frac{F(Z)}{2}\right) + \frac{2}{3}\epsilon (1-\epsilon) \left( \Phi_2(\delta) - \frac{F(Z)}{2} \right) \right]\]

where \(\Phi_i(\delta)\) are the screening functions depending on the screening variable \(\delta\) [eal93].

The value of \(\epsilon\) is sampled using composition and rejection Monte Carlo methods [MC70, eal93, BM60].

After the successful sampling of \(\epsilon\), the process generates the polar angles of the electron with respect to an axis defined along the direction of the parent photon. The electron and the positron are assumed to have a symmetric angular distribution. The energy-angle distribution is given by [Tsa74, Tsa77]:

\[\begin{split}\frac{d \sigma}{dp d\Omega} & = \frac{2 \alpha^2 e^2}{\pi k m^4}\left[ \left(\frac{2x(1-x)}{(1+l)}^2 - \frac{12 lx(1-x)}{(1+l)^4} \right)(Z^2+Z) + \right. \\ &+ \left. \left( \frac{ 2x^2 - 2x+1}{(1+l)^2} + \frac{4lx(1-x)}{(1+l)^4} \right) (X-2Z^2 f((\alpha Z)^2)) \right]\end{split}\]

where \(k\) is the photon energy, \(p\) the momentum and \(E\) the energy of the electron of the \(e^\pm\) pair \(x=E/k\) and \(l = E^2\theta^2/m^2\). The sampling of this cross-section is obtained according to [eal93].

The azimuthal angle \(\phi\) is generated isotropically.

This information together with the momentum conservation is used to calculate the momentum vectors of both decay products and to transform them to the Geant4 coordinate system. The choice of which particle in the pair is the electron/positron is made randomly.

Bibliography

BM60

J.C. Butcher and H. Messel. Nucl. Phys., 20(15):, 1960.

eal93(1,2,3)

René Brun et al. GEANT: Detector Description and Simulation Tool; Oct 1994. CERN Program Library. CERN, Geneva, 1993. Long Writeup W5013. URL: https://cds.cern.ch/record/1082634.

FN78

Richard L. Ford and W. Ralph Nelson. The egs code system: computer programs for the monte carlo simulation of electromagnetic cascade showers (version 3). Technical Report, SLAC, 1978.

MC70

H. Messel and D. Crawford. Electron-Photon shower distribution. Pergamon Press, 1970.

Tsa74

Y. Tsai. Pair production and bremsstrahlung of charged leptons. Rev. Mod. Phys., 46():815, 1974.

Tsa77

Y. Tsai. Pair production and bremsstrahlung of charged leptons. Rev. Mod. Phys., 49():421, 1977.