SetOfHistograms.C

Go to the documentation of this file.

#include "jbltools/sfh/SetOfHistograms.h"
#include "jbltools/sfh/RegH1F.h"
#include "jbltools/sfh/BinningFun.h"

#include <TFile.h>
#include <TF1.h>
#include <TAxis.h>

#include <iostream>
#include <cmath>
#include <cassert>

static const char *ident="@(#)$Id: SetOfHistograms.C,v 1.22 2006/06/30 11:05:00 blist Exp $";

SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  const ROListPoR& hhl, 
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_)
{
}  

SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  Int_t nbinsx, 
                                  Axis_t xlow, 
                                  Axis_t xup,  
                                  const ROListPoR& hhl, 
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_),
  axisbinning (nbinsx, xlow, xup)
{
  initHistos ();
}  


SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  Int_t nbinsx, 
                                  const Double_t* xbins,  
                                  const ROListPoR& hhl, 
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_),
  axisbinning (nbinsx, xbins)
{
  initHistos ();
}  


SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  Int_t nbinsx, 
                                  const Float_t* xbins,  
                                  const ROListPoR& hhl,
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_),
  axisbinning (nbinsx, xbins)
{
  initHistos ();
}  

SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  const Binning& axisbinning_,  
                                  const ROListPoR& hhl, 
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_),
  axisbinning (axisbinning_)
{
  initHistos ();
}  


SetOfHistograms::SetOfHistograms (const char* name_, 
                                  const char* title_, 
                                  TFile& file, 
                                  const ROListPoR& hhl, 
                                  const BinningFunPoR& binning_) 
: SetOfH (name_, title_, hhl, binning_)
{
  BinningFun *binning = getBinning();
  assert (binning);
  for (IndexType i = 0; i < static_cast<IndexType>(binning->getNBins()); i++) {
    const char *hname = genBinName(i);
    new RegH1F (hname, file, *this);
    delete[] hname;
  }
}  


void SetOfHistograms::initHisto (IndexType i)
{
  const char *htitle = genBinTitle(i);
  const char *hname = genBinName(i);
  RegH1F *h = new RegH1F (hname, htitle, axisbinning, 0);
  registerObject (i, h);
  delete[] hname;
  delete[] htitle;
}


RegH1F *SetOfHistograms::getHisto (IndexType i) {
  return dynamic_cast<RegH1F *>(getEntry(i));
}

RegH1F *SetOfHistograms::getHisto (iterator i) {
  return dynamic_cast<RegH1F *>(getEntry(i));
}

const RegH1F *SetOfHistograms::getHisto (const_iterator i) const {
  return dynamic_cast<const RegH1F *>(getEntry(i));
}


SetOfHistograms::~SetOfHistograms () {
}                    

  
RegH1F *SetOfHistograms::getSummary(Option_t* option_, const ROListPoR& hhl) {
  TString opt = option_;
  opt.ToLower();
  
  int number = 99999;
  if (opt.Contains("m")) number = 99998;
  
  int hnamlen = strlen (getName()) + 5 + 1;
  char *hname = new char[hnamlen];
  snprintf (hname, hnamlen, "%s%5d", getName(), number);
  int htitlen = strlen (getTitle()) + 1 + strlen ("Summary") + 1;
  char *htitle = new char[htitlen];
  snprintf (htitle, htitlen, "%s Summary", getTitle());
  RegH1F *result = 0;
  
  BinningFun *binning = getBinning();
  assert (binning);
//   Int_t nbinsx = binning->getNBins();
//   Axis_t xlow = -0.5;
//   Axis_t xup = nbinsx-0.5;
  result = new RegH1F (hname, htitle, *binning, hhl);
  if (result == 0)  {
    std::cerr << "SetOfHistograms::getSummary: Cannot create histo!"
         << std::endl;
  } else {
    for (const_iterator i = begin(); i != end(); ++i) {
      const RegH1F *h = getHisto (i);
      if (h) {
        Int_t ibin = getIndex(i) + 1;
        Stat_t content = 0;
        Stat_t error2 = 0;
        h->getHistInfo (option_, content, error2);
        result->SetBinContent (ibin, content);
//        result->SetBinError2  (ibin, error2);
//  root TH1 has only SetBinError :(
        result->SetBinError  (ibin, sqrt(error2));
      }
    }
  }
  delete[] hname;
  delete[] htitle;
  
  return result;

}
  
RegH1F *SetOfHistograms::getSummary(Option_t* option_, const char *fname,
                       int parnumber, int hnumber, const char *tcomment, const ROListPoR& hhl) {
                       
  if (fname == 0) return 0;
  if (tcomment == 0) tcomment = "";
  
  int hnamlen = strlen (getName()) + 5 + 1;
  char *hname = new char[hnamlen];
  snprintf (hname, hnamlen, "%s%5d", getName(), hnumber);
  int htitlen = strlen (getTitle()) + 20 + strlen(fname) + strlen(tcomment);
  char *htitle = new char[htitlen];
  snprintf (htitle, htitlen, "%s: Par %d of %s %s", getTitle(), parnumber, fname, tcomment);
  RegH1F *result = 0;
  
  BinningFun *binning = getBinning();
  assert (binning);
  result = new RegH1F (hname, htitle, *binning, hhl);
  if (result == 0)  {
    std::cerr << "SetOfHistograms::getSummary: Cannot create histo!"
         << std::endl;
  } else {
    for (const_iterator i = begin(); i != end(); ++i) {
      const RegH1F *h = getHisto (i);
      if (h) {
        TF1 *f = h->GetFunction (fname);
        if (f) {
          double content = f->GetParameter (parnumber);
          double error = f->GetParError (parnumber);
          Int_t ibin = getIndex(i) + 1;
          result->SetBinContent (ibin, content);
          result->SetBinError  (ibin, error);
        }
      }
    }
  }
  delete[] hname;
  delete[] htitle;
  
  return result;

}
 
RegH1F* SetOfHistograms::getSum (const char* name, const char* title, 
                                 const ROListPoR& hhl, int lowbin, int highbin) const {
  IndexType low = (lowbin < 0) ? 0 : lowbin;
  IndexType high = (highbin < lowbin || highbin < 0) ? getHighestEntry() : highbin;
  const_iterator i = lower_bound(low);
  const RegH1F *h = getHisto (i);
  assert (h);
  RegH1F *result = new RegH1F (name, title, *h, hhl);
  if (result) {
    for (++i; i != upper_bound (high); ++i) result->Add (getHisto (i));
  }
  return result;
}

  

---