Changeset 15040 for trunk/psphot/src/psphotSourceSize.c

Timestamp:

Sep 26, 2007, 5:35:47 PM (19 years ago)

Author:

eugene

Message:

various cleanups of I/O operations

File:

• trunk/psphot/src/psphotSourceSize.c (modified) (3 diffs)

trunk/psphot/src/psphotSourceSize.c

@@ -7 +7 @@
 
 bool psphotSourceSize (pmReadout *readout, psArray *sources, psMetadata *recipe) {
-
-    FILE *f = fopen ("srcsize.dat", "w");
 
     // loop over all source
@@ -53 +51 @@
         // P(chisq > chisq_obs; Ndof) = gamma_Q (Ndof/2, chisq/2)
         // Ndof = 4 ? (four measurements, no free parameters)
-        float Ppsf = gsl_sf_gamma_inc_Q (2, 0.5*chisq);
+        // XXX this value is going to be biased low because of systematic errors.
+        // we need to calibrate it somehow
+        source->psfProb = gsl_sf_gamma_inc_Q (2, 0.5*chisq);
 
         float fCR = 0.0;
@@ -87 +87 @@
             nEXT ++;
         }
-        fCR = (nCR > 0) ? fCR / nCR : 0.0;
-        fEXT = (nEXT > 0) ? fEXT / nEXT : 0.0;
-
-        // XXX plot this as a function of flux / magnitude?
-        // XXX overlay outliers on an image?
-        
-        fprintf (f, "%f %f %f   %f %f   %f %f\n", source->peak->xf, source->peak->yf, source->peak->flux,
-                 chisq, Ppsf, fCR, fEXT);
+        source->crNsigma  = (nCR > 0)  ? fCR / nCR : 0.0;
+        source->extNsigma = (nEXT > 0) ? fEXT / nEXT : 0.0;
     }
-    fclose (f);
     return true;
 }