From 237907149fa44e675744c721bd2510763966d24a Mon Sep 17 00:00:00 2001 From: Claire Hebert Date: Wed, 13 Nov 2024 21:13:32 +0000 Subject: [PATCH] change calcSizeResidual from (T_obs-T_psf)/T_psf to (T_obs-T_psf)/T_obs --- .../tools/actions/vector/calcRhoStatistics.py | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/python/lsst/analysis/tools/actions/vector/calcRhoStatistics.py b/python/lsst/analysis/tools/actions/vector/calcRhoStatistics.py index febed1191..22c43d407 100644 --- a/python/lsst/analysis/tools/actions/vector/calcRhoStatistics.py +++ b/python/lsst/analysis/tools/actions/vector/calcRhoStatistics.py @@ -215,17 +215,17 @@ def __call__(self, data: KeyedData, **kwargs) -> KeyedData: calcSizeResidual = FractionalDifference( actionA=CalcMomentSize( - colXx=self.colXx, - colYy=self.colYy, - colXy=self.colXy, - sizeType=self.sizeType, - ), - actionB=CalcMomentSize( colXx=self.colPsfXx, colYy=self.colPsfYy, colXy=self.colPsfXy, sizeType=self.sizeType, ), + actionB=CalcMomentSize( + colXx=self.colXx, + colYy=self.colYy, + colXy=self.colXy, + sizeType=self.sizeType, + ), ) # distortion-type ellipticity has a shear response of 2, so we need to @@ -248,7 +248,7 @@ def __call__(self, data: KeyedData, **kwargs) -> KeyedData: eRes = calcEDiff(data, **kwargs) eRes /= responsitivity # type: ignore e1Res, e2Res = np.real(eRes), np.imag(eRes) - sizeRes = calcSizeResidual(data, **kwargs) + sizeRes = -1 * calcSizeResidual(data, **kwargs) # sign flip residual to T_psf - T_model # Scale the sizeRes by ellipticities e1SizeRes = e1 * sizeRes