Changeset 20581

Timestamp:

Nov 7, 2008, 3:51:01 PM (18 years ago)

Author:

eugene

Message:

require 10 stars per spatial bin, use a soften mag for residual fit

File:

• trunk/psphot/src/psphotApResid.c (modified) (5 diffs)

trunk/psphot/src/psphotApResid.c

@@ -190 +190 @@
     psVector *dMagSys = (psVector *) psBinaryOp (NULL, (void *) dMag, "*", (void *) psScalarAlloc(errorScale, PS_TYPE_F32));
 
-    if (psTraceGetLevel("psphot") >= 5) {
+    if (psTraceGetLevel("psphot") >= 2) {
         FILE *dumpFile = fopen ("apresid.dat", "w");
         for (int i = 0; i < xPos->n; i++) {
@@ -289 +289 @@
             dMo->data.F32[i] = statsM->sampleMean;
             dRo->data.F32[i] = statsS->robustStdev / statsM->sampleMean;
+            fprintf (stderr, "%d (%d) : sys: %f, phot: %f, rat: %f\n", i, j, dSo->data.F32[i], dMo->data.F32[i], dRo->data.F32[i]);
         } else {
             dSo->data.F32[i] = NAN;
@@ -339 +340 @@
         Nx = PS_MAX (1, Nx);
     }
-
-    if (Npsf < 3*Nx*Ny) {
+    
+    // require at least 10 stars per spatial bin
+    if (Npsf < 10*Nx*Ny) {
         return false;
     }
@@ -356 +358 @@
     psf->ApTrend = pmTrend2DAlloc (PM_TREND_MAP, readout->image, Nx, Ny, stats);
 
+    // XXX somewhat arbitrary: soften the errors so the few bright stars do not totally dominate:
+    psVector *dMagSoft = psVectorAlloc (dMag->n, PS_TYPE_F32);
+    for (int i = 0; i < dMag->n; i++) {
+        dMagSoft->data.F32[i] = hypot(dMag->data.F32[i], 0.01);
+    }
+
     // XXX test for errors here
-    pmTrend2DFit (psf->ApTrend, mask, 0xff, xPos, yPos, apResid, dMag);
+    pmTrend2DFit (psf->ApTrend, mask, 0xff, xPos, yPos, apResid, dMagSoft);
     
     // construct the fitted values and the residuals
@@ -373 +381 @@
 
     psFree (stats);
+    psFree (dMagSoft);
     psFree (apResidFit);
     psFree (apResidRes);