Changeset 3015 for trunk/MagicSoft/Mars/manalysis/MHPedestalPixel.cc

Timestamp:

02/04/04 15:03:44 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MHPedestalPixel.cc (modified) (13 diffs)

trunk/MagicSoft/Mars/manalysis/MHPedestalPixel.cc

@@ -48 +48 @@
 const Float_t gkSq2Pi = 2.506628274631;
 const Float_t gkProbLimit = 0.01;
+const Int_t   MHPedestalPixel::gkChargeNbins    = 500 ;
+const Int_t   MHPedestalPixel::gkChargevsNbins  = 1000;
+const Axis_t  MHPedestalPixel::gkChargevsNFirst = -0.5;
+const Axis_t  MHPedestalPixel::gkChargevsNLast  = 999.5;
+const Axis_t  MHPedestalPixel::gkChargeFirst    = -0.5;
+const Axis_t  MHPedestalPixel::gkChargeLast     = 499.5;
 
 // --------------------------------------------------------------------------
@@ -53 +59 @@
 // Default Constructor. 
 //
-MHPedestalPixel::MHPedestalPixel(const char *name, const char *title)
+MHPedestalPixel::MHPedestalPixel()
       : fPixId(-1),
-        fChargeNbins(500),
-        fChargeFirst(-0.5),
-        fChargeLast(499.5),
-        fChargevsNbins(1000),
-        fChargevsNFirst(-0.5),
-        fChargevsNLast(999.5),        
         fGausFit(NULL)
 { 
-
-    fName  = name  ? name  : "MHPedestalPixel";
-    fTitle = title ? title : "Fill the accumulated charges and times of all events and perform fits";
-
-    // Create a large number of bins, later we will rebin
-    fHPedestalCharge = new TH1F("HPedestalCharge","Distribution of Summed FADC Pedestal Slices Pixel ",
-                        fChargeNbins,fChargeFirst,fChargeLast);
-    fHPedestalCharge->SetXTitle("Sum FADC Slices");
-    fHPedestalCharge->SetYTitle("Nr. of events");
-    fHPedestalCharge->Sumw2();
-
-    // We define a reasonable number and later enlarge it if necessary
-    fHPedestalChargevsN = new TH1I("HChargevsN","Sum of Charges vs. Event Number Pixel ",
-                           fChargevsNbins,fChargevsNFirst,fChargevsNLast);
-    fHPedestalChargevsN->SetXTitle("Event Nr.");
-    fHPedestalChargevsN->SetYTitle("Sum of FADC slices");
-
-    fHPedestalCharge->SetDirectory(NULL);
-    fHPedestalChargevsN->SetDirectory(NULL);
-
-    Clear();
+  
+  // Create a large number of bins, later we will rebin
+  fHPedestalCharge = new TH1F("HPedestalCharge","Distribution of Summed FADC Pedestal Slices Pixel ",
+                              gkChargeNbins,gkChargeFirst,gkChargeLast);
+  fHPedestalCharge->SetXTitle("Sum FADC Slices");
+  fHPedestalCharge->SetYTitle("Nr. of events");
+  fHPedestalCharge->Sumw2();
+  
+  // We define a reasonable number and later enlarge it if necessary
+  fHPedestalChargevsN = new TH1I("HChargevsN","Sum of Charges vs. Event Number Pixel ",
+                                 gkChargevsNbins,gkChargevsNFirst,gkChargevsNLast);
+  fHPedestalChargevsN->SetXTitle("Event Nr.");
+  fHPedestalChargevsN->SetYTitle("Sum of FADC slices");
+  
+  fHPedestalCharge->SetDirectory(NULL);
+  fHPedestalChargevsN->SetDirectory(NULL);
+  
+  Clear();
 }
 
@@ -91 +88 @@
   delete fHPedestalCharge;
   delete fHPedestalChargevsN;
-
+  
   if (fGausFit)
     delete fGausFit;
@@ -111 +108 @@
   fChargeNdf       = -1;
 
-  fChargeFirst     = -0.5;
-  fChargeLast      = 499.5;
+  fChargeFirst     = gkChargeFirst;
+  fChargeLast      = gkChargeLast;
 
   if (fGausFit)
     delete fGausFit;
 
+  CLRBIT(fFlags,kFitted);
+  CLRBIT(fFlags,kFitOK);
+  CLRBIT(fFlags,kOscillating);
+
   return;
 }
@@ -160 +161 @@
   fPixId = id;
 
-  TString nameQ = TString(fHPedestalCharge->GetName());
-  nameQ += id;
-  fHPedestalCharge->SetName(nameQ.Data());
-
-  TString nameQvsN  = TString(fHPedestalChargevsN->GetName());
-  nameQvsN += id;
-  fHPedestalChargevsN->SetName(nameQvsN.Data());
-
-  TString titleQ = TString(fHPedestalCharge->GetTitle());
-  titleQ += id;
-  fHPedestalCharge->SetTitle(titleQ.Data());
-
-  TString titleQvsN  = TString(fHPedestalChargevsN->GetTitle());
-  titleQvsN += id;
-  fHPedestalChargevsN->SetTitle(titleQvsN.Data());
-
-}
-
-
-Bool_t MHPedestalPixel::SetupFill(const MParList *plist)
-{
- 
-  Reset();
-
-  return kTRUE;
-}
-
-
-
+  fHPedestalCharge->SetName(Form("%s%d", fHPedestalCharge->GetName(), id));
+  fHPedestalChargevsN->SetName(Form("%s%d", fHPedestalChargevsN->GetName(), id));
+  fHPedestalCharge->SetTitle(Form("%s%d", fHPedestalCharge->GetTitle(), id));
+  fHPedestalChargevsN->SetTitle(Form("%s%d", fHPedestalChargevsN->GetTitle(), id));
+
+}
 
 TObject *MHPedestalPixel::DrawClone(Option_t *option) const
