METAnalyzer.h

Go to the documentation of this file.

#ifndef METAnalyzer_H
#define METAnalyzer_H

#include <memory>
#include <fstream>
#include "TMath.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/EDAnalyzer.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 "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
//
#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/JPTJet.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"

#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"

#include "RecoMET/METAlgorithms/interface/HcalNoiseRBXArray.h"
#include "DataFormats/METReco/interface/BeamHaloSummary.h"
#include "RecoJets/JetProducers/interface/JetIDHelper.h"

#include "DataFormats/MuonReco/interface/MuonMETCorrectionData.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonMETCorrectionData.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlock.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementTrack.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementCluster.h"
#include "DataFormats/ParticleFlowReco/interface/PFCluster.h"

#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/METReco/interface/METCollection.h"
#include "DataFormats/PatCandidates/interface/MET.h"
#include "DataFormats/PatCandidates/interface/Jet.h"

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

#include "DQMOffline/JetMET/interface/JetMETDQMDCSFilter.h"
#include "PhysicsTools/SelectorUtils/interface/JetIDSelectionFunctor.h"
#include "PhysicsTools/SelectorUtils/interface/PFJetIDSelectionFunctor.h"
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include "DataFormats/MuonReco/interface/Muon.h"

#include "CondFormats/L1TObjects/interface/L1GtTriggerMenuFwd.h"
#include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "JetMETCorrections/JetCorrector/interface/JetCorrector.h"
#include "DataFormats/PatCandidates/interface/PackedCandidate.h"

#include <map>
#include <string>

class METAnalyzer : public DQMOneEDAnalyzer<> {
public:
  METAnalyzer(const edm::ParameterSet&);

  ~METAnalyzer() override;

  //  void beginJob(void);
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;

  // Book MonitorElements
  //void bookMESet(std::string);
  //void bookMonitorElement(std::string, bool);

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

  void dqmBeginRun(const edm::Run&, const edm::EventSetup&) override;

  void dqmEndRun(const edm::Run& iRun, const edm::EventSetup& iSetup) override;
  //  void dqmEndRun(const edm::Run& iRun, const edm::EventSetup& iSetup);
  // Fill MonitorElements
  void fillMESet(const edm::Event&,
                 std::string,
                 const reco::MET&,
                 const pat::MET*,
                 const reco::PFMET*,
                 const reco::CaloMET*,
                 const reco::Candidate::PolarLorentzVector&,
                 std::map<std::string, MonitorElement*>&,
                 std::vector<bool>,
                 std::vector<bool>);
  void fillMonitorElement(const edm::Event&,
                          std::string,
                          std::string,
                          const reco::MET&,
                          const pat::MET*,
                          const reco::PFMET*,
                          const reco::CaloMET*,
                          const reco::Candidate::PolarLorentzVector&,
                          std::map<std::string, MonitorElement*>&,
                          bool,
                          bool,
                          std::vector<bool>,
                          std::vector<bool>);
  void makeRatePlot(std::string, double);

  //  bool selectHighPtJetEvent(const edm::Event&);
  //  bool selectLowPtJetEvent(const edm::Event&);
  //  bool selectWElectronEvent(const edm::Event&);
  //  bool selectWMuonEvent(const edm::Event&);

private:
  // Book MonitorElements
  void bookMESet(std::string, DQMStore::IBooker&, std::map<std::string, MonitorElement*>&);
  // Book MonitorElements
  void bookMonitorElement(std::string, DQMStore::IBooker&, std::map<std::string, MonitorElement*>&, bool, bool, bool);

  // ----------member data ---------------------------
  edm::ParameterSet parameters;
  // Switch for verbosity
  int verbose_;

  std::string MetType_;
  std::string FolderName_;

  edm::InputTag metCollectionLabel_;
  edm::InputTag hcalNoiseRBXCollectionTag_;
  edm::InputTag jetCollectionLabel_;
  edm::InputTag vertexTag_;
  edm::InputTag gtTag_;

