Changeset 36757 for trunk/psphot/src

Timestamp:

May 21, 2014, 9:37:09 AM (12 years ago)

Author:

bills

Message:

Changes to full force galaxy shapes analysis

• Put parameters for all extended models into cff file
• in psphotFullForce select models to force and those to skip
• skip sources that are very likely to be stars
• clean up a few "hacky" pieces of the implementation

Location:

trunk/psphot/src

Files:

• psphot.h (modified) (1 diff)
• psphotGalaxyShape.c (modified) (3 diffs)

trunk/psphot/src/psphot.h

@@ -435 +435 @@
     float   NSigma;
     float   clampSN;
-#ifdef notdef
-    float fRmajorMin;
-    float fRmajorMax;
-    float fRmajorDel;
-    float fRminorMin;
-    float fRminorMax;
-    float fRminorDel;
-    float snMinAllTrials;
-    float reduceTrialsFactor;
-    bool  reducedTrials;
-#endif
+    int     extModelType;
 } psphotGalaxyShapeOptions;
 

trunk/psphot/src/psphotGalaxyShape.c

@@ -89 +89 @@
     opt->clampSN = psMetadataLookupF32(&status, recipe, "GALAXY_SHAPES_CLAMP_SN"); 
     psAssert (status, "missing GALAXY_SHAPES_NSIGMA");
+    psString modelTypeToSave = psMetadataLookupStr(&status, recipe, "EXT_MODEL_TYPE_FORCE");
+    if (modelTypeToSave && strcmp(modelTypeToSave, "ALL") && strcmp(modelTypeToSave, "BEST")) {
+        opt->extModelType = pmModelClassGetType(modelTypeToSave);
+    } else {
+        opt->extModelType = pmModelClassGetType("PS_MODEL_SERSIC");
+    }
 
 #ifdef notdef
@@ -251 +257 @@
 
         pmModelType modelType = model->type;
+        model->flags = PM_MODEL_STATUS_NONE;
 
         // we are using fitOptions->mode : be sure this makes sense
@@ -322 +329 @@
 
             // Replace modelEXT with the best model from the first of the model fits if one of them is good
-            if (isfinite(PAR[PM_PAR_I0]) && iModel == 0) {
+            if (isfinite(PAR[PM_PAR_I0]) && modelType == opt->extModelType) {
                 psFree (source->modelEXT);