Changeset 36318 for branches/eam_branches/ipp-20130904/psphot/src/psphotExtendedSourceFits.c

Timestamp:

Nov 28, 2013, 4:09:45 AM (13 years ago)

Author:

eugene

Message:

make sure the convolution scale is passed to all functions that need it; update sersic search a bit

File:

• branches/eam_branches/ipp-20130904/psphot/src/psphotExtendedSourceFits.c (modified) (3 diffs)

branches/eam_branches/ipp-20130904/psphot/src/psphotExtendedSourceFits.c

@@ -120 +120 @@
     }
 
+    float fitNsigmaConv = psMetadataLookupF32 (&status, recipe, "EXT_FIT_NSIGMA_CONV"); // number of sigma for the convolutio
+    if (!status || !isfinite(fitNsigmaConv) || fitNsigmaConv <= 0) {
+        fitNsigmaConv = 5.0;
+    }
+
     bool chisqConvergence = psMetadataLookupBool (&status, recipe, "LMM_FIT_CHISQ_CONVERGENCE"); // Fit tolerance
     if (!status) {
@@ -144 +149 @@
     fitOptions->minTol         = fitMinTol;
     fitOptions->maxTol         = fitMaxTol;
+    fitOptions->nsigma         = fitNsigmaConv;
 
     fitOptions->gainFactorMode   = gainFactorMode;
@@ -656 +662 @@
         if (source->modelEXT->isPCM) {
             // fprintf (stderr, "subtract PCM extended source model for object %d @ %f, %f\n", source->id, source->moments->Mx, source->moments->My);
-            pmPCMCacheModel (source, maskVal, psfSize);
+            pmPCMCacheModel (source, maskVal, psfSize, fitOptions->nsigma);
         } else {
             // fprintf (stderr, "subtract non-PCM extended source model for object %d @ %f, %f\n", source->id, source->moments->Mx, source->moments->My);