  edm::EDGetTokenT<std::vector<reco::Vertex>> vertexToken_;
  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> gtToken_;
  edm::EDGetTokenT<reco::CaloJetCollection> caloJetsToken_;
  edm::EDGetTokenT<reco::PFJetCollection> pfJetsToken_;
  edm::EDGetTokenT<pat::JetCollection> patJetsToken_;
  edm::EDGetTokenT<reco::MuonCollection> MuonsToken_;

  edm::InputTag METFilterMiniAODLabel_;
  edm::EDGetTokenT<edm::TriggerResults> METFilterMiniAODToken_;
  edm::InputTag METFilterMiniAODLabel2_;  //needed for RECO and reRECO differntiation
  edm::EDGetTokenT<edm::TriggerResults> METFilterMiniAODToken2_;

  std::vector<int> miniaodFilterIndex_;
  int miniaodfilterdec;  //if RECO set to 0, if reRECO set to 1, else to -1

  edm::InputTag hbheNoiseFilterResultTag_;
  edm::EDGetTokenT<bool> hbheNoiseFilterResultToken_;
  edm::InputTag hbheNoiseIsoFilterResultTag_;
  edm::EDGetTokenT<bool> hbheIsoNoiseFilterResultToken_;
  edm::InputTag CSCHaloResultTag_;
  edm::EDGetTokenT<bool> CSCHaloResultToken_;
  edm::InputTag CSCHalo2015ResultTag_;
  edm::EDGetTokenT<bool> CSCHalo2015ResultToken_;
  edm::InputTag EcalDeadCellTriggerTag_;
  edm::EDGetTokenT<bool> EcalDeadCellTriggerToken_;
  edm::InputTag eeBadScFilterTag_;
  edm::EDGetTokenT<bool> eeBadScFilterToken_;
  edm::InputTag EcalDeadCellBoundaryTag_;
  edm::EDGetTokenT<bool> EcalDeadCellBoundaryToken_;
  edm::InputTag HcalStripHaloTag_;
  edm::EDGetTokenT<bool> HcalStripHaloToken_;

  edm::EDGetTokenT<pat::METCollection> patMetToken_;
  edm::EDGetTokenT<reco::PFMETCollection> pfMetToken_;
  edm::EDGetTokenT<reco::CaloMETCollection> caloMetToken_;

  edm::InputTag inputJetIDValueMap;
  edm::EDGetTokenT<edm::ValueMap<reco::JetID>> jetID_ValueMapToken_;

  edm::ESGetToken<L1GtTriggerMenu, L1GtTriggerMenuRcd> l1gtTrigMenuToken_;

  JetIDSelectionFunctor jetIDFunctorLoose;
  PFJetIDSelectionFunctor pfjetIDFunctorLoose;

  std::string m_l1algoname_;
  int m_bitAlgTechTrig_;

  double ptThreshold_;

  HLTConfigProvider hltConfig_;
  HLTConfigProvider FilterhltConfig_;
  edm::InputTag triggerResultsLabel_;
  edm::EDGetTokenT<edm::TriggerResults> triggerResultsToken_;

  // list of Jet or MB HLT triggers
  //  std::vector<std::string > HLTPathsJetMBByName_;
  std::vector<std::string> allTriggerNames_;
  std::vector<int> allTriggerDecisions_;

  std::string HBHENoiseStringMiniAOD;
  std::string HBHEIsoNoiseStringMiniAOD;

  edm::EDGetTokenT<reco::JetCorrector> jetCorrectorToken_;

  edm::VParameterSet triggerSelectedSubFolders_;
  std::vector<GenericTriggerEventFlag*> triggerFolderEventFlag_;
  std::vector<std::vector<std::string>> triggerFolderExpr_;
  std::vector<std::string> triggerFolderLabels_;
  std::vector<int> triggerFolderDecisions_;
  //  std::vector<MonitorElement* >           triggerFolderME_;

