psf = PSF.get_psf_object( 'OversampledImagePSF',
x=diff_img.data.shape[1]/2., y=diff_img.data.shape[1]/2.,
oversample_factor=1.,
data=psf_img.data )
I think that x and y should be at the position of the supernova on diff_img, not at the center of the image.
I also think (THINK ABOUT THIS) that we want to round to the nearest integer.