Changeset 20339

Timestamp:

Oct 22, 2008, 5:52:22 PM (18 years ago)

Author:

Paul Price

Message:

Ensure sources are over the S/N limit.

File:

• trunk/psphot/src/psphotMergeSources.c (modified) (2 diffs)

trunk/psphot/src/psphotMergeSources.c

@@ -71 +71 @@
     detections->peaks = psArrayAllocEmpty(100);
 
+    psMetadata *recipe = psMetadataLookupPtr (NULL, config->recipes, PSPHOT_RECIPE);
+    if (!recipe) {
+        psError(PSPHOT_ERR_CONFIG, false, "missing recipe %s", PSPHOT_RECIPE);
+        return false;
+    }
+    float snMin = psMetadataLookupF32(NULL, recipe, "MOMENTS_SN_MIN");
+    if (!isfinite(snMin)) {
+        return NULL;
+    }
+
     for (int i = 0; i < sources->n; i++) {
         pmSource *source = sources->data[i];
         pmModel *model = source->modelPSF;
 
+        if (!isfinite(source->psfMag)) {
+            continue;
+        }
+
+        float flux = powf(10.0, -0.4 * source->psfMag);
         float xpos = model->params->data.F32[PM_PAR_XPOS];
         float ypos = model->params->data.F32[PM_PAR_YPOS];
@@ -81 +96 @@
         peak->xf = xpos;
         peak->yf = ypos;
-        peak->flux = 1.0;
-        peak->SN = INFINITY;
+        peak->flux = flux;
+        peak->SN = snMin + flux;
 
         psArrayAdd (detections->peaks, 100, peak);