Changeset 41637 for branches/eam_branches/relphot.20210521/src/GridOps.c

Timestamp:

Jun 4, 2021, 10:46:36 AM (5 years ago)

Author:

eugene

Message:

add user-control over night, mosaic, image chisq & scatter cuts

File:

• branches/eam_branches/relphot.20210521/src/GridOps.c (modified) (6 diffs)

branches/eam_branches/relphot.20210521/src/GridOps.c

@@ -60 +60 @@
     GridCorr[i] = NULL;
   }
+
+  fprintf (stderr, "Generating grid corrections for %d photcodes\n", NGridCorr);
   return;
 }
@@ -102 +104 @@
     result = GridCorr[i];
     *Nlast = i;
+    break;
   }
   return result;
@@ -127 +130 @@
   Image *images = getimages (&Nimages, NULL);
 
+  int NGridReal = 0;
+
   for (int i = 0; i < Nimages; i++) {
-    int code = images[0].photcode;
+    int code = images[i].photcode;
     myAssert (code >= 0, "oops");
     myAssert (code < NGridCorr, "oops");
@@ -137 +142 @@
 
     if (GridCorr[code]) continue; // already created this one
+
+    NGridReal ++;
 
     ALLOCATE(GridCorr[code], GridCorrectionType, 1); 
@@ -157 +164 @@
   }
   resetMgrid(); // start with values of 0
+  fprintf (stderr, "Init grid corrections for %d photcodes\n", NGridReal);
 }
 
@@ -274 +282 @@
         GridCorr[code]-> Mgrid[ix][iy] = Mgrid;
         GridCorr[code]->dMgrid[ix][iy] = sqrt(r*(Mgrid2 - Mgrid*Mgrid)); // sample stdev
-        fprintf (stderr, "grid code %d, %d x %d : %f +/- %f : %d\n", code, ix, iy,
-                 GridCorr[code]-> Mgrid[ix][iy], GridCorr[code]->dMgrid[ix][iy], GridCorr[code]->nMgrid[ix][iy]);
+        // fprintf (stderr, "grid code %d, %d x %d : %f +/- %f : %d\n", code, ix, iy, GridCorr[code]-> Mgrid[ix][iy], GridCorr[code]->dMgrid[ix][iy], GridCorr[code]->nMgrid[ix][iy]);
       }
     }