JetAnalyzer.h

Go to the documentation of this file.

#ifndef JetAnalyzer_H
#define JetAnalyzer_H

#include <memory>
#include <fstream>
#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/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"

#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
#include "FWCore/Utilities/interface/EDGetToken.h"

#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/JetReco/interface/JPTJet.h"
#include "DataFormats/JetReco/interface/JPTJetCollection.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "JetMETCorrections/JetCorrector/interface/JetCorrector.h"
#include "RecoJets/JetProducers/interface/JetIDHelper.h"
#include "DataFormats/MuonReco/interface/Muon.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 "DQMOffline/JetMET/interface/JetMETDQMDCSFilter.h"

#include "DataFormats/BTauReco/interface/CATopJetTagInfo.h"

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

#include "DataFormats/Scalers/interface/DcsStatus.h"
#include "PhysicsTools/SelectorUtils/interface/JetIDSelectionFunctor.h"
#include "PhysicsTools/SelectorUtils/interface/PFJetIDSelectionFunctor.h"
#include "DataFormats/JetReco/interface/PileupJetIdentifier.h"

#include "CondFormats/L1TObjects/interface/L1GtTriggerMenuFwd.h"
#include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include <map>
#include <string>

//namespace jetAnalysis {
//class TrackPropagatorToCalo;
//class StripSignalOverNoiseCalculator;
//}

class JetAnalyzer : public DQMEDAnalyzer {
public:
  JetAnalyzer(const edm::ParameterSet&);

  ~JetAnalyzer() override;

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

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

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

private:
  // ----------member data ---------------------------
  static bool jetSortingRule(reco::Jet x, reco::Jet y) { return x.pt() > y.pt(); }

  //try to put one collection as start
  edm::InputTag mInputCollection_;
  edm::InputTag theTriggerResultsLabel_;

  std::string jetType_;

  edm::EDGetTokenT<edm::TriggerResults> triggerResultsToken_;
  edm::EDGetTokenT<std::vector<reco::Vertex>> vertexToken_;
  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> gtToken_;
  edm::EDGetTokenT<reco::CaloJetCollection> caloJetsToken_;
  edm::EDGetTokenT<reco::PFJetCollection> pfJetsToken_;

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

  edm::EDGetTokenT<reco::MuonCollection> MuonsToken_;
  edm::EDGetTokenT<pat::JetCollection> patJetsToken_;
  edm::EDGetTokenT<edm::ValueMap<float>> mvaFullPUDiscriminantToken_;
  edm::EDGetTokenT<edm::ValueMap<float>> cutBasedPUDiscriminantToken_;
  edm::EDGetTokenT<edm::ValueMap<int>> cutBasedPUIDToken_;
  edm::EDGetTokenT<edm::ValueMap<int>> mvaPUIDToken_;

  edm::EDGetTokenT<edm::ValueMap<int>> qgMultiplicityToken_;
  edm::EDGetTokenT<edm::ValueMap<float>> qgLikelihoodToken_;
  edm::EDGetTokenT<edm::ValueMap<float>> qgptDToken_;
  edm::EDGetTokenT<edm::ValueMap<float>> qgaxis2Token_;

  //edm::EDGetTokenT<reco::JPTJetCollection>        jptJetsToken_;

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

  //Cleaning parameters
  edm::ParameterSet cleaningParameters_;
  edm::InputTag vertexLabel_;
  edm::InputTag gtLabel_;

  //check later if we need only one set of parameters
  edm::ParameterSet parameters_;

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

  std::string JetIDQuality_;
  std::string JetIDVersion_;
  JetIDSelectionFunctor::Quality_t jetidquality;
  JetIDSelectionFunctor::Version_t jetidversion;
  JetIDSelectionFunctor jetIDFunctor;

  PFJetIDSelectionFunctor::Quality_t pfjetidquality;
  PFJetIDSelectionFunctor::Version_t pfjetidversion;

  PFJetIDSelectionFunctor pfjetIDFunctor;

  std::vector<std::string> folderNames_;

  std::string DirName;

  // Book MonitorElements
  void bookMESetSelection(std::string, DQMStore::IBooker&);
  //void bookMonitorElement(std::string, bool);

  int verbose_;
  //histo binning parameters -> these are PART of ALL analyzers - move it up
  int etaBin_;
  double etaMin_;
  double etaMax_;

  int phiBin_;
  double phiMin_;
  double phiMax_;

  int ptBin_;
  double ptMin_;
  double ptMax_;

  int eBin_;
  double eMin_;
  double eMax_;

  int pBin_;
  double pMin_;
  double pMax_;

  int nbinsPV_;
  double nPVlow_;
  double nPVhigh_;

