EfficiencyAnalyzer Class Reference

#include <EfficiencyAnalyzer.h>

Inheritance diagram for EfficiencyAnalyzer:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	EfficiencyAnalyzer (const edm::ParameterSet &pset)
	~EfficiencyAnalyzer () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
int	_numPV
bool	doPVCheck_
int	etaBin_
double	etaMax_
double	etaMin_
MonitorElement *	h_allProbes_EB_ID_nVtx
MonitorElement *	h_allProbes_EB_ID_pt
MonitorElement *	h_allProbes_EB_pt
MonitorElement *	h_allProbes_EE_ID_nVtx
MonitorElement *	h_allProbes_EE_ID_pt
MonitorElement *	h_allProbes_EE_pt
MonitorElement *	h_allProbes_eta
MonitorElement *	h_allProbes_hp_eta
MonitorElement *	h_allProbes_ID_nVtx
MonitorElement *	h_allProbes_ID_pt
MonitorElement *	h_allProbes_inner_eta
MonitorElement *	h_allProbes_inner_phi
MonitorElement *	h_allProbes_inner_pt
MonitorElement *	h_allProbes_phi
MonitorElement *	h_allProbes_pt
MonitorElement *	h_failProbes_ID_eta
MonitorElement *	h_failProbes_ID_phi
MonitorElement *	h_failProbes_ID_pt
MonitorElement *	h_passProbes_detIsoID_nVtx
MonitorElement *	h_passProbes_detIsoID_pt
MonitorElement *	h_passProbes_EB_detIsoID_nVtx
MonitorElement *	h_passProbes_EB_detIsoID_pt
MonitorElement *	h_passProbes_EB_pfIsodBID_nVtx
MonitorElement *	h_passProbes_EB_pfIsodBID_pt
MonitorElement *	h_passProbes_EB_pfIsoID_nVtx
MonitorElement *	h_passProbes_EB_pfIsoID_pt
MonitorElement *	h_passProbes_EE_detIsoID_nVtx
MonitorElement *	h_passProbes_EE_detIsoID_pt
MonitorElement *	h_passProbes_EE_pfIsodBID_nVtx
MonitorElement *	h_passProbes_EE_pfIsodBID_pt
MonitorElement *	h_passProbes_EE_pfIsoID_nVtx
MonitorElement *	h_passProbes_EE_pfIsoID_pt
MonitorElement *	h_passProbes_ID_EB_pt
MonitorElement *	h_passProbes_ID_EE_pt
MonitorElement *	h_passProbes_ID_eta
MonitorElement *	h_passProbes_ID_hp_eta
MonitorElement *	h_passProbes_ID_inner_eta
MonitorElement *	h_passProbes_ID_inner_phi
MonitorElement *	h_passProbes_ID_inner_pt
MonitorElement *	h_passProbes_ID_phi
MonitorElement *	h_passProbes_ID_pt
MonitorElement *	h_passProbes_pfIsodBID_nVtx
MonitorElement *	h_passProbes_pfIsodBID_pt
MonitorElement *	h_passProbes_pfIsoID_nVtx
MonitorElement *	h_passProbes_pfIsoID_pt
std::string	ID_
std::string	metname
edm::ParameterSet	parameters
int	phiBin_
double	phiMax_
double	phiMin_
int	ptBin_
double	ptMax_
double	ptMin_
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_
std::string	theFolder
edm::EDGetTokenT< edm::View< reco::Muon > >	theMuonCollectionLabel_
edm::EDGetTokenT< reco::TrackCollection >	theTrackCollectionLabel_
edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_
int	vtxBin_
double	vtxMax_
double	vtxMin_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Class EfficiencyAnalyzer

DQM monitoring for dimuon mass

Author: S.Folgueras, A. Calderon

Definition at line 31 of file EfficiencyAnalyzer.h.

Constructor & Destructor Documentation

EfficiencyAnalyzer()

