Changeset 41608 for branches/eam_branches/relphot.20210521/src/MosaicOps.c

Timestamp:

May 22, 2021, 1:44:24 PM (5 years ago)

Author:

eugene

Message:

fix inconsistency between Mgrid and Mflat

File:

• branches/eam_branches/relphot.20210521/src/MosaicOps.c (modified) (3 diffs)

branches/eam_branches/relphot.20210521/src/MosaicOps.c

@@ -1069 +1069 @@
     // apply that offset as well here for this image (in other words, each detection is
     // being compared to the model, excluding the zero point, Mcal.  The model includes
-    // the flat-correction.  NOTE the sign of Mflat (Image.Mcal = Measure.Mcal - Mflat)
+    // the flat-correction.  NOTE the sign of Mflat (Image.Mcal = Measure.Mcal + Mflat).
+    // this was inconsistent pre r41606
 
     float Mflat = getMflat (m, c, catalog);
@@ -1094 +1095 @@
     assert (Nbright >= 0);
 
-    float Moff =  Mcal + Mgrp + Mgrid - Mflat;
+    float Moff =  Mcal + Mgrp + Mgrid + Mflat;
 
     psfStars->alldata-> yVector[N] = MsysPSF - MrelPSF - Moff;
@@ -1205 +1206 @@
     float MsysPSF = PhotSysTiny (&catalog[c].measureT[m], &catalog[c].averageT[n], &catalog[c].secfilt[n*Nsecfilt], MAG_CLASS_PSF);
     
-    float delta = MsysPSF - MrelPSF - Mcal - Mgrid + Mflat;
+    float delta = MsysPSF - MrelPSF - Mcal - Mgrid - Mflat;
 
     int isBad = (catalog[c].measureT[m].dbFlags & MEAS_BAD);