Changeset 37924 for trunk/Ohana/src/relphot

Timestamp:

Feb 20, 2015, 3:57:48 PM (11 years ago)

Author:

eugene

Message:

calculate psfQfMax (and perfect) in diff mode; add to dvopsps diffobj

File:

• trunk/Ohana/src/relphot/src/setMrelCatalog.c (modified) (3 diffs)

trunk/Ohana/src/relphot/src/setMrelCatalog.c

@@ -729 +729 @@
     }
 
+    float psfQfMax = 0.0;
+    float psfQfPerfMax = 0.0;
+
     off_t meas = measureOffset;
 
@@ -773 +776 @@
       float Fcal = MagToFlux(-Mcal);
 
+      // these are only used below if IS_DIFF_DB
+      if (IS_DIFF_DB) {
+        if (measure[k].psfQF     > psfQfMax)     psfQfMax     = measure[k].psfQF;
+        if (measure[k].psfQFperf > psfQfPerfMax) psfQfPerfMax = measure[k].psfQFperf;
+      }
+
       // in the calculations below, 
       // ...list gives the error per measurement, wlist gives the weight
@@ -836 +845 @@
     liststats (Fpsflist, dpsflist, wpsflist, Npsf, psfstats);
 
+    if (IS_DIFF_DB) {
+      // for non DIFF_DB, these are set in setMrelAverageExposure
+      secfilt[Nsec].psfQfMax     = psfQfMax;
+      secfilt[Nsec].psfQfPerfMax = psfQfPerfMax;
+    }
+
     secfilt[Nsec].FpsfWrp  = psfstats->mean;
     secfilt[Nsec].dFpsfWrp = psfstats->error;