METTester.h

Go to the documentation of this file.

#ifndef METTESTER_H
#define METTESTER_H

// author: Mike Schmitt (The University of Florida)
// date: 8/24/2006
// modification: Bobby Scurlock
// date: 03.11.2006
// note: added RMS(METx) vs SumET capability
// modification: Rick Cavanaugh
// date: 05.11.2006
// note: added configuration parameters
// modification: Mike Schmitt
// date: 02.28.2007
// note: code rewrite

// Rewritten by Viola Sordini, Matthias Artur Weber, Robert Schoefbeck Nov./Dez.
// 2013

#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/METReco/interface/GenMETCollection.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METCollection.h"
#include "DataFormats/METReco/interface/METFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonMETCorrectionData.h"
#include "DataFormats/PatCandidates/interface/MET.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "TMath.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

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

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

private:
  std::map<std::string, MonitorElement *> me;

  // Inputs from Configuration File

  edm::InputTag mInputCollection_;
  edm::InputTag inputMETLabel_;
  std::string METType_;

  edm::InputTag inputCaloMETLabel_;

  // Tokens
  edm::EDGetTokenT<std::vector<reco::Vertex>> pvToken_;
  edm::EDGetTokenT<reco::CaloMETCollection> caloMETsToken_;
  edm::EDGetTokenT<reco::PFMETCollection> pfMETsToken_;
  // edm::EDGetTokenT<reco::METCollection> tcMETsToken_;
  edm::EDGetTokenT<reco::GenMETCollection> genMETsToken_;
  edm::EDGetTokenT<reco::GenMETCollection> genMETsTrueToken_;
  edm::EDGetTokenT<reco::GenMETCollection> genMETsCaloToken_;
  edm::EDGetTokenT<pat::METCollection> patMETToken_;

  // Events variables
  MonitorElement *mNvertex;

  // Common variables
  MonitorElement *mMEx;
  MonitorElement *mMEy;
  MonitorElement *mMETSig;
  MonitorElement *mMET;
  MonitorElement *mMETFine;
  MonitorElement *mMET_Nvtx;
  MonitorElement *mMETPhi;
  MonitorElement *mSumET;
  MonitorElement *mMETDifference_GenMETTrue;
  MonitorElement *mMETDeltaPhi_GenMETTrue;
  MonitorElement *mMETDifference_GenMETCalo;
  MonitorElement *mMETDeltaPhi_GenMETCalo;

  // MET Uncertainity Variables
  MonitorElement *mMETUnc_JetResUp;
  MonitorElement *mMETUnc_JetResDown;
  MonitorElement *mMETUnc_JetEnUp;
  MonitorElement *mMETUnc_JetEnDown;
  MonitorElement *mMETUnc_MuonEnUp;
  MonitorElement *mMETUnc_MuonEnDown;
  MonitorElement *mMETUnc_ElectronEnUp;
  MonitorElement *mMETUnc_ElectronEnDown;
  MonitorElement *mMETUnc_TauEnUp;
  MonitorElement *mMETUnc_TauEnDown;
  MonitorElement *mMETUnc_UnclusteredEnUp;
  MonitorElement *mMETUnc_UnclusteredEnDown;
  MonitorElement *mMETUnc_PhotonEnUp;
  MonitorElement *mMETUnc_PhotonEnDown;

  // CaloMET variables

  MonitorElement *mCaloMaxEtInEmTowers;
  MonitorElement *mCaloMaxEtInHadTowers;
  MonitorElement *mCaloEtFractionHadronic;
  MonitorElement *mCaloEmEtFraction;
  MonitorElement *mCaloHadEtInHB;
  MonitorElement *mCaloHadEtInHO;
  MonitorElement *mCaloHadEtInHE;
  MonitorElement *mCaloHadEtInHF;
  MonitorElement *mCaloHadEtInEB;
  MonitorElement *mCaloHadEtInEE;
  MonitorElement *mCaloEmEtInHF;
  MonitorElement *mCaloSETInpHF;
  MonitorElement *mCaloSETInmHF;
  MonitorElement *mCaloEmEtInEE;
  MonitorElement *mCaloEmEtInEB;

  // GenMET variables
  MonitorElement *mNeutralEMEtFraction;
  MonitorElement *mNeutralHadEtFraction;
  MonitorElement *mChargedEMEtFraction;
  MonitorElement *mChargedHadEtFraction;
  MonitorElement *mMuonEtFraction;
  MonitorElement *mInvisibleEtFraction;

  // MET variables

  // PFMET variables
  MonitorElement *mPFphotonEtFraction;
  MonitorElement *mPFphotonEt;
  MonitorElement *mPFneutralHadronEtFraction;
  MonitorElement *mPFneutralHadronEt;
  MonitorElement *mPFelectronEtFraction;
  MonitorElement *mPFelectronEt;
  MonitorElement *mPFchargedHadronEtFraction;
  MonitorElement *mPFchargedHadronEt;
  MonitorElement *mPFmuonEtFraction;
  MonitorElement *mPFmuonEt;
  MonitorElement *mPFHFHadronEtFraction;
  MonitorElement *mPFHFHadronEt;
  MonitorElement *mPFHFEMEtFraction;
  MonitorElement *mPFHFEMEt;

  MonitorElement *mMETDifference_GenMETTrue_MET0to20;
  MonitorElement *mMETDifference_GenMETTrue_MET20to40;
  MonitorElement *mMETDifference_GenMETTrue_MET40to60;
  MonitorElement *mMETDifference_GenMETTrue_MET60to80;
  MonitorElement *mMETDifference_GenMETTrue_MET80to100;
  MonitorElement *mMETDifference_GenMETTrue_MET100to150;
  MonitorElement *mMETDifference_GenMETTrue_MET150to200;
  MonitorElement *mMETDifference_GenMETTrue_MET200to300;
  MonitorElement *mMETDifference_GenMETTrue_MET300to400;
  MonitorElement *mMETDifference_GenMETTrue_MET400to500;
  MonitorElement *mMETDifference_GenMETTrue_MET500;
  // moved into postprocessor
  // MonitorElement* mMETDifference_GenMETTrue_METResolution;

  bool isCaloMET;
  //  bool isCorMET;
  //  bool isTcMET;
  bool isPFMET;
  bool isGenMET;
  bool isMiniAODMET;
};

#endif  // METTESTER_H