EfficiencyAnalyzer::EfficiencyAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 29 of file EfficiencyAnalyzer.cc.

                                                                  {
  parameters = pSet;

  // DATA
  theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> >(parameters.getParameter<edm::InputTag>("MuonCollection"));
  theTrackCollectionLabel_ = consumes<reco::TrackCollection>(parameters.getParameter<edm::InputTag>("TrackCollection"));
  theVertexLabel_ = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
  theBeamSpotLabel_ = mayConsume<reco::BeamSpot>(parameters.getParameter<edm::InputTag>("BeamSpotLabel"));

  //Vertex requirements
  doPVCheck_ = parameters.getParameter<bool>("doPrimaryVertexCheck");

  ptBin_ = parameters.getParameter<int>("ptBin");
  ptMin_ = parameters.getParameter<double>("ptMin");
  ptMax_ = parameters.getParameter<double>("ptMax");

  etaBin_ = parameters.getParameter<int>("etaBin");
  etaMin_ = parameters.getParameter<double>("etaMin");
  etaMax_ = parameters.getParameter<double>("etaMax");

  phiBin_ = parameters.getParameter<int>("phiBin");
  phiMin_ = parameters.getParameter<double>("phiMin");
  phiMax_ = parameters.getParameter<double>("phiMax");

  vtxBin_ = parameters.getParameter<int>("vtxBin");
  vtxMin_ = parameters.getParameter<double>("vtxMin");
  vtxMax_ = parameters.getParameter<double>("vtxMax");

  ID_ = parameters.getParameter<string>("ID");
  theFolder = parameters.getParameter<string>("folder");
}

~EfficiencyAnalyzer()

EfficiencyAnalyzer::~EfficiencyAnalyzer ( )

override

Definition at line 61 of file EfficiencyAnalyzer.cc.

{}

Member Function Documentation

analyze()

void EfficiencyAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

overridevirtual

to be read from output as "generalTracks"

Reimplemented from DQMEDAnalyzer.

Definition at line 158 of file EfficiencyAnalyzer.cc.

