CRPropa is a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles. https://crpropa.github.io/CRPropa3/
I tried to write a source code NeutrinoOutput1D.cpp similar as PhotonOutput1D.cpp (also NeutrinoOutput1D.h) to calculate secondary neutrinos. I changed "if ((pid != 22) and (abs(pid) != 11))" to "if ((abs(pid) != 12) and (abs(pid) != 14))", then recompile the code, but it seems doesn't work.
I tried to write a source code NeutrinoOutput1D.cpp similar as PhotonOutput1D.cpp (also NeutrinoOutput1D.h) to calculate secondary neutrinos. I changed "if ((pid != 22) and (abs(pid) != 11))" to "if ((abs(pid) != 12) and (abs(pid) != 14))", then recompile the code, but it seems doesn't work.