next up previous index
Next: Total cross-sections Up: Method Previous: Method

Continuous energy loss

The integration of (gif) leads to the Berger-Seltzer formulae [,,,,,] :

{dEdx}= {2 πr02mn β2}[ln{2(τ+ 2)(I/m)2}+ F±(τ, Δ)- δ],

where

The functions F± are given by
F+(τ,Δ) = ln(τΔ) -{Δ2τ}[τ+ 2 Δ-{3Δ2y 2}-(Δ- {Δ33}) y2- ({Δ22}- τ{Δ33}+ {Δ44}) y3]

F-(τ,Δ) = -1 -β2+ln[(τ- Δ)Δ]+ {ττ-Δ}+{[{Δ22}+ ( 2τ+1) ln(1- {Δτ}) ]γ2},

where y = 1/(γ+1) .

The density effect correction is calculated as in []: δ=
0 if x<x0

2ln10 + x+C+a(x1-x)m if x0≤x≤x1

2ln10 + x+C if x1<x,

.

where x = ln(γ2-1) / (2ln10)

The quantities n, I and the parameters of the density effect correction (x0, x1, C, a, m ) are computed in the routine GPROBI, and we give the corresponding formulae here. The electron density of the medium, n, can be written as

n =
NAvρ{ZA} for elements

NAvρ{∑iniZiiniAi} for compounds/mixtures,

.

where

N
Avogadro's number;
Zi

atomic number;
Ai

atomic weight;
ρ

density of the material;
ni

proportion by number of the ith element in the material (for a mixture ni= W pi/Ai

where pi the proportion by weight and W is the molecular weight).

The average mean ionisation energy can be calculated as [] [] [] []:

I (GeV) =
16 10-9 Z0.9 for a chemical element

exp[ ∑iniZilnI(Zi) / ∑iniZ i] for a compound or mixture

.

The density effect correction parameters can be computed (for condensed medium []) as
C = 1 + 2ln{I28.8 10-9ρ∑ niZi/ ∑ niAi}

m = 3

xa = {C2ln10}

a = 2ln10 {(xa-x0)(x1- x0)}m

I C x0 x1

<10-7 <3.681 0.2 2

<10-7 ≥3.681 -0.326C-1 2

≥10-7 5.215 0.2 3

≥10-7 ≥5.215 -0.326C-1.5 3



next up previous index
Next: Total cross-sections Up: Method Previous: Method


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