def get_ia(z):
mabs = np.random.normal(-19.3, 0.3)
x1 = np.random.normal(0, 1.0)
c = np.random.normal(0, 0.1)
lc = generate_ia_light_curve(z, mabs, x1, c)
try:
res = get_summary_stats(z, lc, method="iminuit")
assert fit_is_good(*res[0])
except Exception:
return None
return res
