Changeset 38774

Timestamp:

Sep 25, 2015, 12:26:34 PM (11 years ago)

Author:

eugene

Message:

use psfFluxErr instead of apFluxErr as apFluxErr is broken

File:

• trunk/Ohana/src/relphot/src/setMrelCatalog.c (modified) (3 diffs)

trunk/Ohana/src/relphot/src/setMrelCatalog.c

@@ -217 +217 @@
     if (isSetMrelFinal) {
       Map = PhotCat (&measure[k], MAG_CLASS_APER);
-      dMap = MAX (hypot(measure[k].dMap,   code->photomErrSys), MIN_ERROR);
+
+      // NOTE: for PV3, apFluxErr is broken (see pmSourcePhotometry.c:245).  we are going to
+      // use PSF flux error instead here:
+      // dMap = MAX (hypot(measure[k].dMap,   code->photomErrSys), MIN_ERROR);
+      dMap = MAX (hypot(measure[k].dM,   code->photomErrSys), MIN_ERROR);
+
       if (CHECK_VALID_MAG(Map, dMap)) {
         int Nap = results->aperData[Nsec].Nlist;
@@ -638 +643 @@
     secfilt[Nsec].dFkronStk = zpFactor * measure[k].dFluxKron;
     secfilt[Nsec].FapStk    = zpFactor * measure[k].FluxAp;
-    secfilt[Nsec].dFapStk   = zpFactor * measure[k].dFluxAp;
+
+    // NOTE: for PV3, apFluxErr is broken (see pmSourcePhotometry.c:245).  we are going to
+    // use PSF flux error instead here:
+    // secfilt[Nsec].dFapStk   = zpFactor * measure[k].dFluxAp;
+    secfilt[Nsec].dFapStk   = zpFactor * measure[k].dFluxPSF;
 
     // Jy to AB mags
@@ -773 +782 @@
 
     float Fap = PhotFluxCat (&measure[k], MAG_CLASS_APER);
-    float dFap = PhotFluxCatErr (&measure[k], MAG_CLASS_APER);
+
+    // NOTE: for PV3, apFluxErr is broken (see pmSourcePhotometry.c:245).  we are going to
+    // use PSF flux error instead here:
+    // float dFap = PhotFluxCatErr (&measure[k], MAG_CLASS_APER);
+    float dFap = PhotFluxCatErr (&measure[k], MAG_CLASS_PSF);
+
     if (CHECK_VALID_FLUX(Fap, dFap)) {
       dFap = MAX (hypot(dFap, code->photomErrSys*Fap), MIN_ERROR*Fap); // bump up the error by a systematic floor