Subroutines

CALL GDRELA

GDRELA initialises the ionisation energy-loss tables for different materials and particles. The energy binning is in the array ELOW (/GCMULO/) initialised in the routine GPHYSI. The tables are filled with the quantity $dE/dx$ in GeV $cm-1$ , for elements, mixtures and compounds. A temporary table is also filled, containing additional points for integration, to obtain the range-energy table (see [PHYS010]).

GDRELP computes the ionisation energy loss for protons, and the values are stored in the energy loss table. This value is used for other charged hadrons with scaled energy bin.

GDRELM computes the ionisation energy loss for muons, and sums it in the muon table with the energy loss due to the other processes (bremsstrahlung, pair production and nuclear interactions).

The following pointers are used:

JMA = LQ(JMATE-I)

$Ith$ material bank;

JEL1 = LQ(JMA-1)

$dE/dx$ for / (see [PHYS330]);

JEL2 = LQ(JMA-2)

$dE/dx$ for $\mu +/\mu -$ ;

JEL3 = LQ(JMA-3)

$dE/dx$ for protons (used also for other charged particles).

GDRELA is called during initialisation by GPHYSI.

CALL GDRELP (A,Z,DENS,T,DEDX*)

GDRELP computes the ionisation energy loss for protons. It is called by the routine GDRELA. See GDRELX below for the meaning of the parameters.

CALL GDRELM (A,Z,DENS,T,DEDX*)

GDRELM computes the energy loss due to ionisation for muons. It is called by the routine GDRELA. See GDRELX below for the meaning of the parameters.

CALL GDRELX (A,Z,DENS,T,HMASS,DEDX*)

A

( REAL) mass number of the material;

Z

( REAL) atomic number of the material;

DENS

( REAL) density of the material in g $cm-3$ ; it is only used to calculate the density correction in case of high values of T;

T

( REAL) kinetic energy of the particle in GeV;

HMASS

( REAL) mass of the particle;

DEDX

( REAL) energy loss of the particle in GeV $cm2$ $g-1$ .

This routine calculates the ionisation energy loss for a muon or a proton. It is called by GDRELM and GDRELP.

CALL GDRSGA

CDRSGA calculates the total cross-section for the production of $\delta$ -rays by electrons (Möller scattering), positrons (Bhabha scattering) and muons. For hadrons, this value is calculated at tracking time. The mean free path as a function of the energy is stored for every medium. The following pointers are used:

JMA = LQ(JMATE-I)

pointer to the $Ith$ material;

JDRAY = LQ(JMA-11)

$\delta$ -ray cross-section for electrons;

JDRAY+NEK1

$\delta$ -ray cross-section for positrons;

JDRAY+2*NEK1

$\delta$ -ray cross-section for muons.

The routine is called during initialisation by GPHYSI.