Some additional comments

hpc4cmb · Sep 12, 2023 · d34fdb8 · d34fdb8
1 parent 4de8372
commit d34fdb8
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/src/toast/ops/noise_model.py b/src/toast/ops/noise_model.py
@@ -429,8 +429,6 @@ def _get_err_ret(self, psd_unit):
         eret = dict()
         eret["fit_result"] = types.SimpleNamespace()
         eret["fit_result"].success = False
-        # Use a large value, to avoid divide by zero elsewhere
-        # in the code if the detector flags are not examined.
         eret["NET"] = 0.0 * np.sqrt(1.0 * psd_unit)
         eret["fmin"] = 0.0 * u.Hz
         eret["fknee"] = 0.0 * u.Hz
@@ -665,6 +663,9 @@ def _exec(self, data, detectors=None, **kwargs):
             net_std = None
             fknee_mean = None
             fknee_std = None
+            # If the noise model came from fitting, then detectors with a bad
+            # fit are already flagged in addition to NET being set to zero.
+            # This check is just an additional safeguard.
             good_fit = local_net > 0.0
             if obs.comm_col is None:
                 all_net = np.array(local_net[good_fit])
@@ -718,7 +719,7 @@ def _exec(self, data, detectors=None, **kwargs):
                         msg = f"obs {obs.name}, det {det} has fknee "
                         msg += f"{local_fknee[idet]} that is > {self.sigma_fknee} "
                         msg += f"x {fknee_std} from {fknee_mean}"
-                        log.info(msg)
+                        log.debug(msg)
                         obs.detdata[self.det_flags][det, :] |= self.det_flag_mask
                         new_flags[det] = cur_flag | self.det_flag_mask
             obs.update_local_detector_flags(new_flags)