  //  GenericTriggerEventFlag * highPtJetEventFlag_;
  //  GenericTriggerEventFlag * lowPtJetEventFlag_;
  //  GenericTriggerEventFlag * minBiasEventFlag_;
  //  GenericTriggerEventFlag * highMETEventFlag_;
  //  GenericTriggerEventFlag * eleEventFlag_;
  //  GenericTriggerEventFlag * muonEventFlag_;
  //
  //  std::vector<std::string> highPtJetExpr_;
  //  std::vector<std::string> lowPtJetExpr_;
  //  std::vector<std::string> highMETExpr_;
  //  //  std::vector<std::string> lowMETExpr_;
  //  std::vector<std::string> muonExpr_;
  //  std::vector<std::string> elecExpr_;
  //  std::vector<std::string> minbiasExpr_;
  //  MonitorElement* hTriggerName_HighPtJet;
  //  MonitorElement* hTriggerName_LowPtJet;
  //  MonitorElement* hTriggerName_MinBias;
  //  MonitorElement* hTriggerName_HighMET;
  //  //  MonitorElement* hTriggerName_LowMET;
  //  MonitorElement* hTriggerName_Ele;
  //  MonitorElement* hTriggerName_Muon;
  MonitorElement* hMETRate;

  edm::ParameterSet cleaningParameters_;
  std::vector<edm::ParameterSet> diagnosticsParameters_;

  std::string hltPhysDec_;

  int nbinsPV_;
  double nPVMin_;
  double nPVMax_;

  int LSBegin_;
  int LSEnd_;

  bool bypassAllPVChecks_;
  bool bypassAllDCSChecks_;
  bool runcosmics_;
  bool onlyCleaned_;

  //  int trigJetMB_;
  //  int trigHighPtJet_;
  //  int trigLowPtJet_;
  //  int trigMinBias_;
  //  int trigHighMET_;
  //  int trigEle_;
  //  int trigMuon_;
  //  int trigPhysDec_;

  //  double highPtJetThreshold_;
  //  double lowPtJetThreshold_;
  //  double highMETThreshold_;

  int numPV_;
  // Et threshold for MET plots
  //  double etThreshold_;

  // HF calibration factor (in 31X applied by TcProducer)
  //delete altogether not used anymore
  double hfCalibFactor_;  //

  // DCS filter
  JetMETDQMDCSFilter* DCSFilter_;

  std::vector<std::string> folderNames_;
  //
  math::XYZPoint beamSpot_;

  //trigger histos
  // lines commented out have been removed to improve the bin usage of JetMET DQM

  //for all MET types
  bool hTriggerLabelsIsSet_;
  //only in for PF
  //  MonitorElement* meTriggerName_PhysDec;

  MonitorElement* lumisecME;
  MonitorElement* hTrigger;
  //MonitorElement* hNevents;
  MonitorElement* hMEx;
  MonitorElement* hMEy;
  //MonitorElement* hEz;
  MonitorElement* hMETSig;
  MonitorElement* hMET;
  MonitorElement* hMET_2;

  MonitorElement* hMET_HBHENoiseFilter;
  MonitorElement* hMET_CSCTightHaloFilter;
  MonitorElement* hMET_eeBadScFilter;
  MonitorElement* hMET_EcalDeadCellTriggerFilter;
  MonitorElement* hMET_EcalDeadCellBoundaryFilter;
  MonitorElement* hMET_HBHEIsoNoiseFilter;
  MonitorElement* hMET_CSCTightHalo2015Filter;
  MonitorElement* hMET_HcalStripHaloFilter;

  MonitorElement* hMETPhi;
  MonitorElement* hSumET;

  MonitorElement* hMExLS;
  MonitorElement* hMEyLS;

  MonitorElement* hMET_logx;
  MonitorElement* hSumET_logx;

  //CaloMET specific stuff
  MonitorElement* hCaloMETPhi020;
  MonitorElement* hCaloEtFractionHadronic;
  MonitorElement* hCaloEmEtFraction;

  //MonitorElement* hCaloEmEtFraction002;
  //MonitorElement* hCaloEmEtFraction010;
  MonitorElement* hCaloEmEtFraction020;

