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