From 8cc82d01ff60d0a949a449b941242353f0610cee Mon Sep 17 00:00:00 2001 From: John Parejko Date: Fri, 22 Mar 2024 10:31:40 -0700 Subject: [PATCH] SDSS centroid error fallback --- python/lsst/ip/diffim/dipoleFitTask.py | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) diff --git a/python/lsst/ip/diffim/dipoleFitTask.py b/python/lsst/ip/diffim/dipoleFitTask.py index 44941c88..a3b9e92a 100644 --- a/python/lsst/ip/diffim/dipoleFitTask.py +++ b/python/lsst/ip/diffim/dipoleFitTask.py @@ -1026,7 +1026,7 @@ def _setupSchema(self, config, name, schema, metadata): self.centroidKey = measBase.CentroidResultKey.addFields(schema, name, "Dipole centroid position.", - measBase.UncertaintyEnum.NO_UNCERTAINTY) + measBase.UncertaintyEnum.SIGMA_ONLY) self.orientationKey = schema.addField( schema.join(name, "orientation"), type=float, units="deg", @@ -1115,6 +1115,10 @@ def measureDipoles(self, measRecord, exposure, posExp=None, negExp=None): if not self.config.fitAllDiaSources: measRecord[self.centroidKey.getX()] = measRecord["base_SdssCentroid_x"] measRecord[self.centroidKey.getY()] = measRecord["base_SdssCentroid_y"] + self.centroidKey.getCentroidErr().setElement(measRecord, 0, 0, + measRecord["base_SdssCentroid_xErr"]) + self.centroidKey.getCentroidErr().setElement(measRecord, 1, 1, + measRecord["base_SdssCentroid_yErr"]) measRecord[self.flagKey] = measRecord["base_SdssCentroid_flag"] return @@ -1235,6 +1239,10 @@ def fail(self, measRecord, error=None): """ measRecord[self.centroidKey.getX()] = measRecord["base_SdssCentroid_x"] measRecord[self.centroidKey.getY()] = measRecord["base_SdssCentroid_y"] + self.centroidKey.getCentroidErr().setElement(measRecord, 0, 0, + measRecord["base_SdssCentroid_xErr"]) + self.centroidKey.getCentroidErr().setElement(measRecord, 1, 1, + measRecord["base_SdssCentroid_yErr"]) measRecord.set(self.flagKey, True) if error is not None: if error.getFlagBit() == self.FAILURE_EDGE: