Changeset 5029 for trunk/MagicSoft/Mars/mhcalib

Timestamp:

09/15/04 14:27:20 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhcalib

Files:

• MHCalibrationCam.cc (modified) (1 diff)
• MHCalibrationChargeCam.cc (modified) (18 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc

@@ -1226 +1226 @@
       {
         hist.BypassFit();
-        bad.SetUncalibrated( fittyp );
+        if (pix.IsHiGainSaturation())
+          bad.SetUncalibrated( fittyp );
       }
   

trunk/MagicSoft/Mars/mhcalib/MHCalibrationChargeCam.cc

@@ -931 +931 @@
 
   TH1F *h = NULL;
+
+  MCalibrationCam *chargecam = fIntensCam ? fIntensCam->GetCam() : fCam;
+  MBadPixelsCam   *badcam    = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
       
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
@@ -940 +943 @@
         continue;
       
-      MCalibrationChargePix  &pix    = fIntensCam 
-        ? (MCalibrationChargePix&)(*fIntensCam)[i] 
-        : (MCalibrationChargePix&)(*fCam)[i];
+      MCalibrationChargePix  &pix = (MCalibrationChargePix&)(*chargecam)[i] ;
 
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
@@ -952 +953 @@
         }
 
-      MBadPixelsPix &bad = fIntensBad ? (*fIntensBad)[i] : (*fBadPixels)[i];
+      MBadPixelsPix &bad = (*badcam)[i];
 
       h = histhi.GetHGausHist();
@@ -982 +983 @@
         
         MHCalibrationChargePix &histlo = (MHCalibrationChargePix&)(*this)(i);
-        MBadPixelsPix &bad = fIntensBad ? (*fIntensBad)[i] : (*fBadPixels)[i];
+        MBadPixelsPix          &bad    = (*badcam)[i];
         
         if (histlo.IsExcluded())
@@ -1016 +1017 @@
           }
         
-        MCalibrationChargePix &pix = fIntensCam 
-          ? (MCalibrationChargePix&)(*fIntensCam)[i]
-          : (MCalibrationChargePix&)(*fCam)[i];
+        MCalibrationChargePix  &pix    = (MCalibrationChargePix&)(*chargecam)[i] ;
         
         if (pix.IsHiGainSaturation())
@@ -1028 +1027 @@
       
       MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)GetAverageHiGainArea(j);      
-      MCalibrationChargePix  &pix    = fIntensCam 
-        ? (MCalibrationChargePix&)fIntensCam->GetAverageArea(j)
-        : (MCalibrationChargePix&)fCam->GetAverageArea(j);
+      MCalibrationChargePix  &pix    = (MCalibrationChargePix&)chargecam->GetAverageArea(j);
 
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
@@ -1040 +1037 @@
         }
 
-      MBadPixelsPix &bad    = fIntensCam 
-        ? fIntensCam->GetAverageBadArea(j)
-        : fCam->GetAverageBadArea(j);
-
+      MBadPixelsPix &bad = chargecam->GetAverageBadArea(j);
       FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
    }
@@ -1060 +1054 @@
           }
         
-        MCalibrationChargePix  &pix    = fIntensCam
-          ? (MCalibrationChargePix&)fIntensCam->GetAverageArea(j) 
-          : (MCalibrationChargePix&)fCam->GetAverageArea(j) ;
+        MCalibrationChargePix  &pix    = (MCalibrationChargePix&)chargecam->GetAverageArea(j);
 
         if (pix.IsHiGainSaturation())
           {
-            MBadPixelsPix &bad    = fIntensCam
-              ? fIntensCam->GetAverageBadArea(j)  
-              : fCam->GetAverageBadArea(j);      
+            MBadPixelsPix &bad = chargecam->GetAverageBadArea(j);
             FinalizeAbsTimes(histlo, pix, bad, fFirstLoGain, fLastLoGain);
           }
@@ -1078 +1068 @@
       
       MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)GetAverageHiGainSector(j);      
