Changeset 41555 for trunk/Ohana/src/relphot/src/liststats.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/liststats.c (modified) (6 diffs)

trunk/Ohana/src/relphot/src/liststats.c

@@ -26 +26 @@
   stats->Upper80 = NAN;
   stats->Lower20 = NAN;
+  stats->Upper90 = NAN;
+  stats->Lower10 = NAN;
   stats->total   = NAN;
   stats->Nmeas = 0;
@@ -49 +51 @@
   if (N < 1) return (FALSE);
 
+  int N90 = MIN (N-1, 0.9*N);
   int N80 = MIN (N-1, 0.8*N);
   int N20 = MAX (0, 0.2*N);
+  int N10 = MAX (0, 0.1*N);
 
   if (weight) {
@@ -68 +72 @@
   stats[0].Upper80 = value[N80];
   stats[0].Lower20 = value[N20];
+  stats[0].Upper90 = value[N90];
+  stats[0].Lower10 = value[N10];
 
   switch (stats->statmode) {
@@ -116 +122 @@
   WeightedMean |= (stats->statmode == STATS_INNER_WTMEAN);
   WeightedMean |= (stats->statmode == STATS_CHI_INNER_WTMEAN);
+  if (!dvalue) WeightedMean = FALSE;  // warn the user?
 
   /* calculating the per-star offset based on the desired weighting scheme */
@@ -137 +144 @@
         M  += value[i] * weight[i];
         W  +=            weight[i];
-        dM += SQ (weight[i] * dvalue[i]);
+        dM += dvalue ? SQ (weight[i] * dvalue[i]) : 0.0;
         R  += weight[i];
         Nm += 1.0;  
       } 
       Mo  = M / W;
-      dMo = sqrt (dM) / R;
+      dMo = dvalue ? sqrt (dM) / R : NAN;
     }
   } else {
@@ -159 +166 @@
       for (i = ks; i < ke; i++) {
         M   += value[i];
-        dM  += SQ (dvalue[i]);
+        dM  += dvalue ? SQ (dvalue[i]) : 0.0;
         Nm  += 1.0;  
       } 
       Mo = M / Nm;
-      dMo = sqrt (dM) / Nm;
+      dMo = dvalue ? sqrt (dM) / Nm : NAN;
     }
   }