Probability model for peak fluxes of solar proton events