def PSF(x, y, z, h): num_objects = x.size objects = numpy.zeros((num_objects, 5)) for i in range(num_objects): [sx, sy] = multi_fit_c.calcSxSy(wx_params, wy_params, z[i] * 0.001) objects[i, :] = [x[i], y[i], h[i], sx, sy] return objects
def PSFIntegral(z, h): integral = numpy.zeros(z.size) for i in range(z.size): [sx, sy] = multi_fit_c.calcSxSy(wx_params, wy_params, z[i] * 0.001) integral[i] = 2.0 * numpy.pi * h[i] * sx * sy return integral