Changeset 23844 for trunk/psastro/src/psastroMaskUpdates.c

Timestamp:

Apr 13, 2009, 5:18:31 PM (17 years ago)

Author:

eugene

Message:

count ghosts

File:

• trunk/psastro/src/psastroMaskUpdates.c (modified) (4 diffs)

trunk/psastro/src/psastroMaskUpdates.c

@@ -52 +52 @@
         return false;
     }
+    bool COUNT_GHOSTS = psMetadataLookupF32 (&status, recipe, "REFSTAR_COUNT_GHOSTS");
+    double GHOST_MAX_MAG = psMetadataLookupF32 (&status, recipe, "GHOST_MAX_MAG");
+    int nGhosts = 0;
 
     psLogMsg ("psastro", PS_LOG_INFO, "generating a bright-star mask");
@@ -88 +91 @@
     REFSTAR_MASK_SATSPIKE_MAG_MAX += MagOffset;
     REFSTAR_MASK_BLEED_MAG_MAX += MagOffset;
+    GHOST_MAX_MAG += MagOffset;
 
     // select the output mask image :: we mosaic to chip mosaic format
@@ -183 +187 @@
                 for (int i = 0; i < refstars->n; i++) {
                     pmAstromObj *ref = refstars->data[i];
+                    
+                    if (COUNT_GHOSTS) {
+                        if (ref->Mag > GHOST_MAX_MAG) {
+                            nGhosts ++;
+                        }
+                    }
+
                     if (ref->Mag > REFSTAR_MASK_MAX_MAG) continue;
 
@@ -274 +285 @@
     if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE;
 
+    if (COUNT_GHOSTS) {
+        // save nGhosts to update header.
+        psMetadata *updates = psMetadataLookupMetadata (&status, fpa->analysis, "PSASTRO.HEADER");
+        if (!updates) {
+            updates = psMetadataAlloc ();
+            psMetadataAddMetadata (fpa->analysis, PS_LIST_TAIL, "PSASTRO.HEADER",  PS_META_REPLACE, "psastro header stats", updates);
+            psFree (updates);
+        }
+        psMetadataAddS32 (updates, PS_LIST_TAIL, "NGHOSTS", PS_META_REPLACE, "total expected ghosts", nGhosts);
+    }
+
     // deactivate all files
     pmFPAfileActivate (config->files, false, NULL);