  //variables which are present both in
  int jetLoPass_;
  int jetHiPass_;
  int leadJetFlag_;
  double ptThreshold_;
  double ptThresholdUnc_;
  double asymmetryThirdJetCut_;
  double balanceThirdJetCut_;

  //
  int fillJIDPassFrac_;
  std::string m_l1algoname_;
  int m_bitAlgTechTrig_;

  //the histos
  MonitorElement* jetME;

  // --- Used for Data Certification - use for Calo, PF and JPT jets
  MonitorElement* mPt;
  MonitorElement* mPt_1;
  MonitorElement* mPt_2;
  MonitorElement* mPt_3;
  MonitorElement* mEta;
  MonitorElement* mPhi;
  MonitorElement* mPt_uncor;
  MonitorElement* mEta_uncor;
  MonitorElement* mPhi_uncor;
  MonitorElement* mConstituents_uncor;

  MonitorElement* mJetEnergyCorr;
  MonitorElement* mJetEnergyCorrVSEta;
  MonitorElement* mJetEnergyCorrVSPt;

  MonitorElement* mConstituents;
  MonitorElement* mHFrac;
  MonitorElement* mEFrac;
  MonitorElement* mPhiVSEta;

  MonitorElement* mPt_Barrel;
  MonitorElement* mPhi_Barrel;
  MonitorElement* mConstituents_Barrel;
  MonitorElement* mHFrac_Barrel;
  MonitorElement* mEFrac_Barrel;

  MonitorElement* mPt_EndCap;
  MonitorElement* mPhi_EndCap;
  MonitorElement* mConstituents_EndCap;
  MonitorElement* mHFrac_EndCap;
  MonitorElement* mEFrac_EndCap;

  MonitorElement* mPt_Forward;
  MonitorElement* mPhi_Forward;
  MonitorElement* mConstituents_Forward;
  MonitorElement* mHFrac_Forward;
  MonitorElement* mEFrac_Forward;

  MonitorElement* mPt_Barrel_Hi;
  MonitorElement* mPhi_Barrel_Hi;
  MonitorElement* mConstituents_Barrel_Hi;
  MonitorElement* mHFrac_Barrel_Hi;

  MonitorElement* mPt_EndCap_Hi;
  MonitorElement* mPhi_EndCap_Hi;
  MonitorElement* mConstituents_EndCap_Hi;
  MonitorElement* mHFrac_EndCap_Hi;

  MonitorElement* mPt_Forward_Hi;
  MonitorElement* mPhi_Forward_Hi;
  MonitorElement* mConstituents_Forward_Hi;
  MonitorElement* mHFrac_Forward_Hi;

  MonitorElement* mNJets;
  MonitorElement* mDPhi;

  // Leading Jet Parameters
  MonitorElement* mEtaFirst;
  MonitorElement* mPhiFirst;
  MonitorElement* mPtFirst;

  // Events passing the jet triggers
  MonitorElement* mPhi_Lo;
  MonitorElement* mPt_Lo;

  MonitorElement* mEta_Hi;
  MonitorElement* mPhi_Hi;
  MonitorElement* mPt_Hi;

  MonitorElement* mLooseJIDPassFractionVSeta;
  MonitorElement* mLooseJIDPassFractionVSpt;
  MonitorElement* mLooseJIDPassFractionVSptNoHF;

  MonitorElement* mLooseMVAPUJIDPassFractionVSeta;
  MonitorElement* mLooseMVAPUJIDPassFractionVSpt;
  MonitorElement* mMediumMVAPUJIDPassFractionVSeta;
  MonitorElement* mMediumMVAPUJIDPassFractionVSpt;
  MonitorElement* mTightMVAPUJIDPassFractionVSeta;
  MonitorElement* mTightMVAPUJIDPassFractionVSpt;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_Barrel;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_EndCap;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_Forward;
  MonitorElement* mMVAPUJIDDiscriminant_mediumPt_Barrel;
  MonitorElement* mMVAPUJIDDiscriminant_mediumPt_EndCap;
  MonitorElement* mMVAPUJIDDiscriminant_mediumPt_Forward;
  MonitorElement* mMVAPUJIDDiscriminant_highPt_Barrel;
  MonitorElement* mMVAPUJIDDiscriminant_highPt_EndCap;
  MonitorElement* mMVAPUJIDDiscriminant_highPt_Forward;

