Nuclear interactions of 160 MeV protons stopping in copper: a test of Monte Carlo nuclear models

To estimate the influence of nuclear interactions on dose or biological effect one uses Monte Carlo programs which include nuclear models. We introduce an experimental method to check these models at proton therapy energies. We have measured the distribution of charge deposited by 160 MeV protons stopping in a stack of insulated copper plates. A buildup region ahead of the main peak contains approximately 20% of the total charge and is entirely due to charged secondaries from inelastic nuclear interactions. The acceptance for charged secondaries is 100%. Therefore the data are a good benchmark for nuclear models. We have simulated the stack using GEANT with two nuclear models. FLUKA agrees fairly well with the measurement but GHEISHA, designed for much higher energies, does not. The experimental method will work for many other materials, including insulators. Therefore it can also be used for light nuclei.