28.11.2014 Views

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

+ 1 E 2<br />

′ − 1 E 1<br />

′<br />

− C 2 ln E 2E 1<br />

′ ]<br />

E 1 E 2<br />

′<br />

(2.195)<br />

where<br />

and other symbols are the same as in Equation 2.193.<br />

E i = (Ĕi − m)/ ˘T 0 , i = 1, 2 , (2.196)<br />

E ′ i = 1 − E i , i = 1, 2 , (2.197)<br />

The minimum and maximum energies for the scattered electron (by convention, the scattered<br />

electron is taken to be the one which emerges with the lower energy) are A E and ˘T 0<br />

2<br />

+m, respectively.<br />

When these limits are used for Ĕ1 and Ĕ2 in Equation 2.195, we obtain for the total discrete Møller<br />

cross section (which we assume to non-zero only for Ĕ0 > ĔMøller T h<br />

):<br />

˘Σ Møller, Total (Ĕ0) == Equation 2.195 with Ĕ1 = A E and Ĕ2 = ( ˘T 0 /2) + m . (2.198)<br />

PEGS functions AMOLDM, AMOLRM, and AMOLTM evaluate Equations 2.193, 2.195, and 2.198, respectively.<br />

In sampling for the resultant energy, we use the variable E = ˘T / ˘T 0 and obtain<br />

d˘Σ Møller (Ĕ0)<br />

dE<br />

= X 0n2πr 2 0 m<br />

˘T 0 E 0<br />

f(E)g(E) (2.199)<br />

where<br />

E 0 1<br />

f(E) =<br />

1 − 2E 0 E 2 , E ∈ (E 0, 1/2), (2.200)<br />

g(E) = g 1 [1 + g 2 E 2 + r(r − g 3 )], (2.201)<br />

E 0 = T E / ˘T 0 , (2.202)<br />

g 1 = (1 − 2E 0 )/β 2 , (2.203)<br />

g 2 = (γ − 1) 2 /γ 2 , (2.204)<br />

g 3 = (2γ − 1)/γ 2 , (2.205)<br />

r = E/(1 − E) . (2.206)<br />

The sampling procedure is as follows:<br />

1. Compute parameters depending on Ĕ0, but not E: E 0 , g 1 , g 2 , and g 3 .<br />

2. Sample E from f(E) by using<br />

E = T E / [ ˘T 0 − (Ĕ0 − ĔMøller T h<br />

)ζ 1 ]). (2.207)<br />

3. Compute r and the rejection function g(E). If ζ 2 > g(E), reject and return to Step 2.<br />

After the secondary energies have been determined, Equation 2.167 can be used to obtain the<br />

scattering angles, and EGS routine UPHI can be called to select random azimuthal angles and set<br />

up the secondary particles in the usual way.<br />

67

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

Saved successfully!

Ooh no, something went wrong!