Next: Common to All Charged Up: Gamma Incident Previous: Gamma Conversion into an   Contents

Subsections

# Gamma Conversion into a Muon - Anti-muon Pair

The class G4GammaConversionToMuons simulates the process of gamma conversion into muon pairs. Given the photon energy and and of the material in which the photon converts, the probability for the conversions to take place is calculated according to a parameterized total cross section. Next, the sharing of the photon energy between the and is determined. Finally, the directions of the muons are generated. Details of the implementation are given below and can be also found in[1].

## Cross Section and Energy Sharing

In the field of the nucleus, muon pair production on atomic electrons, , has a threshold of . Up to several hundred GeV this process has a much lower cross section than the corresponding process on the nucleus. At higher energies, the cross section on atomic electrons represents a correction of to the total cross section.

For the approximately elastic scattering considered here, momentum, but no energy, is transferred to the nucleon. The photon energy is fully shared by the two muons according to

 (6.39)

or in terms of energy fractions

The differential cross section for electromagnetic pair creation of muons in terms of the energy fractions of the muons is
 (6.40)

where is the charge of the nucleus, is the classical radius of the particles which are pair produced (here muons) and
 (6.41)

where

 (6.42)

These formulae are obtained from the differential cross section for muon bremsstrahlung [2] by means of crossing relations. The formulae take into account the screening of the field of the nucleus by the atomic electrons in the Thomas-Fermi model, as well as the finite size of the nucleus, which is essential for the problem under consideration. The above parameterization gives good results for . The fact that it is approximate close to threshold is of little practical importance. Close to threshold, the cross section is small and the few low energy muons produced will not travel very far. The cross section calculated from Eq.(6.40) is positive for and
 (6.43)

except for very asymmetric pair-production, close to threshold, which can easily be taken care of by explicitly setting whenever .

Note that the differential cross section is symmetric in and and that

where stands for either or . By defining a constant
 (6.44)

the differential cross section Eq.(6.40) can be rewritten as a normalized and symmetric as function of :
 (6.45)

This is shown in Fig.6.1 for several elements and a wide range of photon energies. The asymptotic differential cross section for

is also shown.

## Parameterization of the Total Cross Section

The total cross section is obtained by integration of the differential cross section Eq.(6.40), that is
 (6.46)

is a function of () and () of the element (see Eq.(6.41)). Numerical values of are given in Table6.1.

 W for H W for Be W for Cu W for Pb GeV 1 2.11 1.594 1.3505 5.212 10 19.4 10.85 6.803 43.53 100 191.5 102.3 60.10 332.7 1000 1803 919.3 493.3 1476.1 10000 11427 4671 1824 1028.1 28087 8549 2607 1339.8

Values of the total cross section obtained by numerical integration are listed in Table6.2 for four different elements. Units are in , where .

 , H , Be , Cu , Pb GeV 1 0.01559 0.1515 5.047 30.22 10 0.09720 1.209 49.56 334.6 100 0.1921 2.660 121.7 886.4 1000 0.2873 4.155 197.6 1476 10000 0.3715 5.392 253.7 1880 0.4319 6.108 279.0 2042

Well above threshold, the total cross section rises about linearly in with the slope

 (6.47)

until it saturates due to screening at . Fig.6.2 shows the normalized cross section where
 (6.48)

Numerical values of are listed in Table6.3.

 Element 1/GeV H 0.963169 Be 0.514712 Cu 0.303763 Pb 0.220771

The total cross section can be parameterized as

 (6.49)

with
 (6.50)

and

The threshold behavior in the cross section was found to be well approximated by and the saturation by . The agreement at lower energies is improved using an empirical correction factor, applied to the slope , of the form

where

A comparison of the parameterized cross section with the numerical integration of the exact cross section shows that the accuracy of the parametrization is better than 2%, as seen in Fig.6.3.

## Multi-differential Cross Section and Angular Variables

The angular distributions are based on the multi-differential cross section for lepton pair production in the field of the Coulomb center

 (6.51)

Here
 (6.52)

where
 (6.53)

is the square of the momentum transferred to the target and and are the squares of the components of the vector , which are parallel and perpendicular to the initial photon momentum, respectively. The minimum momentum transfer is .

The muon vectors have the components

 (6.54)

where . The initial photon direction is taken as the -axis. The cross section of Eq.(6.51) does not depend on . Because of azimuthal symmetry, can simply be sampled at random in the interval .

