Changeset 15636

Timestamp:

Nov 16, 2007, 2:37:55 PM (18 years ago)

Author:

eugene

Message:

allow limits to inst mag range

File:

• trunk/psastro/src/psastroConvert.c (modified) (2 diffs)

trunk/psastro/src/psastroConvert.c

@@ -49 +49 @@
     // we are going to select the brighter Nmax subset for astrometry 
     int nMax = psMetadataLookupS32 (&status, recipe, "PSASTRO.MAX.NRAW");
-    if (!status || !nMax) nMax = inStars->n; // 10 is really somewhat absurd as a lower limit
+    if (!status || !nMax) nMax = inStars->n;
+
+    // we are going to select the brighter Nmax subset for astrometry 
+    float iMagMax = psMetadataLookupF32 (&status, recipe, "PSASTRO.MAX.INST.MAG.RAW");
+    float iMagMin = psMetadataLookupF32 (&status, recipe, "PSASTRO.MIN.INST.MAG.RAW");
 
     // choose the first nMax sources
@@ -58 +62 @@
         pmSource *source = sources->data[n];
         if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
+        if ((iMagMax != 0.0) && (source->psfMag > iMagMax)) continue;
+        if ((iMagMin != 0.0) && (source->psfMag < iMagMin)) continue;
         rawStars->data[j] = psMemIncrRefCounter (inStars->data[n]);
         j++;