It turns out I did this in the previous post when I explored using a two dimensional gaussian.
Check here.
I changed the calculation in the signal to noise ratio to:
Signal = The sum of the pixel values within the the effective radius Noise = The RMS of the noise pixels outside 3 * the effective radius