Abstract:
In source-detector radiation transport simulations, pulse height distributions are a useful metric in assessing the effectiveness of nuclear instrumentation. In the area of spectroscopy, pulse height distributions are used to identify an unknown source. It is widely believed that pulse height distributions cannot be created using deterministic methods. This quantity is routinely calculated in Monte Carlo transport codes because each history is individually tracked and the amount of energy deposited by interaction events is easily tallied. We have developed a methodology to calculate the pulse height distribution using deterministic transport codes. The algorithm utilizes the scattered components of scalar flux which are readily calculated by iterating on the scattering source. The approach is applicable in three dimensions but we present here results from an implementation in 1D slab geometry with anisotropic scattering/angular quadrature orders up to P14/S32.
