Changeset 41960 for trunk/Ohana/src/relphot

Timestamp:

Dec 10, 2021, 3:37:47 PM (5 years ago)

Author:

cclin33

Message:

fix average magnitude issue on different photcode

File:

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

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

@@ -179 +179 @@
     // these bits should not be set unless we use them in this pass
     // (note that we can only un-set them in the final pass when we actually have measure, not just measureT
-    if (measure) measure[k].dbFlags &= ~(ID_MEAS_PHOTOM_PSF | ID_MEAS_PHOTOM_APER | ID_MEAS_PHOTOM_KRON);
+    // CCL updated 2021.12.09:  the following line will reset the dbFlags for measure so I comment it out.
+    // if (measure) measure[k].dbFlags &= ~(ID_MEAS_PHOTOM_PSF | ID_MEAS_PHOTOM_APER | ID_MEAS_PHOTOM_KRON);
 
     // skip measurements that do not have a valid photcode (raise exception?)
@@ -322 +323 @@
         results->aperData[Nsec].errlist[Nap] = dMap;
         results->aperData[Nsec].wgtlist[Nap] = modifiedWeight;
-        results->aperData[Nsec].ranking[Nap] = measureRank[k];
+        results->aperData[Nsec].ranking[Nap] = measureRank ? measureRank[k] : 0;
         results->aperData[Nsec].measSeq[Nap] = k;
         results->aperData[Nsec].msklist[Nap] = 0;
@@ -383 +384 @@
   float psfQfPerfMax = 0.0;
 
-  // now calculate the mean stats for the Nsec bands.
-  for (Nsec = 0; Nsec < Nsecfilt; Nsec++) {
+  // now calculate the mean stats for the selected Nsec bands.
+  for (int Ns = 0; Ns < Nphotcodes; Ns++) {
+    int thisCode = photcodes[Ns][0].code;
+    int Nsec = GetPhotcodeNsec(thisCode);
 
     // -preserve-ps1 means keep an existing average PS1 value
@@ -477 +480 @@
       secfilt[Nsec].dMkronChp = kronstats->error; 
       secfilt[Nsec].sMkronChp = kronstats->sigma; 
-      secfilt[Nsec].NusedKron = NrankingKRON;
+      // actually used is Nmeas, NrankingKRON is including masked data
+      secfilt[Nsec].NusedKron = kronstats->Nmeas;
     }
 
@@ -521 +525 @@
         secfilt[Nsec].dMapChp = apstats->error; 
         secfilt[Nsec].sMapChp = apstats->sigma; 
-        secfilt[Nsec].NusedAp = NrankingAPER;
+        // actually used is Nmeas, NrankingAPER is including masked data
+        secfilt[Nsec].NusedAp = apstats->Nmeas;
       }
       int minRankAper = (NrankingAPER > 0) ? results->aperData[Nsec].ranking[0] : 10;