Changeset 30045

Timestamp:

Dec 15, 2010, 11:34:46 AM (15 years ago)

Author:

bills

Message:

use higher sigma limit when finding detections for PPSUB.POS2.SOURCES

File:

• tags/ipp-20101206/ppSub/src/ppSubLoop.c (modified) (1 diff)

tags/ipp-20101206/ppSub/src/ppSubLoop.c

@@ -70 +70 @@
     }
     if (data->forcedPhot2) {
+        // XXX: Change the recipe to use a higher nsigma limit and quit after pass1
+        psMetadata *recipe = psMetadataLookupPtr (NULL, config->recipes, PSPHOT_RECIPE);
+
+        psF32 nsigma_peak_save = psMetadataLookupF32 (NULL, recipe, "PEAKS_NSIGMA_LIMIT");
+        char *breakPt_save =  psMetadataLookupStr (NULL, recipe, "BREAK_POINT");
+
+        // XXX: add recipe entries for this instead of hard coding it
+        psMetadataAddF32(recipe, PS_LIST_TAIL, "PEAKS_NSIGMA_LIMIT", PS_META_REPLACE, "", 25.0);
+        psMetadataAddStr(recipe, PS_LIST_TAIL, "BREAK_POINT", PS_META_REPLACE, "", "PASS1");
+
         bool foundDetections = false;
         if (!ppSubInputDetections(&foundDetections, "PPSUB.POS2.SOURCES", "PPSUB.REF", data)) {
             psError(psErrorCodeLast(), false, "Unable to measure positive detections (2)");
-            return false;
-        }
+            psMetadataAddF32(recipe, PS_LIST_TAIL, "PEAKS_NSIGMA_LIMIT", PS_META_REPLACE, "", nsigma_peak_save);
+            psMetadataAddStr(recipe, PS_LIST_TAIL, "BREAK_POINT", PS_META_REPLACE, "", breakPt_save);
+            return false;
+        }
+        psMetadataAddF32(recipe, PS_LIST_TAIL, "PEAKS_NSIGMA_LIMIT", PS_META_REPLACE, "", nsigma_peak_save);
+        psMetadataAddStr(recipe, PS_LIST_TAIL, "BREAK_POINT", PS_META_REPLACE, "", breakPt_save);
         // if nothing was found, don't bother doing the forced photometry below
         if (!foundDetections) {
-            psWarning ("no sources found in positive image 1, skipping forced photometry");
-            data->forcedPhot1 = false;
+            psWarning ("no sources found in positive image 2, skipping forced photometry");
+            data->forcedPhot2 = false;
         }
     }