  MonitorElement* hCaloHadEtInHB;
  MonitorElement* hCaloHadEtInHO;
  MonitorElement* hCaloHadEtInHE;
  MonitorElement* hCaloHadEtInHF;
  MonitorElement* hCaloEmEtInHF;
  MonitorElement* hCaloEmEtInEE;
  MonitorElement* hCaloEmEtInEB;

  MonitorElement* hCaloEmMEx;
  MonitorElement* hCaloEmMEy;
  //MonitorElement* hCaloEmEz;
  MonitorElement* hCaloEmMET;
  MonitorElement* hCaloEmMETPhi;
  //MonitorElement* hCaloEmSumET;

  MonitorElement* hCaloHaMEx;
  MonitorElement* hCaloHaMEy;
  //MonitorElement* hCaloHaEz;
  MonitorElement* hCaloHaMET;
  MonitorElement* hCaloHaMETPhi;
  //MonitorElement* hCaloHaSumET;

  //now PF only things
  MonitorElement* mePhotonEtFraction;
  MonitorElement* mePhotonEt;
  MonitorElement* meNeutralHadronEtFraction;
  MonitorElement* meNeutralHadronEt;
  MonitorElement* meElectronEt;
  MonitorElement* meChargedHadronEtFraction;
  MonitorElement* meChargedHadronEt;
  MonitorElement* meMuonEt;
  MonitorElement* meHFHadronEtFraction;
  MonitorElement* meHFHadronEt;
  MonitorElement* meHFEMEtFraction;
  MonitorElement* meHFEMEt;
  //MEs where we fill if the previous two bunches are empty (25 ns bunch spacing)
  /*
  MonitorElement* mePhotonEtFraction_BXm2BXm1Empty;
  MonitorElement* meNeutralHadronEtFraction_BXm2BXm1Empty;
  MonitorElement* meChargedHadronEtFraction_BXm2BXm1Empty;
  MonitorElement* meMET_BXm2BXm1Empty;
  MonitorElement* meSumET_BXm2BXm1Empty;

  MonitorElement* meMETPhiChargedHadronsBarrel_BXm2BXm1Empty;
  MonitorElement* meMETPhiChargedHadronsEndcapPlus_BXm2BXm1Empty;
  MonitorElement* meMETPhiChargedHadronsEndcapMinus_BXm2BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsBarrel_BXm2BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsEndcapPlus_BXm2BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsEndcapMinus_BXm2BXm1Empty;
  MonitorElement* meMETPhiPhotonsBarrel_BXm2BXm1Empty;
  MonitorElement* meMETPhiPhotonsEndcapPlus_BXm2BXm1Empty;
  MonitorElement* meMETPhiPhotonsEndcapMinus_BXm2BXm1Empty;
  MonitorElement* meMETPhiHFHadronsPlus_BXm2BXm1Empty;
  MonitorElement* meMETPhiHFHadronsMinus_BXm2BXm1Empty;
  MonitorElement* meMETPhiHFEGammasPlus_BXm2BXm1Empty;
  MonitorElement* meMETPhiHFEGammasMinus_BXm2BXm1Empty;
  */
  //MEs where we fill if the previous bunch is empty (25 ns bunch spacing)
  MonitorElement* mePhotonEtFraction_BXm1Empty;
  MonitorElement* meNeutralHadronEtFraction_BXm1Empty;
  MonitorElement* meChargedHadronEtFraction_BXm1Empty;
  MonitorElement* meMET_BXm1Empty;
  MonitorElement* meSumET_BXm1Empty;

