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

(

const edm::ParameterSet &

pset

)

Definition at line 29 of file EfficiencyAnalyzer.cc.

References etaMax_(), etaMin_(), and ptMin_().

                                                                  {
  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 ( )

override

Definition at line 61 of file EfficiencyAnalyzer.cc.

{}

Member Function Documentation

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 reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), cms::cuda::bs, reco::MuonIsolation::emEt, relativeConstraints::error, edm::Event::getByToken(), reco::MuonIsolation::hadEt, muon::isLooseMuon(), muon::isMediumMuon(), muon::isTightMuon(), edm::HandleBase::isValid(), LogTrace, SiStripPI::max, metname, patZpeak::muons, reco::BeamSpot::position(), position, edm::Handle< T >::product(), reco::BeamSpot::sigmaZ(), reco::MuonIsolation::sumPt, tracks, findQualityFiles::v, and GoodVertex_cfg::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);
      }
    }
  }
}

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, etaMax_(), etaMin_(), ptMin_(), 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

int EfficiencyAnalyzer::_numPV

private

Definition at line 119 of file EfficiencyAnalyzer.h.

bool EfficiencyAnalyzer::doPVCheck_

private

Definition at line 126 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::etaBin_

private

Definition at line 50 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::etaMax_

private

Definition at line 58 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::etaMin_

private

Definition at line 57 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_nVtx

private

Definition at line 108 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_pt

private

Definition at line 105 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_pt

private

Definition at line 96 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_nVtx

private

Definition at line 109 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_pt

private

Definition at line 106 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_pt

private

Definition at line 97 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_eta

private

Definition at line 98 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_hp_eta

private

Definition at line 99 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_nVtx

private

Definition at line 107 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_pt

private

Definition at line 101 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_eta

private

Definition at line 103 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_phi

private

Definition at line 104 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_pt

private

Definition at line 102 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_phi

private

Definition at line 100 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_pt

private

Definition at line 95 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_eta

private

Definition at line 92 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_phi

private

Definition at line 93 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_pt

private

Definition at line 91 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_nVtx

private

Definition at line 84 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_pt

private

Definition at line 78 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_nVtx

private

Definition at line 86 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_pt

private

Definition at line 79 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_nVtx

private

Definition at line 116 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_pt

private

Definition at line 113 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_nVtx

private

Definition at line 88 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_pt

private

Definition at line 82 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_nVtx

private

Definition at line 87 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_pt

private

Definition at line 80 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_nVtx

private

Definition at line 117 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_pt

private

Definition at line 114 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_nVtx

private

Definition at line 89 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_pt

private

Definition at line 83 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EB_pt

private

Definition at line 73 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EE_pt

private

Definition at line 74 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_eta

private

Definition at line 75 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_hp_eta

private

Definition at line 76 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_eta

private

Definition at line 71 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_phi

private

Definition at line 72 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_pt

private

Definition at line 70 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_phi

private

Definition at line 77 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_pt

private

Definition at line 69 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_nVtx

private

Definition at line 115 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_pt

private

Definition at line 112 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_nVtx

private

Definition at line 85 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_pt

private

Definition at line 81 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::ID_

private

Definition at line 67 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::metname

private

Definition at line 47 of file EfficiencyAnalyzer.h.

edm::ParameterSet EfficiencyAnalyzer::parameters

private

Definition at line 44 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::phiBin_

private

Definition at line 51 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::phiMax_

private

Definition at line 61 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::phiMin_

private

Definition at line 60 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::ptBin_

private

Definition at line 52 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::ptMax_

private

Definition at line 55 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::ptMin_

private

Definition at line 54 of file EfficiencyAnalyzer.h.

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

private

Definition at line 128 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::theFolder

private

Definition at line 130 of file EfficiencyAnalyzer.h.

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

private

Definition at line 122 of file EfficiencyAnalyzer.h.

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

private

Definition at line 123 of file EfficiencyAnalyzer.h.

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

private

Definition at line 127 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::vtxBin_

private

Definition at line 63 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::vtxMax_

private

Definition at line 65 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::vtxMin_

private

Definition at line 64 of file EfficiencyAnalyzer.h.