def model(self, height, s, m):
    dx = (self.x - (self.cutout_size / 2))
    dy = (self.y - (self.cutout_size / 2))
    psf = (height * np.exp((- (((dx / (2 * s)) ** 2) + ((dy / (2 * s)) ** 2)))))
    return (psf + m)
def model(self, height, s, m):
    dx = (self.x - (self.cutout_size / 2))
    dy = (self.y - (self.cutout_size / 2))
    psf = (height * np.exp((- (((dx / (2 * s)) ** 2) + ((dy / (2 * s)) ** 2)))))
    return (psf + m)