22.05.2014 Views

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Geant 3.16 GEANT User’s Guide PHYS210<br />

Origin : L. Urbán Submitted: 10.11.84<br />

Revision : Revised: 16.12.93<br />

Documentation :<br />

1 Subroutines<br />

Total cross-section for e+e- pair production by photons<br />

CALL GPRSGA<br />

GPRSGA calculates the total cross-section for the pair production of photons and the direct pair production<br />

of muons in all materials. It tabulates the mean free path, λ =1/Σ (in cm), as a function of the medium and<br />

the energy. For the pair production by photons, it calls the function GPRSGG. The energy bins are the one of<br />

the array ELOW in common /CGMULO/ initialised in the routine GPHYSI. The following pointers are used:<br />

JMA = LQ(JMATE-I)<br />

JPAIR = LQ(JMA-10)<br />

pointer to the I th material<br />

pointer to pair production cross-sections<br />

GPRSGA is called at initialisation time by GPHYSI.<br />

VALUE = GPRSGG (Z,E)<br />

GPRSGG calculates the total cross-section for the pair production of a photon with energy E in material with<br />

atomic number Z. It is called by GPRSGA.<br />

2 Method<br />

We have parameterised the total cross-section of e + e − -pair production by photon as:<br />

σ(Z, E γ )=Z(Z +1)(F 1 (X)+F 2 (X)Z + F 3 (X)/Z) (barn/atom) (1)<br />

where<br />

X = ln(E γ /m),<br />

m = electron mass,<br />

E γ = photon energy<br />

and<br />

F i (X) =<br />

5∑<br />

c n X n<br />

n=0<br />

The parameters c n were taken from a least-square fit to the data [9] and can be found in the DATA statement<br />

in the function GPRSGG(Z,E) which computes the formula (1) (in barn/atom).<br />

This parameterisation gives a good description of the data in the range:<br />

1 ≤ Z ≤ 100<br />

1.5 MeV ≤ E γ ≤ 100 GeV<br />

}<br />

∆ σ<br />

σ<br />

≤ 5% with a mean value of ≈ 2.2% (2)<br />

This mean free path is tabulated at initialisation as a function of the medium and of the energy by the routine<br />

GPRSGA.<br />

211 PHYS210 – 1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!