This implementation works, but is not very elegant. It could probably be made easier using the astropy.coordinates framework at least for some of the work.
We might not want to create a new coordinate frame for every photon to reduce the overhead, but at least for the nominal direction we could do that.
This implementation works, but is not very elegant. It could probably be made easier using the astropy.coordinates framework at least for some of the work. We might not want to create a new coordinate frame for every photon to reduce the overhead, but at least for the nominal direction we could do that.