References cms::cuda::bs, boostedTaus_cff::chargedIso, reco::MuonIsolation::emEt, relativeConstraints::error, reco::MuonIsolation::hadEt, iEvent, muon::isLooseMuon(), muon::isMediumMuon(), muon::isTightMuon(), LogTrace, SiStripPI::max, metname, PDWG_BPHSkim_cff::muons, boostedTaus_cff::neutralIso, gedPhotonSequence_cff::photonIso, position, displacedMuons_cfi::pu, reco::MuonIsolation::sumPt, tracks, findQualityFiles::v, bphysicsOniaDQM_cfi::vertex, and spclusmultinvestigator_cfi::vertexCollection.

                                                                                    {
  LogTrace(metname) << "[EfficiencyAnalyzer] Analyze the mu in different eta regions";
  // ==========================================================
  // BEGIN READ DATA:
  // Muon information
  edm::Handle<edm::View<reco::Muon> > muons;
  iEvent.getByToken(theMuonCollectionLabel_, muons);

  // Tracks information
  edm::Handle<reco::TrackCollection> tracks;
  iEvent.getByToken(theTrackCollectionLabel_, tracks);  

  //Vertex information
  edm::Handle<reco::VertexCollection> vertex;
  iEvent.getByToken(theVertexLabel_, vertex);
  // END READ DATA
  // ==========================================================

  _numPV = 0;
  bool bPrimaryVertex = true;
  if (doPVCheck_) {
    bPrimaryVertex = false;

    if (!vertex.isValid()) {
      LogTrace(metname) << "[EfficiencyAnalyzer] Could not find vertex collection" << std::endl;
      bPrimaryVertex = false;
    }

    if (vertex.isValid()) {
      const reco::VertexCollection& vertexCollection = *(vertex.product());
      int vertex_number = vertexCollection.size();

      reco::VertexCollection::const_iterator v = vertexCollection.begin();
      for (; v != vertexCollection.end(); ++v) {
        double vertex_chi2 = v->normalizedChi2();
        double vertex_ndof = v->ndof();
        bool fakeVtx = v->isFake();
        double vertex_Z = v->z();

        if (!fakeVtx && vertex_number >= 1 && vertex_ndof > 4 && vertex_chi2 < 999 && fabs(vertex_Z) < 24.) {
          bPrimaryVertex = true;
          ++_numPV;
        }
      }
    }
  }

  // =================================================================================
  // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
  reco::Vertex::Point posVtx;
  reco::Vertex::Error errVtx;
  unsigned int theIndexOfThePrimaryVertex = 999.;
  if (vertex.isValid()) {
    for (unsigned int ind = 0; ind < vertex->size(); ++ind) {
      if ((*vertex)[ind].isValid() && !((*vertex)[ind].isFake())) {
        theIndexOfThePrimaryVertex = ind;
        break;
      }
    }
  }

  if (theIndexOfThePrimaryVertex < 100) {
    posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
    errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
  } else {
    LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";

    edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
    iEvent.getByToken(theBeamSpotLabel_, recoBeamSpotHandle);
    reco::BeamSpot bs = *recoBeamSpotHandle;

    posVtx = bs.position();
    errVtx(0, 0) = bs.BeamWidthX();
    errVtx(1, 1) = bs.BeamWidthY();
    errVtx(2, 2) = bs.sigmaZ();
  }

  const reco::Vertex thePrimaryVertex(posVtx, errVtx);
  // ==========================================================

  if (!muons.isValid())
    return;

  // Loop on muon collection
  TLorentzVector Mu1, Mu2;

  bool isMB = false;
  bool isME = false;

  for (edm::View<reco::Muon>::const_iterator muon1 = muons->begin(); muon1 != muons->end(); ++muon1) {
    LogTrace(metname) << "[EfficiencyAnalyzer] loop over first muons" << endl;

    //--- Define combined isolation
    reco::MuonIsolation Iso_muon = muon1->isolationR03();
    float combIso = (Iso_muon.emEt + Iso_muon.hadEt + Iso_muon.sumPt);

    //--- Is Global Muon
    if (!muon1->isGlobalMuon())
      continue;

    // get the track combinig the information from both the Tracker and the Spectrometer
    reco::TrackRef recoCombinedGlbTrack1 = muon1->combinedMuon();
    float muPt1 = recoCombinedGlbTrack1->pt();
    Mu1.SetPxPyPzE(recoCombinedGlbTrack1->px(),
                   recoCombinedGlbTrack1->py(),
                   recoCombinedGlbTrack1->pz(),
                   recoCombinedGlbTrack1->p());

    //--- Define if it is a tight muon
    // Change the Tight muon definition by using the implemented method in: MuonSelectors.cc
    if (ID_ == "Loose" && !muon::isLooseMuon(*muon1))
      continue;
    if (ID_ == "Medium" && !muon::isMediumMuon(*muon1))
      continue;
    if (ID_ == "Tight" && !muon::isTightMuon(*muon1, thePrimaryVertex))
      continue;

    //-- is isolated muon
    if (muPt1 <= 15)
      continue;
    if (combIso / muPt1 > 0.1)
      continue;

    for (edm::View<reco::Muon>::const_iterator muon2 = muons->begin(); muon2 != muons->end(); ++muon2) {
      LogTrace(metname) << "[EfficiencyAnalyzer] loop over second muon" << endl;
      if (muon2 == muon1)
        continue;

      if (muon2->eta() < 1.479)
        isMB = true;
      if (muon2->eta() >= 1.479)
        isME = true;

      //--> should we apply track quality cuts???
      Mu2.SetPxPyPzE(muon2->px(), muon2->py(), muon2->pz(), muon2->p());

      float Minv = (Mu1 + Mu2).M();
      if (!muon2->isTrackerMuon())
        continue;
      if (muon2->pt() < 5)
        continue;
      if ((muon1->charge()) * (muon2->charge()) > 0)
        continue;
      if (Minv < 70 || Minv > 110)
        continue;

      h_allProbes_pt->Fill(muon2->pt());
      h_allProbes_eta->Fill(muon2->eta());
      h_allProbes_phi->Fill(muon2->phi());
      if (muon2->innerTrack()->extra().isAvailable()) {
        h_allProbes_inner_pt->Fill(muon2->innerTrack()->innerMomentum().Rho());
        h_allProbes_inner_eta->Fill(muon2->innerTrack()->innerPosition().Eta());
        h_allProbes_inner_phi->Fill(muon2->innerTrack()->innerPosition().Phi());
      }
      if (isMB)
        h_allProbes_EB_pt->Fill(muon2->pt());
      if (isME)
        h_allProbes_EE_pt->Fill(muon2->pt());
      if (muon2->pt() > 20)
        h_allProbes_hp_eta->Fill(muon2->eta());

      // Probes passing the tight muon criteria
      if (ID_ == "Loose" && !muon::isLooseMuon(*muon2))
        continue;
      if (ID_ == "Medium" && !muon::isMediumMuon(*muon2))
        continue;
      if (ID_ == "Tight" && !muon::isTightMuon(*muon2, thePrimaryVertex))
        continue;

      h_passProbes_ID_pt->Fill(muon2->pt());
      h_passProbes_ID_eta->Fill(muon2->eta());
      h_passProbes_ID_phi->Fill(muon2->phi());
      if (muon2->innerTrack()->extra().isAvailable()) {
        h_passProbes_ID_inner_pt->Fill(muon2->innerTrack()->innerMomentum().Rho());
        h_passProbes_ID_inner_eta->Fill(muon2->innerTrack()->innerPosition().Eta());
        h_passProbes_ID_inner_phi->Fill(muon2->innerTrack()->innerPosition().Phi());
      }

      if (isMB)
        h_passProbes_ID_EB_pt->Fill(muon2->pt());
      if (isME)
        h_passProbes_ID_EE_pt->Fill(muon2->pt());
      if (muon2->pt() > 20)
        h_passProbes_ID_hp_eta->Fill(muon2->eta());

      h_allProbes_ID_pt->Fill(muon2->pt());
      if (isMB)
        h_allProbes_EB_ID_pt->Fill(muon2->pt());
      if (isME)
        h_allProbes_EE_ID_pt->Fill(muon2->pt());

      //------- For PU monitoring -------//
      if (bPrimaryVertex)
        h_allProbes_ID_nVtx->Fill(_numPV);
      if (bPrimaryVertex && isMB)
        h_allProbes_EB_ID_nVtx->Fill(_numPV);
      if (bPrimaryVertex && isME)
        h_allProbes_EE_ID_nVtx->Fill(_numPV);

      //-- Define det relative isolation
      float tkIso = muon2->isolationR03().sumPt;
      float emIso = muon2->isolationR03().emEt;
      float hadIso = muon2->isolationR03().hadEt + muon2->isolationR03().hoEt;
      float relDetIso = (tkIso + emIso + hadIso) / (muon2->pt());

      if (relDetIso < 0.05) {
        h_passProbes_detIsoID_pt->Fill(muon2->pt());
        if (isMB)
          h_passProbes_EB_detIsoID_pt->Fill(muon2->pt());
        if (isME)
          h_passProbes_EE_detIsoID_pt->Fill(muon2->pt());

        if (bPrimaryVertex)
          h_passProbes_detIsoID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isMB)
          h_passProbes_EB_detIsoID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isME)
          h_passProbes_EE_detIsoID_nVtx->Fill(_numPV);
      }

      //-- Define PF relative isolation
      float chargedIso = muon2->pfIsolationR04().sumChargedHadronPt;
      float neutralIso = muon2->pfIsolationR04().sumNeutralHadronEt;
      float photonIso = muon2->pfIsolationR04().sumPhotonEt;
      float relPFIso = (chargedIso + neutralIso + photonIso) / (muon2->pt());

      float pu = muon2->pfIsolationR04().sumPUPt;
      float neutralphotonPUCorrected = std::max(0.0, (neutralIso + photonIso - 0.5 * pu));
      float relPFIsoPUCorrected = (chargedIso + neutralphotonPUCorrected) / (muon2->pt());

      if (relPFIso < 0.12) {
        h_passProbes_pfIsoID_pt->Fill(muon2->pt());
        if (isMB)
          h_passProbes_EB_pfIsoID_pt->Fill(muon2->pt());
        if (isME)
          h_passProbes_EE_pfIsoID_pt->Fill(muon2->pt());

        if (bPrimaryVertex)
          h_passProbes_pfIsoID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isMB)
          h_passProbes_EB_pfIsoID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isME)
          h_passProbes_EE_pfIsoID_nVtx->Fill(_numPV);
      }

      // Apply deltaBeta PU corrections to the PF isolation eficiencies.
      if (relPFIsoPUCorrected < 0.12) {
        h_passProbes_pfIsodBID_pt->Fill(muon2->pt());
        if (isMB)
          h_passProbes_EB_pfIsodBID_pt->Fill(muon2->pt());
        if (isME)
          h_passProbes_EE_pfIsodBID_pt->Fill(muon2->pt());

        if (bPrimaryVertex)
          h_passProbes_pfIsodBID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isMB)
          h_passProbes_EB_pfIsodBID_nVtx->Fill(_numPV);
        if (bPrimaryVertex && isME)
          h_passProbes_EE_pfIsodBID_nVtx->Fill(_numPV);
      }
    }
  }
}

