RegH2F.C

Go to the documentation of this file.

#include "jbltools/sfh/RegH2F.h"

#include "jbltools/sfh/Binning.h"

#include <cassert>

static const char *ident="@(#)$Id: RegH2F.C,v 1.7 2005/07/13 11:35:50 blist Exp $";

RegH2F::RegH2F (const char* name, const char* title, 
                Int_t nbinsx, Axis_t xlow, Axis_t xup,
                Int_t nbinsy, Axis_t ylow, Axis_t yup, 
                const ROListPoR& hhl) 
: TH2F (name, title, nbinsx, xlow, xup, nbinsy, ylow, yup), 
  RegO (hhl) 
{
  if(fSumw2.fN == 0) Sumw2();
}
  
  
RegH2F::RegH2F (const char* name, const char* title, 
                Int_t nbinsx, const Double_t* xbins,
                Int_t nbinsy, Axis_t ylow, Axis_t yup, 
                const ROListPoR& hhl) 
: TH2F (name, title, nbinsx, xbins, nbinsy, ylow, yup), 
  RegO (hhl) 
{
  if(fSumw2.fN == 0) Sumw2();
}
  
  
RegH2F::RegH2F (const char* name, const char* title, 
                Int_t nbinsx, Axis_t xlow, Axis_t xup,
                Int_t nbinsy, const Double_t* ybins, 
                const ROListPoR& hhl) 
: TH2F (name, title, nbinsx, xlow, xup, nbinsy, ybins), 
  RegO (hhl) 
{
  if(fSumw2.fN == 0) Sumw2();
}
  
  
RegH2F::RegH2F (const char* name, const char* title, 
                Int_t nbinsx, const Double_t* xbins,
                Int_t nbinsy, const Double_t* ybins, 
                const ROListPoR& hhl) 
: TH2F (name, title, nbinsx, xbins, nbinsy, ybins), 
  RegO (hhl) 
{
  if(fSumw2.fN == 0) Sumw2();
}

RegH2F::RegH2F (const char* name, const char* title, 
                Int_t nbinsx, const Float_t* xbins,
                Int_t nbinsy, const Float_t* ybins, 
                const ROListPoR& hhl) 
: TH2F (name, title, nbinsx, xbins, nbinsy, ybins), 
  RegO (hhl) 
{
  if(fSumw2.fN == 0) Sumw2();
}

RegH2F::RegH2F (const char* name, const char* title, 
                const Binning& xbinning,  
                const Binning& ybinning,  
                const ROListPoR& hhl) 
: TH2F (name, title, 
        xbinning.getNBins(), xbinning.getLowerEdge(), xbinning.getUpperEdge(), 
        ybinning.getNBins(), ybinning.getLowerEdge(), ybinning.getUpperEdge()), 
  RegO (hhl) 
{  
  // For an equidistant binning, everything is OK already;
  // for a non-equidistant binning, the number of bins is fine,
  // but we have to adjust the bin edges
  if (!xbinning.isEquidistant()) {
    TAxis *xaxis = GetXaxis();
    assert (xaxis);
    assert (xbinning.getEdges());
    xaxis->Set (xbinning.getNBins(), xbinning.getEdges());
  }
  if (xbinning.hasBinLabels()) {
    TAxis *xaxis = GetXaxis();
    assert (xaxis);
    for (int i = 0; i< xbinning.getNBins(); ++i) {
      if (const char *label = xbinning.getBinLabel (i)) {
        xaxis->SetBinLabel (i+1, label);
      }
    }
    LabelsOption ("u", "X");
  }
  if (!ybinning.isEquidistant()) {
    TAxis *yaxis = GetYaxis();
    assert (yaxis);
    assert (ybinning.getEdges());
    yaxis->Set (ybinning.getNBins(), ybinning.getEdges());
  }
  if (ybinning.hasBinLabels()) {
    TAxis *yaxis = GetYaxis();
    assert (yaxis);
    for (int i = 0; i< ybinning.getNBins(); ++i) {
      if (const char *label = ybinning.getBinLabel (i)) {
        yaxis->SetBinLabel (i+1, label);
      }
    }
    LabelsOption ("h", "Y");
  }

  if(fSumw2.fN == 0) Sumw2();
}


RegH2F::~RegH2F() 
{}

---