QGLikelihoodDBReader.cc

Go to the documentation of this file.

#include <memory>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "CondFormats/JetMETObjects/interface/QGLikelihoodObject.h"
#include "JetMETCorrections/Objects/interface/JetCorrectionsRecord.h"
#include "CondFormats/DataRecord/interface/QGLikelihoodRcd.h"

class QGLikelihoodDBReader : public edm::EDAnalyzer{
public:
  explicit QGLikelihoodDBReader(const edm::ParameterSet&);
  ~QGLikelihoodDBReader() override{};
  
private:
  void beginJob() override{};
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void endJob() override{};
 
  std::string mPayloadName;
  bool mCreateTextFile,mPrintScreen;
};


QGLikelihoodDBReader::QGLikelihoodDBReader(const edm::ParameterSet& iConfig){
  mPayloadName    = iConfig.getUntrackedParameter<std::string>("payloadName");
  mPrintScreen    = iConfig.getUntrackedParameter<bool>("printScreen");
  mCreateTextFile = iConfig.getUntrackedParameter<bool>("createTextFile");
}


void QGLikelihoodDBReader::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){
  edm::LogInfo("UserOutput") << "Inspecting QGLikelihood payload with label:" << mPayloadName << std::endl;
  edm::ESHandle<QGLikelihoodObject> QGLParamsColl;
  QGLikelihoodRcd const& rcdhandle = iSetup.get<QGLikelihoodRcd>();
  rcdhandle.get(mPayloadName, QGLParamsColl);

  edm::LogInfo("UserOutput") << "Ranges in which the QGTagger could be applied:"
                 << "  pt: " << QGLParamsColl->qgValidRange.PtMin << " --> " << QGLParamsColl->qgValidRange.PtMax
                 << ", eta: " << QGLParamsColl->qgValidRange.EtaMin << " --> " << QGLParamsColl->qgValidRange.EtaMax
                 << ", rho: " << QGLParamsColl->qgValidRange.RhoMin << " --> " << QGLParamsColl->qgValidRange.RhoMax << std::endl;

  std::vector<QGLikelihoodObject::Entry> const& data = QGLParamsColl->data;
  edm::LogInfo("UserOutput") <<  "There are " << data.size() << " entries (categories with associated PDF):" << std::endl;
  for(auto idata = data.begin(); idata != data.end(); ++idata){    
    int varIndex = idata->category.VarIndex;
    int qgBin = idata->category.QGIndex;
    double etaMin = idata->category.EtaMin;
    double etaMax = idata->category.EtaMax;
    double rhoMin = idata->category.RhoMin;
    double rhoMax = idata->category.RhoMax;
    double ptMin = idata->category.PtMin;
    double ptMax = idata->category.PtMax;

    char buff[1000];
    sprintf(buff, "var=%1d, qg=%1d, ptMin=%8.2f, ptMax=%8.2f, etaMin=%3.1f, etaMax=%3.1f, rhoMin=%6.2f, rhoMax=%6.2f", varIndex, qgBin, ptMin, ptMax, etaMin, etaMax, rhoMin, rhoMax);
    edm::LogVerbatim("UserOutput") << buff << std::endl;
  }
}

DEFINE_FWK_MODULE(QGLikelihoodDBReader);