BTagEntry Class Reference

#include <BTagEntry.h>

Classes
struct	Parameters

Public Types
enum	JetFlavor { FLAV_B =0, FLAV_C =1, FLAV_UDSG =2 }
enum	OperatingPoint { OP_LOOSE =0, OP_MEDIUM =1, OP_TIGHT =2, OP_RESHAPING =3 }

Public Member Functions
	BTagEntry ()
	BTagEntry (const std::string &csvLine, bool validate)
	BTagEntry (const std::string &func, Parameters p)
	BTagEntry (const TF1 *func, Parameters p)
	BTagEntry (const TH1 *histo, Parameters p)
std::string	makeCSVLine () const
	~BTagEntry ()

Static Public Member Functions
static std::string	makeCSVHeader ()
static std::string	trimStr (std::string str)

Public Attributes
std::string	formula
Parameters	params

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Friends
class	boost::serialization::access
template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

BTagEntry

Represents one pt- or discriminator-dependent calibration function.

measurement_type: e.g. comb, ttbar, di-mu, boosted, ... sys_type: e.g. central, plus, minus, plus_JEC, plus_JER, ...

Everything is converted into a function, as it is easiest to store it in a txt or json file.

Definition at line 24 of file BTagEntry.h.

Member Enumeration Documentation

enum BTagEntry::JetFlavor

Enumerator
FLAV_B
FLAV_C
FLAV_UDSG

Definition at line 33 of file BTagEntry.h.

                 {
    FLAV_B=0,
    FLAV_C=1,
    FLAV_UDSG=2,
  };

enum BTagEntry::OperatingPoint

Enumerator
OP_LOOSE
OP_MEDIUM
OP_TIGHT
OP_RESHAPING

Definition at line 27 of file BTagEntry.h.

                      {
    OP_LOOSE=0,
    OP_MEDIUM=1,
    OP_TIGHT=2,
    OP_RESHAPING=3,
  };

Constructor & Destructor Documentation

BTagEntry::BTagEntry ( )

inline

Definition at line 67 of file BTagEntry.h.

References jets_cff::func, trackerHits::histo, AlCaHLTBitMon_ParallelJobs::p, and AlCaHLTBitMon_QueryRunRegistry::string.

{}

BTagEntry::BTagEntry	(	const std::string &	csvLine,
		bool	validate
	)

Definition at line 36 of file BTagEntry.cc.

References begin, end, Exception, connectstrParser::f1, formula, mps_fire::i, params, AlCaHLTBitMon_QueryRunRegistry::string, and trimStr().

                                                            {
  // make tokens
  std::stringstream buff(csvLine);
  std::vector<std::string> vec;
  std::string token;
  while (std::getline(buff, token, ","[0])) {
    token = BTagEntry::trimStr(token);
    if (token.empty()) {
      continue;
    }
    vec.push_back(token);
  }
  if (vec.size() != 11) {
    throw cms::Exception("BTagCalibration")
          << "Invalid csv line; num tokens != 11: "
          << csvLine;
  }

  // clean string values
  char chars[] = " \"\n";
  for (unsigned int i = 0; i < strlen(chars); ++i) {
    vec[1].erase(remove(vec[1].begin(),vec[1].end(),chars[i]),vec[1].end());
    vec[2].erase(remove(vec[2].begin(),vec[2].end(),chars[i]),vec[2].end());
    vec[10].erase(remove(vec[10].begin(),vec[10].end(),chars[i]),vec[10].end());
  }

  // make formula
  formula = vec[10];
  if (validate) {
    TF1 f1("", formula.c_str());  // compile formula to check validity
    if (f1.IsZombie()) {
      throw cms::Exception("BTagCalibration") << "Invalid csv line; formula does not compile: " << csvLine;
    }
  }

  // make parameters
  unsigned op = stoi(vec[0]);
  if (op > 3) {
    throw cms::Exception("BTagCalibration")
          << "Invalid csv line; OperatingPoint > 3: "
          << csvLine;
  }
  unsigned jf = stoi(vec[3]);
  if (jf > 2) {
    throw cms::Exception("BTagCalibration")
          << "Invalid csv line; JetFlavor > 2: "
          << csvLine;
  }
  params = BTagEntry::Parameters(
    BTagEntry::OperatingPoint(op),
    vec[1],
    vec[2],
    BTagEntry::JetFlavor(jf),
    stof(vec[4]),
    stof(vec[5]),
    stof(vec[6]),
    stof(vec[7]),
    stof(vec[8]),
    stof(vec[9])
  );
}

BTagEntry::BTagEntry	(	const std::string &	func,
		BTagEntry::Parameters	p
	)

Definition at line 98 of file BTagEntry.cc.

References Exception, connectstrParser::f1, formula, and jets_cff::func.

                                                                :
  formula(func),
  params(p)
{
  TF1 f1("", formula.c_str());  // compile formula to check validity
  if (f1.IsZombie()) {
    throw cms::Exception("BTagCalibration")
          << "Invalid func string; formula does not compile: "
          << func;
  }
}

