QGLikelihoodCalculator Class Reference

#include <QGLikelihoodCalculator.h>

Public Member Functions
float	computeQGLikelihood (edm::ESHandle< QGLikelihoodObject > &QGLParamsColl, float pt, float eta, float rho, std::vector< float > vars) const
	Compute likelihood for a jet using the QGLikelihoodObject information and a set of variables. More...
	QGLikelihoodCalculator ()
float	systematicSmearing (edm::ESHandle< QGLikelihoodSystematicsObject > &QGLParamsColl, float pt, float eta, float rho, float qgValue, int qgIndex) const
	~QGLikelihoodCalculator ()

Private Member Functions
const QGLikelihoodObject::Entry *	findEntry (std::vector< QGLikelihoodObject::Entry > const &data, float eta, float pt, float rho, int qgIndex, int varIndex) const
	Find matching entry in vector for a given eta, pt, rho, qgIndex and varIndex. More...
bool	isValidRange (float pt, float rho, float eta, const QGLikelihoodCategory &qgValidRange) const
	Check the valid range of this qg tagger. More...
float	smearingFunction (float x0, float a, float b, float min, float max) const
	Return the smeared qgLikelihood value, given input x0 and parameters a, b, min and max. More...

Detailed Description

The QGLikelihoodCalculater calculates the likelihood for a jet It takes information on the valid range of the tool, the binning of the categories, and their PDFs from the QGLikelihoodObject The variables in the vars vector should match with the variables in the QGLikelihoodObject, in which they are identified by the varIndex Authors: andre.nosp@m.a.ca.nosp@m.rlo.m.nosp@m.arin.nosp@m.i@cer.nosp@m.n.ch, tom.c.nosp@m.orne.nosp@m.lis@c.nosp@m.ern..nosp@m.ch, cms-q.nosp@m.g-wo.nosp@m.rking.nosp@m.grou.nosp@m.p@cer.nosp@m.n.ch

Definition at line 13 of file QGLikelihoodCalculator.h.

Constructor & Destructor Documentation

QGLikelihoodCalculator()

QGLikelihoodCalculator::QGLikelihoodCalculator ( )

inline

Definition at line 15 of file QGLikelihoodCalculator.h.

{};

~QGLikelihoodCalculator()

QGLikelihoodCalculator::~QGLikelihoodCalculator ( )

inline

Definition at line 16 of file QGLikelihoodCalculator.h.

{};

Member Function Documentation

computeQGLikelihood()

float QGLikelihoodCalculator::computeQGLikelihood	(	edm::ESHandle< QGLikelihoodObject > &	QGLParamsColl,
		float	pt,
		float	eta,
		float	rho,
		std::vector< float >	vars
	)		const

Compute likelihood for a jet using the QGLikelihoodObject information and a set of variables.

Definition at line 8 of file QGLikelihoodCalculator.cc.

                                                                                                                 {
  if (!isValidRange(pt, rho, eta, QGLParamsColl->qgValidRange))
    return -1;
 
  float Q = 1., G = 1.;
  for (unsigned int varIndex = 0; varIndex < vars.size(); ++varIndex) {
    auto quarkEntry = findEntry(QGLParamsColl->data, eta, pt, rho, 0, varIndex);
    auto gluonEntry = findEntry(QGLParamsColl->data, eta, pt, rho, 1, varIndex);
    if (!quarkEntry || !gluonEntry)
      return -2;
 
    int binQ = quarkEntry->histogram.findBin(vars[varIndex]);
    float Qi = quarkEntry->histogram.binContent(binQ);
 
    int binG = gluonEntry->histogram.findBin(vars[varIndex]);
    float Gi = gluonEntry->histogram.binContent(binG);
 
    Q *= Qi;
    G *= Gi;
  }
 
  if (Q == 0)
    return 0;
  return Q / (Q + G);
}

References QGLikelihoodObject::data, PVValHelper::eta, findEntry(), callgraph::G, isValidRange(), DiDispStaMuonMonitor_cfi::pt, class-composition::Q, and QGLikelihoodObject::qgValidRange.

Referenced by QGTagger::produce().

findEntry()

const QGLikelihoodObject::Entry * QGLikelihoodCalculator::findEntry ( std::vector< QGLikelihoodObject::Entry > const &  data, float  eta, float  pt, float  rho, int  qgIndex, int  varIndex  ) const

private

