Changeset 3014 for trunk/MagicSoft/Mars/manalysis/MPedestalPix.cc

Timestamp:

02/04/04 15:01:19 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MPedestalPix.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/manalysis/MPedestalPix.cc

@@ -40 +40 @@
 using namespace std;
 
-MPedestalPix::MPedestalPix() : fHist(NULL), fFlags(0)
+MPedestalPix::MPedestalPix()
 {
   Clear();
 }
 
-MPedestalPix::~MPedestalPix()
-{
-  if (fHist)
-    delete fHist;
-}
-  
 
 // ------------------------------------------------------------------------
@@ -61 +55 @@
   fPedestalRms = -1.;
 
-  fMean       = -1.;
-  fMeanErr    = -1.;
-  fSigma      = -1.;
-  fProb       = -1.;
-
-  CLRBIT(fFlags,kFitted);
-  CLRBIT(fFlags,kFitValid);  
-
 }
 
@@ -75 +61 @@
   fPedestal = 0.;
   fPedestalRms = 0.;
-
-  fHist = new MHPedestalPixel("MHPedestalPixel","Pedestal Histograms Pixel ");
-
-  if (!fHist)
-    *fLog << warn << dbginf << " Could not create MHCalibrationPixel " << endl;
-
-}
-
-void MPedestalPix::FillHists(const Float_t f)
-{
-
-  if (!fHist)
-    return;
-  
-  fHist->FillCharge(f);
-  fHist->FillChargevsN(f);
 }
 
@@ -105 +75 @@
 }
 
-Bool_t MPedestalPix::IsFitted() const 
-{
-  return TESTBIT(fFlags,kFitted);
-}
-
-Bool_t MPedestalPix::IsFitValid() const 
-{
-  return TESTBIT(fFlags,kFitValid);
-}
-
-// --------------------------------------------------------------------------
-//
-// 1) Return if the charge distribution is already succesfully fitted  
-//    or if the histogram is empty
-// 2) Cut the histograms empty edges
-// 3) Fit the histograms with a Gaussian
-// 4) In case of failure print out the fit results
-// 5) Retrieve the results and store them in this class
-//
-void MPedestalPix::FitCharge()
-{
-  
-  //
-  // 1) Return if the charge distribution is already succesfully fitted  
-  //    or if the histogram is empty
-  //
-  if (fHist->IsFitOK() || fHist->IsEmpty())
-    return;
-
-
-  fHist->CutAllEdges();
-      
-  //
-  // 2) Fit the Lo Gain histograms with a Gaussian
-  //
-  if(fHist->FitCharge())
-    SETBIT(fFlags,kFitted);
-  else
-    CLRBIT(fFlags,kFitted);
-
-  //
-  // 6) Retrieve the results and store them in this class
-  //
-  fMean         = fHist->GetChargeMean();
-  fMeanErr      = fHist->GetChargeMeanErr(); 
-  fSigma        = fHist->GetChargeSigma();
-  fSigmaErr     = fHist->GetChargeSigmaErr();
-  fProb         = fHist->GetChargeProb();
-
-  if (CheckFitValidity())
-    SETBIT(fFlags,kFitValid);
-  else
-    CLRBIT(fFlags,kFitValid);
-
-  return;
-}
-
-Bool_t MPedestalPix::CheckFitValidity()
-{
-  
-  if (fProb < 0.001)
-    return kFALSE;
-  
-  return kTRUE;
-  
-}