modified detector_temp check algorithm

aeri_qc/cooler_checks.py

@@ -87,7 +87,7 @@ def flatten_out(original_flags, pds):
         std_data['mean_std'].append(mean_std)
         std_mean_std = np.std(std_data['mean_std'][lower:upper])
         std_data['std_mean_std'].append(std_mean_std)
-        if std_mean_std < 20:
+        if std_mean_std < 20 or np.isnan(std_mean_std):
             flags[lower:upper] = flags[lower:upper].replace(False, True)
         else:
             flags[lower:upper] = flags[lower:upper].replace(True, False)
@@ -98,7 +98,7 @@ def flatten_out(original_flags, pds):
 
     new_flags = flags & (np.asarray(std_data['mean']) > 0.35)
 
-    return new_flags, std_data
+    return new_flags
 
 @checklist.add_check(depends=['detectorTemp', 'datetime'], affects_calibration=True, description='test for strong periodicity in temperature on the order of 2-10 mins')
 def detector_temp_check(frame, parameters):
@@ -110,12 +110,12 @@ def detector_temp_check(frame, parameters):
     flags_total = {}
     for win_type in 'bartlett','parzen','hann','nuttall':
         flags,pds = flag_detector_temp_oscillations(detector_temp,win_type=win_type)
-        flags, std_data = flatten_out(flags.copy(), pds)
+        flags = flatten_out(flags.copy(), pds)
         flags_total[win_type] = flags
 
-    flags_final = flags_total['bartlett'] | flags_total['parzen']
-    flags_final = flags_final | flags_total['hann']
-    flags_final = flags_final | flags_total['nuttall']
+    flags_final = flags_total['bartlett'] & flags_total['parzen']
+    flags_final = flags_final & flags_total['hann']
+    flags_final = flags_final & flags_total['nuttall']
 
     return flags_final