Find matching entry in vector for a given eta, pt, rho, qgIndex and varIndex.

Definition at line 36 of file QGLikelihoodCalculator.cc.

                                                                                       {
  QGLikelihoodParameters myParameters;
  myParameters.Rho = rho;
  myParameters.Pt = pt;
  myParameters.Eta = fabs(eta);
  myParameters.QGIndex = qgIndex;
  myParameters.VarIndex = varIndex;
 
  auto myDataObject = data.begin();
  while (!(myParameters == myDataObject->category)) {
    ++myDataObject;
    if (myDataObject == data.end()) {
      edm::LogWarning("QGLCategoryNotFound")
          << "Jet passed qgValidRange criteria, but no category found with rho=" << rho << ", pt=" << pt
          << ", eta=" << eta << "\nPlease contact cms-qg-workinggroup@cern.ch" << std::endl;
      return nullptr;
    }
  }
  return &*myDataObject;
}

References data, PVValHelper::eta, hydjet2DefaultParameters_cff::myParameters, and DiDispStaMuonMonitor_cfi::pt.

Referenced by computeQGLikelihood().

isValidRange()

bool QGLikelihoodCalculator::isValidRange ( float  pt, float  rho, float  eta, const QGLikelihoodCategory &  qgValidRange  ) const

private

Check the valid range of this qg tagger.

Definition at line 63 of file QGLikelihoodCalculator.cc.

                                                                                          {
  if (pt < qgValidRange.PtMin)
    return false;
  if (pt > qgValidRange.PtMax)
    return false;
  if (rho < qgValidRange.RhoMin)
    return false;
  if (rho > qgValidRange.RhoMax)
    return false;
  if (fabs(eta) < qgValidRange.EtaMin)
    return false;
  if (fabs(eta) > qgValidRange.EtaMax)
    return false;
  return true;
}

References PVValHelper::eta, QGLikelihoodCategory::EtaMax, QGLikelihoodCategory::EtaMin, DiDispStaMuonMonitor_cfi::pt, QGLikelihoodCategory::PtMax, QGLikelihoodCategory::PtMin, QGLikelihoodCategory::RhoMax, and QGLikelihoodCategory::RhoMin.

Referenced by computeQGLikelihood().

smearingFunction()

float QGLikelihoodCalculator::smearingFunction ( float  x0, float  a, float  b, float  min, float  max  ) const

private

Return the smeared qgLikelihood value, given input x0 and parameters a, b, min and max.

Definition at line 83 of file QGLikelihoodCalculator.cc.

                                                                                                     {
  float x = (x0 - min) / (max - min);
  if (x < 0.)
    x = 0.;
  if (x > 1.)
    x = 1.;
 
  float x1 = tanh(a * atanh(2. * x - 1.) + b) / 2. + .5;
  if (x <= 0.)
    x1 = 0.;
  if (x >= 1.)
    x1 = 1.;
 
  return x1 * (max - min) + min;
}

References a, b, SiStripPI::max, min(), and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by systematicSmearing().

systematicSmearing()

float QGLikelihoodCalculator::systematicSmearing	(	edm::ESHandle< QGLikelihoodSystematicsObject > &	QGLParamsColl,
		float	pt,
		float	eta,
		float	rho,
		float	qgValue,
		int	qgIndex
	)		const

Definition at line 100 of file QGLikelihoodCalculator.cc.

                                                                    {
  if (qgValue < 0 || qgValue > 1)
    return -1.;
 
  QGLikelihoodParameters myParameters;
  myParameters.Rho = rho;
  myParameters.Pt = pt;
  myParameters.Eta = fabs(eta);
  myParameters.QGIndex = qgIndex;
  myParameters.VarIndex = -1;
 
  auto myDataObject = QGLSystematicsColl->data.begin();
  while (!(myParameters == myDataObject->systCategory)) {
    ++myDataObject;
    if (myDataObject == QGLSystematicsColl->data.end())
      return -1;  //Smearing not available in the whole qgValidRange: do not throw warnings or errors
  }
  return smearingFunction(qgValue, myDataObject->a, myDataObject->b, myDataObject->lmin, myDataObject->lmax);
}

References QGLikelihoodSystematicsObject::data, PVValHelper::eta, hydjet2DefaultParameters_cff::myParameters, DiDispStaMuonMonitor_cfi::pt, and smearingFunction().

Referenced by QGTagger::produce().