The neutrino (and \(\gamma \) ) fluxes from the decays of mesons produced in \(p\gamma \)  interactions can be computed as [1]}
\(\frac{{\rm d}N_{i}}{{\rm d}E^{\prime }_i}(E^{\prime }_i)=\int \frac{{\rm d}E^{\prime }_p}{E^{\prime }_p}\,{\rm d}\varepsilon ^{\prime } \, f_p(E^{\prime }_p)f_{\rm ph}(\varepsilon ^{\prime })\Phi _{i}(\eta ,x_i),\) 
