Changeset 30093

Timestamp:

Dec 17, 2010, 9:18:06 AM (15 years ago)

Author:

eugene

Message:

add functions to write out XRAD (radial apertures)

File:

• branches/eam_branches/ipp-20101205/psModules/src/objects/pmSourceIO.c (modified) (13 diffs)

branches/eam_branches/ipp-20101205/psModules/src/objects/pmSourceIO.c

@@ -59 +59 @@
 
 // lookup the EXTNAME values used for table data and image header segments
-static bool sourceExtensions(psString *headname, // Extension name for header
-                             psString *dataname, // Extension name for data
-                             psString *deteffname, // Extension name for detection efficiency
-                             psString *xsrcname, // Extension name for extended sources
-                             psString *xfitname, // Extension name for extended fits
-                             const pmFPAfile *file, // File of interest
-                             const pmFPAview *view // View to level of interest
-                             )
+bool pmSourceIOextnames(psString *headname,    // Extension name for image header
+                        psString *dataname,    // Extension name for PSF table data
+                        psString *deteffname,  // Extension name for detection efficiency
+                        psString *xsrcname,    // Extension name for extended non-parametric measurements
+                        psString *xfitname,    // Extension name for extended fitted measurements
+                        psString *xradname,    // Extension name for radial apertures
+                        const pmFPAfile *file, // File of interest
+                        const pmFPAview *view  // View to level of interest
+    )
 {
     bool status;                        // Status of MD lookup
@@ -87 +88 @@
     }
 
-    // EXTNAME for table data
+    // EXTNAME for PSF table data
     if (dataname) {
         const char *rule = psMetadataLookupStr(&status, menu, "CMF.DATA");
@@ -107 +108 @@
     }
 
-    // EXTNAME for extended source data table
+    // EXTNAME for extended source non-parametric measurements
     if (xsrcname) {
         const char *rule = psMetadataLookupStr(&status, menu, "CMF.XSRC");
@@ -117 +118 @@
     }
 
+    // EXTNAME for extended source fitted measurements
     if (xfitname) {
-        // EXTNAME for extended source data table
         const char *rule = psMetadataLookupStr(&status, menu, "CMF.XFIT");
         if (!rule) {
@@ -125 +126 @@
         }
         *xfitname = pmFPAfileNameFromRule (rule, file, view);
+    }
+
+    // EXTNAME for radial apertures
+    if (xradname) {
+        const char *rule = psMetadataLookupStr(&status, menu, "CMF.XRAD");
+        if (!rule) {
+            psError(PS_ERR_UNKNOWN, true, "missing entry for CMF.XRAD in EXTNAME.RULES in camera.config");
+            return false;
+        }
+        *xradname = pmFPAfileNameFromRule (rule, file, view);
     }
 
@@ -344 +355 @@
             status &= pmSourcesWrite_##TYPE##_XFIT (file->fits, readout, sources, file->header, xfitname); \
         }                                                               \
+        if (xradname) {                                                 \
+            status &= pmSourcesWrite_##TYPE##_XRAD (file->fits, readout, sources, file->header, xradname, recipe); \
+        }                                                               \
     }
 
@@ -449 +463 @@
         bool XSRC_OUTPUT = psMetadataLookupBool(&status, recipe, "EXTENDED_SOURCE_ANALYSIS");
         bool XFIT_OUTPUT = psMetadataLookupBool(&status, recipe, "EXTENDED_SOURCE_FITS");
+        bool XRAD_OUTPUT = psMetadataLookupBool(&status, recipe, "RADIAL_APERTURES");
 
         // define the EXTNAME values for the different data segments:
@@ -456 +471 @@
         psString xsrcname = NULL;
         psString xfitname = NULL;
-        if (!sourceExtensions(&headname, &dataname, &deteffname, 
-                              XSRC_OUTPUT ? &xsrcname : NULL,
-                              XFIT_OUTPUT ? &xfitname : NULL, 
-                              file, view)) {
+        psString xradname = NULL;
+        if (!pmSourceIOextnames(&headname, &dataname, &deteffname, 
+                                XSRC_OUTPUT ? &xsrcname : NULL,
+                                XFIT_OUTPUT ? &xfitname : NULL, 
+                                XRAD_OUTPUT ? &xradname : NULL, 
+                                file, view)) {
             return false;
         }
@@ -539 +556 @@
                 psMetadataAddStr (outhead, PS_LIST_TAIL, "XFITNAME", PS_META_REPLACE, "name of XFIT table extension", xfitname);
             }
-    
+            if (xradname) {
+                psMetadataAddStr (outhead, PS_LIST_TAIL, "XRADNAME", PS_META_REPLACE, "name of XRAD table extension", xradname);
+            }
 
             // these are case-sensitive since the EXTYPE is case-sensitive
@@ -563 +582 @@
         }
 
-
         // write out the detection efficiency TABLE segments
         if (deteffname) {
@@ -583 +601 @@
         psFree (xsrcname);
         psFree (xfitname);
+        psFree (xradname);
         psFree (deteffname);
 
@@ -593 +612 @@
         psFree (xsrcname);
         psFree (xfitname);
+        psFree (xradname);
         psFree (deteffname);
         return false;
@@ -939 +959 @@
         psString dataname = NULL;
         psString deteffname = NULL;
-        if (!sourceExtensions(&headname, &dataname, &deteffname, NULL, NULL, file, view)) {
+        if (!pmSourceIOextnames(&headname, &dataname, &deteffname, NULL, NULL, NULL, file, view)) {
             return false;
         }