Changeset 10203

Timestamp:

Nov 26, 2006, 12:38:36 PM (20 years ago)

Author:

eugene

Message:

fixed chiDBL wrt chiEXT (was not dividing by nDOF)

File:

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

trunk/psphot/src/psphotSourceFits.c

@@ -80 +80 @@
         // if this one failed, skip it
         if (!psphotEvalPSF (blend, model)) {
-            psTrace ("psphot", 5, "failed on blend %d of %ld\n", i, modelSet->n);
+            psTrace ("psphot", 4, "failed on blend %d of %ld\n", i, modelSet->n);
             continue;
         }
@@ -92 +92 @@
         blend->mode &= ~PM_SOURCE_MODE_TEMPSUB;
     }
+    // evaluate the primary object
+    if (!psphotEvalPSF (source, PSF)) {
+        psTrace ("psphot", 4, "failed on blend 0 of %ld\n", modelSet->n);
+        psFree (PSF);
+        psFree (modelSet);
+        psFree (sourceSet);
+        return false;
+    }
     psFree (modelSet);
     psFree (sourceSet);
-
-    // evaluate the primary object
-    if (!psphotEvalPSF (source, PSF)) {
-        psTrace ("psphot", 5, "failed on blend 0 of %ld\n", modelSet->n);
-        psFree (PSF);
-        return false;
-    }
 
     psTrace ("psphot", 5, "fitted primary as PSF\n");
@@ -210 +211 @@
 
     // save chisq for double-star/galaxy comparison
-    chiDBL = ONE->chisqNorm;
+    chiDBL = ONE->chisq / ONE->nDOF;
 
     // correct second model chisqs for flux trend