Distribution fitting optimisation for SPE proton fluence as a function of energy for exponential cut-off power law (green), truncated power law (red) and lognormal distribution (blue). Top left panel: degrees of freedom for optimised solutions. Bottom left panel: (reduced) goodness-of-fit parameters (see text for details). Bottom right panel: resulting thresholds to determine inclusion of SPEs from the REL.