HTDQM.cc

Go to the documentation of this file.

#include "DQMOffline/Trigger/plugins/HTDQM.h"

HTDQM::HTDQM()= default;

HTDQM::~HTDQM()= default;

void HTDQM::initialise(const edm::ParameterSet& iConfig ){

  ht_variable_binning_ = iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("htBinning");
  met_variable_binning_ = iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("metBinning2");
  ht_binning_ = getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("htPSet"));
  ls_binning_ = getHistoLSPSet (iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("htlsPSet"));
  
}

void HTDQM::bookHistograms(DQMStore::IBooker     & ibooker) 
{  
  
  std::string histname, histtitle;

  histname = "ht_variable"; histtitle = "PFHT";
  bookME(ibooker,htME_variableBinning_,histname,histtitle,ht_variable_binning_);
  setMETitle(htME_variableBinning_,"PF HT [GeV]","events / [GeV]");

  histname = "htVsMET"; histtitle = "PFHT vs PFMET";
  bookME(ibooker,htVsMET_,histname,histtitle,met_variable_binning_,ht_variable_binning_);
  setMETitle(htVsMET_,"PF MET [GeV]","PF HT [GeV]");

  histname = "htVsLS"; histtitle = "PFHT vs LS";
  bookME(ibooker,htVsLS_,histname,histtitle,ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,ht_binning_.xmin, ht_binning_.xmax);
  setMETitle(htVsLS_,"LS","PF HT [GeV]");

}

void HTDQM::fillHistograms(const std::vector<reco::PFJet> & htjets,
               const double & met,
               const int & ls,
               const bool passCond){

  // filling histograms (denominator)  
  double htSum = 0;
  for (auto const & htjet : htjets){
    htSum += htjet.pt(); 
  }

  htME_variableBinning_.denominator -> Fill(htSum);

  htVsMET_.denominator -> Fill(met, htSum);
  htVsLS_.denominator -> Fill(ls, htSum);
  
  // applying selection for numerator
  if (passCond){
    // filling histograms (num_genTriggerEventFlag_)  
    htME_variableBinning_.numerator -> Fill(htSum);
    htVsMET_.numerator -> Fill(met, htSum);
    htVsLS_.numerator -> Fill(ls, htSum);
  }

}

void HTDQM::fillHtDescription(edm::ParameterSetDescription & histoPSet){

  edm::ParameterSetDescription htPSet;
  fillHistoPSetDescription(htPSet);
  histoPSet.add<edm::ParameterSetDescription>("htPSet", htPSet);

  std::vector<double> bins = {0.,50.,100.,150.,200.,250.,300.,350.,400.,450.,500.,550.,600.,650.,700.,750.,800.,900.,1000.,1200.,1500.,2000.};
  histoPSet.add<std::vector<double> >("htBinning", bins);

  std::vector<double> metbins = {0.,20.,40.,60.,80.,100.,120.,140.,160.,180.,200.,220.,240.,260.,280.,300.,320.,340.,360.,380.,400.,450.,500.,1000.};
  histoPSet.add<std::vector<double> >("metBinning2", metbins);

  edm::ParameterSetDescription lsPSet;
  fillHistoLSPSetDescription(lsPSet);
  histoPSet.add<edm::ParameterSetDescription>("htlsPSet", lsPSet);

}