Changeset 39515 for trunk

Timestamp:

Apr 6, 2016, 3:05:19 PM (10 years ago)

Author:

eugene

Message:

save numerical values for DEV and EXP model fits; test for INDEX, use DEV and EXP index values if INDEX is missing

File:

• trunk/Ohana/src/addstar/src/loadgalphot_readstars.c (modified) (4 diffs)

trunk/Ohana/src/addstar/src/loadgalphot_readstars.c

@@ -61 +61 @@
   // in xfit, we have strings.  these are listed with their numerical 
   // match in the PHU.  I need to read the list of names and set up a string hash
+
+  int ModelTypeDev = -1;
+  int ModelTypeExp = -1;
 
   int ModelTypeN;
@@ -98 +101 @@
     if (!gfits_scan (&header_xfit, name, "%d", 1, &tmpnumber)) { myAbort ("fail"); }
     ModelTypeNum[i] = tmpnumber;
+    if (!strcasecmp (ModelTypeName[i], "PS_MODEL_EXP")) {
+      ModelTypeExp = ModelTypeNum[i];
+    }
+    if (!strcasecmp (ModelTypeName[i], "PS_MODEL_DEV")) {
+      ModelTypeDev = ModelTypeNum[i];
+    }
   }
 
@@ -117 +126 @@
   GET_COLUMN (xfit, INDEX,          "EXT_PAR_07",    float);
   int Nfit = Nrow;
+
+  float DEV_INDEX_VALUE = 0.5 / 4.0;
+  float EXP_INDEX_VALUE = 0.5 / 1.0;
+
+  // note that INDEX here is defined as 1.0 / (2 \nu), where \nu is the traditional sersic
+  // expression.
 
   myAssert (NcharModel < 256, "max model type name is very long: %d", NcharModel);
@@ -247 +262 @@
     stars[i].galphot.theta     = EXT_THETA[ifit];
     stars[i].galphot.thetaErr  = EXT_THETA_ERR[ifit];
-    stars[i].galphot.index     = INDEX[ifit];
     stars[i].galphot.Npix      = NPIX[i];
     stars[i].galphot.modelType = MODEL_TYPE_gal[i];
+    if (INDEX) {
+      stars[i].galphot.index     = INDEX[ifit];
+    } else {
+      if (MODEL_TYPE_gal[i] == ModelTypeDev) { stars[i].galphot.index = DEV_INDEX_VALUE; }
+      if (MODEL_TYPE_gal[i] == ModelTypeExp) { stars[i].galphot.index = EXP_INDEX_VALUE; }
+    }
     stars[i].galphot.detID     = ID_gal[i];
     stars[i].galphot.photcode  = options->photcode;