Changeset 41555 for trunk/Ohana/src/relphot/src/ZeroPointModes.c

Timestamp:

Apr 24, 2021, 2:28:09 PM (5 years ago)

Author:

eugene

Message:

add function to measure the residual stats for each chip, fix resets and flags on resets; fix update of image mcal from tgroup and mosaic; improved determination of bad nights and bad mosaics; allow tgroups to fit only zero points without airmass terms; clean up the test detection dumps; soften the errors to avoid overweighting extreme outliers with high S/N; more careful sequencing of fitting modes; add manual iteration option

File:

• trunk/Ohana/src/relphot/src/ZeroPointModes.c (modified) (9 diffs)

trunk/Ohana/src/relphot/src/ZeroPointModes.c

@@ -4 +4 @@
 
 void SetZptIteration (int current) {
+  // XXX add 'manual' option here: ask for desired zpt iteration
+  if (MANUAL_ITERATION) {
+    int newValue;
+    fprintf (stdout, "Set Iteration (current = %d): ", CurrentLoop);
+    fscanf (stdin, "%d", &newValue);
+    fprintf (stdout, "Using iteration %d\n", newValue);
+    CurrentLoop = newValue;
+    return;
+  }
+
+  // automatic iteration:
   CurrentLoop = current;
 }
@@ -14 +25 @@
 
   if (TGROUP_ZEROPT && MOSAIC_ZEROPT) {
-    if ((CurrentLoop > 4) && (CurrentLoop <= 8)) {
+    if (CurrentLoop < 6) {
+      TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_ALL;
+      IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_NONE;
+      MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_NONE;
+    }
+    if ((CurrentLoop >= 6) && (CurrentLoop < 12)) {
       // after iterating a few times on the TGroup zero points:
       // * identify the photometric nights
@@ -23 +39 @@
       clean_tgroups(); // do this on every pass to update the status of nights 
     }
-    if ((CurrentLoop > 8) && (CurrentLoop <= 999)) {
+    if ((CurrentLoop >= 12) && (CurrentLoop < 18)) {
       // after iterating a few times on the TGroup & Mosaic zero points:
       // * identify the good / bad Mosaics
@@ -33 +49 @@
       clean_mosaics(); // do this on every pass 
     }
+    if ((CurrentLoop >= 18) && (CurrentLoop <= 999)) {
+      // after iterating a few times on the TGroup & Mosaic zero points:
+      // * identify the good / bad Mosaics
+      // * fit the chips from the bad Mosaics
+      TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_GOOD_NIGHT; // stop fitting the bad nights
+      MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_BAD_NIGHT_GOOD_MOSAIC; // only fit good mosaics on bad nights
+      IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_BAD_NIGHT_BAD_MOSAIC; // only fit bad mosaics on bad nights
+      clean_tgroups(); // do this on every pass 
+      clean_mosaics(); // do this on every pass 
+      clean_images(); // do this on every pass 
+    }
   }
 
   if (TGROUP_ZEROPT && !MOSAIC_ZEROPT) {
+    if (CurrentLoop <= 4) {
+      TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_ALL;
+      IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_NONE;
+      MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_NONE;
+    }
     if ((CurrentLoop > 4) && (CurrentLoop <= 999)) {
       // after iterating a few times on the TGroup zero points:
@@ -48 +80 @@
 
   if (!TGROUP_ZEROPT && MOSAIC_ZEROPT) {
+    if (CurrentLoop <= 4) {
+      TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_NONE;
+      MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_ALL;
+      IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_NONE;
+    }
     if ((CurrentLoop > 4) && (CurrentLoop <= 999)) {
       // after iterating a few times on the Mosaic zero points:
@@ -60 +97 @@
 
   if (!TGROUP_ZEROPT && !MOSAIC_ZEROPT) {
-    if ((CurrentLoop > 4) && (CurrentLoop <= 999)) {
-      // after iterating a few times on the Image zero points:
-      // * ???
-      TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_NONE;
-      MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_NONE;
-      IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_ALL;
-    }
+    // if we are not fitting TGROUP or MOSAIC, just fit all images
+    TGROUP_ZPT_MODE = TGROUP_ZPT_MODE_NONE;
+    MOSAIC_ZPT_MODE = MOSAIC_ZPT_MODE_NONE;
+    IMAGE_ZPT_MODE  = IMAGE_ZPT_MODE_ALL;
   }
 }
@@ -83 +117 @@
         return FALSE;
       case ZPT_MOSAIC:
+        if (CurrentLoop < 9) return TRUE;
+        return FALSE;
+      case ZPT_IMAGES:
+      default:
+        if (CurrentLoop < 15) return TRUE;
+        return FALSE;
+    }
+  }
+
+  if (TGROUP_ZEROPT && !MOSAIC_ZEROPT) {
+    switch (mode) {
+      case ZPT_STARS: 
+        if (CurrentLoop < 3) return TRUE;
+        return FALSE;
+      case ZPT_TGROUP: 
+      case ZPT_MOSAIC:
+        if (CurrentLoop < 3) return TRUE;
+        return FALSE;
+      case ZPT_IMAGES:
+      default:
         if (CurrentLoop < 7) return TRUE;
         return FALSE;
-      case ZPT_IMAGES:
-      default:
-        if (CurrentLoop < 11) return TRUE;
-        return FALSE;
-    }
-  }
-
-  if (TGROUP_ZEROPT && !MOSAIC_ZEROPT) {
+    }
+  }
+
+  if (!TGROUP_ZEROPT && MOSAIC_ZEROPT) {
     switch (mode) {
       case ZPT_STARS: 
@@ -108 +158 @@
   }
 
-  if (!TGROUP_ZEROPT && MOSAIC_ZEROPT) {
-    switch (mode) {
-      case ZPT_STARS: 
-        if (CurrentLoop < 3) return TRUE;
-        return FALSE;
-      case ZPT_TGROUP: 
-      case ZPT_MOSAIC:
-        if (CurrentLoop < 3) return TRUE;
-        return FALSE;
-      case ZPT_IMAGES:
-      default:
-        if (CurrentLoop < 7) return TRUE;
-        return FALSE;
-    }
-  }
-
   if (!TGROUP_ZEROPT && !MOSAIC_ZEROPT) {
     switch (mode) {
@@ -139 +173 @@
   return FALSE;
 }
+
+/*  TGROUP & MOSAIC:
+    0 - TGROUP ONLY & OLS
+    1 - TGROUP ONLY & OLS
+    2 - TGROUP ONLY & OLS
+
+    3 - TGROUP ONLY & IRLS
+    4 - TGROUP ONLY & IRLS
+    5 - TGROUP ONLY & IRLS
+
+    6 - TGROUP + MOSAIC & OLS MOSAIC
+    7 - TGROUP + MOSAIC & OLS MOSAIC
+    8 - TGROUP + MOSAIC & OLS MOSAIC
+
+    9 - TGROUP + MOSAIC & IRLS MOSAIC
+   10 - TGROUP + MOSAIC & IRLS MOSAIC
+   11 - TGROUP + MOSAIC & IRLS MOSAIC
+
+   12 - TGROUP + MOSAIC + IMAGE & OLS IMAGE
+   13 - TGROUP + MOSAIC + IMAGE & OLS IMAGE
+   14 - TGROUP + MOSAIC + IMAGE & OLS IMAGE
+
+   15 - TGROUP + MOSAIC + IMAGE & IRLS IMAGE
+   16 - TGROUP + MOSAIC + IMAGE & IRLS IMAGE
+   17 - TGROUP + MOSAIC + IMAGE & IRLS IMAGE
+
+*/
+