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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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_
MuonServiceProxy *	theService
edm::EDGetTokenT< reco::TrackCollection >	theTrackCollectionLabel_
edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_
int	vtxBin_
double	vtxMax_
double	vtxMin_

Detailed Description

Class EfficiencyAnalyzer

DQM monitoring for dimuon mass

Author: S.Folgueras, A. Calderon

Definition at line 36 of file EfficiencyAnalyzer.h.

Constructor & Destructor Documentation

EfficiencyAnalyzer::EfficiencyAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 29 of file EfficiencyAnalyzer.cc.

References MuonServiceProxy_cff::MuonServiceProxy.

                                                                 {
  parameters = pSet;
  
  theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));

  // 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 63 of file EfficiencyAnalyzer.cc.

                                        {
  delete theService;
}

Member Function Documentation

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

override

to be read from output as "generalTracks"

Definition at line 131 of file EfficiencyAnalyzer.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), taus_cff::chargedIso, reco::MuonIsolation::emEt, relativeConstraints::error, edm::Event::getByToken(), reco::MuonIsolation::hadEt, muon::isLooseMuon(), muon::isMediumMuon(), muon::isTightMuon(), edm::HandleBase::isValid(), LogTrace, SiStripPI::max, metname, extraflags_cff::muons, taus_cff::neutralIso, gedPhotonSequence_cff::photonIso, reco::BeamSpot::position(), position, edm::Handle< T >::product(), muons2muons_cfi::pu, reco::BeamSpot::sigmaZ(), reco::MuonIsolation::sumPt, l1t::tracks, findQualityFiles::v, and particleFlowSuperClusterECAL_cfi::vertexCollection.

                                                                                   {

  LogTrace(metname)<<"[EfficiencyAnalyzer] Analyze the mu in different eta regions";
  theService->update(iSetup);
  // ==========================================================
  // 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 &    )

override

Definition at line 67 of file EfficiencyAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), gather_cfg::cout, and DQMStore::IBooker::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 127 of file EfficiencyAnalyzer.h.

bool EfficiencyAnalyzer::doPVCheck_

private

Definition at line 134 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::etaBin_

private

Definition at line 57 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::etaMax_

private

Definition at line 65 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::etaMin_

private

Definition at line 64 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_nVtx

private

Definition at line 115 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_ID_pt

private

Definition at line 112 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EB_pt

private

Definition at line 103 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_nVtx

private

Definition at line 116 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_ID_pt

private

Definition at line 113 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_EE_pt

private

Definition at line 104 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_eta

private

Definition at line 105 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_hp_eta

private

Definition at line 106 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_nVtx

private

Definition at line 114 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_ID_pt

private

Definition at line 108 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_eta

private

Definition at line 110 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_phi

private

Definition at line 111 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_inner_pt

private

Definition at line 109 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_phi

private

Definition at line 107 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_allProbes_pt

private

Definition at line 102 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_eta

private

Definition at line 99 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_phi

private

Definition at line 100 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_failProbes_ID_pt

private

Definition at line 98 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_nVtx

private

Definition at line 91 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_detIsoID_pt

private

Definition at line 85 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_nVtx

private

Definition at line 93 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_detIsoID_pt

private

Definition at line 86 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_nVtx

private

Definition at line 124 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsodBID_pt

private

Definition at line 121 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_nVtx

private

Definition at line 95 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EB_pfIsoID_pt

private

Definition at line 89 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_nVtx

private

Definition at line 94 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_detIsoID_pt

private

Definition at line 87 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_nVtx

private

Definition at line 125 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsodBID_pt

private

Definition at line 122 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_nVtx

private

Definition at line 96 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_EE_pfIsoID_pt

private

Definition at line 90 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EB_pt

private

Definition at line 80 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_EE_pt

private

Definition at line 81 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_eta

private

Definition at line 82 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_hp_eta

private

Definition at line 83 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_eta

private

Definition at line 78 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_phi

private

Definition at line 79 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_inner_pt

private

Definition at line 77 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_phi

private

Definition at line 84 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_ID_pt

private

Definition at line 76 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_nVtx

private

Definition at line 123 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsodBID_pt

private

Definition at line 120 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_nVtx

private

Definition at line 92 of file EfficiencyAnalyzer.h.

MonitorElement* EfficiencyAnalyzer::h_passProbes_pfIsoID_pt

private

Definition at line 88 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::ID_

private

Definition at line 74 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::metname

private

Definition at line 54 of file EfficiencyAnalyzer.h.

edm::ParameterSet EfficiencyAnalyzer::parameters

private

Definition at line 50 of file EfficiencyAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int EfficiencyAnalyzer::phiBin_

private

Definition at line 58 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::phiMax_

private

Definition at line 68 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::phiMin_

private

Definition at line 67 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::ptBin_

private

Definition at line 59 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::ptMax_

private

Definition at line 62 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::ptMin_

private

Definition at line 61 of file EfficiencyAnalyzer.h.

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

private

Definition at line 136 of file EfficiencyAnalyzer.h.

std::string EfficiencyAnalyzer::theFolder

private

Definition at line 138 of file EfficiencyAnalyzer.h.

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

private

Definition at line 130 of file EfficiencyAnalyzer.h.

MuonServiceProxy* EfficiencyAnalyzer::theService

private

Definition at line 51 of file EfficiencyAnalyzer.h.

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

private

Definition at line 131 of file EfficiencyAnalyzer.h.

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

private

Definition at line 135 of file EfficiencyAnalyzer.h.

int EfficiencyAnalyzer::vtxBin_

private

Definition at line 70 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::vtxMax_

private

Definition at line 72 of file EfficiencyAnalyzer.h.

double EfficiencyAnalyzer::vtxMin_

private

Definition at line 71 of file EfficiencyAnalyzer.h.