Changeset 36860 for trunk/psModules/src/objects/pmSourceOutputs.c

Timestamp:

Jun 13, 2014, 8:47:47 AM (12 years ago)

Author:

eugene

Message:

add psfCore to the outputs structure, assign only for PS1_V1 and QGAUSS models; add apNpixels to source & set with pmSourcePhotometryAper; use new pmModel.class.functions; add AP_NPIX,PSF_CORE,PSF_FWHM_MAJ,PSF_FWHM_MIN to output CMF (PS1_V5+, PS1_DV4+, PS1_SV3+), AP_MAG_RADIUS was S32 on read (how did this work??)

File:

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

trunk/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;
     }