-      MCalibrationChargePix  &pix    = fIntensCam
-        ? (MCalibrationChargePix&)fIntensCam->GetAverageSector(j)  
-        : (MCalibrationChargePix&)fCam->GetAverageSector(j);
+      MCalibrationChargePix  &pix    = (MCalibrationChargePix&)chargecam->GetAverageSector(j);
 
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
@@ -1090 +1078 @@
         }
 
-      MBadPixelsPix &bad = fIntensCam 
-        ? fIntensCam->GetAverageBadSector(j)
-        : fCam->GetAverageBadSector(j);      
-
+      MBadPixelsPix &bad = chargecam->GetAverageBadSector(j);
       FinalizeAbsTimes(histhi, pix, bad, fFirstHiGain, fLastHiGain);
     }
@@ -1102 +1087 @@
         
         MHCalibrationChargePix &histlo = (MHCalibrationChargePix&)GetAverageLoGainSector(j);      
-        MBadPixelsPix          &bad    = fIntensCam
-          ? fIntensCam->GetAverageBadSector(j)
-          : fCam->GetAverageBadSector(j);        
+        MBadPixelsPix          &bad    = chargecam->GetAverageBadSector(j);
         
         if (histlo.GetSaturated() > fNumLoGainSaturationLimit*histlo.GetHGausHist()->GetEntries())
@@ -1115 +1098 @@
           }
         
-        MCalibrationChargePix  &pix    = fIntensCam 
-          ? (MCalibrationChargePix&)fIntensCam->GetAverageSector(j)
-          : (MCalibrationChargePix&)fCam->GetAverageSector(j);
+        MCalibrationChargePix  &pix    = (MCalibrationChargePix&)chargecam->GetAverageSector(j);
 
         if (pix.IsHiGainSaturation())
@@ -1126 +1107 @@
   // Perform the fitting for the High Gain (done in MHCalibrationCam)
   //
-  FitHiGainArrays(fIntensCam ? (MCalibrationCam&)(*fIntensCam->GetCam()) : (MCalibrationCam&)(*fCam),
-                  fIntensBad ?                   (*fIntensBad->GetCam()) : *fBadPixels,
+  FitHiGainArrays(*chargecam, *badcam, 
                   MBadPixelsPix::kHiGainNotFitted,
                   MBadPixelsPix::kHiGainOscillating);
@@ -1135 +1115 @@
   //
   if (IsLoGain())
-    FitLoGainArrays(fIntensCam ? (MCalibrationCam&)(*fIntensCam->GetCam()) : (MCalibrationCam&)(*fCam),
-                    fIntensBad ?                   (*fIntensBad->GetCam()) : *fBadPixels,
+    FitLoGainArrays(*chargecam, *badcam, 
                     MBadPixelsPix::kLoGainNotFitted,
                     MBadPixelsPix::kLoGainOscillating);
@@ -1195 +1174 @@
 {
 
-  MBadPixelsCam *badcam = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
+  MBadPixelsCam   *badcam    = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
+  MCalibrationCam *chargecam = fIntensCam ? fIntensCam->GetCam() : fCam;
       
   for (Int_t i=0; i<badcam->GetSize(); i++)
@@ -1201 +1181 @@
       
       MBadPixelsPix    &bad    = (*badcam)[i];
-      MCalibrationPix  &pix    = fIntensCam ? (*fIntensCam)[i] : (*fCam)[i];
+      MCalibrationPix  &pix    = (*chargecam)[i];
 
       if (bad.IsUncalibrated( MBadPixelsPix::kHiGainNotFitted ))
@@ -1210 +1190 @@
         if (pix.IsHiGainSaturation())
           bad.SetUnsuitable(   MBadPixelsPix::kUnreliableRun    );
- 
+
       if (bad.IsUncalibrated( MBadPixelsPix::kLoGainSaturation ))
           bad.SetUnsuitable(   MBadPixelsPix::kUnsuitableRun    );
@@ -1224 +1204 @@
         }
     }
+
+
+      
 }