Multiple scattering calculations for proton beams: Comparison of results from the general-purpose Monte Carlo codes penh, fluka and topas

Abstract

Purpose: To test the multiple Coulomb scattering theories implemented in the Monte Carlo simulation codes PENH, FLUKA and TOPAS. Methods: Simulations with the three codes of proton beams with initial energies between 100 and 220 MeV impinging normally on slabs of 14 different compositions and various thicknesses were performed. The simulated angular distributions of transmitted protons are very approximately Gaussian with a characteristic angle θ0, which measures the spread of the distribution. The characteristic angles resulting from the simulations are compared to experimental data available from two experiments. The degree of agreement is quantified by the χ2 statistic. The 151 cases considered are analyzed by grouping the data in various ways (all together, by experiment, according to the energy of the proton beam or the slab material). Results: In general, PENH produced the better description of the experimental data. If all data are included, the χ2 values were 3.7±0.6 for PENH, 18.7±0.4 for FLUKA and 7.4±1.5 for TOPAS. The ranges in the values of χ2 obtained for the various data groups are [0.1±0.2;10.7±4.1] for PENH, [0.10±0.03;221.1±7.8] for FLUKA, and [0.2±0.3;46.2±26.4] for TOPAS. The minimum and maximum values in these ranges occur, in the three codes, for the Zn and the brass slabs, respectively, with 158.6MeV protons. Conclusions: The three codes provide a reasonably accurate description of multiple scattering distributions of protons transmitted through material slabs. Brass is the material whose experimental data are worst described by the three codes analyzed.