#include <MuonPFAnalyzer.h>

Inheritance diagram for MuonPFAnalyzer:

Public Types
typedef std::vector< RecoGenPair >	RecoGenCollection

typedef std::pair< const reco::Muon , const reco::GenParticle >	RecoGenPair

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 >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	MuonPFAnalyzer (const edm::ParameterSet &)
	Constructor. More...

	~MuonPFAnalyzer () override
	Destructor. More...

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	bookHistos (DQMStore::IBooker &, const std::string &)

float	combRelIso (const reco::Muon *muon)

float	fDeltaPhi (float phi1, float phi2)

void	fillInRange (MonitorElement *plot, int nAxis, double x, double y=0)

MonitorElement *	getPlot (const std::string &group, const std::string &type)

const reco::Vertex	getPrimaryVertex (edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)

int	muonTrackType (const reco::Muon *muon, bool usePF)

void	recoToGenMatch (edm::Handle< reco::MuonCollection > &reco, edm::Handle< reco::GenParticleCollection > &gen)

void	setCodeLabels (MonitorElement *plot, int nAxis)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_

std::string	theFolder

edm::EDGetTokenT< reco::GenParticleCollection >	theGenLabel_

double	theHighPtTh

double	theIsoCut

std::vector< std::string >	theMuonKinds

std::map< std::string, std::map< std::string, MonitorElement * > >	thePlots

RecoGenCollection	theRecoGen

double	theRecoGenR

edm::EDGetTokenT< reco::MuonCollection >	theRecoLabel_

bool	theRunOnMC

edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_

Additional Inherited Members
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

Definition at line 30 of file MuonPFAnalyzer.h.

Member Typedef Documentation

◆ RecoGenCollection

typedef std::vector<RecoGenPair> MuonPFAnalyzer::RecoGenCollection

Definition at line 33 of file MuonPFAnalyzer.h.

◆ RecoGenPair

typedef std::pair<const reco::Muon *, const reco::GenParticle *> MuonPFAnalyzer::RecoGenPair

Definition at line 32 of file MuonPFAnalyzer.h.

Constructor & Destructor Documentation

◆ MuonPFAnalyzer()

MuonPFAnalyzer::MuonPFAnalyzer ( const edm::ParameterSet & pSet )

explicit

Constructor.

Definition at line 39 of file MuonPFAnalyzer.cc.

                                                        {
   LogTrace("MuonPFAnalyzer") << "[MuonPFAnalyzer] Initializing configuration from parameterset.\n";
  
   theGenLabel_ = consumes<GenParticleCollection>(pSet.getParameter<InputTag>("inputTagGenParticles"));
   theRecoLabel_ = consumes<MuonCollection>(pSet.getParameter<InputTag>("inputTagMuonReco"));
   theVertexLabel_ = consumes<VertexCollection>(pSet.getParameter<InputTag>("inputTagVertex"));
   theBeamSpotLabel_ = consumes<BeamSpot>(pSet.getParameter<InputTag>("inputTagBeamSpot"));
  
   theHighPtTh = pSet.getParameter<double>("highPtThreshold");
   theRecoGenR = pSet.getParameter<double>("recoGenDeltaR");
   theIsoCut = pSet.getParameter<double>("relCombIsoCut");
   theRunOnMC = pSet.getParameter<bool>("runOnMC");
  
   theFolder = pSet.getParameter<string>("folder");
  
   theMuonKinds.push_back("");          // all TUNEP/PF muons
   theMuonKinds.push_back("Tight");     // tight TUNEP/PF muons
   theMuonKinds.push_back("TightIso");  // tight/iso TUNEP/PF muons
 }

References edm::ParameterSet::getParameter(), and LogTrace.

◆ ~MuonPFAnalyzer()

MuonPFAnalyzer::~MuonPFAnalyzer ( )

override

Destructor.

Definition at line 59 of file MuonPFAnalyzer.cc.

59 { LogTrace("MuonPFAnalyzer") << "[MuonPFAnalyzer] Destructor called.\n"; }

References LogTrace.

Member Function Documentation

◆ analyze()

void MuonPFAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	context
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 77 of file MuonPFAnalyzer.cc.

                                                                           {
   Handle<reco::MuonCollection> muons;
   event.getByToken(theRecoLabel_, muons);
  
   Handle<GenParticleCollection> genMuons;
   event.getByToken(theGenLabel_, genMuons);
  
   Handle<BeamSpot> beamSpot;
   event.getByToken(theBeamSpotLabel_, beamSpot);
  
   Handle<VertexCollection> vertex;
   event.getByToken(theVertexLabel_, vertex);
  
   const Vertex primaryVertex = getPrimaryVertex(vertex, beamSpot);
  
   recoToGenMatch(muons, genMuons);
  
   RecoGenCollection::const_iterator recoGenIt = theRecoGen.begin();
   RecoGenCollection::const_iterator recoGenEnd = theRecoGen.end();
  
   for (; recoGenIt != recoGenEnd; ++recoGenIt) {
     const Muon *muon = recoGenIt->first;
     TrackRef tunePTrack = muon->tunePMuonBestTrack();
  
     const GenParticle *genMuon = recoGenIt->second;
  
     vector<string>::const_iterator kindIt = theMuonKinds.begin();
     vector<string>::const_iterator kindEnd = theMuonKinds.end();
  
     for (; kindIt != kindEnd; ++kindIt) {
       const string &kind = (*kindIt);
  
       if (kind.find("Tight") != string::npos && !muon::isTightMuon((*muon), primaryVertex))
         continue;
  
       if (kind.find("Iso") != string::npos && combRelIso(muon) > theIsoCut)
         continue;
  
       if (theRunOnMC && genMuon && !muon->innerTrack().isNull())  // has matched gen muon
       {
         if (!tunePTrack.isNull()) {
           string group = "TUNEP" + kind;
  
           float pt = tunePTrack->pt();
           float phi = tunePTrack->phi();
           float eta = tunePTrack->eta();
  
           float genPt = genMuon->pt();
           float genPhi = genMuon->p4().phi();
           float genEta = genMuon->p4().eta();
  
           float dPtOverPt = (pt / genPt) - 1;
  
           if (pt < theHighPtTh) {
             fillInRange(getPlot(group, "code"), 1, muonTrackType(muon, false));
             fillInRange(getPlot(group, "deltaPtOverPt"), 1, dPtOverPt);
           } else {
             fillInRange(getPlot(group, "codeHighPt"), 1, muonTrackType(muon, false));
             fillInRange(getPlot(group, "deltaPtOverPtHighPt"), 1, dPtOverPt);
           }
  
           fillInRange(getPlot(group, "deltaPt"), 1, (pt - genPt));
           fillInRange(getPlot(group, "deltaPhi"), 1, fDeltaPhi(genPhi, phi));
           fillInRange(getPlot(group, "deltaEta"), 1, genEta - eta);
         }
  
         if (muon->isPFMuon()) {
           string group = "PF" + kind;
  
           // Assumes that default in muon is PF
           float pt = muon->pt();
           float phi = muon->p4().phi();
           float eta = muon->p4().eta();
  
           float genPt = genMuon->pt();
           float genPhi = genMuon->p4().phi();
           float genEta = genMuon->p4().eta();
  
           float dPtOverPt = (pt / genPt) - 1;
  
           if (pt < theHighPtTh) {
             fillInRange(getPlot(group, "code"), 1, muonTrackType(muon, true));
             fillInRange(getPlot(group, "deltaPtOverPt"), 1, dPtOverPt);
           } else {
             fillInRange(getPlot(group, "codeHighPt"), 1, muonTrackType(muon, true));
             fillInRange(getPlot(group, "deltaPtOverPtHighPt"), 1, dPtOverPt);
           }
  
           fillInRange(getPlot(group, "deltaPt"), 1, pt - genPt);
           fillInRange(getPlot(group, "deltaPhi"), 1, fDeltaPhi(genPhi, phi));
           fillInRange(getPlot(group, "deltaEta"), 1, genEta - eta);
         }
       }
  
       if (muon->isPFMuon() && !tunePTrack.isNull() && !muon->innerTrack().isNull())  // Compare PF with TuneP + Tracker
       {                                                                              // No gen matching needed
  
         string group = "PFvsTUNEP" + kind;
  
         float pt = tunePTrack->pt();
         float phi = tunePTrack->phi();
         float eta = tunePTrack->eta();
  
         // Assumes that default in muon is PF
         float pfPt = muon->pt();
         float pfPhi = muon->p4().phi();
         float pfEta = muon->p4().eta();
         float dPtOverPt = (pfPt / pt) - 1;  // TUNEP vs PF pt used as denum.
  
         if (pt < theHighPtTh) {
           fillInRange(getPlot(group, "code"), 2, muonTrackType(muon, false), muonTrackType(muon, true));
           fillInRange(getPlot(group, "deltaPtOverPt"), 1, dPtOverPt);
         } else {
           fillInRange(getPlot(group, "codeHighPt"), 2, muonTrackType(muon, false), muonTrackType(muon, true));
           fillInRange(getPlot(group, "deltaPtOverPtHighPt"), 1, dPtOverPt);
         }
  
         fillInRange(getPlot(group, "deltaPt"), 1, pfPt - pt);
         fillInRange(getPlot(group, "deltaPhi"), 1, fDeltaPhi(pfPhi, phi));
         fillInRange(getPlot(group, "deltaEta"), 1, pfEta - eta);
  
         if (theRunOnMC && genMuon)  // has a matched gen muon
  
         {
           float genPt = genMuon->pt();
           float dPtOverPtGen = (pt / genPt) - 1;
           float dPtOverPtGenPF = (pfPt / genPt) - 1;
  
           if (pt < theHighPtTh) {
             fillInRange(getPlot(group, "deltaPtOverPtPFvsTUNEP"), 2, dPtOverPtGen, dPtOverPtGenPF);
           } else {
             fillInRange(getPlot(group, "deltaPtOverPtHighPtPFvsTUNEP"), 2, dPtOverPtGen, dPtOverPtGenPF);
           }
         }
       }
     }
   }
 }