  MonitorElement* mLooseCutPUJIDPassFractionVSeta;
  MonitorElement* mLooseCutPUJIDPassFractionVSpt;
  MonitorElement* mMediumCutPUJIDPassFractionVSeta;
  MonitorElement* mMediumCutPUJIDPassFractionVSpt;
  MonitorElement* mTightCutPUJIDPassFractionVSeta;
  MonitorElement* mTightCutPUJIDPassFractionVSpt;
  MonitorElement* mCutPUJIDDiscriminant_lowPt_Barrel;
  MonitorElement* mCutPUJIDDiscriminant_lowPt_EndCap;
  MonitorElement* mCutPUJIDDiscriminant_lowPt_Forward;
  MonitorElement* mCutPUJIDDiscriminant_mediumPt_Barrel;
  MonitorElement* mCutPUJIDDiscriminant_mediumPt_EndCap;
  MonitorElement* mCutPUJIDDiscriminant_mediumPt_Forward;
  MonitorElement* mCutPUJIDDiscriminant_highPt_Barrel;
  MonitorElement* mCutPUJIDDiscriminant_highPt_EndCap;
  MonitorElement* mCutPUJIDDiscriminant_highPt_Forward;

  //dijet analysis quantities
  MonitorElement* mDijetBalance;
  MonitorElement* mDijetAsymmetry;

  // NPV profiles
  //----------------------------------------------------------------------------
  MonitorElement* mNJets_profile;
  MonitorElement* mPt_profile;
  MonitorElement* mEta_profile;
  MonitorElement* mPhi_profile;
  MonitorElement* mConstituents_profile;
  MonitorElement* mHFrac_profile;
  MonitorElement* mEFrac_profile;

  bool hltInitialized_;
  bool bypassAllPVChecks_;

  HLTConfigProvider hltConfig_;
  std::string processname_;

  //MonitorElement* hltpathME;
  MonitorElement* cleanupME;
  MonitorElement* verticesME;

  GenericTriggerEventFlag* highPtJetEventFlag_;
  GenericTriggerEventFlag* lowPtJetEventFlag_;

  std::vector<std::string> highPtJetExpr_;
  std::vector<std::string> lowPtJetExpr_;

  bool jetCleaningFlag_;
  bool filljetsubstruc_;
  double pt_min_boosted_;

  bool runcosmics_;

  //  bool energycorrected;

  // CaloJet specific
  MonitorElement* mMaxEInEmTowers;
  MonitorElement* mMaxEInHadTowers;
  MonitorElement* mHadEnergyInHO;
  MonitorElement* mHadEnergyInHB;
  MonitorElement* mHadEnergyInHF;
  MonitorElement* mHadEnergyInHE;
  MonitorElement* mEmEnergyInEB;
  MonitorElement* mEmEnergyInEE;
  MonitorElement* mEmEnergyInHF;
  MonitorElement* mN90Hits;
  MonitorElement* mfHPD;
  MonitorElement* mfRBX;
  MonitorElement* mresEMF;
  MonitorElement* mEMF;

  //now define PFJet only flags
  MonitorElement* mCHFrac_lowPt_Barrel;
  MonitorElement* mNHFrac_lowPt_Barrel;
  MonitorElement* mPhFrac_lowPt_Barrel;
  MonitorElement* mCHFrac_mediumPt_Barrel;
  MonitorElement* mNHFrac_mediumPt_Barrel;
  MonitorElement* mPhFrac_mediumPt_Barrel;
  MonitorElement* mCHFrac_highPt_Barrel;
  MonitorElement* mNHFrac_highPt_Barrel;
  MonitorElement* mPhFrac_highPt_Barrel;
  MonitorElement* mCHEn_lowPt_Barrel;
  MonitorElement* mNHEn_lowPt_Barrel;
  MonitorElement* mPhEn_lowPt_Barrel;
  MonitorElement* mElEn_lowPt_Barrel;
  MonitorElement* mMuEn_lowPt_Barrel;
  MonitorElement* mCHEn_mediumPt_Barrel;
  MonitorElement* mNHEn_mediumPt_Barrel;
  MonitorElement* mPhEn_mediumPt_Barrel;
  MonitorElement* mElEn_mediumPt_Barrel;
  MonitorElement* mMuEn_mediumPt_Barrel;
  MonitorElement* mCHEn_highPt_Barrel;
  MonitorElement* mNHEn_highPt_Barrel;
  MonitorElement* mPhEn_highPt_Barrel;
  MonitorElement* mElEn_highPt_Barrel;
  MonitorElement* mMuEn_highPt_Barrel;
  MonitorElement* mChMultiplicity_lowPt_Barrel;
  MonitorElement* mNeutMultiplicity_lowPt_Barrel;
  MonitorElement* mMuMultiplicity_lowPt_Barrel;
  MonitorElement* mChMultiplicity_mediumPt_Barrel;
  MonitorElement* mNeutMultiplicity_mediumPt_Barrel;
  MonitorElement* mMuMultiplicity_mediumPt_Barrel;
  MonitorElement* mChMultiplicity_highPt_Barrel;
  MonitorElement* mNeutMultiplicity_highPt_Barrel;
  MonitorElement* mMuMultiplicity_highPt_Barrel;

