DM-40556: Add LATISS color terms and turn on by default in initial calibration.

config/latiss/calibrate.py

@@ -1,10 +1,23 @@
+import os.path
+
 from lsst.meas.astrom import FitAffineWcsTask
+from lsst.meas.algorithms import ColorLimit
+
+ObsConfigDir = os.path.dirname(__file__)
 
 # Configure the photometry to use atlas_refcat2.
-config.connections.photoRefCat = 'atlas_refcat2_20220201'
+config.connections.photoRefCat = "atlas_refcat2_20220201"
+config.photoRefObjLoader.load(os.path.join(ObsConfigDir, "filterMap.py"))
+
 config.photoCal.match.referenceSelection.magLimit.fluxField = "r_flux"
+colors = config.photoCal.match.referenceSelection.colorLimits
+colors["g-r"] = ColorLimit(primary="g_flux", secondary="r_flux", minimum=0.4, maximum=2.0)
+
+config.photoCal.applyColorTerms = True
+config.photoCal.photoCatName="atlas_refcat2_20220201"
+config.photoCal.colorterms.load(os.path.join(ObsConfigDir, "colorterms.py"))
 
-# We often have very few sources due to smaller aperture so use affine task. 
+# We often have very few sources due to smaller aperture so use affine task.
 config.astrometry.wcsFitter.retarget(FitAffineWcsTask)
 
 # Note that the following two config values were validated on data taken in

config/latiss/colorterms.py

@@ -0,0 +1,37 @@
+from lsst.pipe.tasks.colorterms import Colorterm, ColortermDict
+
+
+config.data = {
+    "atlas_refcat2*": ColortermDict(data={
+        "SDSSg_65mm~empty": Colorterm(
+            primary="g",
+            secondary="r",
+            c0=-0.09034144345111599,
+            c1=0.1710923238086337,
+            c2=-0.038260355621929296,
+        ),
+        "SDSSr_65mm~empty": Colorterm(
+            primary="r",
+            secondary="i",
+            c0=0.0073632488906825045,
+            c1=-0.026620900037027242,
+            c2=-0.03203533692013322,
+        ),
+        "SDSSi_65mm~empty": Colorterm(
+            primary="i",
+            secondary="r",
+            c0=0.016940180565664747,
+            c1=0.0610018330811135,
+            c2=-0.0722575356707918,
+        ),
+        # The following two are blank until we have data to measure them.
+        "SDSSz_65mm~empty": Colorterm(
+            primary="z",
+            secondary="z",
+        ),
+        "SDSSy_65mm~empty": Colorterm(
+            primary="y",
+            secondary="y",
+        ),
+    }),
+}

config/latiss/filterMap.py

@@ -0,0 +1,12 @@
+# Mapping of camera filter name: reference catalog filter name
+# This file is appropriate for ATLAS Refcat2.
+
+for source, target in [
+        ("SDSSg_65mm~empty", "g"),
+        ("SDSSr_65mm~empty", "r"),
+        ("SDSSi_65mm~empty", "i"),
+        ("SDSSz_65mm~empty", "z"),
+        # ATLAS Refcat2 does not have y band.
+        ("SDSSy_65mm~empty", "z"),
+    ]:
+    config.filterMap[source] = target