Four separate geometries were simulated in MCNPX. The first was a bare HPGe crystal (with the aluminum housing) 10 cm away from a mixed Cs-137 and Co-60 point source. The second (side scattering) was with 1 cm thick lead cylinder encasing the air between the source and the detector. The third (all scatter) was with a 1 cm thick lead cylinder on the sides, and a 1 cm thick slab below the source. The final was with a 2" (5.08 cm) thick (all scatter thick bottom) slab below the source and a 1 cm thick side scattering cylinder. 1 cm of lead was chosen because 1.07 cm is the half-thickness for lead (Co-60 source).
![]() |
| X-Z profile of side scattering |
![]() |
| X-Z profile of back scattering and side scatter (all scattering) |
This also provides cute little examples of solid angle calculations. The simplest fraction of the solid angle the detector subtends is the ratio of the area of the detector to the area of the sphere of radius the distance from the source to the detector is. For example since the simulated detector is a cylinder along the z-axis with a radius of 2.54 cm, and is 10 cm away from the point source, as is shown below.
\begin{align}
\eta &= \frac{A_{detector}}{A_{Sphere}} \newline
&= \frac{\pi(2.54cm)^2}{4\pi(10cm)^2} \newline
&= 0.016
\end{align}
A more accurate solid angle fraction can be calculated as follows:
\begin{align}
\eta &= \frac{2\pi\left(1-\cos(\theta)\right)}{4\pi}\newline
&=\frac{2\pi\left(1-\cos \left(\tan^{-1}\frac{2.54}{10}\right)\right)}{4\pi} \newline
&= 0.015
\end{align}
Okay, so we do get around a 30% increase in the incident flux with reflection. This was higher than I expected. In addition, we don't take that much of a hit on our detector response, as shown below. I used a gaussian energy broadening on the pulse height tally in order to accurate simulate the HPGe response, which is something that I haven't done before and was is really nifty. What could be super nifty would be to compare this to a measurement, but I don't think that is gonna happen.



