Changeset 6712 for branches/rel10_ifa/psModules/src/objects/pmSourceIO_RAW.c

Timestamp:

Mar 27, 2006, 4:14:56 PM (20 years ago)

Author:

magnier

Message:

substantial work on the pmFPAfile I/O concepts

File:

• branches/rel10_ifa/psModules/src/objects/pmSourceIO_RAW.c (modified) (9 diffs)

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

@@ +1 @@
+/** @file  pmSourceIO.c
+ *
+ *  @author EAM, IfA
+ *
+ *  @version $Revision: 1.1.2.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-03-28 02:14:56 $
+ *
+ *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
+ *
+ */
+
+#include <stdio.h>
+#include <math.h>
+#include <string.h>
+#include "pslib.h"
+#include "psLine.h"
+#include "psEllipse.h"
+#include "pmHDU.h"
+#include "pmFPA.h"
+#include "pmPeaks.h"
+#include "pmMoments.h"
+#include "pmGrowthCurve.h"
+#include "pmModel.h"
+#include "pmPSF.h"
+#include "pmSource.h"
+#include "pmModelGroup.h"
+#include "pmSourcePhotometry.h"
+#include "pmSourceIO.h"
+
+/***** Text Output Methods *****/
+bool pmSourcesWriteRAW (psArray *sources, char *filename)
+{
+
+    char *name = (char *) psAlloc (strlen(filename) + 10);
+
+    sprintf (name, "%s.psf.dat", filename);
+    pmSourcesWritePSFs (sources, name);
+
+    sprintf (name, "%s.ext.dat", filename);
+    pmSourcesWriteEXTs (sources, name);
+
+    sprintf (name, "%s.nul.dat", filename);
+    pmSourcesWriteNULLs (sources, name);
+
+    sprintf (name, "%s.mnt.dat", filename);
+    pmMomentsWriteText (sources, name);
+
+    psFree (name);
+    return true;
+}
+
 // write the PSF sources to an output file
-bool pmModelWritePSFs (psArray *sources, char *filename)
-{
-
-    double dPos, dMag;
+bool pmSourcesWritePSFs (psArray *sources, char *filename)
+{
+
+    double dPos;
     int i, j;
     FILE *f;
@@ -11 +62 @@
     f = fopen (filename, "w");
     if (f == NULL) {
-        psLogMsg ("pmModelWritePSFs", 3, "can't open output file for moments%s\n", filename);
+        psLogMsg (__func__, 3, "can't open output file for PSFs: %s\n", filename);
         return false;
     }
@@ -18 +69 @@
     for (i = 0; i < sources->n; i++) {
         pmSource *source = (pmSource *) sources->data[i];
-        if (source->type != PM_SOURCE_STAR)
+        if (source->type != PM_SOURCE_TYPE_STAR)
             continue;
         model = source->modelPSF;
@@ -29 +80 @@
         // dPos is positional error, dMag is mag error
         dPos = hypot (dPAR[2], dPAR[3]);
-        dMag = dPAR[1] / PAR[1];
 
         fprintf (f, "%7.1f %7.1f  %7.1f %8.4f  %7.4f %7.4f  ",
-                 PAR[2], PAR[3], PAR[0], source->fitMag, dMag, dPos);
+                 PAR[2], PAR[3], PAR[0], source->psfMag, source->errMag, dPos);
 
         for (j = 4; j < model->params->n; j++) {
@@ -54 +104 @@
 
 // dump the sources to an output file
-bool pmModelWriteEXTs (psArray *sources, char *filename)
-{
-
-    double dPos, dMag;
+bool pmSourcesWriteEXTs (psArray *sources, char *filename)
+{
+
+    double dPos;
     int i, j;
     FILE *f;
@@ -73 +123 @@
         pmSource *source = (pmSource *) sources->data[i];
 
-        if (source->type != PM_SOURCE_EXTENDED)
+        if (source->type != PM_SOURCE_TYPE_EXTENDED)
             continue;
         model = pmSourceMagnitudes (source, NULL, 0.0);
@@ -85 +135 @@
         // XXX these are hardwired for SGAUSS
         dPos = hypot ((dPAR[4] / PAR[4]), (dPAR[5] / PAR[5]));
-        dMag = dPAR[1] / PAR[1];
 
         fprintf (f, "%7.1f %7.1f  %7.1f %8.4f  %7.4f %7.4f  ",
-                 PAR[2], PAR[3], PAR[0], source->fitMag, dMag, dPos);
+                 PAR[2], PAR[3], PAR[0], source->extMag, source->errMag, dPos);
 
         for (j = 4; j < model->params->n; j++) {
@@ -111 +160 @@
 
 // dump the sources to an output file
-bool pmModelWriteNULLs (psArray *sources, char *filename)
+bool pmSourcesWriteNULLs (psArray *sources, char *filename)
 {
 
@@ -132 +181 @@
 
         // skip these sources (in PSF or EXT)
-        if (source->type == PM_SOURCE_STAR)
-            continue;
-        if (source->type == PM_SOURCE_EXTENDED)
+        if (source->type == PM_SOURCE_TYPE_STAR)
+            continue;
+        if (source->type == PM_SOURCE_TYPE_EXTENDED)
             continue;