28.11.2014 Views

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

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.

defined by<br />

δ ij = 1 if i = j,<br />

0 otherwise. (2.42)<br />

Using this notation, we start with the following formulas for the bremsstrahlung and pair production<br />

differential cross sections: 5 .<br />

dσ Brem (Z, Ĕ0, ˘k)<br />

= A′ (Z, Ĕ0)r0 2 αZ(Z + ξ(Z))<br />

{ d˘k<br />

˘k<br />

(<br />

× 1 + (Ĕ/Ĕ0) 2) [ φ 1 (δ) − 4 ]<br />

3 ln Z − (4f c(Z) if Ĕ 0 > 50, 0)<br />

and<br />

− 2 3<br />

( ) [ Ĕ<br />

φ 2 (δ) − 4 ] }<br />

Ĕ 0<br />

3 ln Z − (4 f c(Z) if Ĕ 0 > 50, 0)<br />

(2.43)<br />

dσ P air (Z, ˘k, Ĕ + )<br />

= A′ p (Z, ˘k)r0 2 αZ(Z + ξ(Z))<br />

dĔ+<br />

˘k<br />

{ 3<br />

(Ĕ2<br />

× + + −) [<br />

Ĕ2 φ 1 (δ) − 4 ]<br />

3 ln Z − (4f c(Z) if ˘k > 50, 0)<br />

+ 2 [<br />

3 Ĕ+Ĕ− φ 2 (δ) − 4 ]}<br />

3 ln Z − (4 f c(Z) if ˘k > 50, 0)<br />

(2.44)<br />

where<br />

and<br />

∆ =<br />

=<br />

˘km<br />

2Ĕ0Ĕ<br />

˘km<br />

2Ĕ+Ĕ−<br />

δ = 136 Z −1/3 2∆ (2.45)<br />

(for bremsstrahlung) (2.46)<br />

(for pair production) . (2.47)<br />

To avoid confusion it should be noted that our δ is the same as the δ of Butcher and Messel[39] but<br />

we use φ i (δ) to denote their f i (δ). Rossi[141] and Koch and Motz[91] use a variable γ = 100<br />

136δ. Also,<br />

note that our φ i (δ) has the same value as the φ i (γ) of Koch and Motz[91] (e.g., see their Figure 1)<br />

provided “our δ” = 136<br />

100 times “their γ.” For arbitrary screening, φ 1 and φ 2 are given by<br />

φ 1 (δ) = 4<br />

φ 2 (δ) = 4<br />

∫ 1<br />

∆<br />

∫ 1<br />

∆<br />

(q − ∆) 2 [1 − F (q, Z)] 2 dq<br />

q 3 + 4 + 4 ln Z , (2.48)<br />

3<br />

[<br />

( ]<br />

q<br />

q 3 − 6∆ 2 q ln + 3∆<br />

∆) 2 q − 4∆ 3<br />

× [1 − F (q, Z)] 2 dq<br />

q 4 + 10 3 + 4 ln Z (2.49)<br />

3<br />

5 As described earlier, we have adopted the notation that cupped energy variables (e.g., Ĕ) are in MeV and cupped<br />

distance variables (e.g., ˘t) are in radiation lengths. Uncupped variables are in CGS units and constants are defined<br />

explicitly<br />

39

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

Saved successfully!

Ooh no, something went wrong!