Changeset 20054 for trunk/ppSub/src/ppSubReadout.c

Timestamp:

Oct 10, 2008, 3:10:59 PM (18 years ago)

Author:

Paul Price

Message:

Adding additional QA values.

File:

• trunk/ppSub/src/ppSubReadout.c (modified) (3 diffs)

trunk/ppSub/src/ppSubReadout.c

@@ -285 +285 @@
         // Statistics on the matching
         if (stats) {
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MODE);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_STAMPS);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_DEV_MEAN);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_DEV_RMS);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_NORM);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_BGDIFF);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MX);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MY);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MXX);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MXY);
+            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MYY);
+
             psMetadataAddF32(stats, PS_LIST_TAIL, "MATCH_TIME", 0, "Time to match PSFs",
                              psTimerClear("PPSUB_MATCH"));
-            psMetadataAddS32(stats, PS_LIST_TAIL, "SUBTRACTION.NUM", 0, "Number of stamps",
-                             kernels->numStamps);
-            psMetadataAddF32(stats, PS_LIST_TAIL, "SUBTRACTION.MEAN", 0, "Mean stamp deviation",
-                             kernels->mean);
-            psMetadataAddF32(stats, PS_LIST_TAIL, "SUBTRACTION.RMS", 0, "RMS stamp deviation",
-                             kernels->rms);
-            metadataCopySingle(stats, outRO->analysis, PM_SUBTRACTION_ANALYSIS_MODE);
         }
     }
@@ -539 +544 @@
             psArray *sources = psMetadataLookupPtr(NULL, photRO->analysis, "PSPHOT.SOURCES"); // Sources
             psMetadataAddS32(stats, PS_LIST_TAIL, "NUM_SOURCES", 0, "Number of sources detected", sources->n);
-            psMetadataAddF32(stats, PS_LIST_TAIL, "MATCH_TIME", 0, "Time to do photometry",
+            psMetadataAddF32(stats, PS_LIST_TAIL, "PHOT_TIME", 0, "Time to do photometry",
                              psTimerClear("PPSUB_PHOT"));
         }
@@ -589 +594 @@
     }
 
+    // Significance image
+    {
+        psImage *sig = (psImage*)psBinaryOp(NULL, outRO->image, "*", outRO->image);
+        psBinaryOp(sig, sig, "/", outRO->weight);
+        psFits *fits = psFitsOpen("significance.fits", "w");
+        psFitsWriteImage(fits, NULL, sig, 0, NULL);
+        psFitsClose(fits);
+        psFree(sig);
+    }
+
     psFree(outRO);