SFSetOfHistograms.C

Go to the documentation of this file.

#include "jbltools/sfh/SFSetOfHistograms.h"
#include "jbltools/sfh/BaseCutBase.h"
#include "jbltools/sfh/FloatFunBase.h"
#include "jbltools/sfh/BinningFun.h"
#include "jbltools/sfh/ROList.h"
#include "jbltools/sfh/FillIterator.h"

// system headers
#include <iostream>

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

SFSetOfHistograms::SFSetOfHistograms (const char* name_, 
                                      const char* title_, 
                                      const ROListPoR& hhl,  
                                      const FloatFunPoR& xfun_, 
                                      const BaseCutPoR& cut_, 
                                      const FloatFunPoR& wfun_, 
                                      const BinningFunPoR& binning_,  
                                      const FillIteratorPoR& iter_) 
:  SetOfHistograms  (name_, title_, hhl, binning_),
   xfun (xfun_.pff), cut(cut_.pbc), wfun (wfun_.pff), iter(iter_.pfi)
{
  checkIterators(name_);
}              
    
SFSetOfHistograms::SFSetOfHistograms (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 BinningFunPoR& binning_,  
                                      const FillIteratorPoR& iter_) 
:  SetOfHistograms  (name_, title_, nbinsx, xlow, xup, hhl, binning_),
   xfun (xfun_.pff), cut(cut_.pbc), wfun (wfun_.pff), iter(iter_.pfi)
{
  checkIterators(name_);
}              
    
SFSetOfHistograms::SFSetOfHistograms (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_, 
                                      const BinningFunPoR& binning_,  
                                      const FillIteratorPoR& iter_) 
:  SetOfHistograms  (name_, title_, nbinsx, xbins, hhl, binning_),
   xfun (xfun_.pff), cut(cut_.pbc), wfun (wfun_.pff), iter(iter_.pfi)
{
  checkIterators(name_);
}              
    
SFSetOfHistograms::SFSetOfHistograms (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_, 
                                      const BinningFunPoR& binning_,  
                                      const FillIteratorPoR& iter_) 
:  SetOfHistograms  (name_, title_, nbinsx, xbins, hhl, binning_),
   xfun (xfun_.pff), cut(cut_.pbc), wfun (wfun_.pff), iter(iter_.pfi)
{
  checkIterators(name_);
}              
      
SFSetOfHistograms::SFSetOfHistograms (const char* name_,  
                                      const char* title_, 
                                      const Binning& axisbinning,  
                                      const ROListPoR& hhl,  
                                      const FloatFunPoR& xfun_,  
                                      const BaseCutPoR& cut_,  
                                      const FloatFunPoR& wfun_, 
                                      const BinningFunPoR& binning_,  
                                      const FillIteratorPoR& iter_) 
:  SetOfHistograms  (name_, title_, axisbinning, hhl, binning_),
   xfun (xfun_.pff), cut(cut_.pbc), wfun (wfun_.pff), iter(iter_.pfi)
{
  checkIterators(name_);
}              
      
    
SFSetOfHistograms::~SFSetOfHistograms () {}

void SFSetOfHistograms::Fill () {
  const BinningFun *binning = getBinning();
  assert (binning);
  if (!xfun) return;
//   if (!iter || iter->reset()) {
//     do {
//       if (!cut || (*cut)()) {
//         Float_FF w = wfun ? (*wfun)() : 1.;
//         int ibinlow, ibinup;
//         if (binning->isExclusive()) {
//           ibinup =  (*binning)();
//           ibinlow = ibinup > 0 ? ibinup : 0;
//         }
//         else {
//           ibinlow = 0; 
//           ibinup = binning->getNBins()-1; 
//         }
//         for (int ibin = ibinlow; ibin <= ibinup; ++ibin) {
//           if ((*binning)(ibin)) {
//             if (TH1F *th1f=dynamic_cast<TH1F *>(getEntry(ibin))) th1f->Fill ((*xfun)(), w);
//           } 
//         } 
//       } 
//     } while (iter && iter->next());  
//   }

  if (binning->isExclusive()) {
    if (!iter || iter->reset()) {
      do {
        if (!cut || (*cut)()) {
          Float_FF w = wfun ? (*wfun)() : 1.;
          int ibin = (*binning)();
          if (ibin >= 0) {
            if (TH1F *th1f=dynamic_cast<TH1F *>(getEntry(ibin))) th1f->Fill ((*xfun)(), w);
          } 
        } 
      } while (iter && iter->next());  
    }
  }
  else {
    if (!iter || iter->reset()) {
      do {
        if (!cut || (*cut)()) {
          Float_FF w = wfun ? (*wfun)() : 1.;
          Float_FF x = (*xfun)();
          int ibinup = binning->getNBins(); 
          for (int ibin = 0; ibin < ibinup; ++ibin) {
            if ((*binning)(ibin)) {
              if (TH1F *th1f=dynamic_cast<TH1F *>(getEntry(ibin))) th1f->Fill (x, w);
            } 
          } 
        } 
      } while (iter && iter->next());  
    }
  }
}

template<class Fun1, class Fun2>
void SFSetOfHistograms::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 << "SFSetOfHistograms: Iterators of " << fun1name
              << " and " << fun2name
              << " are different for histogram "
              << name << "!" << std::endl;
  }
  assert (!fun1 || !fun1->getIterator() || !fun2  || !fun2->getIterator()  || fun1->getIterator() == fun2->getIterator());
}                       
template<class Fun>
void SFSetOfHistograms::checkOneIterator (const char *name, 
                                          const char *funname, 
                                                Fun  *fun) {
  if (fun && fun->getIterator()) {
    if (iter) {
      if (fun->getIterator() != iter) {
        std::cerr << "SFSetOfHistograms: Iterator '" 
                  << fun->getIterator()->getName() 
                  << "' of " << funname 
                  << " '" << fun->getName() << "'"
                  << " and iter '" << iter->getName() 
                  << "'are different for histogram "
                  << name << "!\n";
      }
    }
    else {
#if(DEBUG>=2)      
      std::cerr << "SFSetOfHistograms: 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 SFSetOfHistograms::checkIterators(const char *name) {
  const BinningFun *binning = getBinning();
  assert (binning);
  checkTwoIterators (name, "xfun", xfun, "cut",     cut);
  checkTwoIterators (name, "xfun", xfun, "wfun",    wfun);
  checkTwoIterators (name, "xfun", xfun, "binning", binning);
  checkTwoIterators (name, "wfun", wfun, "cut",     cut);
  checkTwoIterators (name, "wfun", wfun, "binning", binning);
  checkTwoIterators (name, "cut",  cut,  "binning", binning);

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

---