References pwdgSkimBPark_cfi::beamSpot, PVValHelper::eta, runTauDisplay::genMuons, PFJetFilter_cfi::genPt, watchdog::group, edm::Ref< C, T, F >::isNull(), muon::isTightMuon(), CalibrationSummaryClient_cfi::kind, PDWG_BPHSkim_cff::muons, reco::LeafCandidate::p4(), BeamMonitor_cff::primaryVertex, DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), and bphysicsOniaDQM_cfi::vertex.

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 62 of file MuonPFAnalyzer.cc.

                                                                                                    {
   if (theRunOnMC) {
     bookHistos(ibooker, "PF");
     bookHistos(ibooker, "PFTight");
     bookHistos(ibooker, "PFTightIso");
     bookHistos(ibooker, "TUNEP");
     bookHistos(ibooker, "TUNEPTight");
     bookHistos(ibooker, "TUNEPTightIso");
   }
  
   bookHistos(ibooker, "PFvsTUNEP");
   bookHistos(ibooker, "PFvsTUNEPTight");
   bookHistos(ibooker, "PFvsTUNEPTightIso");
 }

References bookHistos().

◆ bookHistos()

void MuonPFAnalyzer::bookHistos	(	DQMStore::IBooker &	,
		const std::string &
	)

private

Definition at line 216 of file MuonPFAnalyzer.cc.

                                                                              {
   LogTrace("MuonPFAnalyzer") << "[MuonPFAnalyzer] Booking histos for group :" << group << "\n";
  
   ibooker.setCurrentFolder(string(theFolder) + group);
  
   bool isPFvsTUNEP = group.find("PFvsTUNEP") != string::npos;
  
   string hName;
  
   hName = "deltaPtOverPt" + group;
   thePlots[group]["deltaPtOverPt"] = ibooker.book1D(hName.c_str(), hName.c_str(), 101, -1.01, 1.01);
  
   hName = "deltaPtOverPtHighPt" + group;
   thePlots[group]["deltaPtOverPtHighPt"] = ibooker.book1D(hName.c_str(), hName.c_str(), 101, -1.01, 1.01);
  
   hName = "deltaPt" + group;
   thePlots[group]["deltaPt"] = ibooker.book1D(hName.c_str(), hName.c_str(), 201., -10.25, 10.25);
  
   hName = "deltaPhi" + group;
   thePlots[group]["deltaPhi"] = ibooker.book1D(hName.c_str(), hName.c_str(), 51., 0, .0102);
  
   hName = "deltaEta" + group;
   thePlots[group]["deltaEta"] = ibooker.book1D(hName.c_str(), hName.c_str(), 101., -.00505, .00505);
  
   if (isPFvsTUNEP) {
     hName = "code" + group;
     MonitorElement *plot = ibooker.book2D(hName.c_str(), hName.c_str(), 7, -.5, 6.5, 7, -.5, 6.5);
     thePlots[group]["code"] = plot;
     setCodeLabels(plot, 1);
     setCodeLabels(plot, 2);
  
     hName = "codeHighPt" + group;
     plot = ibooker.book2D(hName.c_str(), hName.c_str(), 7, -.5, 6.5, 7, -.5, 6.5);
     thePlots[group]["codeHighPt"] = plot;
     setCodeLabels(plot, 1);
     setCodeLabels(plot, 2);
  
     if (theRunOnMC) {
       hName = "deltaPtOverPtPFvsTUNEP" + group;
       thePlots[group]["deltaPtOverPtPFvsTUNEP"] =
           ibooker.book2D(hName.c_str(), hName.c_str(), 101, -1.01, 1.01, 101, -1.01, 1.01);
  
       hName = "deltaPtOverPtHighPtPFvsTUNEP" + group;
       thePlots[group]["deltaPtOverPtHighPtPFvsTUNEP"] =
           ibooker.book2D(hName.c_str(), hName.c_str(), 101, -1.01, 1.01, 101, -1.01, 1.01);
     }
   } else {
     hName = "code" + group;
     MonitorElement *plot = ibooker.book1D(hName.c_str(), hName.c_str(), 7, -.5, 6.5);
     thePlots[group]["code"] = plot;
     setCodeLabels(plot, 1);
  
     hName = "codeHighPt" + group;
     plot = ibooker.book1D(hName.c_str(), hName.c_str(), 7, -.5, 6.5);
     thePlots[group]["codeHighPt"] = plot;
     setCodeLabels(plot, 1);
   }
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), watchdog::group, LogTrace, plotFactory::plot, and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ combRelIso()