  MonitorElement* mCHFracVSpT_Barrel;
  MonitorElement* mNHFracVSpT_Barrel;
  MonitorElement* mPhFracVSpT_Barrel;
  MonitorElement* mCHFracVSpT_EndCap;
  MonitorElement* mNHFracVSpT_EndCap;
  MonitorElement* mPhFracVSpT_EndCap;
  MonitorElement* mHFHFracVSpT_Forward;
  MonitorElement* mHFEFracVSpT_Forward;

  MonitorElement* mCHFracVSeta_lowPt;
  MonitorElement* mNHFracVSeta_lowPt;
  MonitorElement* mPhFracVSeta_lowPt;
  MonitorElement* mCHFracVSeta_mediumPt;
  MonitorElement* mNHFracVSeta_mediumPt;
  MonitorElement* mPhFracVSeta_mediumPt;
  MonitorElement* mCHFracVSeta_highPt;
  MonitorElement* mNHFracVSeta_highPt;
  MonitorElement* mPhFracVSeta_highPt;

  MonitorElement* mCHFrac_lowPt_EndCap;
  MonitorElement* mNHFrac_lowPt_EndCap;
  MonitorElement* mPhFrac_lowPt_EndCap;
  MonitorElement* mCHFrac_mediumPt_EndCap;
  MonitorElement* mNHFrac_mediumPt_EndCap;
  MonitorElement* mPhFrac_mediumPt_EndCap;
  MonitorElement* mCHFrac_highPt_EndCap;
  MonitorElement* mNHFrac_highPt_EndCap;
  MonitorElement* mPhFrac_highPt_EndCap;

  MonitorElement* mCHEn_lowPt_EndCap;
  MonitorElement* mNHEn_lowPt_EndCap;
  MonitorElement* mPhEn_lowPt_EndCap;
  MonitorElement* mElEn_lowPt_EndCap;
  MonitorElement* mMuEn_lowPt_EndCap;
  MonitorElement* mCHEn_mediumPt_EndCap;
  MonitorElement* mNHEn_mediumPt_EndCap;
  MonitorElement* mPhEn_mediumPt_EndCap;
  MonitorElement* mElEn_mediumPt_EndCap;
  MonitorElement* mMuEn_mediumPt_EndCap;
  MonitorElement* mCHEn_highPt_EndCap;
  MonitorElement* mNHEn_highPt_EndCap;
  MonitorElement* mPhEn_highPt_EndCap;
  MonitorElement* mElEn_highPt_EndCap;
  MonitorElement* mMuEn_highPt_EndCap;
  MonitorElement* mMass_lowPt_Barrel;
  MonitorElement* mMass_lowPt_EndCap;
  MonitorElement* mMass_lowPt_Forward;
  MonitorElement* mMass_mediumPt_Barrel;
  MonitorElement* mMass_mediumPt_EndCap;
  MonitorElement* mMass_mediumPt_Forward;
  MonitorElement* mMass_highPt_Barrel;
  MonitorElement* mMass_highPt_EndCap;
  MonitorElement* mMass_highPt_Forward;

  MonitorElement* mChMultiplicity_lowPt_EndCap;
  MonitorElement* mNeutMultiplicity_lowPt_EndCap;
  MonitorElement* mMuMultiplicity_lowPt_EndCap;
  MonitorElement* mChMultiplicity_mediumPt_EndCap;
  MonitorElement* mNeutMultiplicity_mediumPt_EndCap;
  MonitorElement* mMuMultiplicity_mediumPt_EndCap;
  MonitorElement* mChMultiplicity_highPt_EndCap;
  MonitorElement* mNeutMultiplicity_highPt_EndCap;
  MonitorElement* mMuMultiplicity_highPt_EndCap;

  MonitorElement* mHFEFrac_lowPt_Forward;
  MonitorElement* mHFHFrac_lowPt_Forward;
  MonitorElement* mHFEFrac_mediumPt_Forward;
  MonitorElement* mHFHFrac_mediumPt_Forward;
  MonitorElement* mHFEFrac_highPt_Forward;
  MonitorElement* mHFHFrac_highPt_Forward;
  MonitorElement* mHFEEn_lowPt_Forward;
  MonitorElement* mHFHEn_lowPt_Forward;
  MonitorElement* mHFEEn_mediumPt_Forward;
  MonitorElement* mHFHEn_mediumPt_Forward;
  MonitorElement* mHFEEn_highPt_Forward;
  MonitorElement* mHFHEn_highPt_Forward;
  MonitorElement* mNeutMultiplicity_lowPt_Forward;
  MonitorElement* mNeutMultiplicity_mediumPt_Forward;
  MonitorElement* mNeutMultiplicity_highPt_Forward;

