The optimal balance between quality and efficiency in proton radiography imaging technique at various proton beam energies: A Monte Carlo study.

Proton radiography is a novel imaging modality that allows direct measurement of the proton energy loss in various tissues. Currently, due to the conversion of so-called Hounsfield units from X-ray Computed Tomography (CT) into relative proton stopping powers (RPSP), the uncertainties of RPSP are 3-5% or higher, which need to be minimized down to 1% to make the proton treatment plans more accurate. In this work, we simulated a proton radiography system, with position-sensitive detectors (PSDs) and a residual energy detector (RED). The simulations were built using Geant4, a Monte Carlo simulation toolkit. A phantom, consisting of several materials was placed between the PSDs of various Water Equivalent Thicknesses (WET), corresponding to an ideal detector, a gaseous detector, silicon and plastic scintillator detectors. The energy loss radiograph and the scattering angle distributions of the protons were studied for proton beam energies of 150MeV, 190MeV and 230MeV. To improve the image quality deteriorated by the multiple Coulomb scattering (MCS), protons with small angles were selected. Two ways of calculating a scattering angle were considered using the proton's direction and position. A scattering angle cut of 8.7mrad was applied giving an optimal balance between quality and efficiency of the radiographic image. For the three proton beam energies, the number of protons used in image reconstruction with the direction method was half the number of protons kept using the position method.