float MuonPFAnalyzer::combRelIso ( const reco::Muon * muon )

inlineprivate

Definition at line 292 of file MuonPFAnalyzer.cc.

                                                             {
   MuonIsolation iso = muon->isolationR03();
   float combRelIso = (iso.emEt + iso.hadEt + iso.sumPt) / muon->pt();
  
   return combRelIso;
 }

References reco::MuonIsolation::emEt, reco::MuonIsolation::hadEt, and reco::MuonIsolation::sumPt.

◆ fDeltaPhi()

float MuonPFAnalyzer::fDeltaPhi	(	float	phi1,
		float	phi2
	)

inlineprivate

Definition at line 299 of file MuonPFAnalyzer.cc.

                                                              {
   float fPhiDiff = fabs(acos(cos(phi1 - phi2)));
   return fPhiDiff;
 }

References funct::cos().

◆ fillInRange()

void MuonPFAnalyzer::fillInRange	(	MonitorElement *	plot,
		int	nAxis,
		double	x,
		double	y = `0`
	)

private

Definition at line 328 of file MuonPFAnalyzer.cc.

                                                                                     {
   TH1 *histo = plot->getTH1();
  
   TAxis *axis[2] = {nullptr, nullptr};
   axis[0] = histo->GetXaxis();
   if (nAxis == 2)
     axis[1] = histo->GetYaxis();
  
   double value[2] = {0, 0};
   value[0] = x;
   value[1] = y;
  
   for (int i = 0; i < nAxis; ++i) {
     double min = axis[i]->GetXmin();
     double max = axis[i]->GetXmax();
  
     if (value[i] <= min)
       value[i] = axis[i]->GetBinCenter(1);
  
     if (value[i] >= max)
       value[i] = axis[i]->GetBinCenter(axis[i]->GetNbins());
   }
  
   if (nAxis == 2)
     plot->Fill(value[0], value[1]);
   else
     plot->Fill(value[0]);
 }

References timingPdfMaker::histo, mps_fire::i, edm::max(), edm::min(), and plotFactory::plot.

◆ getPlot()

MuonPFAnalyzer::MonitorElement * MuonPFAnalyzer::getPlot	(	const std::string &	group,
		const std::string &	type
	)

private

