#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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, 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 &
	)

overridevirtual

Implements DQMEDAnalyzer.

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.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation