corrected definition of smoothcorr statistics

desihub · Nov 17, 2024 · e36a655 · e36a655
1 parent 61b5841
commit e36a655
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/py/fastspecfit/continuum.py b/py/fastspecfit/continuum.py
@@ -1526,7 +1526,7 @@ def do_fit_full(objflam, specflux):
         for icam, (ss, ee) in enumerate(data['camerapix']):
             I = ((specflux[ss:ee] != 0.) & (specivar[ss:ee] != 0.) & (smoothcontinuum[ss:ee] != 0.))
             if np.count_nonzero(I) > 3: # require three good pixels to compute the mean
-                smoothstats[icam] = np.mean(1. - smoothcontinuum[ss:ee][I] / specflux[ss:ee][I])
+                smoothstats[icam] = median(smoothcontinuum[ss:ee][I] / specflux[ss:ee][I])
 
     log.debug(f'Deriving the smooth continuum took {time.time()-t0:.2f} seconds.')