Monte Carlo simulations for the evaluation of various influence factors on projections in computed tomography

Abstract

This paper presents Monte Carlo simulations considering all stages of the creation process of two-dimensional projections in a computed tomography (CT) device: excitation of angle dependent X-ray spectra within the X-ray tube using results from a previous study [Chyba et al. (2008). Powder Diffr. 23, 150–153]; interaction of these X-rays and secondary photoelectrons with a simple inhomogeneous sample; and interaction of X-rays and photoelectrons with the components (thin layers) of a matrix scintillation detector. The simulations were carried out by using custom software running on up to 50 nodes of a computer cluster. Comparative calculations were also made by using the software package MCNP [Booth et al. (2003). MCNP—A general Monte Carlo N-particle transport code, Report LAUR 03-1987, Los Alamos National Laboratory, Los Alamos, NM]. Tube spectra were calculated with algorithms proposed by Ebel [(2006). Adv. X-Ray Anal. 49, 267–273]. Measurements for the chosen setup made with an available CT device were in relatively good agreement with calculated results. It was shown that good knowledge of the tube spectra is of importance, but most differences between resulting projections and measurements are caused by uncertainties concerning detector response due to light yield of the scintillator and to internal scattering effects within the thin detector layers which lead to spreading of a detected point signal within the detector matrix into neighboring matrix elements.