SFH1F.C

Go to the documentation of this file.

#include "jbltools/sfh/SFH1F.h"
#include "jbltools/sfh/FloatFunBase.h"
#include "jbltools/sfh/BaseCutBase.h"
#include "jbltools/sfh/FillIterator.h"

// system headers
#include <iostream>

static const char *ident="@(#)$Id: SFH1F.C,v 1.16 2005/08/23 12:05:44 blist Exp $";

SFH1F::SFH1F (const char* name, 
              const char* title, 
              Int_t nbinsx, Axis_t xlow, Axis_t xup,
              const ROListPoR& hhl, 
              const FloatFunPoR& xfun_, 
              const BaseCutPoR& cut_, 
              const FloatFunPoR& wfun_, 
              const FillIteratorPoR& iter_) 
: RegH1F (name, title, nbinsx, xlow, xup, hhl), 
  xfun(xfun_.pff), cut(cut_.pbc), wfun(wfun_.pff), iter(iter_.pfi) 
{
  checkIterators(name);
}
  
  
SFH1F::SFH1F (const char* name, 
              const char* title, 
              Int_t nbinsx, const Float_t* xbins,
              const ROListPoR& hhl,
              const FloatFunPoR& xfun_, 
              const BaseCutPoR& cut_, 
              const FloatFunPoR& wfun_, 
              FillIteratorPoR iter_) 
: RegH1F (name, title, nbinsx, xbins, hhl), 
  xfun(xfun_.pff), cut(cut_.pbc), wfun(wfun_.pff), iter(iter_.pfi) 
{
  checkIterators(name);
}

SFH1F::SFH1F (const char* name, 
              const char* title, 
              Int_t nbinsx, const Double_t* xbins,
              const ROListPoR& hhl, 
              const FloatFunPoR& xfun_, 
              const BaseCutPoR& cut_, 
              const FloatFunPoR& wfun_, 
              FillIteratorPoR iter_) 
: RegH1F (name, title, nbinsx, xbins, hhl), 
  xfun(xfun_.pff), cut(cut_.pbc), wfun(wfun_.pff), iter(iter_.pfi) 
{
  checkIterators(name);
}

SFH1F::SFH1F (const char* name, 
              const char* title, 
              const Binning& binning,
              const ROListPoR& hhl, 
              const FloatFunPoR& xfun_, 
              const BaseCutPoR& cut_, 
              const FloatFunPoR& wfun_, 
              FillIteratorPoR iter_) 
: RegH1F (name, title, binning, hhl), 
  xfun(xfun_.pff), cut(cut_.pbc), wfun(wfun_.pff), iter(iter_.pfi) 
{
  checkIterators(name);
}


SFH1F::SFH1F (const char* name, 
              const char* title, 
              const TH1F& source, 
              const ROListPoR& hhl,
              const FloatFunPoR& xfun_,          
              const BaseCutPoR& cut_, 
              const FloatFunPoR& wfun_, 
              FillIteratorPoR iter_) 
: RegH1F (name, title, source, hhl), 
  xfun(xfun_.pff), cut(cut_.pbc), wfun(wfun_.pff), iter(iter_.pfi) 
{
  if (name) SetName (name);
  if (title) SetTitle (title);
  // we want something self-filling => erease original histogram content!
  for (int i = 0; i < this->GetNbinsX(); ++i) this->SetBinContent (i, 0.);
}

SFH1F::~SFH1F() {}
  
void SFH1F::Fill() {
  if (!xfun) return;
  if (!iter || iter->reset()) {
    do {
      if (!cut || (*cut)()) {
        Float_FF w = wfun ? (*wfun)() : 1.;
        this->TH1F::Fill ((*xfun)(), w);
      } 
    } while (iter && iter->next());  
  }
}

template<class Fun1, class Fun2>
void SFH1F::checkTwoIterators(const char *name, 
                              const char *fun1name, Fun1 *fun1, 
                              const char *fun2name, Fun2 *fun2) {
  if (fun1 && fun1->getIterator() && 
      fun2  && fun2->getIterator() &&
      fun1->getIterator() != fun2->getIterator()) {
    std::cerr << "SFH1F: Iterators '" << fun1->getIterator()->getName() 
              << "' of " << fun1name
              << " '" << fun1->getName() << "'"
              << " and '" << fun2->getIterator()->getName() << "' of " << fun2name
              << " '" << fun2->getName() << "'"
              << " are different for histogram "
              << name << "!\n";
  }
  assert (!fun1 || !fun1->getIterator() || !fun2  || !fun2->getIterator()  || fun1->getIterator() == fun2->getIterator());
}                       
template<class Fun>
void SFH1F::checkOneIterator (const char *name, 
                              const char *funname, 
                                    Fun  *fun) {
  if (fun && fun->getIterator()) {
    if (iter) {
      if (fun->getIterator() != iter) {
        std::cerr << "SFH1F: Iterator '" 
                  << fun->getIterator()->getName() 
                  << "' of " << funname 
                  << " '" << fun->getName() << "'"
                  << " and iter '" << iter->getName() 
                  << "'are different for histogram "
                  << name << "!\n";
      }
    }
    else {
#if(DEBUG>=2)      
      std::cerr << "SFH1F: Iterator '" 
                << fun->getIterator()->getName() << "' of " 
                << funname
                << " '" << fun->getName() << "'"
                << " used as iterator for histogram "
                << name << ".\n";
#endif
      iter = const_cast<FillIterator *>(fun->getIterator());
    }
  }
  assert (!fun || !fun->getIterator() || fun->getIterator() == iter);
}                       

void SFH1F::checkIterators(const char *name) {
  checkTwoIterators (name, "xfun", xfun, "cut",     cut);
  checkTwoIterators (name, "xfun", xfun, "wfun",    wfun);
  checkTwoIterators (name, "wfun", wfun, "cut",     cut);

  checkOneIterator (name, "xfun",    xfun);
  checkOneIterator (name, "wfun",    wfun);
  checkOneIterator (name, "cut",     cut);
}

---