An important step in assessing the quality of an image reconstruction algorithm is the simulation of the medical imaging process. For that purpose, the patient's anatomical structure is substituted in general by more or less simple geometrical objects, as, e.g., the Shepp-Logan phantom. Furthermore, the attenuation of the human body and thus the resulting detector image (e.g., the sinogram in CT) is often computed by integrating the attenuation coefficient along various rays without considering the contribution of scattered photons in the detector signal. We therefore decided to improve the simulation by using an existing Monte Carlo code (EGSnrc) to model the transport of numerous photons from the x-ray tube through the body to the detector. The deflection of photons and creation of secondary particles in scattering events occurs naturally in this program, but can also be avoided artificially. Besides the improved simulation of the irradiation process, this allows us to quantify the amount of scattered radiation in the detector image. The patient is represented by a so-called voxel phantom, which is based on tomographic image data of a real person, adopted to represent the ICRP Reference Man. Our improved modeling process is being applied to determine the amount of scatter radiation in helical multi-slice CT of the thorax compared to a planned circular CT with large flat panel detectors. The new reconstruction algorithm OPED (orthogonal polynomial expansion on disc), developed at GSF and the University of Oregon, might reduce the scatter radiation considerably.