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 "DataFormats/Common/interface/Handle.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METCollection.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/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/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METFwd.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 "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>
#include "DataFormats/PatCandidates/interface/MET.h"
#include "TMath.h"


//class METTester: public edm::EDAnalyzer {
class METTester: public DQMEDAnalyzer {
public:

  explicit METTester(const edm::ParameterSet&);

  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  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;
  //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