Definition at line 275 of file MuonPFAnalyzer.cc.

                                                                                              {
   map<string, map<string, MonitorElement *> >::iterator groupIt = thePlots.find(group);
   if (groupIt == thePlots.end()) {
     LogTrace("MuonPFAnalyzer") << "[MuonPFAnalyzer] GROUP : " << group << " is not a valid plot group. Returning 0.\n";
     return nullptr;
   }
  
   map<string, MonitorElement *>::iterator typeIt = groupIt->second.find(type);
   if (typeIt == groupIt->second.end()) {
     LogTrace("MuonPFAnalyzer") << "[MuonPFAnalyzer] TYPE : " << type << " is not a valid type for GROUP : " << group
                                << ". Returning 0.\n";
     return nullptr;
   }
  
   return typeIt->second;
 }

References watchdog::group, and LogTrace.

◆ getPrimaryVertex()

const reco::Vertex MuonPFAnalyzer::getPrimaryVertex	(	edm::Handle< reco::VertexCollection > &	vertex,
		edm::Handle< reco::BeamSpot > &	beamSpot
	)

private

Definition at line 416 of file MuonPFAnalyzer.cc.

                                                                                                               {
   Vertex::Point posVtx;
   Vertex::Error errVtx;
  
   bool hasPrimaryVertex = false;
  
   if (vertex.isValid()) {
     vector<Vertex>::const_iterator vertexIt = vertex->begin();
     vector<Vertex>::const_iterator vertexEnd = vertex->end();
  
     for (; vertexIt != vertexEnd; ++vertexIt) {
       if (vertexIt->isValid() && !vertexIt->isFake()) {
         posVtx = vertexIt->position();
         errVtx = vertexIt->error();
         hasPrimaryVertex = true;
         break;
       }
     }
   }
  
   if (!hasPrimaryVertex) {
     LogInfo("MuonPFAnalyzer") << "[MuonPFAnalyzer] PrimaryVertex not found, use BeamSpot position instead.\n";
  
     posVtx = beamSpot->position();
     errVtx(0, 0) = beamSpot->BeamWidthX();
     errVtx(1, 1) = beamSpot->BeamWidthY();
     errVtx(2, 2) = beamSpot->sigmaZ();
   }
  
   const Vertex primaryVertex(posVtx, errVtx);
  
   return primaryVertex;
 }

References pwdgSkimBPark_cfi::beamSpot, BeamMonitor_cff::primaryVertex, and bphysicsOniaDQM_cfi::vertex.

◆ muonTrackType()

int MuonPFAnalyzer::muonTrackType	(	const reco::Muon *	muon,
		bool	usePF
	)

private

Definition at line 357 of file MuonPFAnalyzer.cc.

                                                               {
   switch (usePF ? muon->muonBestTrackType() : muon->tunePMuonBestTrackType()) {
     case Muon::InnerTrack:
       return 0;
     case Muon::OuterTrack:
       return 1;
     case Muon::CombinedTrack:
       return 2;
     case Muon::TPFMS:
       return 3;
     case Muon::Picky:
       return 4;
     case Muon::DYT:
       return 5;
     case Muon::None:
       return 6;
   }
  
   return 6;
 }

References None.

◆ recoToGenMatch()

void MuonPFAnalyzer::recoToGenMatch	(	edm::Handle< reco::MuonCollection > &	reco,
		edm::Handle< reco::GenParticleCollection > &	gen
	)

private

Definition at line 378 of file MuonPFAnalyzer.cc.

                                                                                                       {
   theRecoGen.clear();
  
   if (muons.isValid()) {
     MuonCollection::const_iterator muonIt = muons->begin();
     MuonCollection::const_iterator muonEnd = muons->end();
  
     for (; muonIt != muonEnd; ++muonIt) {
       float bestDR = 999.;
       const GenParticle *bestGen = nullptr;
  
       if (theRunOnMC && gens.isValid()) {
         GenParticleCollection::const_iterator genIt = gens->begin();
         GenParticleCollection::const_iterator genEnd = gens->end();
  
         for (; genIt != genEnd; ++genIt) {
           if (abs(genIt->pdgId()) == 13) {
             float muonPhi = muonIt->phi();
             float muonEta = muonIt->eta();
  
             float genPhi = genIt->phi();
             float genEta = genIt->eta();
  
             float dR = deltaR(muonEta, muonPhi, genEta, genPhi);
  
             if (dR < theRecoGenR && dR < bestDR) {
               bestDR = dR;
               bestGen = &(*genIt);
             }
           }
         }
       }
  
       theRecoGen.push_back(RecoGenPair(&(*muonIt), bestGen));
     }
   }
 }

References funct::abs(), PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, edm::HandleBase::isValid(), HLTObjectMonitor_cfi::muonEta, HLTObjectMonitor_cfi::muonPhi, and PDWG_BPHSkim_cff::muons.

◆ setCodeLabels()

void MuonPFAnalyzer::setCodeLabels	(	MonitorElement *	plot,
		int	nAxis
	)

private

Definition at line 304 of file MuonPFAnalyzer.cc.

                                                                   {
   TAxis *axis = nullptr;
  
   TH1 *histo = plot->getTH1();
   if (!histo)
     return;
  
   if (nAxis == 1)
     axis = histo->GetXaxis();
   else if (nAxis == 2)
     axis = histo->GetYaxis();
  
   if (!axis)
     return;
  
   axis->SetBinLabel(1, "Inner Track");
   axis->SetBinLabel(2, "Outer Track");
   axis->SetBinLabel(3, "Combined");
   axis->SetBinLabel(4, "TPFMS");
   axis->SetBinLabel(5, "Picky");
   axis->SetBinLabel(6, "DYT");
   axis->SetBinLabel(7, "None");
 }

References timingPdfMaker::histo, and plotFactory::plot.

Member Data Documentation

◆ theBeamSpotLabel_

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

private

Definition at line 76 of file MuonPFAnalyzer.h.

◆ theFolder

std::string MuonPFAnalyzer::theFolder

private

Definition at line 89 of file MuonPFAnalyzer.h.

◆ theGenLabel_

edm::EDGetTokenT<reco::GenParticleCollection> MuonPFAnalyzer::theGenLabel_

private

Definition at line 73 of file MuonPFAnalyzer.h.

◆ theHighPtTh

double MuonPFAnalyzer::theHighPtTh

private

Definition at line 83 of file MuonPFAnalyzer.h.

◆ theIsoCut

double MuonPFAnalyzer::theIsoCut

private

Definition at line 85 of file MuonPFAnalyzer.h.

◆ theMuonKinds

std::vector<std::string> MuonPFAnalyzer::theMuonKinds

private

Definition at line 78 of file MuonPFAnalyzer.h.

◆ thePlots

std::map<std::string, std::map<std::string, MonitorElement *> > MuonPFAnalyzer::thePlots

private

Definition at line 80 of file MuonPFAnalyzer.h.

◆ theRecoGen

RecoGenCollection MuonPFAnalyzer::theRecoGen

private

Definition at line 81 of file MuonPFAnalyzer.h.

◆ theRecoGenR

double MuonPFAnalyzer::theRecoGenR

private

Definition at line 84 of file MuonPFAnalyzer.h.

◆ theRecoLabel_

edm::EDGetTokenT<reco::MuonCollection> MuonPFAnalyzer::theRecoLabel_

private

Definition at line 74 of file MuonPFAnalyzer.h.

◆ theRunOnMC

bool MuonPFAnalyzer::theRunOnMC

private

Definition at line 87 of file MuonPFAnalyzer.h.

◆ theVertexLabel_

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

private

Definition at line 75 of file MuonPFAnalyzer.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ RecoGenCollection

◆ RecoGenPair

Constructor & Destructor Documentation

◆ MuonPFAnalyzer()

◆ ~MuonPFAnalyzer()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ bookHistos()

◆ combRelIso()

◆ fDeltaPhi()

◆ fillInRange()

◆ getPlot()

◆ getPrimaryVertex()

◆ muonTrackType()

◆ recoToGenMatch()

◆ setCodeLabels()

Member Data Documentation

◆ theBeamSpotLabel_

◆ theFolder

◆ theGenLabel_

◆ theHighPtTh

◆ theIsoCut

◆ theMuonKinds

◆ thePlots

◆ theRecoGen

◆ theRecoGenR

◆ theRecoLabel_

◆ theRunOnMC

◆ theVertexLabel_