The pion production interaction rates are evaluated without interpolation. Instead the next lower tabulated value is used. This can lead to significant errors.
Example: for a proton of 102.3 EeV the rate for 93.2603 EeV is used, resulting in a rate that is ~25% too low.
Energy [EeV] Interaction rate [1/Mpc]
93.2603 0.0272672
102.353 0.0367994
DavidWalz
commented
10 years ago
The pion production interaction rates are evaluated without interpolation. Instead the next lower tabulated value is used. This can lead to significant errors. Example: for a proton of 102.3 EeV the rate for 93.2603 EeV is used, resulting in a rate that is ~25% too low. Energy [EeV] Interaction rate [1/Mpc] 93.2603 0.0272672 102.353 0.0367994
i_in = min( (int) (log10(lRefEnergy / lpInt->E_pionprod(0)) / lpInt->dEtabPion()), 99); rate = (i_in > 0) ? _fChargeNumber * scaleFactor * lpInt->LossRateProton(i_in) : 0;