bookHistograms()

void EfficiencyAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 63 of file EfficiencyAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), gather_cfg::cout, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                {
  ibooker.cd();
  ibooker.setCurrentFolder(theFolder + ID_);

  h_allProbes_pt = ibooker.book1D("allProbes_pt", "All Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_inner_pt = ibooker.book1D("allProbes_inner_pt", "All Probes inner Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_inner_eta = ibooker.book1D("allProbes_inner_eta", "All Probes inner eta", etaBin_, etaMin_, etaMax_);
  h_allProbes_inner_phi = ibooker.book1D("allProbes_inner_phi", "All Probes inner phi", phiBin_, phiMin_, phiMax_);
  h_allProbes_EB_pt = ibooker.book1D("allProbes_EB_pt", "Barrel: all Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_EE_pt = ibooker.book1D("allProbes_EE_pt", "Endcap: all Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_eta = ibooker.book1D("allProbes_eta", "All Probes Eta", etaBin_, etaMin_, etaMax_);
  h_allProbes_hp_eta = ibooker.book1D("allProbes_hp_eta", "High Pt all Probes Eta", etaBin_, etaMin_, etaMax_);
  h_allProbes_phi = ibooker.book1D("allProbes_phi", "All Probes Phi", phiBin_, phiMin_, phiMax_);

  h_allProbes_ID_pt = ibooker.book1D("allProbes_ID_pt", "All ID Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_EB_ID_pt = ibooker.book1D("allProbes_EB_ID_pt", "Barrel: all ID Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_EE_ID_pt = ibooker.book1D("allProbes_EE_ID_pt", "Endcap: all ID Probes Pt", ptBin_, ptMin_, ptMax_);
  h_allProbes_ID_nVtx = ibooker.book1D("allProbes_ID_nVtx", "All Probes (ID) nVtx", vtxBin_, vtxMin_, vtxMax_);
  h_allProbes_EB_ID_nVtx =
      ibooker.book1D("allProbes_EB_ID_nVtx", "Barrel: All Probes (ID) nVtx", vtxBin_, vtxMin_, vtxMax_);
  h_allProbes_EE_ID_nVtx =
      ibooker.book1D("allProbes_EE_ID_nVtx", "Endcap: All Probes (ID) nVtx", vtxBin_, vtxMin_, vtxMax_);

  h_passProbes_ID_pt = ibooker.book1D("passProbes_ID_pt", "ID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_ID_inner_pt =
      ibooker.book1D("passProbes_ID_inner_pt", "ID Passing Probes inner Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_ID_inner_eta =
      ibooker.book1D("passProbes_ID_inner_eta", "ID Passing Probes inner eta", etaBin_, etaMin_, etaMax_);
  h_passProbes_ID_inner_phi =
      ibooker.book1D("passProbes_ID_inner_phi", "ID Passing Probes inner phi", phiBin_, phiMin_, phiMax_);
  h_passProbes_ID_EB_pt = ibooker.book1D("passProbes_ID_EB_pt", "Barrel: ID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_ID_EE_pt = ibooker.book1D("passProbes_ID_EE_pt", "Endcap: ID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_ID_eta = ibooker.book1D("passProbes_ID_eta", "ID Passing Probes #eta", etaBin_, etaMin_, etaMax_);
  h_passProbes_ID_hp_eta =
      ibooker.book1D("passProbes_ID_hp_eta", "High Pt ID Passing Probes #eta", etaBin_, etaMin_, etaMax_);
  h_passProbes_ID_phi = ibooker.book1D("passProbes_ID_phi", "ID Passing Probes #phi", phiBin_, phiMin_, phiMax_);

  h_passProbes_detIsoID_pt =
      ibooker.book1D("passProbes_detIsoID_pt", "detIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_EB_detIsoID_pt =
      ibooker.book1D("passProbes_EB_detIsoID_pt", "Barrel: detIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_EE_detIsoID_pt =
      ibooker.book1D("passProbes_EE_detIsoID_pt", "Endcap: detIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);

  h_passProbes_pfIsoID_pt =
      ibooker.book1D("passProbes_pfIsoID_pt", "pfIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_EB_pfIsoID_pt =
      ibooker.book1D("passProbes_EB_pfIsoID_pt", "Barrel: pfIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);
  h_passProbes_EE_pfIsoID_pt =
      ibooker.book1D("passProbes_EE_pfIsoID_pt", "Endcap: pfIsoID Passing Probes Pt", ptBin_, ptMin_, ptMax_);

  h_passProbes_detIsoID_nVtx =
      ibooker.book1D("passProbes_detIsoID_nVtx", "detIsoID Passing Probes nVtx (R03)", vtxBin_, vtxMin_, vtxMax_);
  h_passProbes_pfIsoID_nVtx =
      ibooker.book1D("passProbes_pfIsoID_nVtx", "pfIsoID Passing Probes nVtx (R04)", vtxBin_, vtxMin_, vtxMax_);
  h_passProbes_EB_detIsoID_nVtx = ibooker.book1D(
      "passProbes_EB_detIsoID_nVtx", "Barrel: detIsoID Passing Probes nVtx (R03)", vtxBin_, vtxMin_, vtxMax_);
  h_passProbes_EE_detIsoID_nVtx = ibooker.book1D(
      "passProbes_EE_detIsoID_nVtx", "Endcap: detIsoID Passing Probes nVtx (R03)", vtxBin_, vtxMin_, vtxMax_);
  h_passProbes_EB_pfIsoID_nVtx = ibooker.book1D(
      "passProbes_EB_pfIsoID_nVtx", "Barrel: pfIsoID Passing Probes nVtx (R04)", vtxBin_, vtxMin_, vtxMax_);
  h_passProbes_EE_pfIsoID_nVtx = ibooker.book1D(
      "passProbes_EE_pfIsoID_nVtx", "Endcap: pfIsoID Passing Probes nVtx (R04)", vtxBin_, vtxMin_, vtxMax_);

  // Apply deltaBeta PU corrections to the PF isolation eficiencies.

  h_passProbes_pfIsodBID_pt = ibooker.book1D(
      "passProbes_pfIsodBID_pt", "pfIsoID Passing Probes Pt (deltaB PU correction)", ptBin_, ptMin_, ptMax_);
  h_passProbes_EB_pfIsodBID_pt = ibooker.book1D(
      "passProbes_EB_pfIsodBID_pt", "Barrel: pfIsoID Passing Probes Pt (deltaB PU correction)", ptBin_, ptMin_, ptMax_);
  h_passProbes_EE_pfIsodBID_pt = ibooker.book1D(
      "passProbes_EE_pfIsodBID_pt", "Endcap: pfIsoID Passing Probes Pt (deltaB PU correction)", ptBin_, ptMin_, ptMax_);
  h_passProbes_pfIsodBID_nVtx = ibooker.book1D("passProbes_pfIsodBID_nVtx",
                                               "pfIsoID Passing Probes nVtx (R04) (deltaB PU correction)",
                                               vtxBin_,
                                               vtxMin_,
                                               vtxMax_);
  h_passProbes_EB_pfIsodBID_nVtx = ibooker.book1D("passProbes_EB_pfIsodBID_nVtx",
                                                  "Barrel: pfIsoID Passing Probes nVtx (R04) (deltaB PU correction)",
                                                  vtxBin_,
                                                  vtxMin_,
                                                  vtxMax_);
  h_passProbes_EE_pfIsodBID_nVtx = ibooker.book1D("passProbes_EE_pfIsodBID_nVtx",
                                                  "Endcap: pfIsoID Passing Probes nVtx (R04) (deltaB PU correction)",
                                                  vtxBin_,
                                                  vtxMin_,
                                                  vtxMax_);

#ifdef DEBUG
  cout << "[EfficiencyAnalyzer] Parameters initialization DONE" << endl;
#endif
}

Member Data Documentation

_numPV

int EfficiencyAnalyzer::_numPV

private

Definition at line 119 of file EfficiencyAnalyzer.h.

doPVCheck_

bool EfficiencyAnalyzer::doPVCheck_

private

Definition at line 126 of file EfficiencyAnalyzer.h.

etaBin_

int EfficiencyAnalyzer::etaBin_

private

Definition at line 50 of file EfficiencyAnalyzer.h.

etaMax_

double EfficiencyAnalyzer::etaMax_

private

Definition at line 58 of file EfficiencyAnalyzer.h.

etaMin_

double EfficiencyAnalyzer::etaMin_

private

Definition at line 57 of file EfficiencyAnalyzer.h.

h_allProbes_EB_ID_nVtx

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_nVtx

private

Definition at line 108 of file EfficiencyAnalyzer.h.

h_allProbes_EB_ID_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_pt

private

Definition at line 105 of file EfficiencyAnalyzer.h.

h_allProbes_EB_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_pt

private

Definition at line 96 of file EfficiencyAnalyzer.h.

h_allProbes_EE_ID_nVtx

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_nVtx

private

Definition at line 109 of file EfficiencyAnalyzer.h.

h_allProbes_EE_ID_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_pt

private

Definition at line 106 of file EfficiencyAnalyzer.h.

h_allProbes_EE_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_pt

private

Definition at line 97 of file EfficiencyAnalyzer.h.

h_allProbes_eta

MonitorElement* EfficiencyAnalyzer::h_allProbes_eta

private

Definition at line 98 of file EfficiencyAnalyzer.h.

h_allProbes_hp_eta

MonitorElement* EfficiencyAnalyzer::h_allProbes_hp_eta

private

Definition at line 99 of file EfficiencyAnalyzer.h.

h_allProbes_ID_nVtx

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_nVtx

private

Definition at line 107 of file EfficiencyAnalyzer.h.

h_allProbes_ID_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_pt

private

Definition at line 101 of file EfficiencyAnalyzer.h.

h_allProbes_inner_eta

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_eta

private

Definition at line 103 of file EfficiencyAnalyzer.h.

h_allProbes_inner_phi

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_phi

private

Definition at line 104 of file EfficiencyAnalyzer.h.

h_allProbes_inner_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_pt

private

Definition at line 102 of file EfficiencyAnalyzer.h.

h_allProbes_phi

MonitorElement* EfficiencyAnalyzer::h_allProbes_phi

private

Definition at line 100 of file EfficiencyAnalyzer.h.

h_allProbes_pt

MonitorElement* EfficiencyAnalyzer::h_allProbes_pt

private

Definition at line 95 of file EfficiencyAnalyzer.h.

h_failProbes_ID_eta

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_eta

private

Definition at line 92 of file EfficiencyAnalyzer.h.

h_failProbes_ID_phi

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_phi

private

Definition at line 93 of file EfficiencyAnalyzer.h.

h_failProbes_ID_pt

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_pt

private

Definition at line 91 of file EfficiencyAnalyzer.h.

h_passProbes_detIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_nVtx

private

Definition at line 84 of file EfficiencyAnalyzer.h.

h_passProbes_detIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_pt

private

Definition at line 78 of file EfficiencyAnalyzer.h.

h_passProbes_EB_detIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_nVtx

private

Definition at line 86 of file EfficiencyAnalyzer.h.

h_passProbes_EB_detIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_pt

private

Definition at line 79 of file EfficiencyAnalyzer.h.

h_passProbes_EB_pfIsodBID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_nVtx

private

Definition at line 116 of file EfficiencyAnalyzer.h.

h_passProbes_EB_pfIsodBID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_pt

private

Definition at line 113 of file EfficiencyAnalyzer.h.

h_passProbes_EB_pfIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_nVtx

private

Definition at line 88 of file EfficiencyAnalyzer.h.

h_passProbes_EB_pfIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_pt

private

Definition at line 82 of file EfficiencyAnalyzer.h.

h_passProbes_EE_detIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_nVtx

private

Definition at line 87 of file EfficiencyAnalyzer.h.

h_passProbes_EE_detIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_pt

private

Definition at line 80 of file EfficiencyAnalyzer.h.

h_passProbes_EE_pfIsodBID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_nVtx

private

Definition at line 117 of file EfficiencyAnalyzer.h.

h_passProbes_EE_pfIsodBID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_pt

private

Definition at line 114 of file EfficiencyAnalyzer.h.

h_passProbes_EE_pfIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_nVtx

private

Definition at line 89 of file EfficiencyAnalyzer.h.

h_passProbes_EE_pfIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_pt

private

Definition at line 83 of file EfficiencyAnalyzer.h.

h_passProbes_ID_EB_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EB_pt

private

Definition at line 73 of file EfficiencyAnalyzer.h.

h_passProbes_ID_EE_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EE_pt

private

Definition at line 74 of file EfficiencyAnalyzer.h.

h_passProbes_ID_eta

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_eta

private

Definition at line 75 of file EfficiencyAnalyzer.h.

h_passProbes_ID_hp_eta

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_hp_eta

private

Definition at line 76 of file EfficiencyAnalyzer.h.

h_passProbes_ID_inner_eta

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_eta

private

Definition at line 71 of file EfficiencyAnalyzer.h.

h_passProbes_ID_inner_phi

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_phi

private

Definition at line 72 of file EfficiencyAnalyzer.h.

h_passProbes_ID_inner_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_pt

private

Definition at line 70 of file EfficiencyAnalyzer.h.

h_passProbes_ID_phi

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_phi

private

Definition at line 77 of file EfficiencyAnalyzer.h.

h_passProbes_ID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_pt

private

Definition at line 69 of file EfficiencyAnalyzer.h.

h_passProbes_pfIsodBID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_nVtx

private

Definition at line 115 of file EfficiencyAnalyzer.h.

h_passProbes_pfIsodBID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_pt

private

Definition at line 112 of file EfficiencyAnalyzer.h.

h_passProbes_pfIsoID_nVtx

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_nVtx

private

Definition at line 85 of file EfficiencyAnalyzer.h.

h_passProbes_pfIsoID_pt

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_pt

private

Definition at line 81 of file EfficiencyAnalyzer.h.

ID_

std::string EfficiencyAnalyzer::ID_

private

Definition at line 67 of file EfficiencyAnalyzer.h.

metname

std::string EfficiencyAnalyzer::metname

private

Definition at line 47 of file EfficiencyAnalyzer.h.

parameters

edm::ParameterSet EfficiencyAnalyzer::parameters

private

Definition at line 44 of file EfficiencyAnalyzer.h.

phiBin_

int EfficiencyAnalyzer::phiBin_

private

Definition at line 51 of file EfficiencyAnalyzer.h.

phiMax_

double EfficiencyAnalyzer::phiMax_

private

Definition at line 61 of file EfficiencyAnalyzer.h.

phiMin_

double EfficiencyAnalyzer::phiMin_

private

Definition at line 60 of file EfficiencyAnalyzer.h.

ptBin_

int EfficiencyAnalyzer::ptBin_

private

Definition at line 52 of file EfficiencyAnalyzer.h.

ptMax_

double EfficiencyAnalyzer::ptMax_

private

Definition at line 55 of file EfficiencyAnalyzer.h.

ptMin_

double EfficiencyAnalyzer::ptMin_

private

Definition at line 54 of file EfficiencyAnalyzer.h.

theBeamSpotLabel_

edm::EDGetTokenT<reco::BeamSpot> EfficiencyAnalyzer::theBeamSpotLabel_

private

Definition at line 128 of file EfficiencyAnalyzer.h.

theFolder

std::string EfficiencyAnalyzer::theFolder

private

Definition at line 130 of file EfficiencyAnalyzer.h.

theMuonCollectionLabel_

edm::EDGetTokenT<edm::View<reco::Muon> > EfficiencyAnalyzer::theMuonCollectionLabel_

private

Definition at line 122 of file EfficiencyAnalyzer.h.

theTrackCollectionLabel_

edm::EDGetTokenT<reco::TrackCollection> EfficiencyAnalyzer::theTrackCollectionLabel_

private

Definition at line 123 of file EfficiencyAnalyzer.h.

theVertexLabel_

edm::EDGetTokenT<reco::VertexCollection> EfficiencyAnalyzer::theVertexLabel_

private

Definition at line 127 of file EfficiencyAnalyzer.h.

vtxBin_

int EfficiencyAnalyzer::vtxBin_

private

Definition at line 63 of file EfficiencyAnalyzer.h.

vtxMax_

double EfficiencyAnalyzer::vtxMax_

private

Definition at line 65 of file EfficiencyAnalyzer.h.

vtxMin_

double EfficiencyAnalyzer::vtxMin_

private

Definition at line 64 of file EfficiencyAnalyzer.h.