  MonitorElement* meMETPhiChargedHadronsBarrel_BXm1Empty;
  MonitorElement* meMETPhiChargedHadronsEndcapPlus_BXm1Empty;
  MonitorElement* meMETPhiChargedHadronsEndcapMinus_BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsBarrel_BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsEndcapPlus_BXm1Empty;
  MonitorElement* meMETPhiNeutralHadronsEndcapMinus_BXm1Empty;
  MonitorElement* meMETPhiPhotonsBarrel_BXm1Empty;
  MonitorElement* meMETPhiPhotonsEndcapPlus_BXm1Empty;
  MonitorElement* meMETPhiPhotonsEndcapMinus_BXm1Empty;
  MonitorElement* meMETPhiHFHadronsPlus_BXm1Empty;
  MonitorElement* meMETPhiHFHadronsMinus_BXm1Empty;
  MonitorElement* meMETPhiHFEGammasPlus_BXm1Empty;
  MonitorElement* meMETPhiHFEGammasMinus_BXm1Empty;

  //MEs where we fill if the previous bunch is filled (25 ns bunch spacing)
  MonitorElement* mePhotonEtFraction_BXm1Filled;
  MonitorElement* meNeutralHadronEtFraction_BXm1Filled;
  MonitorElement* meChargedHadronEtFraction_BXm1Filled;
  MonitorElement* meMET_BXm1Filled;
  MonitorElement* meSumET_BXm1Filled;

  MonitorElement* meMETPhiChargedHadronsBarrel_BXm1Filled;
  MonitorElement* meMETPhiChargedHadronsEndcapPlus_BXm1Filled;
  MonitorElement* meMETPhiChargedHadronsEndcapMinus_BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsBarrel_BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsEndcapPlus_BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsEndcapMinus_BXm1Filled;
  MonitorElement* meMETPhiPhotonsBarrel_BXm1Filled;
  MonitorElement* meMETPhiPhotonsEndcapPlus_BXm1Filled;
  MonitorElement* meMETPhiPhotonsEndcapMinus_BXm1Filled;
  MonitorElement* meMETPhiHFHadronsPlus_BXm1Filled;
  MonitorElement* meMETPhiHFHadronsMinus_BXm1Filled;
  MonitorElement* meMETPhiHFEGammasPlus_BXm1Filled;
  MonitorElement* meMETPhiHFEGammasMinus_BXm1Filled;

  //MEs where we fill if two previous bunches are filled (25 ns bunch spacing)
  /*
  MonitorElement* meChargedHadronEtFraction_BXm2BXm1Filled;
  MonitorElement* mePhotonEtFraction_BXm2BXm1Filled;
  MonitorElement* meNeutralHadronEtFraction_BXm2BXm1Filled;
  MonitorElement* meMET_BXm2BXm1Filled;
  MonitorElement* meSumET_BXm2BXm1Filled;
  */
  MonitorElement* meCHF_Barrel;
  MonitorElement* meCHF_EndcapPlus;
  MonitorElement* meCHF_EndcapMinus;
  MonitorElement* meCHF_Barrel_BXm1Empty;
  MonitorElement* meCHF_EndcapPlus_BXm1Empty;
  MonitorElement* meCHF_EndcapMinus_BXm1Empty;
  //MonitorElement* meCHF_Barrel_BXm2BXm1Empty;
  //MonitorElement* meCHF_EndcapPlus_BXm2BXm1Empty;
  //MonitorElement* meCHF_EndcapMinus_BXm2BXm1Empty;
  MonitorElement* meCHF_Barrel_BXm1Filled;
  MonitorElement* meCHF_EndcapPlus_BXm1Filled;
  MonitorElement* meCHF_EndcapMinus_BXm1Filled;
  //MonitorElement* meCHF_Barrel_BXm2BXm1Filled;
  //MonitorElement* meCHF_EndcapPlus_BXm2BXm1Filled;
  //MonitorElement* meCHF_EndcapMinus_BXm2BXm1Filled;

  MonitorElement* meNHF_Barrel;
  MonitorElement* meNHF_EndcapPlus;
  MonitorElement* meNHF_EndcapMinus;
  MonitorElement* meNHF_Barrel_BXm1Empty;
  MonitorElement* meNHF_EndcapPlus_BXm1Empty;
  MonitorElement* meNHF_EndcapMinus_BXm1Empty;
  //MonitorElement* meNHF_Barrel_BXm2BXm1Empty;
  //MonitorElement* meNHF_EndcapPlus_BXm2BXm1Empty;
  //MonitorElement* meNHF_EndcapMinus_BXm2BXm1Empty;
  MonitorElement* meNHF_Barrel_BXm1Filled;
  MonitorElement* meNHF_EndcapPlus_BXm1Filled;
  MonitorElement* meNHF_EndcapMinus_BXm1Filled;
  //MonitorElement* meNHF_Barrel_BXm2BXm1Filled;
  //MonitorElement* meNHF_EndcapPlus_BXm2BXm1Filled;
  //MonitorElement* meNHF_EndcapMinus_BXm2BXm1Filled;

  MonitorElement* mePhF_Barrel;
  MonitorElement* mePhF_EndcapPlus;
  MonitorElement* mePhF_EndcapMinus;
  MonitorElement* mePhF_Barrel_BXm1Empty;
  MonitorElement* mePhF_EndcapPlus_BXm1Empty;
  MonitorElement* mePhF_EndcapMinus_BXm1Empty;
  //MonitorElement* mePhF_Barrel_BXm2BXm1Empty;
  //MonitorElement* mePhF_EndcapPlus_BXm2BXm1Empty;
  //MonitorElement* mePhF_EndcapMinus_BXm2BXm1Empty;
  MonitorElement* mePhF_Barrel_BXm1Filled;
  MonitorElement* mePhF_EndcapPlus_BXm1Filled;
  MonitorElement* mePhF_EndcapMinus_BXm1Filled;
  //MonitorElement* mePhF_Barrel_BXm2BXm1Filled;
  //MonitorElement* mePhF_EndcapPlus_BXm2BXm1Filled;
  //MonitorElement* mePhF_EndcapMinus_BXm2BXm1Filled;

  MonitorElement* meHFHadF_Plus;
  MonitorElement* meHFHadF_Minus;
  MonitorElement* meHFHadF_Plus_BXm1Empty;
  MonitorElement* meHFHadF_Minus_BXm1Empty;
  //MonitorElement* meHFHadF_Plus_BXm2BXm1Empty;
  //MonitorElement* meHFHadF_Minus_BXm2BXm1Empty;
  MonitorElement* meHFHadF_Plus_BXm1Filled;
  MonitorElement* meHFHadF_Minus_BXm1Filled;
  //MonitorElement* meHFHadF_Plus_BXm2BXm1Filled;
  //MonitorElement* meHFHadF_Minus_BXm2BXm1Filled;

  MonitorElement* meHFEMF_Plus;
  MonitorElement* meHFEMF_Minus;
  MonitorElement* meHFEMF_Plus_BXm1Empty;
  MonitorElement* meHFEMF_Minus_BXm1Empty;
  //MonitorElement* meHFEMF_Plus_BXm2BXm1Empty;
  //MonitorElement* meHFEMF_Minus_BXm2BXm1Empty;
  MonitorElement* meHFEMF_Plus_BXm1Filled;
  MonitorElement* meHFEMF_Minus_BXm1Filled;
  //MonitorElement* meHFEMF_Plus_BXm2BXm1Filled;
  //MonitorElement* meHFEMF_Minus_BXm2BXm1Filled;
  /*
  MonitorElement* meMETPhiChargedHadronsBarrel_BXm2BXm1Filled;
  MonitorElement* meMETPhiChargedHadronsEndcapPlus_BXm2BXm1Filled;
  MonitorElement* meMETPhiChargedHadronsEndcapMinus_BXm2BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsBarrel_BXm2BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsEndcapPlus_BXm2BXm1Filled;
  MonitorElement* meMETPhiNeutralHadronsEndcapMinus_BXm2BXm1Filled;
  MonitorElement* meMETPhiPhotonsBarrel_BXm2BXm1Filled;
  MonitorElement* meMETPhiPhotonsEndcapPlus_BXm2BXm1Filled;
  MonitorElement* meMETPhiPhotonsEndcapMinus_BXm2BXm1Filled;
  MonitorElement* meMETPhiHFHadronsPlus_BXm2BXm1Filled;
  MonitorElement* meMETPhiHFHadronsMinus_BXm2BXm1Filled;
  MonitorElement* meMETPhiHFEGammasPlus_BXm2BXm1Filled;
  MonitorElement* meMETPhiHFEGammasMinus_BXm2BXm1Filled;
  */
  double ptMinCand_;

  //PFcandidate maps
  std::vector<MonitorElement*> profilePFCand_x_, profilePFCand_y_;
  std::vector<std::string> profilePFCand_x_name_, profilePFCand_y_name_;
  std::vector<double> etaMinPFCand_, etaMaxPFCand_, MExPFCand_, MEyPFCand_;
  std::vector<int> typePFCand_, countsPFCand_, etaNBinsPFCand_;

  MonitorElement* meMETPhiChargedHadronsBarrel;
  MonitorElement* meMETPhiChargedHadronsEndcapPlus;
  MonitorElement* meMETPhiChargedHadronsEndcapMinus;
  MonitorElement* meMETPhiNeutralHadronsBarrel;
  MonitorElement* meMETPhiNeutralHadronsEndcapPlus;
  MonitorElement* meMETPhiNeutralHadronsEndcapMinus;
  MonitorElement* meMETPhiPhotonsBarrel;
  MonitorElement* meMETPhiPhotonsEndcapPlus;
  MonitorElement* meMETPhiPhotonsEndcapMinus;
  MonitorElement* meMETPhiHFHadronsPlus;
  MonitorElement* meMETPhiHFHadronsMinus;
  MonitorElement* meMETPhiHFEGammasPlus;
  MonitorElement* meMETPhiHFEGammasMinus;

  edm::EDGetTokenT<std::vector<reco::PFCandidate>> pflowToken_;
  edm::EDGetTokenT<std::vector<pat::PackedCandidate>> pflowPackedToken_;

  // NPV profiles -->
  //----------------------------------------------------------------------------
  MonitorElement* meMEx_profile;
  MonitorElement* meMEy_profile;
  MonitorElement* meMET_profile;
  MonitorElement* meSumET_profile;

  MonitorElement* mePhotonEtFraction_profile;
  MonitorElement* mePhotonEt_profile;
  MonitorElement* meNeutralHadronEtFraction_profile;
  MonitorElement* meNeutralHadronEt_profile;
  MonitorElement* meChargedHadronEtFraction_profile;
  MonitorElement* meChargedHadronEt_profile;
  MonitorElement* meHFHadronEtFraction_profile;
  MonitorElement* meHFHadronEt_profile;
  MonitorElement* meHFEMEtFraction_profile;
  MonitorElement* meHFEMEt_profile;

  MonitorElement* meZJets_u_par;
  MonitorElement* meZJets_u_par_ZPt_0_15;
  MonitorElement* meZJets_u_par_ZPt_15_30;
  MonitorElement* meZJets_u_par_ZPt_30_55;
  MonitorElement* meZJets_u_par_ZPt_55_75;
  MonitorElement* meZJets_u_par_ZPt_75_150;
  MonitorElement* meZJets_u_par_ZPt_150_290;
  MonitorElement* meZJets_u_par_ZPt_290;

  MonitorElement* meZJets_u_perp;
  MonitorElement* meZJets_u_perp_ZPt_0_15;
  MonitorElement* meZJets_u_perp_ZPt_15_30;
  MonitorElement* meZJets_u_perp_ZPt_30_55;
  MonitorElement* meZJets_u_perp_ZPt_55_75;
  MonitorElement* meZJets_u_perp_ZPt_75_150;
  MonitorElement* meZJets_u_perp_ZPt_150_290;
  MonitorElement* meZJets_u_perp_ZPt_290;

  std::map<std::string, MonitorElement*> map_dijet_MEs;

  bool isCaloMet_;
  bool isPFMet_;
  bool isMiniAODMet_;

  bool fill_met_high_level_histo;
  bool fillCandidateMap_histos;
};
#endif