Changeset 3642 for trunk/MagicSoft/Mars/manalysis/MHPedestalPix.cc

Timestamp:

04/04/04 17:50:30 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/manalysis/MHPedestalPix.cc

@@ -27 +27 @@
 //  MHPedestalPix
 //
-//  Histogram class for pedestal analysis. Holds the histogrammed pedestals, 
-//  derives from MHGausEvents, perform Fourier analysis
+//  Histogram class for pedestal analysis. 
+//  Stores and fits the pedestals taken from MPedestalPix on an event-by-event
+//  basis. The results are re-normalized to a value per slice with the formulae:
+//
+// - Mean Pedestal        / slice = Mean Pedestal        / Number slices
+// - Mean Pedestal Error  / slice = Mean Pedestal Error  / Number slices
+// - Sigma Pedestal       / slice = Sigma Pedestal       / Sqrt (Number slices)
+// - Sigma Pedestal Error / slice = Sigma Pedestal Error / Sqrt (Number slices)
+// 
+//  Derives from MHGausEvents, fits the pedestals to a Gaussian and performs
+//  a Fourier analysis.
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -46 +55 @@
 // Default Constructor. 
 //
+// Sets: 
+// - the default number for fNbins        (fgChargeNbins)
+// - the default number for fFirst        (fgChargeFirst)
+// - the default number for fLast         (fgChargeLast)
+//
+// - the default name of the  fHGausHist ("HPedestalCharge")
+// - the default title of the fHGausHist ("Distribution of Summed FADC Pedestal slices Pixel ")
+// - the default x-axis title for fHGausHist ("Sum FADC Slices")
+// - the default y-axis title for fHGausHist ("Nr. of events")
+// - TH1::Sumw2() for fHGausHist
+//
+// Initializes:
+// - fNSlices to 1
+//
 MHPedestalPix::MHPedestalPix(const char *name, const char *title) 
-    : fPixId(-1)
+    : fNSlices(1)
 { 
 
@@ -53 +76 @@
   fTitle = title ? title : "Histogrammed Pedestal events";
 
-  SetChargeNbins();
-  SetChargeFirst();
-  SetChargeLast();
+  SetNbins( fgChargeNbins );
+  SetFirst( fgChargeFirst );
+  SetLast(  fgChargeLast  );
 
   // Create a large number of bins, later we will rebin
@@ -66 +89 @@
 }
 
-MHPedestalPix::~MHPedestalPix()
-{
-}
-
-void MHPedestalPix::Clear(Option_t *o)
+// --------------------------------------------------------------------------
+// 
+// If mean and sigma have not yet been set, returns.
+//
+// Renormalizes the pedestal fit results by the following formulae:
+//
+// - Mean Pedestal        / slice = Mean Pedestal        / Number slices
+// - Mean Pedestal Error  / slice = Mean Pedestal Error  / Number slices
+// - Sigma Pedestal       / slice = Sigma Pedestal       / Sqrt (Number slices)
+// - Sigma Pedestal Error / slice = Sigma Pedestal Error / Sqrt (Number slices)
+// 
+void MHPedestalPix::Renorm()
 {
 
-  fPixId = -1;
-  MHGausEvents::Clear();
-  return;
-}
-
-
-
-void MHPedestalPix::InitBins()
-{
-
-  fHGausHist.SetBins(fChargeNbins,fChargeFirst,fChargeLast);
-
-}
-
-
-void MHPedestalPix::ChangeHistId(Int_t id)
-{
-
-  fPixId = id;
-
-  fHGausHist.SetName(Form("%s%d", fHGausHist.GetName(), id));
-  fHGausHist.SetTitle(Form("%s%d", fHGausHist.GetTitle(), id));
-}
-
-
-void MHPedestalPix::Renorm(const Float_t nslices)
-{
-
-  if (!IsGausFitOK())
+  if (fMean == fMeanErr == fSigma == fSigmaErr == 0.)
     return;
 
-  Float_t sqslices = TMath::Sqrt(nslices);
+  //
+  // One never knows...
+  //
+  if (fNSlices <= 0)
+    return;
   
-  SetMean(GetMean()/nslices);
+  const Float_t sqslices = TMath::Sqrt(fNSlices);
+  
+  SetMean     ( GetMean()    / fNSlices  );
   //
   // Mean error goes with PedestalRMS/Sqrt(entries) -> scale with slices
   // 
-  SetMeanErr(GetMeanErr()/nslices);
+  SetMeanErr  ( GetMeanErr() / fNSlices  );
   //
   // Sigma goes like PedestalRMS -> scale with sqrt(slices)    
   //
-  SetSigma(GetSigma()/sqslices);
+  SetSigma    ( GetSigma()   / sqslices  );
   //
-  // Sigma error goes like PedestalRMS/2.(entries) -> scale with slices
+  // Sigma error goes like PedestalRMS/2.(entries) -> scale with sqrt(slices)
   //
-  SetSigmaErr(GetSigmaErr()/nslices);
+  SetSigmaErr ( GetSigmaErr() / sqslices );
   
 }