BTagCalibrationReader Class Reference

#include <BTagCalibrationReader.h>

Classes
struct	TmpEntry

Public Member Functions
	BTagCalibrationReader ()
	BTagCalibrationReader (const BTagCalibration *c, BTagEntry::OperatingPoint op, std::string measurementType="comb", std::string sysType="central")
double	eval (BTagEntry::JetFlavor jf, float eta, float pt, float discr=0.) const
	~BTagCalibrationReader ()

Protected Member Functions
void	setupTmpData (const BTagCalibration *c)

Protected Attributes
BTagEntry::Parameters	params
std::map< BTagEntry::JetFlavor, std::vector< TmpEntry > >	tmpData_
std::vector< bool >	useAbsEta

Detailed Description

BTagCalibrationReader

Helper class to pull out a specific set of BTagEntry's out of a BTagCalibration. TF1 functions are set up at initialization time.

Definition at line 20 of file BTagCalibrationReader.h.

Constructor & Destructor Documentation

BTagCalibrationReader::BTagCalibrationReader ( )

inline

Definition at line 23 of file BTagCalibrationReader.h.

{}

BTagCalibrationReader::BTagCalibrationReader	(	const BTagCalibration *	c,
		BTagEntry::OperatingPoint	op,
		std::string	measurementType = "comb",
		std::string	sysType = "central"
	)

Definition at line 3 of file BTagCalibrationReader.cc.

References setupTmpData().

                                                               :
  params(BTagEntry::Parameters(op, measurementType, sysType)),
  useAbsEta(true)
{
  setupTmpData(c);
}

BTagCalibrationReader::~BTagCalibrationReader ( )

inline

Definition at line 28 of file BTagCalibrationReader.h.

{}

Member Function Documentation

double BTagCalibrationReader::eval	(	BTagEntry::JetFlavor	jf,
		float	eta,
		float	pt,
		float	discr = 0.
	)		const

Definition at line 13 of file BTagCalibrationReader.cc.

References BTagCalibrationReader::TmpEntry::discrMax, BTagCalibrationReader::TmpEntry::discrMin, alignCSCRings::e, python.tagInventory::entries, eta, BTagCalibrationReader::TmpEntry::etaMax, BTagCalibrationReader::TmpEntry::etaMin, BTagCalibrationReader::TmpEntry::func, i, BTagEntry::OP_RESHAPING, BTagEntry::Parameters::operatingPoint, params, BTagCalibrationReader::TmpEntry::ptMax, BTagCalibrationReader::TmpEntry::ptMin, tmpData_, and useAbsEta.

{
  bool use_discr = (params.operatingPoint == BTagEntry::OP_RESHAPING);
  if (useAbsEta[jf] && eta < 0) {
    eta = -eta;
  }

  // search linearly through eta, pt and discr ranges and eval
  // future: find some clever data structure based on intervals
  const std::vector<TmpEntry> &entries = tmpData_.at(jf);
  for (unsigned i=0; i<entries.size(); ++i) {
    const BTagCalibrationReader::TmpEntry &e = entries.at(i);
    if (
      e.etaMin <= eta && eta < e.etaMax                   // find eta
      && e.ptMin <= pt && pt < e.ptMax                    // check pt
    ){
      if (use_discr) {                                    // discr. reshaping?
        if (e.discrMin <= discr && discr < e.discrMax) {  // check discr
          return e.func.Eval(discr);
        }
      } else {
        return e.func.Eval(pt);
      }
    }
  }

  return 0.;  // default value
}

void BTagCalibrationReader::setupTmpData ( const BTagCalibration *  c )

protected

Definition at line 45 of file BTagCalibrationReader.cc.

References BTagCalibrationReader::TmpEntry::discrMax, BTagEntry::Parameters::discrMax, BTagCalibrationReader::TmpEntry::discrMin, BTagEntry::Parameters::discrMin, python.tagInventory::entries, BTagCalibrationReader::TmpEntry::etaMax, BTagEntry::Parameters::etaMax, BTagCalibrationReader::TmpEntry::etaMin, BTagEntry::Parameters::etaMin, BTagEntry::formula, BTagCalibrationReader::TmpEntry::func, BTagCalibration::getEntries(), i, BTagEntry::Parameters::jetFlavor, BTagEntry::OP_RESHAPING, BTagEntry::Parameters::operatingPoint, params, BTagEntry::params, BTagCalibrationReader::TmpEntry::ptMax, BTagEntry::Parameters::ptMax, BTagCalibrationReader::TmpEntry::ptMin, BTagEntry::Parameters::ptMin, tmpData_, and useAbsEta.

Referenced by BTagCalibrationReader().

{
  useAbsEta = std::vector<bool>(4, true);
  const std::vector<BTagEntry> &entries = c->getEntries(params);
  for (unsigned i=0; i<entries.size(); ++i) {
    const BTagEntry &be = entries[i];
    BTagCalibrationReader::TmpEntry te;
    te.etaMin = be.params.etaMin;
    te.etaMax = be.params.etaMax;
    te.ptMin = be.params.ptMin;
    te.ptMax = be.params.ptMax;
    te.discrMin = be.params.discrMin;
    te.discrMax = be.params.discrMax;

    if (params.operatingPoint == BTagEntry::OP_RESHAPING) {
      te.func = TF1("", be.formula.c_str(),
                    be.params.discrMin, be.params.discrMax);
    } else {
      te.func = TF1("", be.formula.c_str(),
                    be.params.ptMin, be.params.ptMax);
    }

    tmpData_[be.params.jetFlavor].push_back(te);
    if (te.etaMin < 0) {
      useAbsEta[be.params.jetFlavor] = false;
    }
  }
}

Member Data Documentation

BTagEntry::Parameters BTagCalibrationReader::params

protected

Definition at line 47 of file BTagCalibrationReader.h.

Referenced by eval(), and setupTmpData().

std::map<BTagEntry::JetFlavor, std::vector<TmpEntry> > BTagCalibrationReader::tmpData_

protected

Definition at line 48 of file BTagCalibrationReader.h.

Referenced by eval(), and setupTmpData().

std::vector<bool> BTagCalibrationReader::useAbsEta

protected

Definition at line 49 of file BTagCalibrationReader.h.

Referenced by eval(), and setupTmpData().