  MonitorElement* mChargedHadronEnergy;
  MonitorElement* mNeutralHadronEnergy;
  MonitorElement* mChargedEmEnergy;
  MonitorElement* mChargedMuEnergy;
  MonitorElement* mNeutralEmEnergy;
  MonitorElement* mChargedMultiplicity;
  MonitorElement* mNeutralMultiplicity;
  MonitorElement* mMuonMultiplicity;

  //it is there for ak4PFCHS
  MonitorElement* mAxis2_lowPt_Barrel;
  MonitorElement* mpTD_lowPt_Barrel;
  MonitorElement* mMultiplicityQG_lowPt_Barrel;
  MonitorElement* mqgLikelihood_lowPt_Barrel;
  MonitorElement* mAxis2_mediumPt_Barrel;
  MonitorElement* mpTD_mediumPt_Barrel;
  MonitorElement* mMultiplicityQG_mediumPt_Barrel;
  MonitorElement* mqgLikelihood_mediumPt_Barrel;
  MonitorElement* mAxis2_highPt_Barrel;
  MonitorElement* mpTD_highPt_Barrel;
  MonitorElement* mMultiplicityQG_highPt_Barrel;
  MonitorElement* mqgLikelihood_highPt_Barrel;

  MonitorElement* mAxis2_lowPt_EndCap;
  MonitorElement* mpTD_lowPt_EndCap;
  MonitorElement* mMultiplicityQG_lowPt_EndCap;
  MonitorElement* mqgLikelihood_lowPt_EndCap;
  MonitorElement* mAxis2_mediumPt_EndCap;
  MonitorElement* mpTD_mediumPt_EndCap;
  MonitorElement* mMultiplicityQG_mediumPt_EndCap;
  MonitorElement* mqgLikelihood_mediumPt_EndCap;
  MonitorElement* mAxis2_highPt_EndCap;
  MonitorElement* mpTD_highPt_EndCap;
  MonitorElement* mMultiplicityQG_highPt_EndCap;
  MonitorElement* mqgLikelihood_highPt_EndCap;

  MonitorElement* mAxis2_lowPt_Forward;
  MonitorElement* mpTD_lowPt_Forward;
  MonitorElement* mMultiplicityQG_lowPt_Forward;
  MonitorElement* mqgLikelihood_lowPt_Forward;
  MonitorElement* mAxis2_mediumPt_Forward;
  MonitorElement* mpTD_mediumPt_Forward;
  MonitorElement* mMultiplicityQG_mediumPt_Forward;
  MonitorElement* mqgLikelihood_mediumPt_Forward;
  MonitorElement* mAxis2_highPt_Forward;
  MonitorElement* mpTD_highPt_Forward;
  MonitorElement* mMultiplicityQG_highPt_Forward;
  MonitorElement* mqgLikelihood_highPt_Forward;

  //new Plots with Res./ Eff. as function of neutral, charged &  em fraction

  MonitorElement* mNeutralFraction;
  MonitorElement* mNeutralFraction2;

  MonitorElement* mEEffNeutralFraction;
  MonitorElement* mEEffChargedFraction;
  MonitorElement* mEResNeutralFraction;
  MonitorElement* mEResChargedFraction;
  MonitorElement* nEEff;
  //PF specific NPV profiles
  MonitorElement* mChargedHadronEnergy_profile;
  MonitorElement* mNeutralHadronEnergy_profile;
  MonitorElement* mChargedEmEnergy_profile;
  MonitorElement* mChargedMuEnergy_profile;
  MonitorElement* mNeutralEmEnergy_profile;
  MonitorElement* mChargedMultiplicity_profile;
  MonitorElement* mNeutralMultiplicity_profile;
  MonitorElement* mMuonMultiplicity_profile;

  //Monitor Elements for special selections
  //for special selections
  MonitorElement* mCHFrac;
  MonitorElement* mNHFrac;
  MonitorElement* mPhFrac;
  MonitorElement* mHFEMFrac;
  MonitorElement* mHFHFrac;
  MonitorElement* mCHFrac_profile;
  MonitorElement* mNHFrac_profile;
  MonitorElement* mPhFrac_profile;
  MonitorElement* mHFEMFrac_profile;
  MonitorElement* mHFHFrac_profile;

