Double differential distributions of neutrons from 100 - 250 MeV protons

Your figures above patently contradict this statement.

This is because, contrary to the paper authors, you did not set to blackhole the concrete shell outside the scoring regions, and thereby you got the contribution of the neutrons reflected by it.

Also, note that the unit you indicated on the vertical axis is wrong, since the plotted quantity is not expressed in MeV^{-1}, because of the multiplication by the energy implied by the lethargy representation.