Changeset 42267

Timestamp:

Aug 23, 2022, 11:35:52 AM (4 years ago)

Author:

eugene

Message:

add option to use Map instead of Mkron to test for star/galaxy

Location:

branches/eam_branches/ipp-20220316/Ohana/src/relphot

Files:

• include/relphot.h (modified) (1 diff)
• src/args.c (modified) (1 diff)
• src/bcatalog.c (modified) (3 diffs)

branches/eam_branches/ipp-20220316/Ohana/src/relphot/include/relphot.h

@@ -470 +470 @@
 int    PRESERVE_PS1;
 int    REQUIRE_PSFFIT;
+int    USE_APER_FOR_STARGAL;
 int    UPDATE;
 int    SAVE_IMAGE_UPDATES;

branches/eam_branches/ipp-20220316/Ohana/src/relphot/src/args.c

@@ -247 +247 @@
     REQUIRE_PSFFIT = TRUE;
   }
+  USE_APER_FOR_STARGAL = FALSE;
+  if ((N = get_argument (argc, argv, "-use-aper-for-stargal"))) {
+    remove_argument (N, &argc, argv);
+    USE_APER_FOR_STARGAL = TRUE;
+  }
 
   UPDATE = FALSE;

branches/eam_branches/ipp-20220316/Ohana/src/relphot/src/bcatalog.c

@@ -82 +82 @@
       if (DophotSelect && ((catalog[0].measure[offset].photFlags >> 16) != DophotValue)) { Ndophot ++; continue; }
 
+      float Maper = USE_APER_FOR_STARGAL ? catalog[0].measure[offset].Map : catalog[0].measure[offset].Mkron;
+
       // skip garbage measurements
       if (isnan(catalog[0].measure[offset].psfQF)     || (catalog[0].measure[offset].psfQF < 0.95))     { Npsfqf ++; continue; }
       if (isnan(catalog[0].measure[offset].psfQFperf) || (catalog[0].measure[offset].psfQFperf < 0.95)) { Npsfqf ++; continue; }
       if (isnan(catalog[0].measure[offset].M)) { Nnan ++; continue; }
-      if (isnan(catalog[0].measure[offset].Mkron)) { Nnan ++; continue; }
+      if (isnan(Maper)) { Nnan ++; continue; }
 
       // require 0x01 in photFlags (fitted with a PSF) -- add to the test data
       if (REQUIRE_PSFFIT && (catalog[0].measure[offset].photFlags & 0x01) == 0) { Nbad ++; continue; }
 
-      // very loose cut on PSF - Kron 
-      float Mkp = catalog[0].measure[offset].M - catalog[0].measure[offset].Mkron;
+      // very loose cut on PSF - APER mag (Map or Mkron) -- this lightly rejects galaxies & other oddities
+      float Mkp = catalog[0].measure[offset].M - Maper;
       if (fabs(Mkp) > 1.0) { Nbad ++; continue; }
 
@@ -101 +103 @@
       if (catalog[0].measure[offset].photFlags & code->photomPoorMask) { Npoor++; continue;}
 
-      // check for galaxies (XXX skip for now)
-      if (FALSE && !isnan(catalog[0].measure[offset].Map)) {
-        if (catalog[0].measure[offset].M - catalog[0].measure[offset].Map > 0.15) {
-          nEXT ++;
-        } else {
-          nPSF ++;
-        }
+      // weak test for galaxies
+      if (Mkp > 0.5) {
+        nEXT ++;
+      } else {
+        nPSF ++;
       }
 
@@ -145 +145 @@
 
     // skip object if it is likely to be a galaxy
-    if (FALSE && (nEXT >= nPSF)) {
+    if (nEXT >= nPSF) {
       Nmeasure -= Nm;
       Ngalaxy ++;