Changeset 38675 for branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/UpdateObjects.c

Timestamp:

Aug 11, 2015, 9:00:49 AM (11 years ago)

Author:

eugene

Message:

trying to get 2mass correction right

File:

• branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/UpdateObjects.c (modified) (3 diffs)

branches/eam_branches/ipp-20150625/Ohana/src/relastro/src/UpdateObjects.c

@@ -20 +20 @@
   FitStats *fitStatsChips = FitStatsInit (NmeasureMax, N_BOOTSTRAP_SAMPLES);
   FitStats *fitStatsStack = FitStatsInit (NmeasureMax, N_BOOTSTRAP_SAMPLES);
+
+  AstromErrorSetLoop (Nloop);
 
   int i;
@@ -464 +466 @@
   // I've already allocated fit->points (and fit->sample) with space for fit->NpointsAlloc entries
 
+  int has2MASS = FALSE;
+
   int Npoints = fit->Npoints = 0;
   FitAstromPoint *points = fit->points;
@@ -555 +559 @@
     Npoints++;
 
+    if ((measure[0].photcode >= 2011) && (measure[0].photcode <= 2013)) {
+      has2MASS = TRUE;
+    }
+
     myAssert (Npoints <= fit->NpointsAlloc, "oops");
   } // loop over measurements : average[0].Nmeasure 
+
+  // XXX flag measurements from stars with 2MASS 
+  for (k = 0; k < average[0].Nmeasure; k++) {
+    // reset the bit to note that a detection was used (or not)
+    if (has2MASS) {
+      measure[k].dbFlags &= ~ID_MEAS_POOR_PHOTOM;
+    } else {
+      measure[k].dbFlags |=  ID_MEAS_POOR_PHOTOM;
+    }
+  }
 
   fit->Npoints = Npoints;