Changeset 13064 for trunk/psModules/src/objects/pmSourceIO_PS1_DEV_0.c

Timestamp:

Apr 27, 2007, 12:14:08 PM (19 years ago)

Author:

magnier

Message:

converting the PSF model fits to the polariation terms (replacing Sx,Sy,Sxy)

File:

• trunk/psModules/src/objects/pmSourceIO_PS1_DEV_0.c (modified) (5 diffs)

trunk/psModules/src/objects/pmSourceIO_PS1_DEV_0.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-04-26 01:20:29 $
+ *  @version $Revision: 1.9 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-04-27 22:14:08 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -53 +53 @@
     int i;
     psF32 *PAR, *dPAR;
-    psEllipseShape shape;
     psEllipseAxes axes;
     psF32 xPos, yPos;
@@ -75 +74 @@
             yErr = dPAR[PM_PAR_YPOS];
 
-            shape.sx  = PAR[PM_PAR_SXX];
-            shape.sy  = PAR[PM_PAR_SYY];
-            shape.sxy = PAR[PM_PAR_SXY];
-            if ((shape.sx == 0) || (shape.sy == 0)) {
-                axes.major = 0.0;
-                axes.minor = 0.0;
-                axes.theta = 0.0;
-            } else {
-                axes = psEllipseShapeToAxes (shape, 20.0);
-            }
+            axes = pmPSF_ModelToAxes (PAR, 20.0);
         } else {
             xPos = source->peak->xf;
@@ -159 +149 @@
     bool status;
     psF32 *PAR, *dPAR;
-    psEllipseShape shape;
     psEllipseAxes axes;
 
@@ -203 +192 @@
         source->errMag    = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_ERR");
 
-        bool skip = false;
-        skip |= (axes.major <= 0);
-        skip |= (axes.minor <= 0);
-        if (!skip) {
-            shape = psEllipseAxesToShape (axes);
-            PAR[PM_PAR_SXX] = shape.sx;
-            PAR[PM_PAR_SYY] = shape.sy;
-            PAR[PM_PAR_SXY] = shape.sxy;
-        } else {
-            PAR[PM_PAR_SXX] = 0.0;
-            PAR[PM_PAR_SYY] = 0.0;
-            PAR[PM_PAR_SXY] = 0.0;
-        }
+        pmPSF_AxesToModel (PAR, axes);
 
         float lflux = psMetadataLookupF32 (&status, row, "PEAK_FLUX_AS_MAG");