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 ( )

inline

Definition at line 16 of file QGLikelihoodCalculator.h.

{};

QGLikelihoodCalculator::~QGLikelihoodCalculator ( )

inline

Definition at line 17 of file QGLikelihoodCalculator.h.

References a, b, computeQGLikelihood(), data, PVValHelper::eta, findEntry(), isValidRange(), SiStripPI::max, min(), EnergyCorrector::pt, smearingFunction(), and systematicSmearing().

{};

Member Function Documentation

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.

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

Referenced by QGTagger::produce(), and ~QGLikelihoodCalculator().

                                                                                                                                                            {
  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);
}

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 34 of file QGLikelihoodCalculator.cc.

References QGLikelihoodParameters::Eta, QGLikelihoodParameters::Pt, EnergyCorrector::pt, QGLikelihoodParameters::QGIndex, QGLikelihoodParameters::Rho, and QGLikelihoodParameters::VarIndex.

Referenced by computeQGLikelihood(), and ~QGLikelihoodCalculator().

                                                                                                                                                                               {
  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;
}

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 56 of file QGLikelihoodCalculator.cc.

References QGLikelihoodCategory::EtaMax, QGLikelihoodCategory::EtaMin, QGLikelihoodCategory::PtMax, QGLikelihoodCategory::PtMin, QGLikelihoodCategory::RhoMax, and QGLikelihoodCategory::RhoMin.

Referenced by computeQGLikelihood(), and ~QGLikelihoodCalculator().

                                                                                                                       {
  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;
}

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 68 of file QGLikelihoodCalculator.cc.

References min(), and globals_cff::x1.

Referenced by systematicSmearing(), and ~QGLikelihoodCalculator().

                                                                                                  {
  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;
}

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

Definition at line 81 of file QGLikelihoodCalculator.cc.

References QGLikelihoodSystematicsObject::data, QGLikelihoodParameters::Eta, QGLikelihoodParameters::Pt, EnergyCorrector::pt, QGLikelihoodParameters::QGIndex, QGLikelihoodParameters::Rho, smearingFunction(), and QGLikelihoodParameters::VarIndex.

Referenced by QGTagger::produce(), and ~QGLikelihoodCalculator().

                                                                                                                                                                                {
  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);
}