Index: trunk/psphot/src/psphotSourceSize.c
===================================================================
--- trunk/psphot/src/psphotSourceSize.c	(revision 13086)
+++ trunk/psphot/src/psphotSourceSize.c	(revision 15040)
@@ -7,6 +7,4 @@
 
 bool psphotSourceSize (pmReadout *readout, psArray *sources, psMetadata *recipe) {
-
-    FILE *f = fopen ("srcsize.dat", "w");
 
     // loop over all source
@@ -53,5 +51,7 @@
 	// 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,14 +87,7 @@
 	    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;
 }