MuonPFAnalyzer Class Reference

#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 >
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

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
	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 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.

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

                                                       {
  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
}

~MuonPFAnalyzer()

MuonPFAnalyzer::~MuonPFAnalyzer ( )

override

Destructor.

Definition at line 59 of file MuonPFAnalyzer.cc.

References LogTrace.

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

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.

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

                                                                          {
  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);
          }
        }
      }
    }
  }
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 62 of file MuonPFAnalyzer.cc.

References bookHistos().

                                                                                                   {
  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");
}

bookHistos()

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

private

Definition at line 216 of file MuonPFAnalyzer.cc.

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

                                                                             {
  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);
  }
}

combRelIso()

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

inlineprivate

Definition at line 292 of file MuonPFAnalyzer.cc.

References genparticles_cff::iso.

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

  return combRelIso;
}

fDeltaPhi()

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

inlineprivate

Definition at line 299 of file MuonPFAnalyzer.cc.

References funct::cos().

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

fillInRange()

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

private

Definition at line 328 of file MuonPFAnalyzer.cc.

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

                                                                                    {
  TH1 *histo = plot->getTH1();

  // Avoid LLVM analyzer warning
  assert(nAxis == 1 || nAxis == 2);

  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]);
}

getPlot()

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

private

Definition at line 275 of file MuonPFAnalyzer.cc.

References watchdog::group, and LogTrace.

                                                                                             {
  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;
}

getPrimaryVertex()

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

private

Definition at line 420 of file MuonPFAnalyzer.cc.

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

                                                                                                              {
  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;
}

muonTrackType()

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

private

Definition at line 361 of file MuonPFAnalyzer.cc.

References isotrackNtuplerCondor::None.

                                                              {
  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;
}

recoToGenMatch()

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

private

Definition at line 382 of file MuonPFAnalyzer.cc.

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

                                                                                                      {
  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));
    }
  }
}

setCodeLabels()

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

private

Definition at line 304 of file MuonPFAnalyzer.cc.

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

                                                                  {
  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");
}

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.