Neutral Pion production for use with nucdens python package https://pypi.org/project/nucdens/
Currently the code only works for the one body case. The two body case is running but does not yet give physical results.
Initial implimentation based off of this paper Lenkewitz 2011: Neutral pion photoproduction off 3H and 3He in chiral perturbation theory https://arxiv.org/abs/1103.3400
With use of the density formalism as described in this paper. Griesshammer et. al. Scattering Observables from One- and Two-Body Densities: Formalism and Application to γ3He Scattering https://arxiv.org/abs/2005.12207