Closed itaboada closed 3 years ago
Looks like we could pretty easily interface with gammapy for the attenuations:
I don't know of any programs that convert numu to gamma-ray fluxes, though. If anyone has any suggestions, I'm sure we could figure out a way to make it all work.
It turns out that Legend.py simulates gamma rays, not neutrinos. To use Legend.py to simulate neutrinos, we will need the proposed utility to calculate fluxes based on p-p and p-gamma processes.
This should be done for pp and p-gamma interactions. It could be a module that is called as an option itself in FIRESONG, which would output gamma-ray flux or as a standalone program that converts numu fluxes to gamma ray fluxes. Should be easy to code.