Changeset 2525 for trunk/MagicSoft/Mars/mhist

Timestamp:

11/18/03 12:13:51 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MH.cc (modified) (1 diff)
• MH.h (modified) (1 diff)
• MHCalibrationBlindPixel.cc (added)
• MHCalibrationBlindPixel.h (added)
• MHCalibrationConfig.h (added)
• MHCalibrationPINDiode.cc (added)
• MHCalibrationPINDiode.h (added)
• MHCalibrationPixel.cc (added)
• MHCalibrationPixel.h (added)

trunk/MagicSoft/Mars/mhist/MH.cc

@@ -882 +882 @@
 }
 
+// --------------------------------------------------------------------------
+//
+//  Cuts the bins containing only zeros at the edges.
+//
+//  A new number of bins can be defined with nbins != 0
+//        In the case of nbins == 0, no rebinning will take place
+//
+//  Returns the new (real) number of bins
+//
+Int_t MH::CutEdges(TH1 *h, Int_t nbins)
+{
+
+  TAxis* axe = h->GetXaxis();
+
+  Int_t min1   = axe->GetFirst();
+  Int_t max1   = axe->GetLast();
+  Int_t range1 = max1-min1;
+
+  // 
+  // Check for useless zeros 
+  // 
+  if (range1 == 0)
+    return 0;
+
+  Int_t min2 = 0;
+  Int_t max2 = 0;
+
+  for (int i=min1; i<=max1; i++)
+    {
+      Double_t x = h->GetBinContent(i);
+      if (x != 0.)
+        {
+        min2 = i;
+        break;
+        }
+    }
+  
+  //
+  // If the histogram consists of zeros only
+  //
+  if (min2 == max1) 
+    return 0;
+
+  for (int i=max1; i>=min2; i--)
+    {
+      Double_t x = h->GetBinContent(i);      
+      if (x != 0.)
+        {
+        max2 = i;
+        break;
+        }
+    }
+
+  //
+  // Check for rebinning
+  //
+  if (nbins < 1) 
+    {
+      axe->SetRange(min2,max2);
+      return axe->GetLast()-axe->GetFirst();
+    }
+  
+  //
+  // Appying TAxis->SetRange before ReBin does not work ...
+  // But this workaround helps quite fine 
+  //
+  Axis_t min = h->GetBinLowEdge(min2);
+  Axis_t max = h->GetBinLowEdge(max2)+h->GetBinWidth(max2);
+  
+  Int_t ngroup = (int)((max2-min2)*h->GetNbinsX()/nbins/(max1-min1)); 
+
+  if (ngroup > 1)
+    h->Rebin(ngroup);
+
+  axe->SetRangeUser(min,max);
+
+  return axe->GetLast()-axe->GetFirst();
+}
+
 void MH::ProjectionX(TH1D &dest, const TH2 &src, Int_t firstybin, Int_t lastybin)
 {

trunk/MagicSoft/Mars/mhist/MH.h

@@ -88 +88 @@
     static void FindGoodLimits(Int_t nbins, Int_t &newbins, Double_t &xmin, Double_t &xmax, Bool_t isInteger);
     static Double_t GetMinimumGT(const TH1 &h, Double_t gt=0);
-
+    static Int_t CutEdges(TH1 *h, Int_t nbins);
+    
     static void ProjectionX(TH1D &dest, const TH2 &src, Int_t firstybin=-1, Int_t lastybin=9999);
     static void ProjectionY(TH1D &dest, const TH2 &src, Int_t firstxbin=-1, Int_t lastxbin=9999);