Implement higher-order electroweak corrections in νpropa
Implement next-to-leading-order electroweak corrections, or suitable analytic approximations to those corrections, for the inclusive and differential neutrino-interaction cross sections in the νpropa simulation framework.
References
These results indicate that the NLO EW corrections become significant at extremely large centre of mass energies (reaching $-50\%$). We note that the numerical evaluation of the NLO virtual corrections is rather time consuming, and therefore directly evaluating the squared amplitude when sampling the differential cross-section for the various interaction channels does not seem tractable for the full propagation framework. Possible future strategies could be to pre-tabulate the NLO predictions for the inclusive and/or differential cross section. Alternatively, as the dominant numerical contribution of these correction arises from large Sudakov logarithms, it may also be feasible to include analytic approximations for the various interaction channels (at both the differential and inclusive levels)---see, for example, Ref.. The inclusion of these effects as part of the plug-in is left for future work.