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

Timestamp:

Nov 13, 2008, 1:50:19 PM (17 years ago)

Author:

eugene

Message:

make region files optional

File:

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

trunk/psastro/src/psastroMaskUpdates.c

@@ -41 +41 @@
     bool REFSTAR_MASK                      = psMetadataLookupBool (&status, recipe, "REFSTAR_MASK");
     if (!REFSTAR_MASK) return true;
+
+    bool REFSTAR_MASK_REGIONS              = psMetadataLookupBool (&status, recipe, "REFSTAR_MASK_REGIONS");
 
     double REFSTAR_MASK_MAX_MAG            = psMetadataLookupF32 (&status, recipe, "REFSTAR_MASK_MAX_MAG");
@@ -112 +114 @@
         if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE;
 
-        char *filename = NULL;
-        char *chipname = psMetadataLookupStr (&status, chip->concepts, "CHIP.NAME");
-        psStringAppend (&filename, "refstars.mask.%s.dat", chipname);
-        FILE *f = fopen (filename, "w");
-        if (!f) {
+        // text region files for testing
+        FILE *f = NULL;
+        if (REFSTAR_MASK_REGIONS) {
+          char *filename = NULL;
+          char *chipname = psMetadataLookupStr (&status, chip->concepts, "CHIP.NAME");
+          psStringAppend (&filename, "refstars.mask.%s.dat", chipname);
+          FILE *f = fopen (filename, "w");
+          if (!f) {
             psWarning ("cannot create refstar mask file %s\n", filename);
             continue;
-        }
-        psFree (filename);
+          }
+          psFree (filename);
+        }
 
         pmChip *refChip  = pmFPAviewThisChip (view, refMask->fpa);
@@ -176 +182 @@
                     // CIRCLE around the stars (scaled by magnitude)
                     float radius = REFSTAR_MASK_SATSTAR_MAG_SLOPE * (REFSTAR_MASK_SATSTAR_MAG_MAX - ref->Mag);
-                    fprintf (f, "CIRCLE %f %f  %f %f\n", ref->chip->x, ref->chip->y, radius, radius);
+
+                    if (REFSTAR_MASK_REGIONS) {
+                      fprintf (f, "CIRCLE %f %f  %f %f\n", ref->chip->x, ref->chip->y, radius, radius);
+                    }
 
                     // XXX for now, assume cell binning is 1x1 relative to chip
@@ -192 +201 @@
                         float Theta = POSANGLE - REFSTAR_MASK_SATSTAR_POS_ZERO + theta;
 
-                        // lower side 
-                        x0 = ref->chip->x + spikeWidth*sin(Theta);
-                        y0 = ref->chip->y - spikeWidth*cos(Theta);
-                        x1 = ref->chip->x + spikeLength*cos(Theta) + spikeWidth*sin(Theta);
-                        y1 = ref->chip->y + spikeLength*sin(Theta) - spikeWidth*cos(Theta);
-                        dx = x1 - x0;
-                        dy = y1 - y0;
-                        fprintf (f, "LINE %f %f  %f %f\n", x0, y0, dx, dy);
-
-                        // upper side 
-                        x0 = ref->chip->x - spikeWidth*sin(Theta);
-                        y0 = ref->chip->y + spikeWidth*cos(Theta);
-                        x1 = ref->chip->x + spikeLength*cos(Theta) - spikeWidth*sin(Theta);
-                        y1 = ref->chip->y + spikeLength*sin(Theta) + spikeWidth*cos(Theta);
-                        dx = x1 - x0;
-                        dy = y1 - y0;
-                        fprintf (f, "LINE %f %f  %f %f\n", x0, y0, dx, dy);
+                        if (REFSTAR_MASK_REGIONS) {
+                          // lower side 
+                          x0 = ref->chip->x + spikeWidth*sin(Theta);
+                          y0 = ref->chip->y - spikeWidth*cos(Theta);
+                          x1 = ref->chip->x + spikeLength*cos(Theta) + spikeWidth*sin(Theta);
+                          y1 = ref->chip->y + spikeLength*sin(Theta) - spikeWidth*cos(Theta);
+                          dx = x1 - x0;
+                          dy = y1 - y0;
+
+                          fprintf (f, "LINE %f %f  %f %f\n", x0, y0, dx, dy);
+
+                          // upper side 
+                          x0 = ref->chip->x - spikeWidth*sin(Theta);
+                          y0 = ref->chip->y + spikeWidth*cos(Theta);
+                          x1 = ref->chip->x + spikeLength*cos(Theta) - spikeWidth*sin(Theta);
+                          y1 = ref->chip->y + spikeLength*sin(Theta) + spikeWidth*cos(Theta);
+                          dx = x1 - x0;
+                          dy = y1 - y0;
+                          fprintf (f, "LINE %f %f  %f %f\n", x0, y0, dx, dy);
+                        }
 
                         if (readoutMask) {
@@ -219 +231 @@
 
                     // LINE for boundaries of the bleed lines
-                    fprintf (f, "LINE %f %f  %f %f\n", ref->chip->x, ref->chip->y, 0.0, -100.0);
+                    if (REFSTAR_MASK_REGIONS) {
+                      fprintf (f, "LINE %f %f  %f %f\n", ref->chip->x, ref->chip->y, 0.0, -100.0);
+                    }
+
                     if (readoutMask && refCell) {
                         float xCell = 0.0;
@@ -253 +268 @@
 
         if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE;
-        fclose (f);
+        if (REFSTAR_MASK_REGIONS) {
+          fclose (f);
+        }
     }
     if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE;