HTMonitor.h

Go to the documentation of this file.

#ifndef DQMOffline_Trigger_HTMonitor_h
#define DQMOffline_Trigger_HTMonitor_h

#include <string>
#include <vector>
#include <map>

#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>

#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/Registry.h"

#include "CommonTools/Utils/interface/StringCutObjectSelector.h"

//DataFormats
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"

#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/JetReco/interface/JetCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"

#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"

class GenericTriggerEventFlag;
enum quant { HT, MJJ, SOFTDROP };

//
// class declaration
//

class HTMonitor : public DQMEDAnalyzer {
public:
  struct MEHTbinning {
    unsigned nbins;
    double xmin;
    double xmax;
  };

  struct HTME {
    MonitorElement* numerator = nullptr;
    MonitorElement* denominator = nullptr;
  };

  HTMonitor(const edm::ParameterSet&);
  ~HTMonitor() override;
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
  static void fillHistoPSetDescription(edm::ParameterSetDescription& pset);
  static void fillHistoLSPSetDescription(edm::ParameterSetDescription& pset);

protected:
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
  void bookME(DQMStore::IBooker&,
              HTME& me,
              const std::string& histname,
              const std::string& histtitle,
              int nbins,
              double xmin,
              double xmax);
  void bookME(DQMStore::IBooker&,
              HTME& me,
              const std::string& histname,
              const std::string& histtitle,
              const std::vector<double>& binningX);
  void bookME(DQMStore::IBooker&,
              HTME& me,
              const std::string& histname,
              const std::string& histtitle,
              int nbinsX,
              double xmin,
              double xmax,
              double ymin,
              double ymax);
  void bookME(DQMStore::IBooker&,
              HTME& me,
              const std::string& histname,
              const std::string& histtitle,
              int nbinsX,
              double xmin,
              double xmax,
              int nbinsY,
              double ymin,
              double ymax);
  void bookME(DQMStore::IBooker&,
              HTME& me,
              const std::string& histname,
              const std::string& histtitle,
              const std::vector<double>& binningX,
              const std::vector<double>& binningY);
  void setHTitle(HTME& me, const std::string& titleX, const std::string& titleY);

  void analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) override;

private:
  static MEHTbinning getHistoPSet(const edm::ParameterSet& pset);
  static MEHTbinning getHistoLSPSet(const edm::ParameterSet& pset);

  std::string folderName_;
  std::string histoSuffix_;
  quant quantity_;

  edm::InputTag metInputTag_;
  edm::InputTag jetInputTag_;
  edm::InputTag eleInputTag_;
  edm::InputTag muoInputTag_;
  edm::InputTag vtxInputTag_;

  edm::EDGetTokenT<reco::PFMETCollection> metToken_;
  edm::EDGetTokenT<reco::JetView> jetToken_;
  edm::EDGetTokenT<reco::GsfElectronCollection> eleToken_;
  edm::EDGetTokenT<reco::MuonCollection> muoToken_;
  edm::EDGetTokenT<reco::VertexCollection> vtxToken_;

  std::vector<double> ht_variable_binning_;
  MEHTbinning ht_binning_;
  MEHTbinning ls_binning_;

  HTME qME_variableBinning_;
  HTME htVsLS_;
  HTME deltaphimetj1ME_;
  HTME deltaphij1j2ME_;

  std::unique_ptr<GenericTriggerEventFlag> num_genTriggerEventFlag_;
  std::unique_ptr<GenericTriggerEventFlag> den_genTriggerEventFlag_;

  StringCutObjectSelector<reco::MET, true> metSelection_;
  StringCutObjectSelector<reco::Jet, true> jetSelection_;
  StringCutObjectSelector<reco::GsfElectron, true> eleSelection_;
  StringCutObjectSelector<reco::Muon, true> muoSelection_;
  StringCutObjectSelector<reco::Jet, true> jetSelection_HT_;
  unsigned njets_;
  unsigned nelectrons_;
  unsigned nmuons_;
  double dEtaCut_;

  static constexpr double MAXedge_PHI = 3.2;
  static constexpr int Nbin_PHI = 64;
  static constexpr MEHTbinning phi_binning_{Nbin_PHI, -MAXedge_PHI, MAXedge_PHI};

  std::vector<bool> warningPrinted4token_;
};

#endif  // HTMONITOR_H