Proton LET discrepancies in CR39

Hello, I’m hoping someone has a suggestion for how I can investigate the discrepancies in proton LET in CR39 I am recording in my simulations. The plot attached is taken from the following paper Caresana, M., Ferrarini, M., Fuerstner, M., & Mayer, S. (2012). Determination of LET in PADC detectors through track parameter measurement. Nucl. Instrum. Methods Phys. Res. A, 683, 8–15. Redirecting and is what I’m using to benchmark my results. I attached a plot I generated of proton energy vs LET from simulating the irradiation of a CR39 film with the neutron spectrum generated by a 15 MV photon beam at a defined point in the linac bunker. The overall shape seems good to me, but it appears to almost be shifted entirely to the right since key points don’t correspond. I’m using QGSP_BIC_HP physics list and g4em-standard_opt4. I’d appreciate any input!

Geant4 Version: v11.1.3
Operating System: macOS Sonoma 14.3.1