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.

Final<br />

Direction<br />

x<br />

φ<br />

Θ<br />

t( (1 − ζ) K1)<br />

∆x<br />

t(K1)<br />

t(ζ K1)<br />

Θ<br />

∆y<br />

z<br />

y<br />

s<br />

Figure 2.10: Modified random hinge transport mechanics schematic.<br />

Thus, instead of the hinge consisting of steps ζt and (1 − ζ)t, it consists of steps of the distances<br />

corresponding to the accumulation of scattering strengths ζK 1 (t) and (1 − ζ)K 1 (t). In this way,<br />

the total scattering strength K 1 (t) over the step is preserved, and the average location of the hinge<br />

will correspond to that location which preserves the first moments of spatial displacements. The<br />

modified random hinge transport mechanics of <strong>EGS5</strong>, illustrated in Figure 2.10, is therefore capable<br />

of simulating accurately the average final positions of electrons moving long pathlengths through<br />

materials, even when energy loss occurs continuously along the track. The length of an electron<br />

step is therefore limited only by the accuracy of the multiple scattering p.d.f., any need to higher<br />

order spatial moments (as would be required for some differential tallies), and the accuracy of the<br />

method used to compute the integral in Equation 2.372 along the paths between hinge points.<br />

Values of the scattering powers G 1 for electrons and positrons as a function of energy are<br />

computed in PEGS in a new function G1E. For kinetic energies less than 100 MeV, values are taken<br />

from cubic spline fits of data provided in the DCSLIB package described in the previous section.<br />

Above 100 MeV, the Wentzel-shaped scattering cross section used by Molière is assumed, in which<br />

case Equation 2.367 for l = 1 can be integrated analytically to yield<br />

G 1 (E) = N [ ( )]<br />

aρ 2πZ 2 e 4 π 2 + χ 2 a<br />

M p 2 v 2 ln<br />

χ 2 −<br />

π2<br />

a π 2 + χ 2 , (2.373)<br />

a<br />

where χ 2 a is taken from Equation 2.322. Discontinuities at 100 MeV between values of G 1 computed<br />

from the partial-wave cross section data of DCSLIB and from the Wentzel cross section above are<br />

usually small (on the order of fractions of 1 percent), and hence are smoothed only crudely.<br />

Since the Goudsmit-Saunderson multiple scattering distributions are pre-computed in PEGS,<br />

100

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

Saved successfully!

Ooh no, something went wrong!