Changeset 33435 for trunk/psModules

Timestamp:

Mar 7, 2012, 11:28:04 AM (14 years ago)

Author:

bills

Message:

When choosing sources for stamps, check that the apMag is finite before comparing
it to the psfMag. If it is not finite check whether apMagRaw is finite.
If it is compare it to apMag instead. If neither apMag or apMagRaw is finite
reject the source.
This was the cause of the psphotStack failure to find stamps problem

File:

• trunk/psModules/src/imcombine/pmSubtractionStamps.c (modified) (1 diff)

trunk/psModules/src/imcombine/pmSubtractionStamps.c

@@ -1208 +1208 @@
         // XXX this is somewhat arbitrary...
         if (source->psfMagErr > 0.05) continue;
-        if (fabs(source->psfMag - source->apMag) > 0.5) continue;
+        if (isfinite(source->apMag)) {
+            if (fabs(source->psfMag - source->apMag) > 0.5) continue;
+        } else if (isfinite(source->apMagRaw)) {
+            if (fabs(source->psfMag - source->apMagRaw) > 0.5) continue;
+        } else {
+            // XXX: Should we carry on or drop this source?
+            // drop it for now
+            continue;
+        }
 
         if (source->modelPSF) {