Changeset 41662 for trunk/Ohana/src/relphot/src/TGroupOps.c

Timestamp:

Jun 14, 2021, 11:02:20 AM (5 years ago)

Author:

eugene

Message:

add code to rationalize the zero points on each pass

File:

• trunk/Ohana/src/relphot/src/TGroupOps.c (modified) (1 diff)

trunk/Ohana/src/relphot/src/TGroupOps.c

@@ -1365 +1365 @@
 }
 
+double get_median_zpt_tgroups (short photcode) {
+
+  double *mlist;
+
+  if (!TGROUP_ZEROPT) return NAN;
+
+  ALLOCATE (mlist, double, NtgroupTimes);
+
+  int N = 0;
+  for (int i = 0; i < NtgroupTimes; i++) {
+    TGroup *tgroup = tgroupTimes[i][0].byCode;
+    for (int j = 0; j < tgroupTimes[i][0].nCode; j++) {
+      if (tgroup[j].photcode != photcode) continue;
+      if (!(tgroup[j].flags & ID_IMAGE_TGROUP_PHOTCAL)) continue;
+      mlist[N] = tgroup[j].McalPSF;
+      N++;
+    }
+  }
+
+  if (N == 0) {
+    free (mlist);
+    return NAN;
+  }
+
+  StatType stats;
+  liststats_setmode (&stats, "MEAN");
+
+  liststats (mlist, NULL, NULL, N, &stats);
+  free (mlist);
+
+  return stats.median; 
+}
+
+void set_median_zpt_tgroups (short photcode, double zpt) {
+
+  if (!TGROUP_ZEROPT) return;
+  if (!isfinite(zpt)) return; // do not break the zero points
+
+  for (int i = 0; i < NtgroupTimes; i++) {
+    TGroup *tgroup = tgroupTimes[i][0].byCode;
+    for (int j = 0; j < tgroupTimes[i][0].nCode; j++) {
+      if (tgroup[j].photcode != photcode) continue;
+      tgroup[j].McalPSF -= zpt;
+      tgroup[j].McalAPER -= zpt;
+    }
+  }
+  return; 
+}
+
 void plot_tgroup_fields (Catalog *catalog) {