Changeset 6556 for branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c

Timestamp:

Mar 8, 2006, 5:14:23 PM (20 years ago)

Author:

magnier

Message:

major rework of objects code

File:

• branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c (modified) (8 diffs)

branches/rel10_ifa/psModules/src/objects/pmFPAviewWriteObjects.c

@@ -1 @@
-# include "pmFPAviewWriteObjects.h"
+/** @file  pmSourcePhotometry.c
+ *
+ *  @author EAM, IfA; GLG, MHPCC
+ *
+ *  @version $Revision: 1.1.2.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-09 03:14:23 $
+ *
+ *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ *
+ */
+
+#include <stdio.h>
+#include <math.h>
+#include <string.h>
+#include "pslib.h"
+#include "pmHDU.h"
+#include "pmFPA.h"
+#include "pmFPAview.h"
+#include "pmFPAfile.h"
+
+#include "pmPeaks.h"
+#include "pmMoments.h"
+#include "pmModel.h"
+#include "pmSource.h"
+
+#include "pmSourceIO.h"
+#include "pmFPAviewObjectsIO.h"
 
 // Given a FITS file pointer, read the table of object data
@@ -42 +68 @@
 
 // read in all chip-level Objects files for this FPA
-pmFPAWriteObjects (pmFPA *fpa, pmFPAview *view, pmFPAfile *file)
+bool pmFPAWriteObjects (pmFPA *fpa, pmFPAview *view, pmFPAfile *file)
 {
 
     for (int i = 0; i < fpa->chips->n; i++) {
 
-        pmChip *chip = fpa->data[i];
+        pmChip *chip = fpa->chips->data[i];
         pmChipWriteObjects (chip, view, file);
     }
@@ -54 +80 @@
 
 // read in all cell-level Objects files for this chip
-pmChipWriteObjects (pmChip *chip, pmFPAview *view, pmFPAfile *file)
+bool pmChipWriteObjects (pmChip *chip, pmFPAview *view, pmFPAfile *file)
 {
 
     for (int i = 0; i < chip->cells->n; i++) {
 
-        pmCell *cell = chip->data[i];
+        pmCell *cell = chip->cells->data[i];
         pmCellWriteObjects (cell, view, file);
     }
@@ -66 +92 @@
 
 // read in all readout-level Objects files for this cell
-pmCellWriteObjects (pmCell *cell, pmFPAview *view, pmFPAfile *file)
+bool pmCellWriteObjects (pmCell *cell, pmFPAview *view, pmFPAfile *file)
 {
 
     for (int i = 0; i < cell->readouts->n; i++) {
 
-        pmReadout *readout = cell->data[i];
+        pmReadout *readout = cell->readouts->data[i];
         pmReadoutWriteObjects (readout, view, file);
     }
@@ -78 +104 @@
 
 // read in all readout-level Objects files for this cell
-pmReadoutWriteObjects (pmReadout *readout, pmFPAview *view, pmFPAfile *file)
+bool pmReadoutWriteObjects (pmReadout *readout, pmFPAview *view, pmFPAfile *file)
 {
+
+    bool status;
+    char *filename;
+    char *dataname;
+    char *headname;
+    pmHDU *hdu;
 
     psArray *sources = psMetadataLookupPtr (&status, readout->analysis, "PSPHOT.SOURCES");
@@ -85 +117 @@
     switch (file->type) {
     case PM_FPA_FILE_OBJ:
-        filename = pmConfigNameFromRule (file->filerule, view);
-        pmSourcesWriteOBJ (filename, sources);
+        filename = pmFPAviewNameFromRule (file->filerule, view);
+        pmSourcesWriteOBJ (sources, filename);
         break;
 
     case PM_FPA_FILE_SX:
-        filename = pmConfigNameFromRule (file->filerule, view);
-        pmSourcesWriteSX (filename, sources);
+        filename = pmFPAviewNameFromRule (file->filerule, view);
+        pmSourcesWriteSX (sources, filename);
         break;
 
     case PM_FPA_FILE_CMP:
         // a SPLIT format : only one header and object table per file
-        pmHDU *hdu = pmFPAviewThisHDU (view);
-        filename = pmConfigNameFromRule (file->filerule, view);
-        pmSourcesWriteCMP (filename, sources, hdu->header);
+        hdu = pmFPAviewThisHDU (view);
+        filename = pmFPAviewNameFromRule (file->filerule, view);
+        pmSourcesWriteCMP (sources, filename, hdu->header);
         break;
 
     case PM_FPA_FILE_CMF:
         // write header, if not yet written?
-        pmHDU *hdu = pmFPAviewThisHDU (view);
+        hdu = pmFPAviewThisHDU (view);
 
         // XXX is this needed? is it automatically added?
@@ -111 +143 @@
         if (hdu->header) {
             psMetadataItem *mdi = NULL;
-            headname = pmConfigNameFromRule (file->extxtra, view);
+            headname = pmFPAviewNameFromRule (file->extxtra, view);
 
             // set NAXIS to 0 (we don't write out the data array)
@@ -119 +151 @@
             mdi->type = PS_DATA_S32;
 
-            psFitsWriteHeader (file->fits, hdu->header, headname);
+            if (strcasecmp (headname, "PHU")) {
+                psMetadataAdd (hdu->header, PS_LIST_HEAD, "EXTNAME", PS_DATA_STRING, "extension name", headname);
+            }
+            psFitsWriteHeader (hdu->header, file->fits);
         }
-        dataname = pmConfigNameFromRule (file->extrule, view);
+        dataname = pmFPAviewNameFromRule (file->extrule, view);
         pmSourcesWriteCMF (file->fits, sources, hdu->header, dataname);
+        break;
+
+    default:
+        fprintf (stderr, "warning: type mismatch\n");
         break;
     }