  JetMETDQMDCSFilter* DCSFilterForJetMonitoring_;
  JetMETDQMDCSFilter* DCSFilterForDCSMonitoring_;
  /*
  MonitorElement* mePhFracBarrel_BXm2BXm1Empty;
  MonitorElement* meNHFracBarrel_BXm2BXm1Empty;
  MonitorElement* meCHFracBarrel_BXm2BXm1Empty;
  MonitorElement* mePtBarrel_BXm2BXm1Empty;
  MonitorElement* mePhFracEndCapMinus_BXm2BXm1Empty;
  MonitorElement* meNHFracEndCapMinus_BXm2BXm1Empty;
  MonitorElement* meCHFracEndCapMinus_BXm2BXm1Empty;
  MonitorElement* mePtEndCapMinus_BXm2BXm1Empty;
  MonitorElement* mePhFracEndCapPlus_BXm2BXm1Empty;
  MonitorElement* meNHFracEndCapPlus_BXm2BXm1Empty;
  MonitorElement* meCHFracEndCapPlus_BXm2BXm1Empty;
  MonitorElement* mePtEndCapPlus_BXm2BXm1Empty;
  MonitorElement* meHFHFracMinus_BXm2BXm1Empty;
  MonitorElement* meHFEMFracMinus_BXm2BXm1Empty;
  MonitorElement* mePtForwardMinus_BXm2BXm1Empty;
  MonitorElement* meHFHFracPlus_BXm2BXm1Empty;
  MonitorElement* meHFEMFracPlus_BXm2BXm1Empty;
  MonitorElement* mePtForwardPlus_BXm2BXm1Empty;
  MonitorElement* meEta_BXm2BXm1Empty;
  */
  MonitorElement* mePhFracBarrel_BXm1Empty;
  MonitorElement* meNHFracBarrel_BXm1Empty;
  MonitorElement* meCHFracBarrel_BXm1Empty;
  MonitorElement* mePtBarrel_BXm1Empty;
  MonitorElement* mePhFracEndCapMinus_BXm1Empty;
  MonitorElement* meNHFracEndCapMinus_BXm1Empty;
  MonitorElement* meCHFracEndCapMinus_BXm1Empty;
  MonitorElement* mePtEndCapMinus_BXm1Empty;
  MonitorElement* mePhFracEndCapPlus_BXm1Empty;
  MonitorElement* meNHFracEndCapPlus_BXm1Empty;
  MonitorElement* meCHFracEndCapPlus_BXm1Empty;
  MonitorElement* mePtEndCapPlus_BXm1Empty;
  MonitorElement* meHFHFracMinus_BXm1Empty;
  MonitorElement* meHFEMFracMinus_BXm1Empty;
  MonitorElement* mePtForwardMinus_BXm1Empty;
  MonitorElement* meHFHFracPlus_BXm1Empty;
  MonitorElement* meHFEMFracPlus_BXm1Empty;
  MonitorElement* mePtForwardPlus_BXm1Empty;
  MonitorElement* meEta_BXm1Empty;
  /*
  MonitorElement* mePhFracBarrel_BXm2BXm1Filled;
  MonitorElement* meNHFracBarrel_BXm2BXm1Filled;
  MonitorElement* meCHFracBarrel_BXm2BXm1Filled;
  MonitorElement* mePtBarrel_BXm2BXm1Filled;
  MonitorElement* mePhFracEndCapMinus_BXm2BXm1Filled;
  MonitorElement* meNHFracEndCapMinus_BXm2BXm1Filled;
  MonitorElement* meCHFracEndCapMinus_BXm2BXm1Filled;
  MonitorElement* mePtEndCapMinus_BXm2BXm1Filled;
  MonitorElement* mePhFracEndCapPlus_BXm2BXm1Filled;
  MonitorElement* meNHFracEndCapPlus_BXm2BXm1Filled;
  MonitorElement* meCHFracEndCapPlus_BXm2BXm1Filled;
  MonitorElement* mePtEndCapPlus_BXm2BXm1Filled;
  MonitorElement* meHFHFracMinus_BXm2BXm1Filled;
  MonitorElement* meHFEMFracMinus_BXm2BXm1Filled;
  MonitorElement* mePtForwardMinus_BXm2BXm1Filled;
  MonitorElement* meHFHFracPlus_BXm2BXm1Filled;
  MonitorElement* meHFEMFracPlus_BXm2BXm1Filled;
  MonitorElement* mePtForwardPlus_BXm2BXm1Filled;
  MonitorElement* meEta_BXm2BXm1Filled;
  */
  MonitorElement* mePhFracBarrel_BXm1Filled;
  MonitorElement* meNHFracBarrel_BXm1Filled;
  MonitorElement* meCHFracBarrel_BXm1Filled;
  MonitorElement* mePtBarrel_BXm1Filled;
  MonitorElement* mePhFracEndCapMinus_BXm1Filled;
  MonitorElement* meNHFracEndCapMinus_BXm1Filled;
  MonitorElement* meCHFracEndCapMinus_BXm1Filled;
  MonitorElement* mePtEndCapMinus_BXm1Filled;
  MonitorElement* mePhFracEndCapPlus_BXm1Filled;
  MonitorElement* meNHFracEndCapPlus_BXm1Filled;
  MonitorElement* meCHFracEndCapPlus_BXm1Filled;
  MonitorElement* mePtEndCapPlus_BXm1Filled;
  MonitorElement* meHFHFracMinus_BXm1Filled;
  MonitorElement* meHFEMFracMinus_BXm1Filled;
  MonitorElement* mePtForwardMinus_BXm1Filled;
  MonitorElement* meHFHFracPlus_BXm1Filled;
  MonitorElement* meHFEMFracPlus_BXm1Filled;
  MonitorElement* mePtForwardPlus_BXm1Filled;
  MonitorElement* meEta_BXm1Filled;

