EfficiencyPlotter.h

Go to the documentation of this file.

#ifndef EfficiencyPlotter_H
#define EfficiencyPlotter_H


#include "FWCore/Framework/interface/Frameworkfwd.h"
#include <FWCore/Framework/interface/EDAnalyzer.h>
#include "DataFormats/Common/interface/Handle.h"
#include <FWCore/Framework/interface/ESHandle.h>
#include <FWCore/Framework/interface/Event.h>
#include <FWCore/Framework/interface/MakerMacros.h>
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <FWCore/Framework/interface/LuminosityBlock.h>

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/Run.h"

#include <memory>
#include <iostream>
#include <string>
#include <vector>
#include <map>
#include "TH1F.h"

class EfficiencyPlotter: public edm::EDAnalyzer{

public:

  EfficiencyPlotter(const edm::ParameterSet& ps);
  
  virtual ~EfficiencyPlotter();

protected:

  void beginJob(void);

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

  void endJob();

  void beginRun(edm::Run const& run, edm::EventSetup const& eSetup);
  void beginLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& context) ;

  void endLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& c);
  void endRun(edm::Run const& run, edm::EventSetup const& eSetup);


private:

  // counters
  int nevents;
  unsigned int nLumiSegs;
  int prescaleFactor;
  int run;
  // Switch for verbosity
  std::string metname;

  DQMStore* theDbe;
  edm::ParameterSet parameters;

   //histo binning parameters
  int etaBin;
  double etaMin;
  double etaMax;

  int ptBin;
  double ptMin;
  double ptMax;

  int phiBin;
  double phiMin;
  double phiMax;
  
  int vtxBin;
  double vtxMin;
  double vtxMax;

  // efficiency histograms
  MonitorElement* h_eff_pt_TightMu;
  MonitorElement* h_eff_pt_EB_TightMu;
  MonitorElement* h_eff_pt_EE_TightMu;
  MonitorElement* h_eff_eta_TightMu;
  MonitorElement* h_eff_hp_eta_TightMu;
  MonitorElement* h_eff_phi_TightMu;
  MonitorElement* h_eff_pt_detIsoTightMu;
  MonitorElement* h_eff_pt_EB_detIsoTightMu;
  MonitorElement* h_eff_pt_EE_detIsoTightMu;
  MonitorElement* h_eff_pt_pfIsoTightMu;
  MonitorElement* h_eff_pt_EB_pfIsoTightMu;
  MonitorElement* h_eff_pt_EE_pfIsoTightMu;

  MonitorElement* h_eff_vtx_detIsoTightMu;
  MonitorElement* h_eff_vtx_pfIsoTightMu;
  MonitorElement* h_eff_vtx_EB_detIsoTightMu;
  MonitorElement* h_eff_vtx_EE_detIsoTightMu;
  MonitorElement* h_eff_vtx_EB_pfIsoTightMu;
  MonitorElement* h_eff_vtx_EE_pfIsoTightMu;

  MonitorElement* h_eff_pt_pfIsodBTightMu;
  MonitorElement* h_eff_pt_EB_pfIsodBTightMu;
  MonitorElement* h_eff_pt_EE_pfIsodBTightMu;

  MonitorElement* h_eff_vtx_pfIsodBTightMu;
  MonitorElement* h_eff_vtx_EB_pfIsodBTightMu;
  MonitorElement* h_eff_vtx_EE_pfIsodBTightMu;

};

#endif