Reject diaSource on diffraction spikes of bright stars

Author: isullivan

python/lsst/ip/diffim/detectAndMeasure.py

@@ -301,6 +301,8 @@ class DetectAndMeasureConfig(pipeBase.PipelineTaskConfig,
                  "base_PixelFlags_flag_nodataCenter",
                  "base_PixelFlags_flag_saturatedCenter",
                  "base_PixelFlags_flag_saturated_templateCenter",
+                 "base_PixelFlags_flag_spikeCenter",
+                 "base_PixelFlags_flag_spike_templateCenter",
                  ),
     )
     clearMaskPlanes = lsst.pex.config.ListField(
@@ -424,9 +426,11 @@ def setDefaults(self):
 
         # Keep track of which footprints contain streaks
         self.measurement.plugins["base_PixelFlags"].masksFpAnywhere = [
-            "STREAK", "INJECTED", "INJECTED_TEMPLATE", "HIGH_VARIANCE", "SATURATED_TEMPLATE"]
+            "STREAK", "INJECTED", "INJECTED_TEMPLATE", "HIGH_VARIANCE", "SATURATED_TEMPLATE",
+            "SPIKE", "SPIKE_TEMPLATE"]
         self.measurement.plugins["base_PixelFlags"].masksFpCenter = [
-            "STREAK", "INJECTED", "INJECTED_TEMPLATE", "HIGH_VARIANCE", "SATURATED_TEMPLATE"]
+            "STREAK", "INJECTED", "INJECTED_TEMPLATE", "HIGH_VARIANCE", "SATURATED_TEMPLATE",
+            "SPIKE", "SPIKE_TEMPLATE"]
         self.skySources.avoidMask = ["DETECTED", "DETECTED_NEGATIVE", "BAD", "NO_DATA", "EDGE"]
 
     def validate(self):

python/lsst/ip/diffim/subtractImages.py

@@ -284,7 +284,7 @@ class AlardLuptonSubtractBaseConfig(lsst.pex.config.Config):
     )
     renameTemplateMask = lsst.pex.config.ListField(
         dtype=str,
-        default=("SAT", "INJECTED", "INJECTED_CORE",),
+        default=("SAT", "INJECTED", "INJECTED_CORE", "SPIKE"),
         doc="Mask planes from the template to propagate to the image difference"
         "with '_TEMPLATE' appended to the name."
     )