Changeset 6571 for trunk/psphot/src/psphotMagnitudes.c

Timestamp:

Mar 11, 2006, 5:27:13 PM (20 years ago)

Author:

eugene

Message:

major revisions to use the new pmFPAfile,view I/O handling

File:

• trunk/psphot/src/psphotMagnitudes.c (modified) (2 diffs)

trunk/psphot/src/psphotMagnitudes.c

@@ -1 +1 @@
 # include "psphot.h"
-
-// XXX EAM : the apMag should only be calculated for the brighter sources?
-// XXX EAM : SN limit set by user?
-// XXX EAM : masked region should be (optionally) elliptical
-pmModel *pmSourceMagnitudes (pmSource *source, pmPSF *psf, float apRadius) {
-
-    int status;
-    bool isPSF;
-    float x, y;
-    float rflux;
-    float radius;
-    pmModel *model;
-
-    switch (source->type) {
-      case PM_SOURCE_STAR:
-        model = source->modelPSF;
-        if (model == NULL) return NULL;
-        radius = (apRadius > 0) ? apRadius : model->radius;
-        isPSF = true;
-        break;
-
-      case PM_SOURCE_EXTENDED:
-        model = source->modelEXT;
-        if (model == NULL) return NULL;
-        radius = model->radius;
-        isPSF = false;
-        break;
-            
-      default:
-        return NULL;
-    }
-
-    x = model->params->data.F32[2];
-    y = model->params->data.F32[3];
-
-    // replace source flux
-    pmSourceAddModel (source->pixels, source->mask, model, false, false);
-
-    // set aperture mask circle of PSF_FIT_RADIUS 
-    psImageKeepCircle (source->mask, x, y, radius, "OR", PSPHOT_MASK_MARKED);
-
-    // measure object photometry
-    status = pmSourcePhotometry (&source->fitMag, &source->apMag, model, source->pixels, source->mask);
-
-    // for PSFs, correct both apMag and fitMag to same system, consistent with infinite flux star in aperture RADIUS
-    if (isPSF && (psf != NULL)) {
-      if (psf->growth != NULL) {
-          source->apMag += pmGrowthCurveCorrect (psf->growth, model->radius);
-      }
-
-      rflux   = pow (10.0, 0.4*source->fitMag);
-      source->apMag  -= PS_SQR(model->radius)*rflux * psf->skyBias + psf->skySat / rflux;
-      source->fitMag += psPolynomial4DEval (psf->ApTrend, x, y, 0.0, 0.0);
-    }
-
-    // unmask aperture
-    psImageKeepCircle (source->mask, x, y, radius, "AND", ~PSPHOT_MASK_MARKED);
-
-    // subtract object, leave local sky
-    pmSourceSubModel (source->pixels, source->mask, model, false, false);
-
-    if (!status) return NULL;
-    return model;
-}
-
-/*
-aprMag' - fitMag = flux*skySat + r^2*rflux*skyBias + ApTrend(x,y)
-(aprMag - flux*skySat - r^2*rflux*skyBias) - fitMAg = ApTrend(x,y)
-(aprMag - flux*skySat - r^2*rflux*skyBias) = fitMAg + ApTrend(x,y)
-
-*/
 
 bool psphotMagnitudes (psMetadata *config, psArray *sources, pmPSF *psf) {
@@ -88 +17 @@
 }
 
+// XXX should this function use RADIUS or should it