Changeset 3960 for trunk/MagicSoft/Mars/mcalib/MCalibrationRelTimeCalc.cc

Timestamp:

05/04/04 16:53:25 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationRelTimeCalc.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationRelTimeCalc.cc

@@ -264 +264 @@
   const UInt_t nareas   = fGeom->GetNumAreas();
 
-  Float_t lowlim      [nareas];
-  Float_t upplim      [nareas];
-  Float_t areasum     [nareas];
-  Float_t areasum2    [nareas];
-  Int_t   numareavalid[nareas];
+  Float_t lowlim       [nareas];
+  Float_t upplim       [nareas];
+  Float_t areasum      [nareas];
+  Float_t areasum2     [nareas];
+  Int_t   numareavalid [nareas];
+  Int_t   useunreliable[nareas];
 
   memset(lowlim        ,0, nareas   * sizeof(Float_t));
@@ -275 +276 @@
   memset(areasum2      ,0, nareas   * sizeof(Float_t));
   memset(numareavalid  ,0, nareas   * sizeof(Int_t  ));
-
+  memset(useunreliable ,0, nareas   * sizeof(Int_t  ));
+
+  //
+  // Apero loop: Count number of unreliable pixels:
+  //     
+  for (UInt_t i=0; i<npixels; i++)
+    {
+      MBadPixelsPix &bad = (*fBadPixels)[i];      
+      const Int_t  aidx  = (*fGeom)[i].GetAidx();
+
+      if (bad.IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
+        continue;
+
+      if (bad.IsUnsuitable(MBadPixelsPix::kUnreliableRun))
+        continue;
+
+      numareavalid[aidx] ++;
+    }
+  
+  for (UInt_t aidx=0; aidx<nareas; aidx++)
+    if (numareavalid[aidx] < 100)
+      useunreliable[aidx] = 1;
+
+  memset(numareavalid  ,0, nareas   * sizeof(Int_t  ));
   //
   // First loop: Get mean time resolution the RMS
@@ -292 +316 @@
         continue;
 
+      const Int_t   aidx  = (*fGeom)[i].GetAidx();
+
+      if (!useunreliable[aidx])
+        if (bad.IsUnsuitable(MBadPixelsPix::kUnreliableRun))
+          continue;
+
       const Float_t res   = pix.GetTimePrecision();
-      const Int_t   aidx  = (*fGeom)[i].GetAidx();
 
       areasum     [aidx] += res;
@@ -316 +345 @@
       lowlim  [aidx]  = areasum [aidx] - fRelTimeRelErrLimit*areasum2[aidx];
       upplim  [aidx]  = areasum [aidx] + fRelTimeRelErrLimit*areasum2[aidx];
-
+      
+      *fLog << endl;
       *fLog << inf << GetDescriptor() << ": Limits for acceptance of time resolution: ["
-            << Form("%4.2f%s%4.2f",lowlim[aidx],",",upplim[aidx]) << "] in area index " << aidx << endl;
-    }
-  
+            << Form("%4.2f%s%4.2f",lowlim[aidx],",",upplim[aidx]) << "] in area index " << aidx;
+    }
+  *fLog << endl;