@@ -223 +201 @@
   gROOT->SetSelectedPad(NULL);
   
-  TCanvas *c = MakeDefCanvas(this,600,900); 
+  TCanvas *c = MH::MakeDefCanvas(this,600,900); 
   
   c->Divide(1,2);
@@ -233 +211 @@
   if (fHPedestalCharge->Integral() > 0)
     gPad->SetLogy(1);
-  else
-    gPad->SetLogy(0);
 
   fHPedestalCharge->Draw(opt);
@@ -243 +219 @@
   if (fGausFit)
     {
-      if (IsFitOK())
-        fGausFit->SetLineColor(kGreen);          
-      else
-        fGausFit->SetLineColor(kRed);
-      
+      fGausFit->SetLineColor(IsFitOK() ? kGreen : kRed);
       fGausFit->Draw("same");
     }
@@ -265 +237 @@
 
 
+// --------------------------------------------------------------------------
+//
+// 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
+//
 Bool_t MHPedestalPixel::FitCharge(Option_t *option)
 {
@@ -282 +263 @@
   //
   const Stat_t   entries     = fHPedestalCharge->Integral();
-  const Double_t area_guess  = entries/gkSq2Pi;
   const Double_t mu_guess    = fHPedestalCharge->GetBinCenter(fHPedestalCharge->GetMaximumBin());
   const Double_t sigma_guess = mu_guess/15.;
-
-  TString name = TString("GausFit");
-  name += fPixId;
-
-  fGausFit = new TF1(name.Data(),"gaus",rmin,rmax);
+  const Double_t area_guess  = entries/gkSq2Pi/sigma_guess;
+
+  fGausFit = new TF1(Form("%s%d","GausFit ",fPixId),"gaus",rmin,rmax);
 
   if (!fGausFit) 
@@ -312 +290 @@
     {
 
-      Axis_t rtry = fGausFit->GetParameter(1) - 3.0*fGausFit->GetParameter(2);
+      Axis_t rtry = fGausFit->GetParameter(1) - 2.0*fGausFit->GetParameter(2);
       rmin        = (rtry < rmin ? rmin : rtry);
-      rmax        = fGausFit->GetParameter(1) + 3.0*fGausFit->GetParameter(2);
+      rmax        = fGausFit->GetParameter(1) + 2.0*fGausFit->GetParameter(2);
       fGausFit->SetRange(rmin,rmax);  
 
@@ -353 +331 @@
 {
 
-  Int_t nbins = 50;
-
-  CutEdges(fHPedestalCharge,nbins);
+  Int_t nbins = 30;
+
+  MH::CutEdges(fHPedestalCharge,nbins);
 
   fChargeFirst = fHPedestalCharge->GetBinLowEdge(fHPedestalCharge->GetXaxis()->GetFirst());
   fChargeLast  = fHPedestalCharge->GetBinLowEdge(fHPedestalCharge->GetXaxis()->GetLast())
-                      +fHPedestalCharge->GetBinWidth(0);
-
-  CutEdges(fHPedestalChargevsN,0);
+                +fHPedestalCharge->GetBinWidth(0);
+
+  MH::CutEdges(fHPedestalChargevsN,0);
 
 }
@@ -370 +348 @@
   *fLog << all << "Pedestal Fits Pixel: "  << fPixId  << endl;
 
-  if (TESTBIT(fFlags,kFitted))
-    {
-      
-      *fLog << all << "Results of the Summed Charges Fit: "                      << endl;
-      *fLog << all << "Chisquare: "        << fChargeChisquare                   << endl;
-      *fLog << all << "DoF: "              << fChargeNdf                         << endl;
-      *fLog << all << "Probability: "      << fChargeProb                        << endl;
-      *fLog << all << "Results of fit: "   ;
-        if (TESTBIT(fFlags,kFitOK))
-          *fLog << inf << "OK" << endl;
-        else 
-          *fLog << err << "NOT OK" << endl;          
-      *fLog << all                                                               << endl;
-    }
-  else
-    {
-      *fLog << all << "Pedestal Histogram has not yet been fitted"               << endl;      
-      *fLog << all                                                               << endl;      
-    }
+  if (!TESTBIT(fFlags,kFitted))
+    {
+      gLog << "Pedestal Histogram has not yet been fitted" << endl;
+      return;
+    }
+  
+  gLog << "Results of the Summed Charges Fit: "   << endl;
+  gLog << "Chisquare:      " << fChargeChisquare  << endl;
+  gLog << "DoF:            " << fChargeNdf        << endl;
+  gLog << "Probability:    " << fChargeProb       << endl;
+  gLog << "Results of fit: " << (TESTBIT(fFlags,kFitOK) ? "Ok.": "Not OK!") << endl;
     
 }