  //miniaod specific variables, especially for substructure
  MonitorElement* mSoftDropMass;
  MonitorElement* mPrunedMass;
  MonitorElement* mTrimmedMass;
  MonitorElement* mFilteredMass;
  MonitorElement* mtau2_over_tau1;
  MonitorElement* mtau3_over_tau2;
  MonitorElement* mCATopTag_topMass;
  MonitorElement* mCATopTag_minMass;
  MonitorElement* mCATopTag_nSubJets;

  MonitorElement* mnSubJetsCMSTopTag;
  MonitorElement* mSubJet1_CMSTopTag_pt;
  MonitorElement* mSubJet1_CMSTopTag_eta;
  MonitorElement* mSubJet1_CMSTopTag_phi;
  MonitorElement* mSubJet1_CMSTopTag_mass;
  MonitorElement* mSubJet2_CMSTopTag_pt;
  MonitorElement* mSubJet2_CMSTopTag_eta;
  MonitorElement* mSubJet2_CMSTopTag_phi;
  MonitorElement* mSubJet2_CMSTopTag_mass;
  MonitorElement* mSubJet3_CMSTopTag_pt;
  MonitorElement* mSubJet3_CMSTopTag_eta;
  MonitorElement* mSubJet3_CMSTopTag_phi;
  MonitorElement* mSubJet3_CMSTopTag_mass;
  MonitorElement* mSubJet4_CMSTopTag_pt;
  MonitorElement* mSubJet4_CMSTopTag_eta;
  MonitorElement* mSubJet4_CMSTopTag_phi;
  MonitorElement* mSubJet4_CMSTopTag_mass;

  MonitorElement* mnSubJetsSoftDrop;
  MonitorElement* mSubJet1_SoftDrop_pt;
  MonitorElement* mSubJet1_SoftDrop_eta;
  MonitorElement* mSubJet1_SoftDrop_phi;
  MonitorElement* mSubJet1_SoftDrop_mass;
  MonitorElement* mSubJet2_SoftDrop_pt;
  MonitorElement* mSubJet2_SoftDrop_eta;
  MonitorElement* mSubJet2_SoftDrop_phi;
  MonitorElement* mSubJet2_SoftDrop_mass;

  //miniaod specific variables, especially for substructure for a boosted regime
  MonitorElement* mSoftDropMass_boosted;
  MonitorElement* mPrunedMass_boosted;
  MonitorElement* mTrimmedMass_boosted;
  MonitorElement* mFilteredMass_boosted;
  MonitorElement* mtau2_over_tau1_boosted;
  MonitorElement* mtau3_over_tau2_boosted;
  MonitorElement* mCATopTag_topMass_boosted;
  MonitorElement* mCATopTag_minMass_boosted;
  MonitorElement* mCATopTag_nSubJets_boosted;

  MonitorElement* mnSubJetsCMSTopTag_boosted;
  MonitorElement* mSubJet1_CMSTopTag_pt_boosted;
  MonitorElement* mSubJet1_CMSTopTag_eta_boosted;
  MonitorElement* mSubJet1_CMSTopTag_phi_boosted;
  MonitorElement* mSubJet1_CMSTopTag_mass_boosted;
  MonitorElement* mSubJet2_CMSTopTag_pt_boosted;
  MonitorElement* mSubJet2_CMSTopTag_eta_boosted;
  MonitorElement* mSubJet2_CMSTopTag_phi_boosted;
  MonitorElement* mSubJet2_CMSTopTag_mass_boosted;
  MonitorElement* mSubJet3_CMSTopTag_pt_boosted;
  MonitorElement* mSubJet3_CMSTopTag_eta_boosted;
  MonitorElement* mSubJet3_CMSTopTag_phi_boosted;
  MonitorElement* mSubJet3_CMSTopTag_mass_boosted;
  MonitorElement* mSubJet4_CMSTopTag_pt_boosted;
  MonitorElement* mSubJet4_CMSTopTag_eta_boosted;
  MonitorElement* mSubJet4_CMSTopTag_phi_boosted;
  MonitorElement* mSubJet4_CMSTopTag_mass_boosted;

