fix the mask in the CHISQ_TOT calculation

segasai · Nov 4, 2024 · 9fc9d96 · 9fc9d96
1 parent 03660f3
commit 9fc9d96
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/py/rvspecfit/spec_fit.py b/py/rvspecfit/spec_fit.py
@@ -701,9 +701,9 @@ def get_chisq(specdata,
 
             cur_deviation = ((curmodel - curdata.spec) / curdata.espec)
             if curdata.badmask is not None:
-                cur_mask = curdata.badmask
+                cur_mask = ~curdata.badmask
             else:
-                cur_mask = np.zeros(len(cur_deviation), dtype=bool)
+                cur_mask = np.ones(len(cur_deviation), dtype=bool)
             cur_true_chisq = np.sum(cur_deviation[cur_mask]**2)
             chisq_array.append(cur_true_chisq)
             cur_npix = cur_mask.sum()