Changeset 6556 for branches/rel10_ifa/psModules/src/objects/pmSourceIO_SX.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/pmSourceIO_SX.c (modified) (4 diffs)

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

@@ -1 @@
-# include "pmSource.h"
+/** @file  pmSourceIO.c
+ *
+ *  @author EAM, IfA
+ *
+ *  @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 "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 "pmSourceIO.h"
 
 // elixir-mode / sextractor-style output list with fixed line width
@@ -6 +32 @@
 
     psF32 *PAR, *dPAR;
-    float dmag;
+    EllipseShape shape;
+    EllipseAxes axes;
 
     psLine *line = psLineAlloc (110);  // 110 is sextractor line length
@@ -19 +46 @@
     for (int i = 0; i < sources->n; i++) {
         pmSource *source = (pmSource *) sources->data[i];
-        pmModel *model = pmModelSelect (source);
+        pmModel *model = pmSourceSelectModel (source);
         if (model == NULL)
             continue;
@@ -28 +55 @@
         // pmSourceSextractType (source, &type, &flags);
 
-        dmag = dPAR[1] / PAR[1];
+        shape.sx  = PAR[4];
+        shape.sy  = PAR[5];
+        shape.sxy = PAR[6];
+        axes = EllipseShapeToAxes (shape);
 
         psLineInit (line);
-        psLineAdd (line, "%5.2f", 0.0); // should be type
+        psLineAdd (line, "%5.2f",  0.0); // should be type
         psLineAdd (line, "%11.3f", PAR[2]);
         psLineAdd (line, "%11.3f", PAR[3]);
-        psLineAdd (line, "%9.4f", source->fitMag);
-        psLineAdd (line, "%9.4f", dmag);
+        psLineAdd (line, "%9.4f",  source->psfMag);
+        psLineAdd (line, "%9.4f",  source->errMag);
         psLineAdd (line, "%13.4f", PAR[0]);
-        psLineAdd (line, "%9.2f", 0.0); // should be FWHMx
-        psLineAdd (line, "%9.2f", 0.0); // should be FWHMy
-        psLineAdd (line, "%6.1f", 0.0); // should be Theta
-        psLineAdd (line, "%9.4f", 99.999); // should be MAG_ISO
-        psLineAdd (line, "%9.4f", source->apMag);
-        psLineAdd (line, "%4d\n", 0); // should be flags
+        psLineAdd (line, "%9.2f",  axes.major);
+        psLineAdd (line, "%9.2f",  axes.minor);
+        psLineAdd (line, "%6.1f",  axes.theta);
+        psLineAdd (line, "%9.4f",  source->extMag);
+        psLineAdd (line, "%9.4f",  source->apMag);
+        psLineAdd (line, "%4d\n",  0); // should be flags
         fwrite (line->line, 1, line->Nline, f);
     }