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

Private Member Functions
const QGLikelihoodObject::Entry *	findEntry (std::vector< QGLikelihoodObject::Entry > const &data, float eta, float pt, float rho, int qgIndex, int varIndex)
	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)
	Check the valid range of this qg tagger. More...
float	smearingFunction (float x0, float a, float b, float min, float max)
	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.

{};

Member Function Documentation

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

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

Definition at line 8 of file QGLikelihoodCalculator.cc.

References delta, geometryDiff::epsilon, findEntry(), callgraph::G, and isValidRange().

Referenced by QGTagger::produce().

                                                                                                                                                      {
  if(!isValidRange(pt, rho, eta, QGLParamsColl->qgValidRange)) return -1;

  float Q=1., G=1.;
  for(unsigned int varIndex = 0; varIndex < vars.size(); ++varIndex){

    auto qgEntry = findEntry(QGLParamsColl->data, eta, pt, rho, 0, varIndex); 
    if(!qgEntry) return -1; 
    float Qi = qgEntry->histogram.binContent(qgEntry->histogram.findBin(vars[varIndex]));
    float mQ = qgEntry->mean;

    qgEntry = findEntry(QGLParamsColl->data, eta, pt, rho, 1, varIndex); 
    if(!qgEntry) return -1;
    float Gi = qgEntry->histogram.binContent(qgEntry->histogram.findBin(vars[varIndex]));
    float mG = qgEntry->mean;

    float epsilon=0;
    float delta=0.000001;
    if(Qi <= epsilon && Gi <= epsilon){
      if(mQ>mG){
    if(vars[varIndex] > mQ){ Qi = 1-delta; Gi = delta;}
    else if(vars[varIndex] < mG){ Qi = delta; Gi = 1-delta;}
      }
      else if(mQ<mG){
    if(vars[varIndex]<mQ) { Qi = 1-delta; Gi = delta;}
    else if(vars[varIndex]>mG){Qi = delta;Gi = 1-delta;}
      }
    } 
    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  )

private

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

Definition at line 46 of file QGLikelihoodCalculator.cc.

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

Referenced by computeQGLikelihood().

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

private

Check the valid range of this qg tagger.

Definition at line 68 of file QGLikelihoodCalculator.cc.

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

Referenced by computeQGLikelihood().

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

private

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

Definition at line 80 of file QGLikelihoodCalculator.cc.

References min(), and x.

Referenced by systematicSmearing().

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

Definition at line 93 of file QGLikelihoodCalculator.cc.

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

Referenced by QGTagger::produce().

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