Commit b14f5681 authored by Médéric Boquien's avatar Médéric Boquien

Now that upper limits are stored as positive values, we do not need to take...

Now that upper limits are stored as positive values, we do not need to take the opposite when correcting the χ².
parent 79188fa0
......@@ -273,14 +273,12 @@ def compute_chi2(models, obs, corr_dz, wz, lim_flag):
model = models.flux[band][wz]
chi2 -= 2. * np.log(.5 *
(1. + erf(((obs.flux_ul[band] -
model * scaling) /
(-np.sqrt(2.)*obs_error)))))
model * scaling) / (np.sqrt(2.)*obs_error)))))
for band, obs_error in obs.extprop_ul_err.items():
model = models.extprop[band][wz]
chi2 -= 2. * np.log(.5 *
(1. + erf(((obs.extprop_ul[band] -
model * scaling) /
(-np.sqrt(2.)*obs_error)))))
model * scaling) / (np.sqrt(2.)*obs_error)))))
return chi2, scaling
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment