My group is currently attempting to simulate 3-photon annihilations to mimic that of ortho-Positronium self annihilation decay.
We have concluded that fully simulating o-Ps formation and then its annihilation is perhaps too complicated. Instead we want to simply simulate the 3-photon kinematics after annihilation. However, we are having trouble figuring out how to properly and randomly select the momentum vectors and energy of each photon such that the net momentum is 0 and the net energy sums to 1022 keV.
Could someone point us to a physics model that does this well or perhaps a clever method that we could implement ourselves with particle generator such that the 3 gammas meet the specifications I mentioned above?
We have tried to use the G4eplusTo3GammaOKCModel() and while this does allow for 3-photon annihilations, the result is not what we want, neither momentum or energy seems to be conserved. We also ensured to set annihilations to occur at rest, yet this model still allowed for 3-photon decays that had a net momentum. This may be because this model was intended for a specific purpose, not o-Ps self annihilations.
Thank you, I am happy to provide any additional information!