Eq.(6.51) is too complicated for efficient Monte Carlo generation. To simplify, the cross section is rewritten to be symmetric in , using a new variable and small parameters , where and . When higher powers in small parameters are dropped, the differential cross section in terms of becomes

 (6.55)

where, in this approximation,

For Monte Carlo generation, it is convenient to replace () by the polar coordinates () with and . Integrating Eq. 6.55 over and using symbolically where yields
 (6.56)

Integration with logarithmic accuracy over gives
 (6.57)

Within the logarithmic accuracy, can be replaced by , so that
 (6.58)

Making the substitution , gives
 (6.59)

Atomic screening and the finite nuclear radius may be taken into account by multiplying the differential cross section determined by Eq.(6.56) with the factor
 (6.60)

where and are atomic and nuclear form factors. Please note that after integrating Eq. 6.56 over , the -dependence is lost.

## Procedure for the Generation of Muon - Anti-muon Pairs

Given the photon energy and and of the material in which the converts, the probability for the conversions to take place is calculated according to the parametrized total cross section Eq.(6.49). The next step, determining how the photon energy is shared between the and , is done by generating according to Eq.(6.40). The directions of the muons are then generated via the auxilliary variables . In more detail, the final state is generated by the following five steps, in which are random numbers with a flat distribution in the interval [0,1]. The generation proceeds as follows.

1) Sampling of the positive muon energy .
This is done using the rejection technique. is first sampled from a flat distribution within kinematic limits using

and then brought to the shape of Eq.(6.40) by keeping all which satisfy

Here is the maximum value of , obtained for symmetric pair production at . About 60% of the events are kept in this step. Results of a Monte Carlo generation of are illustrated in Fig.6.4. The shape of the histograms agrees with the differential cross section illustrated in Fig.6.1.

2) Generate .
The distribution in is obtained from Eq.(6.59) as
 (6.61)

with form factors taken into account by
 (6.62)

In the interval considered, the function will always be bounded from above by

For small and large , approaches unity, as shown in Fig.6.5.

The Monte Carlo generation is done using the rejection technique. About 70% of the generated numbers are kept in this step. Generated -distributions are shown in Fig.6.6.

3) Generate by the rejection technique using generated in the previous step for the frequency distribution

 (6.63)

The maximum of is
 (6.64)

Generated distributions in are shown in Fig.6.7.

4) Generate .
The distribution in has the form
 (6.65)

where
 (6.66)

and
 (6.67)

The distribution is obtained by a direct transformation applied to uniform random numbers according to
 (6.68)

where
 (6.69)

Generated distributions of are shown in Fig.6.8

5) Calculate and from with

 (6.70)

according to
 (6.71)

The muon vectors can now be constructed from Eq.(6.54), where is chosen randomly between 0 and . Fig.6.9 shows distributions of at different photon energies (in beryllium). The spectra peak around as expected.

The most probable values are . In the small angle approximation used here, the values of and can in principle be any positive value from 0 to . In the simulation, this may lead (with a very small probability, of the order of ) to unphysical events in which or is greater than . To avoid this, a limiting angle is introduced, and the angular sampling repeated, whenever .

Figs.6.10,6.11 and 6.12 show distributions of the simulated angular characteristics of muon pairs in comparison with results of exact calculations. The latter were obtained by means of numerical integration of the squared matrix elements with respective nuclear and atomic form factors. All these calculations were made for iron, with and . As seen from Fig.6.10, wide angle pairs (at low values of the argument in the figure) are suppressed in comparison with the Coulomb center approximation. This is due to the influence of the finite nuclear size which is comparable to the inverse mass of the muon. Typical angles of particle emission are of the order of (Fig.6.11). Fig.6.12 illustrates the influence of the momentum transferred to the target on the angular characteristics of the produced pair. In the frame of the often used model which neglects target recoil, the pair particles would be symmetric in transverse momenta, and coplanar with the initial photon.

## Status of this document

28.05.02 created by H.Burkhardt.
01.12.02 re-worded by D.H. Wright

## Bibliography

1. H. Burkhardt, S. Kelner, and R. Kokoulin, Monte Carlo Generator for Muon Pair Production''. CERN-SL-2002-016 (AP) and CLIC Note 511, May 2002.

2. S. R. Kelner, R. P. Kokoulin, and A. A. Petrukhin, About cross section for high energy muon bremsstrahlung,''. Moscow Phys. Eng. Inst. 024-95, 1995. 31pp.

Next: Common to All Charged Up: Gamma Incident Previous: Gamma Conversion into an   Contents