QGLikelihoodObject.h

Go to the documentation of this file.

#ifndef QGLikelihoodObject_h
#define QGLikelihoodObject_h

#include "CondFormats/PhysicsToolsObjects/interface/Histogram.h"
#include "CondFormats/Serialization/interface/Serializable.h"
#include <vector>

struct QGLikelihoodCategory {
  float RhoMin, RhoMax, PtMin, PtMax, EtaMin, EtaMax;
  int QGIndex, VarIndex;
  COND_SERIALIZABLE;
};

struct QGLikelihoodParameters {
  float Rho, Pt, Eta;
  int QGIndex, VarIndex;
};

struct QGLikelihoodObject {
  typedef PhysicsTools::Calibration::HistogramF Histogram;

  struct Entry {
    QGLikelihoodCategory category;
    Histogram histogram;
    float mean;
    COND_SERIALIZABLE;
  };

  QGLikelihoodCategory qgValidRange;
  std::vector<Entry> data;
  COND_SERIALIZABLE;
};

struct QGLikelihoodSystematicsObject {
  struct Entry {
    QGLikelihoodCategory systCategory;
    float a, b, lmin, lmax;
    COND_SERIALIZABLE;
  };
  std::vector<Entry> data;
  COND_SERIALIZABLE;
};

inline bool operator==(const QGLikelihoodParameters& lhs, const QGLikelihoodCategory& rhs) {
  if (lhs.QGIndex != rhs.QGIndex)
    return false;
  if (lhs.VarIndex != rhs.VarIndex)
    return false;
  if (lhs.Eta < rhs.EtaMin)
    return false;
  if (lhs.Eta > rhs.EtaMax)
    return false;
  if (lhs.Rho < rhs.RhoMin)
    return false;
  if (lhs.Rho > rhs.RhoMax)
    return false;
  if (lhs.Pt < rhs.PtMin)
    return false;
  if (lhs.Pt > rhs.PtMax)
    return false;
  return true;
}

#endif