Changeset 36884 for branches/eam_branches/ipp-20140610/psModules/src/objects/pmSourceOutputs.c

Timestamp:

Jun 17, 2014, 5:58:38 AM (12 years ago)

Author:

eugene

Message:

merge in changes from trunk (up through pv3 tag)

Location:

branches/eam_branches/ipp-20140610

Files:

• . (modified) (1 prop)
• psModules/src/objects/pmSourceOutputs.c (modified) (3 diffs)

branches/eam_branches/ipp-20140610/psModules/src/objects/pmSourceOutputs.c

@@ -26 +26 @@
 #include "pmMoments.h"
 #include "pmModelFuncs.h"
+#include "pmModelClass.h"
 #include "pmModel.h"
 #include "pmModelUtils.h"
-#include "pmModelClass.h"
 #include "pmSourceMasks.h"
 #include "pmSourceExtendedPars.h"
@@ -108 +108 @@
         }
         if (isfinite(PAR[PM_PAR_SXX]) && isfinite(PAR[PM_PAR_SXY]) && isfinite(PAR[PM_PAR_SYY])) {
-            axes = pmPSF_ModelToAxes (PAR, model->type);
+            axes = pmPSF_ModelToAxes (PAR, model->class->useReff);
             outputs->psfMajor = axes.major;
             outputs->psfMinor = axes.minor;
             outputs->psfTheta = axes.theta*PS_DEG_RAD;
+
+            // some models (PS1_V1, QGAUSS) have an extra 'core' parameter
+            outputs->psfCore = NAN;
+            if (model->type == pmModelClassGetType ("PS_MODEL_PS1_V1")) {
+                outputs->psfCore = PAR[PM_PAR_7];
+            }
+            if (model->type == pmModelClassGetType ("PS_MODEL_QGAUSS")) {
+                outputs->psfCore = PAR[PM_PAR_7];
+            }
+
+            outputs->psfMajorFWHM = model->class->modelSetFWHM(model->params, axes.major);
+            outputs->psfMinorFWHM = model->class->modelSetFWHM(model->params, axes.minor);
         } else {
             outputs->psfMajor = NAN;
             outputs->psfMinor = NAN;
             outputs->psfTheta = NAN;
+            outputs->psfCore = NAN;
         }
         outputs->chisq = model->chisq;
         outputs->nDOF = model->nDOF;
         outputs->nPix = model->nPix;
-        outputs->apRadius = source->apRadius;
     } else {
         bool useMoments = pmSourcePositionUseMoments(source);
@@ -138 +150 @@
         outputs->psfMinor = NAN;
         outputs->psfTheta = NAN;
+        outputs->psfCore = NAN;
         outputs->chisq = NAN;
         outputs->nDOF = 0;
         outputs->nPix = 0;
-        outputs->apRadius = NAN;
     }