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.

tally, as in EGS4. Depositing energy continually in this manner, rather than only in discrete,<br />

large, chunks at the hinge points, obviates requiring very small values of DEINITIAL to attain<br />

small statistical uncertainty in most energy deposition problems. It does introduce some artifacts,<br />

however, in that while translating a particle to an energy hinge point, energy is being deposited<br />

into the medium without being decremented from the particle. Similarly, along the track segment<br />

between the hinge and the end of the energy transport step (the part described by DERESID), the<br />

full transport step energy loss will have been decremented from the electron, but not all of that will<br />

have been deposited via continuous stopping loss until the end of the transport step. Thus energy<br />

is not strictly conserved, in that a particle’s total energy plus the energy of its daughter particles<br />

plus the energy that it has deposited do not sum to the particle’s initial energy except at the exact<br />

transitions between hinge steps (which are not stopping points on the translation steps). This has<br />

some interesting consequences. For instance, if an electron escapes the problem geometry before<br />

reaching its next energy hinge, a hinge must be imposed at the boundary to and the amount of<br />

energy deposited prior to the escape of the particle much be decremented from the particle before<br />

it is tallied as an escaped particle. Similarly, if an electron escapes during the residual part of a<br />

previous hinge (i.e., before the full DERESID has been deposited), its kinetic energy must adjusted<br />

upward to account for the fact that the full hinge energy has already been decremented from the<br />

particle, but a portion of it has yet to be deposited, because the end of the hinge step had yet to<br />

be reached. See Figure 2.14 for a schematic illustrating this problem.<br />

Likewise, when a hard collision occurs, a particle will be somewhere in the midst of some<br />

transport step, on either the residual side of the previous energy hinge or the initial side of the<br />

upcoming hinge, and so the kinetic energy available for the collision must be adjusted to account<br />

for the energy deposited during the translation step. If the electron has yet to reach the end of<br />

the previous hinge still (i.e., it’s still on the residual side), we must adjust the particle kinetic<br />

energy by increasing it to account for energy decremented but not deposited before proceeding the<br />

collision analysis. And, as with the case of escaping particles, if the electron has passed the end of<br />

the previous hinge step and is on the initial part of the next hinge step, we must impose a hinge<br />

and decrement the already deposited energy from the electron. These adjustments are necessary<br />

not only to preserve energy deposition, but also to determine the kinetic energy available for the<br />

hard collision. In terminating hinges before hard collisions, we adjust the electron energy exactly<br />

as was done for escaping electrons, by setting E = E − DEINITIAL + DENSTEP.<br />

A final interesting problem occurs at final electron energy hinge. Because the electron energy<br />

is set to be equivalent to ECUT at the final hinge point, which is reached before the electron has<br />

been translated through the DERESID portion of its transport step, some provision must be made<br />

for tracking electrons during the residual parts of their final hinge steps. This requirement is<br />

complicated somewhat in that such electrons will have total energy which is exactly equal to the<br />

problem cutoff energy. In the present <strong>EGS5</strong> algorithm, when the last energy hinge point is reached,<br />

the electron energy is set to ECUT, DENSTEP is set to DERESID, and DEINITIAL is set to zero (DERESID<br />

is then set to zero as well). The multiple scattering hinge step at this point is set to infinity, as<br />

is the distance to the next hard collision. The particle is thus transported linearly through the<br />

distance corresponding to DERESID, depositing the appropriate CSDA energy along the way.<br />

106

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

Saved successfully!

Ooh no, something went wrong!