Changeset 5860 for trunk/MagicSoft/Mars/mcalib

Timestamp:

01/16/05 17:37:14 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrateData.cc

@@ -143 +143 @@
 //
 MCalibrateData::MCalibrateData(CalibrationMode_t calmode,const char *name, const char *title) 
-    : fGeomCam(NULL),   fBadPixels(NULL), fCalibrations(NULL),
+    : fGeomCam(NULL),   fBadPixels(NULL), fCalibrations(NULL), fCalibUpdate(NULL),
       fQEs(NULL), fSignals(NULL), fCerPhotEvt(NULL), 
       fPedestalFlag(kNo), fSignalType(kPhot), fRenormFactor(1.)
@@ -232 +232 @@
             return kFALSE;
         }
+
+        SetCalibUpdate( fCalibrations );
     }
 
@@ -433 +435 @@
 Bool_t MCalibrateData::UpdateConversionFactors()
 {
+
+  *fLog << inf << GetDescriptor() 
+        << ": Updating Conversion Factors... " << endl;
+
     //
     // For the moment, we use only a dummy zenith for the calibration:
@@ -452 +458 @@
         Float_t calibQEVar    = 0.;
 
+        Float_t calibUpdate   = 1.;
+
         if(fCalibrationMode!=kNone)
           {
-            MCalibrationChargePix &pix = (MCalibrationChargePix&)(*fCalibrations)[pixidx];
-            
-            hiloconv   = pix.GetConversionHiLo   ();
-            hiloconverr= pix.GetConversionHiLoErr();
-            
             if ((*fBadPixels)[pixidx].IsUnsuitable())
               {
@@ -465 +468 @@
               }
             
+            MCalibrationChargePix &pix   = (MCalibrationChargePix&)(*fCalibrations)[pixidx];
+            MCalibrationChargePix &avpix = (MCalibrationChargePix&)fCalibrations->GetAverageArea(0);
+            
+            hiloconv   = pix.GetConversionHiLo   ();
+            hiloconverr= pix.GetConversionHiLoErr();
+            
             calibConv    = pix.GetMeanConvFADC2Phe();
             calibConvVar = pix.GetMeanConvFADC2PheVar();
@@ -471 +480 @@
             MCalibrationQEPix &qe = (MCalibrationQEPix&) (*fQEs)[pixidx];
 
+            if (fCalibUpdate)
+              {
+                MCalibrationChargePix &upix = (MCalibrationChargePix&)(*fCalibUpdate)[pixidx];
+                //
+                // Correct for the possible change in amplification of the individual pixels chain
+                // 
+                const Float_t pixmean = upix.GetConvertedMean();
+                calibUpdate = (pixmean == 0.) ? 1. : pix.GetConvertedMean() / pixmean;
+                //
+                // Correct for global shifts in light emission
+                // 
+                MCalibrationChargePix &ugpix = (MCalibrationChargePix&)fCalibUpdate->GetAverageArea(0);
+                MBadPixelsPix         &ubad  = (MBadPixelsPix&)        fCalibUpdate->GetAverageBadArea(0);
+                
+                //                if (ubad.IsUnsuitable())
+                //                  *fLog << err << "Average Area is unsuitable!!!!" << endl;
+
+                const Float_t globmean = avpix.GetConvertedMean();
+                calibUpdate = (globmean == 0.) ? 1. : ugpix.GetConvertedMean() / globmean;
+              }
+
             switch(fCalibrationMode)
               {
               case kFlatCharge:
                 {
-                  MCalibrationChargePix &avpix = (MCalibrationChargePix&)fCalibrations->GetAverageArea(0);
-                  calibConv    = avpix.GetMean() / (pix.GetMean() * fGeomCam->GetPixRatio(pixidx));
+                  calibConv    = avpix.GetConvertedMean() 
+                              / (pix.GetConvertedMean() * fGeomCam->GetPixRatio(pixidx));
                   calibConvVar = (avpix.GetMeanRelVar() + pix.GetMeanRelVar()) * calibConv * calibConv;
                   if (pix.IsFFactorMethodValid())
@@ -531 +561 @@
                 hiloconv    = 1.;
                 hiloconverr = 0.;
+                calibUpdate = 1.;
                 break;
               } /* switch calibration mode */
@@ -551 +582 @@
           }
 
-        calibConv *= fRenormFactor;
+        calibConv *= fRenormFactor * calibUpdate;
 
         fHiLoConv     [pixidx] = hiloconv;