DataCertificationJetMET.h

Go to the documentation of this file.

#ifndef DATACERTIFICATIONJETMET_H
#define DATACERTIFICATIONJETMET_H
 
// author: Kenichi Hatakeyama (Rockefeller U.)
 
// system include files
#include <memory>
#include <cstdio>
#include <cmath>
#include <sstream>
 
// user include files
#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 "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "DQMServices/Core/interface/DQMEDHarvester.h"
//
// class decleration
//
 
class DataCertificationJetMET : public DQMEDHarvester {
public:
  explicit DataCertificationJetMET(const edm::ParameterSet&);
  ~DataCertificationJetMET() override;
 
private:
  void dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter&) override;
 
  MonitorElement* reportSummary;
  MonitorElement* CertificationSummary;
  MonitorElement* reportSummaryMap;
  MonitorElement* CertificationSummaryMap;
 
  // ----------member data ---------------------------
 
  edm::ParameterSet conf_;
  edm::Service<TFileService> fs_;
  int verbose_;
  bool InMemory_;
  bool isData;
  std::string metFolder;
  std::string jetAlgo;
 
  edm::InputTag inputMETLabelRECO_;
  edm::InputTag inputMETLabelRECOUncleaned_;
  edm::InputTag inputMETLabelMiniAOD_;
  edm::InputTag inputJetLabelRECO_;
  edm::InputTag inputJetLabelMiniAOD_;
 
  int nbinsPV_;
  double nPVMin_;
  double nPVMax_;
 
  int etaBin_;
  double etaMin_;
  double etaMax_;
 
  int ptBin_;
  double ptMin_;
  double ptMax_;
 
  std::string folderName;
 
  bool caloJetMeanTest;
  bool caloJetKSTest;
  bool pfJetMeanTest;
  bool pfJetKSTest;
  bool jptJetMeanTest;
  bool jptJetKSTest;
  bool caloMETMeanTest;
  bool caloMETKSTest;
  bool pfMETMeanTest;
  bool pfMETKSTest;
  bool tcMETMeanTest;
  bool tcMETKSTest;
 
  bool jetTests[5][2];  //one for each type of jet certification/test type
  bool metTests[5][2];  //one for each type of met certification/test type
 
  bool isHI;
 
  //MET: filter efficiencies, started from uncleaned directories
  MonitorElement* mMET_EffHBHENoiseFilter;
  MonitorElement* mMET_EffCSCTightHaloFilter;
  MonitorElement* mMET_EffeeBadScFilter;
  MonitorElement* mMET_EffEcalDeadCellTriggerFilter;
  MonitorElement* mMET_EffEcalDeadCellBoundaryFilter;
  MonitorElement* mMET_EffHBHEIsoNoiseFilter;
  MonitorElement* mMET_EffCSCTightHalo2015Filter;
  MonitorElement* mMET_EffHcalStripHaloFilter;
 
  //MET: RECO vs MiniAOD histos
  MonitorElement* mMET_MiniAOD_over_Reco;
  MonitorElement* mMEy_MiniAOD_over_Reco;
  MonitorElement* mSumET_MiniAOD_over_Reco;
  MonitorElement* mMETPhi_MiniAOD_over_Reco;
  MonitorElement* mMET_logx_MiniAOD_over_Reco;
  MonitorElement* mSumET_logx_MiniAOD_over_Reco;
  MonitorElement* mChargedHadronEtFraction_MiniAOD_over_Reco;
  MonitorElement* mNeutralHadronEtFraction_MiniAOD_over_Reco;
  MonitorElement* mPhotonEtFraction_MiniAOD_over_Reco;
  MonitorElement* mHFHadronEtFraction_MiniAOD_over_Reco;
  MonitorElement* mHFEMEtFraction_MiniAOD_over_Reco;
  MonitorElement* mMET_nVtx_profile_MiniAOD_over_Reco;
  MonitorElement* mSumET_nVtx_profile_MiniAOD_over_Reco;
  MonitorElement* mChargedHadronEtFraction_nVtx_profile_MiniAOD_over_Reco;
  MonitorElement* mNeutralHadronEtFraction_nVtx_profile_MiniAOD_over_Reco;
  MonitorElement* mPhotonEtFraction_nVtx_profile_MiniAOD_over_Reco;
 
  //Jets: RECO vs MiniAOD histos
  MonitorElement* mPt_MiniAOD_over_Reco;
  MonitorElement* mEta_MiniAOD_over_Reco;
  MonitorElement* mPhi_MiniAOD_over_Reco;
  MonitorElement* mNjets_MiniAOD_over_Reco;
  MonitorElement* mPt_uncor_MiniAOD_over_Reco;
  MonitorElement* mEta_uncor_MiniAOD_over_Reco;
  MonitorElement* mPhi_uncor_MiniAOD_over_Reco;
  MonitorElement* mJetEnergyCorr_MiniAOD_over_Reco;
  MonitorElement* mJetEnergyCorrVSeta_MiniAOD_over_Reco;
  MonitorElement* mDPhi_MiniAOD_over_Reco;
  MonitorElement* mLooseJIDPassFractionVSeta_MiniAOD_over_Reco;
  MonitorElement* mPt_Barrel_MiniAOD_over_Reco;
  MonitorElement* mPt_EndCap_MiniAOD_over_Reco;
  MonitorElement* mPt_Forward_MiniAOD_over_Reco;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_Barrel_MiniAOD_over_Reco;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_EndCap_MiniAOD_over_Reco;
  MonitorElement* mMVAPUJIDDiscriminant_lowPt_Forward_MiniAOD_over_Reco;
  MonitorElement* mMVAPUJIDDiscriminant_mediumPt_EndCap_MiniAOD_over_Reco;
  MonitorElement* mMVAPUJIDDiscriminant_highPt_Barrel_MiniAOD_over_Reco;
  MonitorElement* mCHFracVSpT_Barrel_MiniAOD_over_Reco;
  MonitorElement* mNHFracVSpT_EndCap_MiniAOD_over_Reco;
  MonitorElement* mPhFracVSpT_Barrel_MiniAOD_over_Reco;
  MonitorElement* mHFHFracVSpT_Forward_MiniAOD_over_Reco;
  MonitorElement* mHFEFracVSpT_Forward_MiniAOD_over_Reco;
  MonitorElement* mCHFrac_MiniAOD_over_Reco;
  MonitorElement* mNHFrac_MiniAOD_over_Reco;
  MonitorElement* mPhFrac_MiniAOD_over_Reco;
  MonitorElement* mChargedMultiplicity_MiniAOD_over_Reco;
  MonitorElement* mNeutralMultiplicity_MiniAOD_over_Reco;
  MonitorElement* mMuonMultiplicity_MiniAOD_over_Reco;
  MonitorElement* mNeutralFraction_MiniAOD_over_Reco;
};
 
#endif