next up previous index
Next: PHYS441 Simulation of Up: Method Previous: Method

Parameterisation of energy loss and total cross-section

The cross-sections from [] have been used to calculate data points for (gif) and (gif). These have been parameterised as:

σ(Z,T,kc) = Z[Z+ ξσ(1+ γlnZ)][ln{kmaxkc}]αFσ(Z,X,Y) in barn

Elossbrem(Z,T,kc) = Z[Z+ξl(1+δlnZ)][ln{kcE}]βFl(Z,X,Y) in GeV barn

where kmax is the maximum possible value of the photon energy. The functions Fi(Z,X,Y) (i=σ,l ) are polynomials:
Fi(Z,X,Y) = (C1+C2X+...+C6X5)+(C7+C8X+...+C12X5)Y

+ ...+ (C31+C32X+...+C36X5)Y5

+ Z[(C37+C38X+...+C40X3)+(C41+C42X+ ...+C44X3)Y

+ ...+ (C48+C49X+...+C52X3)Y3)]

A least-squares fit has been performed on more than 2000 points for Z = 1, 6, 13, 26, 50, 82, 92 and 1 GeV ≤T ≤ 10 TeV and 10 keV ≤kc≤T . The resulting values of ξσ,l , γ , α ,Ci , ξl , δ and β

can be found in the DATA statements within the functions GBRSGM and GBRELM which compute formulae (gif) and (gif) respectively.

The accuracy of the fit can be estimated as:

The contribution of the bremsstrahlung to the total energy loss of the muons is less than 1% for T ≤ 5 GeV.

When kc≥kmax , a parameterisation different from (gif) can be used for the total muon energy loss due to bremsstrahlung:
Elossbrem(Z,T) = Elossbrem(Z,T,k=kmax)

= Z(Z+1) kmax[d1+(d2X+d3Y) + (d4X2+d5XY+d6Y2)

+ ...+(d22X6+d23X5Y+...+d28Y6)]

where Y=Z1/3 . The accuracy of the formula (gif) for 1 ≤Z ≤100

is rather good:

{ΔElossbremElossbrem}=
≤1.5 if T > 1  GeV

≤1 if T > 5  GeV

.

PHYS441



next up previous index
Next: PHYS441 Simulation of Up: Method Previous: Method


Janne Saarela
Mon Apr 3 12:46:29 METDST 1995