  MonitorElement* mnSubJetsSoftDrop_boosted;
  MonitorElement* mSubJet1_SoftDrop_pt_boosted;
  MonitorElement* mSubJet1_SoftDrop_eta_boosted;
  MonitorElement* mSubJet1_SoftDrop_phi_boosted;
  MonitorElement* mSubJet1_SoftDrop_mass_boosted;
  MonitorElement* mSubJet2_SoftDrop_pt_boosted;
  MonitorElement* mSubJet2_SoftDrop_eta_boosted;
  MonitorElement* mSubJet2_SoftDrop_phi_boosted;
  MonitorElement* mSubJet2_SoftDrop_mass_boosted;

  //miniaod only variables
  MonitorElement* mPt_CaloJet;
  MonitorElement* mEMF_CaloJet;
  MonitorElement* mMass_Barrel;
  MonitorElement* mMass_EndCap;
  MonitorElement* mMass_Forward;

  //now ZJets plots
  MonitorElement* mDPhiZJet;
  MonitorElement* mZMass;
  MonitorElement* mZJetAsymmetry;
  MonitorElement* mJetZBalance_lowZPt_J_Barrel;
  MonitorElement* mJetZBalance_mediumZPt_J_Barrel;
  MonitorElement* mJetZBalance_highZPt_J_Barrel;
  MonitorElement* mJetZBalance_lowZPt_J_EndCap;
  MonitorElement* mJetZBalance_mediumZPt_J_EndCap;
  MonitorElement* mJetZBalance_highZPt_J_EndCap;
  MonitorElement* mJetZBalance_lowZPt_J_Forward;
  MonitorElement* mJetZBalance_mediumZPt_J_Forward;
  MonitorElement* mJetZBalance_highZPt_J_Forward;
  MonitorElement* mJ1Pt_over_ZPt_J_Barrel;
  MonitorElement* mJ1Pt_over_ZPt_J_EndCap;
  MonitorElement* mJ1Pt_over_ZPt_J_Forward;
  MonitorElement* mMPF_J_Barrel;
  MonitorElement* mMPF_J_EndCap;
  MonitorElement* mMPF_J_Forward;
  MonitorElement* mJ1Pt_over_ZPt_lowZPt_J_Barrel;
  MonitorElement* mJ1Pt_over_ZPt_mediumZPt_J_Barrel;
  MonitorElement* mJ1Pt_over_ZPt_highZPt_J_Barrel;
  MonitorElement* mJ1Pt_over_ZPt_lowZPt_J_EndCap;
  MonitorElement* mJ1Pt_over_ZPt_mediumZPt_J_EndCap;
  MonitorElement* mJ1Pt_over_ZPt_highZPt_J_EndCap;
  MonitorElement* mJ1Pt_over_ZPt_lowZPt_J_Forward;
  MonitorElement* mJ1Pt_over_ZPt_mediumZPt_J_Forward;
  MonitorElement* mJ1Pt_over_ZPt_highZPt_J_Forward;
  MonitorElement* mMPF_lowZPt_J_Barrel;
  MonitorElement* mMPF_mediumZPt_J_Barrel;
  MonitorElement* mMPF_highZPt_J_Barrel;
  MonitorElement* mMPF_lowZPt_J_EndCap;
  MonitorElement* mMPF_mediumZPt_J_EndCap;
  MonitorElement* mMPF_highZPt_J_EndCap;
  MonitorElement* mMPF_lowZPt_J_Forward;
  MonitorElement* mMPF_mediumZPt_J_Forward;
  MonitorElement* mMPF_highZPt_J_Forward;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_30_55_J_Barrel;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_55_75_J_Barrel;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_75_150_J_Barrel;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_150_290_J_Barrel;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_290_J_Barrel;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_30_55_J_EndCap;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_55_75_J_EndCap;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_75_150_J_EndCap;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_150_290_J_EndCap;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_290_J_EndCap;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_30_55_J_Forward;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_55_100_J_Forward;
  MonitorElement* mDeltaPt_Z_j1_over_ZPt_100_J_Forward;

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

  bool isCaloJet_;
  bool isPFJet_;
  bool isMiniAODJet_;

  bool fill_jet_high_level_histo;

  bool fill_CHS_histos;
};
#endif