BTagEntry::BTagEntry	(	const TF1 *	func,
		BTagEntry::Parameters	p
	)

Definition at line 110 of file BTagEntry.cc.

References Exception.

                                                          :
  formula(std::string(func->GetExpFormula("p").Data())),
  params(p)
{
  if (func->IsZombie()) {
    throw cms::Exception("BTagCalibration")
          << "Invalid TF1 function; function is zombie: "
          << func->GetName();
  }
}

BTagEntry::BTagEntry	(	const TH1 *	histo,
		BTagEntry::Parameters	p
	)

Definition at line 185 of file BTagEntry.cc.

References BTagEntry::Parameters::discrMax, BTagEntry::Parameters::discrMin, Exception, connectstrParser::f1, formula, pileupCalc::nbins, OP_RESHAPING, BTagEntry::Parameters::operatingPoint, params, BTagEntry::Parameters::ptMax, BTagEntry::Parameters::ptMin, th1ToFormulaBinTree(), and th1ToFormulaLin().

                                                          :
  params(p)
{
  int nbins = hist->GetNbinsX();
  TAxis const* axis = hist->GetXaxis();

  // overwrite bounds with histo values
  if (params.operatingPoint == BTagEntry::OP_RESHAPING) {
    params.discrMin = axis->GetBinLowEdge(1);
    params.discrMax = axis->GetBinUpEdge(nbins);
  } else {
    params.ptMin = axis->GetBinLowEdge(1);
    params.ptMax = axis->GetBinUpEdge(nbins);
  }

  // balanced full binary tree height = ceil(log(2*n_leaves)/log(2))
  // breakes even around 10, but lower values are more propable in pt-spectrum
  if (nbins < 15) {
    formula = th1ToFormulaLin(hist);
  } else {
    formula = th1ToFormulaBinTree(hist);
  }

  // compile formula to check validity
  TF1 f1("", formula.c_str());
  if (f1.IsZombie()) {
    throw cms::Exception("BTagCalibration")
          << "Invalid histogram; formula does not compile (>150 bins?): "
          << hist->GetName();
  }
}

BTagEntry::~BTagEntry ( )

inline

Definition at line 72 of file BTagEntry.h.

References makeCSVHeader(), makeCSVLine(), str, AlCaHLTBitMon_QueryRunRegistry::string, and trimStr().

{}

Member Function Documentation

std::string BTagEntry::makeCSVHeader ( )

static

Definition at line 217 of file BTagEntry.cc.

Referenced by BTagCalibration::makeCSV(), and ~BTagEntry().

{
  return "OperatingPoint, "
         "measurementType, "
         "sysType, "
         "jetFlavor, "
         "etaMin, "
         "etaMax, "
         "ptMin, "
         "ptMax, "
         "discrMin, "
         "discrMax, "
         "formula \n";
}

std::string BTagEntry::makeCSVLine

(

)

const

Definition at line 232 of file BTagEntry.cc.

References BTagEntry::Parameters::discrMax, BTagEntry::Parameters::discrMin, BTagEntry::Parameters::etaMax, BTagEntry::Parameters::etaMin, formula, BTagEntry::Parameters::jetFlavor, BTagEntry::Parameters::measurementType, BTagEntry::Parameters::operatingPoint, params, BTagEntry::Parameters::ptMax, BTagEntry::Parameters::ptMin, and BTagEntry::Parameters::sysType.

Referenced by ~BTagEntry().

{
  std::stringstream buff;
  buff << params.operatingPoint
       << ", " << params.measurementType
       << ", " << params.sysType
       << ", " << params.jetFlavor
       << ", " << params.etaMin
       << ", " << params.etaMax
       << ", " << params.ptMin
       << ", " << params.ptMax
       << ", " << params.discrMin
       << ", " << params.discrMax
       << ", \"" << formula
       << "\" \n";
  return buff.str();
}

template<class Archive >

void BTagEntry::serialize ( Archive &  ar, const unsigned int  version  )

private

std::string BTagEntry::trimStr ( std::string  str )

static

Definition at line 250 of file BTagEntry.cc.

References MillePedeFileConverter_cfg::e, and alignCSCRings::s.

Referenced by BTagEntry(), BTagCalibration::readCSV(), and ~BTagEntry().

                                          {
  size_t s = str.find_first_not_of(" \n\r\t");
  size_t e = str.find_last_not_of (" \n\r\t");

  if((std::string::npos == s) || (std::string::npos == e))
    return "";
  else
    return str.substr(s, e-s+1);
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 81 of file BTagEntry.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 81 of file BTagEntry.h.

Member Data Documentation

std::string BTagEntry::formula

Definition at line 78 of file BTagEntry.h.

Referenced by BTagEntry(), and makeCSVLine().

Parameters BTagEntry::params

Definition at line 79 of file BTagEntry.h.

Referenced by BTagCalibration::addEntry(), BTagEntry(), and makeCSVLine().