Changeset 37881 for trunk/psphot

Timestamp:

Jan 29, 2015, 2:15:11 PM (11 years ago)

Author:

bills

Message:

drop models with major axis > 100 pixels.

Location:

trunk/psphot/src

Files:

• psphotExtendedSourceFits.c (modified) (1 diff)
• psphotGalaxyParams.c (modified) (4 diffs)

trunk/psphot/src/psphotExtendedSourceFits.c

@@ -534 +534 @@
               Nconvolve ++;
               if (!(modelFit->flags & badModel)) {
+                  float *PAR = modelFit->params->data.F32;
+                  psEllipseAxes axes = pmPSF_ModelToAxes (PAR, modelFit->class->useReff);
+                  if (axes.major >= 100) {
+                      Nfail ++;
+                      source->mode2 |= PM_SOURCE_MODE2_EXT_FITS_FAIL;
+                      psFree(modelFit);
+                      continue;
+                  }
                   NconvolvePass ++;
                   source->mode |= PM_SOURCE_MODE_EXTENDED_FIT;

trunk/psphot/src/psphotGalaxyParams.c

@@ -152 +152 @@
 
 // set this to 0 to run without threading
-# if (0)            
+# if (1)            
             if (!psThreadJobAddPending(job)) {
                 psError(PS_ERR_UNKNOWN, false, "Unable to guess model.");
@@ -270 +270 @@
     badModel |= PM_MODEL_PCM_FAIL_GUESS;
 
+    badModel |= PM_MODEL_STATUS_LIMITS;
+
     pmModelType sersicModelType = pmModelClassGetType("PS_MODEL_SERSIC");
 
@@ -334 +336 @@
         if (!sersicModel) SKIP (Nskip8);
         if (sersicModel->flags & badModel) SKIP (Nskip8);
+
+        // check out the model params and skip those that are too large
+        psF32 *PAR = sersicModel->params->data.F32;
+        psF32 Xo = PAR[PM_PAR_XPOS];
+        psF32 Yo = PAR[PM_PAR_YPOS];
+        psF32 Io = PAR[PM_PAR_I0];
+        //
+        // The half light radius to the major axis of the model
+        psEllipseAxes axes = pmPSF_ModelToAxes (PAR, sersicModel->class->useReff);
+        if (axes.major >= 100) SKIP (Nskip8);
+
+        source->extpars->ghalfLightRadius = axes.major;
         Next ++;
 
@@ -353 +367 @@
             // do we need this?
             source->mode |= PM_SOURCE_MODE_NONLINEAR_FIT;
-#ifdef notdef
-#endif
-        }
-
-
-        psF32 *PAR = sersicModel->params->data.F32;
-        psF32 Xo = PAR[PM_PAR_XPOS];
-        psF32 Yo = PAR[PM_PAR_YPOS];
-        psF32 Io = PAR[PM_PAR_I0];
-        // For now: set the half light radius to the major axis of the model
-        psEllipseAxes axes = pmPSF_ModelToAxes (PAR, sersicModel->class->useReff);
-        source->extpars->ghalfLightRadius = axes.major;
+        }
+
+
         // and we will integrate out to 2 * half light radius
         psF32 R = 2 * axes.major;