EtaPtBin.cc

Go to the documentation of this file.

#include "DQMOffline/RecoB/interface/EtaPtBin.h"

#include <algorithm>
#include <sstream>

EtaPtBin::EtaPtBin(const bool& etaActive_,
                   const double& etaMin_,
                   const double& etaMax_,
                   const bool& ptActive_,
                   const double& ptMin_,
                   const double& ptMax_)
    : etaActive(etaActive_), etaMin(etaMin_), etaMax(etaMax_), ptActive(ptActive_), ptMin(ptMin_), ptMax(ptMax_) {
  descriptionString = buildDescriptionString(etaActive, etaMin, etaMax, ptActive, ptMin, ptMax);
}

std::string EtaPtBin::buildDescriptionString(const bool& etaActive_,
                                             const double& etaMin_,
                                             const double& etaMax_,
                                             const bool& ptActive_,
                                             const double& ptMin_,
                                             const double& ptMax_) {
  // create string only from the active parts
  std::stringstream stream("");

  if (etaActive_) {
    stream << "_ETA_" << etaMin_ << "-" << etaMax_;
  }

  if (ptActive_) {
    stream << "_PT_" << ptMin_ << "-" << ptMax_;
  }
  if (!(etaActive_ || ptActive_))
    stream << "_GLOBAL";

  std::string descr(stream.str());
  // remove blanks which are introduced when adding doubles
  std::remove(descr.begin(), descr.end(), ' ');
  std::replace(descr.begin(), descr.end(), '.', 'v');

  return descr;
}

bool EtaPtBin::inBin(const reco::Jet& jet, const double jec) const { return inBin(jet.eta(), jet.pt() * jec); }

// bool EtaPtBin::inBin(const BTagMCTools::JetFlavour & jetFlavour) const
// {
//   return inBin(jetFlavour.underlyingParton4Vec().Eta(),
//         jetFlavour.underlyingParton4Vec().Pt());
// }

bool EtaPtBin::inBin(const double& eta, const double& pt) const {
  if (etaActive) {
    if (fabs(eta) < etaMin)
      return false;
    if (fabs(eta) > etaMax)
      return false;
  }

  if (ptActive) {
    if (pt < ptMin)
      return false;
    if (pt > ptMax)
      return false;
  }

  return true;
}