Adding gamma flux
The gamma fluxes described in https://arxiv.org/abs/0807.4151v1 can be used to calculate the neutrino emission spectrum.
This is also implemented in https://git.km3net.de/common/flux (see https://git.km3net.de/common/flux/-/blob/master/nufromgamma.hh) but the ROOT dependency is still there, so we cannot easily use that. For now, I propose to transcribe the formulas in the paper to Python (preferably using Numba's JIT power) so that it's maximum performance. Later if flux.so
is a thing, we can exchange it easily.
I had a quick look and implemented the formulas 30-33 in Python. Here is an example of some of the mentioned gamma fluxes:
Edited by Tamas Gal