Changeset 33405 for branches/eam_branches/ipp-20111122/Ohana/src/relphot/src/StarOps.c

Timestamp:

Mar 4, 2012, 5:34:35 PM (14 years ago)

Author:

eugene

Message:

remove a bunch of test print statements; use new ubercalDist field in secfilt

File:

• branches/eam_branches/ipp-20111122/Ohana/src/relphot/src/StarOps.c (modified) (7 diffs)

branches/eam_branches/ipp-20111122/Ohana/src/relphot/src/StarOps.c

@@ -56 +56 @@
   int Nsecfilt = GetPhotcodeNsecfilt ();
   Nfew = Nsys = Nbad = Ncal = Nmos = Ngrid = 0;
-
-  FILE *f = fopen ("few.v1.dat", "w");
 
   for (i = 0; i < Ncatalog; i++) {
@@ -133 +131 @@
         // when performing the grid analysis, STAR_TOOFEW will be set to 1;
         if (N <= STAR_TOOFEW) { /* too few measurements */
-          fprintf (f, "%10.6f %10.6f %d %d %d\n", catalog[i].averageT[j].R, catalog[i].averageT[j].D, catalog[i].measureT[catalog[i].averageT[j].measureOffset].imageID, N, STAR_TOOFEW); 
+          // fprintf (f, "%10.6f %10.6f %d %d %d\n", catalog[i].averageT[j].R, catalog[i].averageT[j].D, catalog[i].measureT[catalog[i].averageT[j].measureOffset].imageID, N, STAR_TOOFEW); 
           catalog[i].secfilt[Nsecfilt*j+Nsec].flags |= ID_STAR_FEW;
           Nfew ++;
@@ -150 +148 @@
     }
   }
-  fclose (f);
   fprintf (stderr, "%d stars marked having too few measurements (Nbad: %d, Ncal: %d, Nmos: %d, Ngrid: %d, Nsys: %d)\n", Nfew, Nbad, Ncal, Nmos, Ngrid, Nsys);
 
@@ -288 +285 @@
         catalog[i].secfilt[Nsecfilt*j+Nsec].Nused = stats.Nmeas;
 
-        catalog[i].secfilt[Nsecfilt*j+Nsec].M_80 = stats.Upper80; // M_80 == ubercalDist
-        catalog[i].secfilt[Nsecfilt*j+Nsec].M_20 = stats.Lower20; // M_80 == ubercalDist
-        catalog[i].secfilt[Nsecfilt*j+Nsec].M_80 = minUbercalDist; // M_80 == ubercalDist
+        catalog[i].secfilt[Nsecfilt*j+Nsec].M_80 = stats.Upper80;
+        catalog[i].secfilt[Nsecfilt*j+Nsec].M_20 = stats.Lower20;
+        catalog[i].secfilt[Nsecfilt*j+Nsec].ubercalDist = minUbercalDist;
       }
     }
@@ -577 +574 @@
   TOOFEW = MAX (5, STAR_TOOFEW);
 
-  FILE *f = final ? fopen ("poor.dat", "a") : fopen ("poor.dat", "w");
-
   Ndel = Nave = 0;
   Ncal = Nmos = Ngrid = Nfew = 0;
@@ -675 +670 @@
           float mySigma = hypot (stats.sigma, dlist[k]);
           if (fabs (list[k] - stats.median) > NSIGMA_REJECT*mySigma) {
-            fprintf (f, "%10.6f %10.6f %d %f %f %f %f 0x%08x\n", catalog[i].averageT[j].R, catalog[i].averageT[j].D, catalog[i].measureT[ilist[k]].imageID, list[k], stats.median, stats.sigma, mySigma, catalog[i].measureT[ilist[k]].photFlags); 
+            // fprintf (f, "%10.6f %10.6f %d %f %f %f %f 0x%08x\n", catalog[i].averageT[j].R, catalog[i].averageT[j].D, catalog[i].measureT[ilist[k]].imageID, list[k], stats.median, stats.sigma, mySigma, catalog[i].measureT[ilist[k]].photFlags); 
             catalog[i].measureT[ilist[k]].dbFlags |= ID_MEAS_POOR_PHOTOM;
             if (final) {
@@ -692 +687 @@
     }
   }
-  fclose (f);
   initstats (STATMODE);
   if (VERBOSE) fprintf (stderr, OFF_T_FMT" measures marked poor, "OFF_T_FMT" total\n", Ndel, Nave);