L1TPhase2MuonOffline Class Reference

#include <L1TPhase2MuonOffline.h>

Inheritance diagram for L1TPhase2MuonOffline:

Public Types
enum	EffType { kEffPt, kEffPhi, kEffEta, kEffTypes }
enum	EtaRegion {   kEtaRegionAll, kEtaRegionBmtf, kEtaRegionOmtf, kEtaRegionEmtf,   kNEtaRegions }
enum	MuType { kSAMuon, kTkMuon, kNMuTypes }
enum	QualLevel { kQualOpen, kQualDouble, kQualSingle, kNQualLevels }
enum	ResType {   kResPt, kRes1OverPt, kResQOverPt, kResPhi,   kResEta, kResCh, kNResTypes }
enum	VarType {   kPt, kEta, kPhi, kIso,   kQual, kZ0, kD0, kNVarTypes }
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
	L1TPhase2MuonOffline (const edm::ParameterSet &ps)
	~L1TPhase2MuonOffline () override=default
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
	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

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	bookControlHistos (DQMStore::IBooker &, MuType type)
void	bookEfficiencyHistos (DQMStore::IBooker &ibooker, MuType type)
void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &iSetup) override
void	bookResolutionHistos (DQMStore::IBooker &ibooker, MuType type)
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override
void	fillControlHistos ()
void	fillEfficiencyHistos ()
void	fillResolutionHistos ()
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
std::vector< float >	getHistBinsEff (EffType eff)
std::tuple< int, double, double >	getHistBinsRes (ResType res)
void	getMuonGmtPairs (edm::Handle< l1t::MuonBxCollection > &gmtCands)
void	matchMuonsToGen (std::vector< const reco::GenParticle *> genmus)

Private Attributes
MonitorElement *	controlHistos_ [kNMuTypes][kNVarTypes]
std::vector< std::pair< int, QualLevel > >	cuts_
const std::vector< edm::ParameterSet >	cutsVPSet_
MonitorElement *	efficiencyDen_ [kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes]
MonitorElement *	efficiencyNum_ [kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes]
const std::vector< EffType >	effTypes_
const std::vector< double >	effVsEtaBins_
const std::vector< double >	effVsPhiBins_
const std::vector< double >	effVsPtBins_
std::map< EtaRegion, std::string >	etaNames_
const std::vector< EtaRegion >	etaRegions_
edm::Handle< std::vector< reco::GenParticle > >	genparticles_
edm::EDGetTokenT< std::vector< reco::GenParticle > >	genParticleToken_
edm::EDGetTokenT< l1t::SAMuonCollection >	gmtMuonToken_
edm::Handle< l1t::SAMuonCollection >	gmtSAMuon_
std::vector< GenMuonGMTPair >	gmtSAMuonPairs_
edm::Handle< l1t::TrackerMuonCollection >	gmtTkMuon_
std::vector< GenMuonGMTPair >	gmtTkMuonPairs_
edm::EDGetTokenT< l1t::TrackerMuonCollection >	gmtTkMuonToken_
const std::string	histFolder_
const float	lsb_d0 = Phase2L1GMT::LSBSAd0
const float	lsb_eta = Phase2L1GMT::LSBeta
const float	lsb_phi = Phase2L1GMT::LSBphi
const float	lsb_pt = Phase2L1GMT::LSBpt
const float	lsb_z0 = Phase2L1GMT::LSBSAz0
const double	maxGmtMuonDR_
std::map< MuType, std::string >	muonNames_
const std::vector< MuType >	muonTypes_
const std::vector< QualLevel >	qualLevels_
std::map< QualLevel, std::string >	qualNames_
std::map< ResType, std::string >	resLabels_
std::map< ResType, std::string >	resNames_
MonitorElement *	resolutionHistos_ [kNMuTypes][kNEtaRegions][kNQualLevels][kNResTypes]
const std::vector< ResType >	resTypes_

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 43 of file L1TPhase2MuonOffline.h.

Member Enumeration Documentation

EffType

enum L1TPhase2MuonOffline::EffType

Enumerator
kEffPt
kEffPhi
kEffEta
kEffTypes

Definition at line 50 of file L1TPhase2MuonOffline.h.

{ kEffPt, kEffPhi, kEffEta, kEffTypes };

EtaRegion

enum L1TPhase2MuonOffline::EtaRegion

Enumerator
kEtaRegionAll
kEtaRegionBmtf
kEtaRegionOmtf
kEtaRegionEmtf
kNEtaRegions

Definition at line 52 of file L1TPhase2MuonOffline.h.

{ kEtaRegionAll, kEtaRegionBmtf, kEtaRegionOmtf, kEtaRegionEmtf, kNEtaRegions };

MuType

enum L1TPhase2MuonOffline::MuType

Enumerator
kSAMuon
kTkMuon
kNMuTypes

Definition at line 48 of file L1TPhase2MuonOffline.h.

{ kSAMuon, kTkMuon, kNMuTypes };

QualLevel

enum L1TPhase2MuonOffline::QualLevel

Enumerator
kQualOpen
kQualDouble
kQualSingle
kNQualLevels

Definition at line 53 of file L1TPhase2MuonOffline.h.

{ kQualOpen, kQualDouble, kQualSingle, kNQualLevels };

ResType

enum L1TPhase2MuonOffline::ResType

Enumerator
kResPt
kRes1OverPt
kResQOverPt
kResPhi
kResEta
kResCh
kNResTypes

Definition at line 51 of file L1TPhase2MuonOffline.h.

{ kResPt, kRes1OverPt, kResQOverPt, kResPhi, kResEta, kResCh, kNResTypes };

VarType

enum L1TPhase2MuonOffline::VarType

Enumerator
kPt
kEta
kPhi
kIso
kQual
kZ0
kD0
kNVarTypes

Definition at line 49 of file L1TPhase2MuonOffline.h.

{ kPt, kEta, kPhi, kIso, kQual, kZ0, kD0, kNVarTypes };

Constructor & Destructor Documentation

L1TPhase2MuonOffline()

L1TPhase2MuonOffline::L1TPhase2MuonOffline

(

const edm::ParameterSet &

ps

)

Definition at line 94 of file L1TPhase2MuonOffline.cc.

References kSAMuon, and kTkMuon.

    : gmtMuonToken_(consumes<l1t::SAMuonCollection>(ps.getParameter<edm::InputTag>("gmtMuonToken"))),
      gmtTkMuonToken_(consumes<l1t::TrackerMuonCollection>(ps.getParameter<edm::InputTag>("gmtTkMuonToken"))),
      genParticleToken_(
          consumes<std::vector<reco::GenParticle>>(ps.getUntrackedParameter<edm::InputTag>("genParticlesInputTag"))),
      muonTypes_({kSAMuon, kTkMuon}),
      effTypes_({kEffPt, kEffPhi, kEffEta}),
      resTypes_({kResPt, kResQOverPt, kResPhi, kResEta}),
      etaRegions_({kEtaRegionAll, kEtaRegionBmtf, kEtaRegionOmtf, kEtaRegionEmtf}),
      qualLevels_({kQualOpen, kQualDouble, kQualSingle}),
      resNames_({{kResPt, "pt"},
                 {kRes1OverPt, "1overpt"},
                 {kResQOverPt, "qoverpt"},
                 {kResPhi, "phi"},
                 {kResEta, "eta"},
                 {kResCh, "charge"}}),
      resLabels_({{kResPt, "(p_{T}^{L1} - p_{T}^{reco})/p_{T}^{reco}"},
                  {kRes1OverPt, "(p_{T}^{reco} - p_{T}^{L1})/p_{T}^{L1}"},
                  {kResQOverPt, "(q^{L1}*q^{reco}*p_{T}^{reco} - p_{T}^{L1})/p_{T}^{L1}"},
                  {kResPhi, "#phi_{L1} - #phi_{reco}"},
                  {kResEta, "#eta_{L1} - #eta_{reco}"},
                  {kResCh, "charge^{L1} - charge^{reco}"}}),
      etaNames_({{kEtaRegionAll, "etaMin0_etaMax2p4"},
                 {kEtaRegionBmtf, "etaMin0_etaMax0p83"},
                 {kEtaRegionOmtf, "etaMin0p83_etaMax1p24"},
                 {kEtaRegionEmtf, "etaMin1p24_etaMax2p4"}}),
      qualNames_({{kQualOpen, "qualOpen"}, {kQualDouble, "qualDouble"}, {kQualSingle, "qualSingle"}}),
      muonNames_({{kSAMuon, "SAMuon"}, {kTkMuon, "TkMuon"}}),
      histFolder_(ps.getUntrackedParameter<string>("histFolder")),
      cutsVPSet_(ps.getUntrackedParameter<std::vector<edm::ParameterSet>>("cuts")),
      effVsPtBins_(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPtBins")),
      effVsPhiBins_(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPhiBins")),
      effVsEtaBins_(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsEtaBins")),
      maxGmtMuonDR_(ps.getUntrackedParameter<double>("maxDR")) {
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::L1TPhase2MuonOffline()" << endl;

  for (const auto& c : cutsVPSet_) {
    const auto qCut = c.getUntrackedParameter<int>("qualCut");
    QualLevel qLevel = kQualOpen;
    if (qCut > 11) {
      qLevel = kQualSingle;
    } else if (qCut > 7) {
      qLevel = kQualDouble;
    } else if (qCut > 3) {
      qLevel = kQualOpen;
    }
    cuts_.emplace_back(std::make_pair(c.getUntrackedParameter<int>("ptCut"), qLevel));
  }
}

~L1TPhase2MuonOffline()

L1TPhase2MuonOffline::~L1TPhase2MuonOffline ( )

overridedefault

Member Function Documentation

analyze()

void L1TPhase2MuonOffline::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 164 of file L1TPhase2MuonOffline.cc.

References funct::abs(), fillControlHistos(), fillEfficiencyHistos(), fillResolutionHistos(), genparticles_, genParticleToken_, gmtMuonToken_, gmtSAMuon_, gmtTkMuon_, gmtTkMuonToken_, iEvent, and matchMuonsToGen().

                                                                                    {
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::analyze() " << endl;

  // COLLECT GEN MUONS
  iEvent.getByToken(genParticleToken_, genparticles_);

  std::vector<const reco::GenParticle*> genmus;
  for (const reco::GenParticle& gen : *genparticles_) {
    if (std::abs(gen.pdgId()) != 13)
      continue;
    genmus.push_back(&gen);
  }
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::analyze() N of genmus: " << genmus.size() << endl;

  // Collect both muon collection:
  iEvent.getByToken(gmtMuonToken_, gmtSAMuon_);
  iEvent.getByToken(gmtTkMuonToken_, gmtTkMuon_);

  // Fill Control histograms
  edm::LogInfo("L1TPhase2MuonOffline") << "Fill Control histograms for GMT Muons" << endl;
  fillControlHistos();

  // Match each muon to a gen muon, if possible.
  if (genmus.empty())
    return;
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::analyze() calling matchMuonsToGen() " << endl;
  matchMuonsToGen(genmus);

  // Fill efficiency and resolution once, matching has been done...
  fillEfficiencyHistos();
  fillResolutionHistos();
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::analyze() Computation finished" << endl;
}

bookControlHistos()

void L1TPhase2MuonOffline::bookControlHistos ( DQMStore::IBooker &  ibooker, MuType  type  )

protected

Definition at line 199 of file L1TPhase2MuonOffline.cc.

References dqm::implementation::IBooker::book1D(), controlHistos_, histFolder_, kD0, kEta, kIso, kPhi, kPt, kQual, kZ0, muonNames_, and dqm::implementation::NavigatorBase::setCurrentFolder().

Referenced by bookHistograms().

                                                                                    {
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::bookControlHistos()" << endl;

  ibooker.setCurrentFolder(histFolder_ + "/" + muonNames_[mutype] + "/control_variables");

  controlHistos_[mutype][kPt] = ibooker.book1D(muonNames_[mutype] + "Pt", "MuonPt; p_{T}", 50, 0., 100.);
  controlHistos_[mutype][kPhi] = ibooker.book1D(muonNames_[mutype] + "Phi", "MuonPhi; #phi", 66, -3.3, 3.3);
  controlHistos_[mutype][kEta] = ibooker.book1D(muonNames_[mutype] + "Eta", "MuonEta; #eta", 50, -2.5, 2.5);
  controlHistos_[mutype][kIso] = ibooker.book1D(muonNames_[mutype] + "Iso", "MuonIso; RelIso", 50, 0, 1.0);
  controlHistos_[mutype][kQual] = ibooker.book1D(muonNames_[mutype] + "Qual", "MuonQual; Quality", 15, 0.5, 15.5);
  controlHistos_[mutype][kZ0] = ibooker.book1D(muonNames_[mutype] + "Z0", "MuonZ0; Z_{0}", 50, 0, 50.0);
  controlHistos_[mutype][kD0] = ibooker.book1D(muonNames_[mutype] + "D0", "MuonD0; D_{0}", 50, 0, 200.);
}

bookEfficiencyHistos()

void L1TPhase2MuonOffline::bookEfficiencyHistos ( DQMStore::IBooker &  ibooker, MuType  type  )

protected

Definition at line 213 of file L1TPhase2MuonOffline.cc.

References dqm::implementation::IBooker::book1D(), efficiencyDen_, efficiencyNum_, PVValHelper::eta, etaNames_, etaRegions_, getHistBinsEff(), histFolder_, kEffEta, kEffPhi, kEffPt, muonNames_, submitPVResolutionJobs::q, qualLevels_, qualNames_, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms().

                                                                                       {
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::bookEfficiencyHistos()" << endl;

  ibooker.setCurrentFolder(histFolder_ + "/" + muonNames_[mutype] + "/nums_and_dens");

  std::string histoname = "";
  for (const auto eta : etaRegions_) {
    for (const auto q : qualLevels_) {
      histoname = "Eff_" + muonNames_[mutype] + "_" + etaNames_[eta] + "_" + qualNames_[q];

      auto histBins = getHistBinsEff(kEffPt);
      efficiencyNum_[mutype][eta][q][kEffPt] =
          ibooker.book1D(histoname + "_Pt_Num", "MuonPt; p_{T} ;", histBins.size() - 1, &histBins[0]);
      efficiencyDen_[mutype][eta][q][kEffPt] =
          ibooker.book1D(histoname + "_Pt_Den", "MuonPt; p_{T} ;", histBins.size() - 1, &histBins[0]);

      histBins = getHistBinsEff(kEffEta);
      efficiencyNum_[mutype][eta][q][kEffEta] =
          ibooker.book1D(histoname + "_Eta_Num", "MuonEta; #eta ;", histBins.size() - 1, &histBins[0]);
      efficiencyDen_[mutype][eta][q][kEffEta] =
          ibooker.book1D(histoname + "_Eta_Den", "MuonEta; #eta ;", histBins.size() - 1, &histBins[0]);

      histBins = getHistBinsEff(kEffPhi);
      efficiencyNum_[mutype][eta][q][kEffPhi] =
          ibooker.book1D(histoname + "_Phi_Num", "MuonPhi; #phi ;", histBins.size() - 1, &histBins[0]);
      efficiencyDen_[mutype][eta][q][kEffPhi] =
          ibooker.book1D(histoname + "_Phi_Den", "MuonPhi; #phi ;", histBins.size() - 1, &histBins[0]);
    }
  }
}

bookHistograms()

void L1TPhase2MuonOffline::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  run, const edm::EventSetup &  iSetup  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 150 of file L1TPhase2MuonOffline.cc.

References bookControlHistos(), bookEfficiencyHistos(), bookResolutionHistos(), and muonTypes_.

                                                                       {
  edm::LogInfo("L1TPhase2MuonOFfline") << "L1TPhase2MuonOffline::bookHistograms" << endl;

  //book histos
  for (const auto mutype : muonTypes_) {
    bookControlHistos(ibooker, mutype);
    bookEfficiencyHistos(ibooker, mutype);
    bookResolutionHistos(ibooker, mutype);
  }
}

bookResolutionHistos()

void L1TPhase2MuonOffline::bookResolutionHistos ( DQMStore::IBooker &  ibooker, MuType  type  )

protected

Definition at line 244 of file L1TPhase2MuonOffline.cc.

References dqm::implementation::IBooker::book1D(), PVValHelper::eta, etaNames_, etaRegions_, getHistBinsRes(), histFolder_, muonNames_, LaserClient_cfi::nbins, submitPVResolutionJobs::q, qualLevels_, qualNames_, resLabels_, resNames_, resolutionHistos_, resTypes_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, ALCARECOEcalPhiSym_cff::var, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

Referenced by bookHistograms().

                                                                                       {
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::bookResolutionHistos()" << endl;

  ibooker.setCurrentFolder(histFolder_ + "/" + muonNames_[mutype] + "/resolution");
  std::string histoname = "";
  for (const auto eta : etaRegions_) {
    for (const auto q : qualLevels_) {
      for (const auto var : resTypes_) {
        histoname = "Res_" + muonNames_[mutype] + "_" + etaNames_[eta] + "_" + qualNames_[q] + "_" + resNames_[var];
        auto nbins = std::get<0>(getHistBinsRes(var));
        auto xmin = std::get<1>(getHistBinsRes(var));
        auto xmax = std::get<2>(getHistBinsRes(var));
        resolutionHistos_[mutype][eta][q][var] =
            ibooker.book1D(histoname, resNames_[var] + ";" + resLabels_[var], nbins, xmin, xmax);
      }
    }
  }
}

dqmBeginRun()

void L1TPhase2MuonOffline::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  iSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 145 of file L1TPhase2MuonOffline.cc.

                                                                                     {
  edm::LogInfo("L1TPhase2MuonOFfline") << "L1TPhase2MuonOffline::dqmBeginRun" << endl;
}

fillControlHistos()

void L1TPhase2MuonOffline::fillControlHistos ( )

protected

Definition at line 264 of file L1TPhase2MuonOffline.cc.

References controlHistos_, dqm::impl::MonitorElement::Fill(), gmtSAMuon_, gmtTkMuon_, kD0, kEta, kIso, kPhi, kPt, kQual, kSAMuon, kTkMuon, kZ0, lsb_d0, lsb_eta, lsb_phi, lsb_pt, and lsb_z0.

Referenced by analyze().

                                             {
  for (auto& muIt : *gmtSAMuon_) {
    controlHistos_[kSAMuon][kPt]->Fill(lsb_pt * muIt.hwPt());
    controlHistos_[kSAMuon][kPhi]->Fill(lsb_phi * muIt.hwPhi());
    controlHistos_[kSAMuon][kEta]->Fill(lsb_eta * muIt.hwEta());
    controlHistos_[kSAMuon][kIso]->Fill(muIt.hwIso());
    controlHistos_[kSAMuon][kQual]->Fill(muIt.hwQual());
    controlHistos_[kSAMuon][kZ0]->Fill(lsb_z0 * muIt.hwZ0());
    controlHistos_[kSAMuon][kD0]->Fill(lsb_d0 * muIt.hwD0());
  }

  for (auto& muIt : *gmtTkMuon_) {
    controlHistos_[kTkMuon][kPt]->Fill(lsb_pt * muIt.hwPt());
    controlHistos_[kTkMuon][kPhi]->Fill(lsb_phi * muIt.hwPhi());
    controlHistos_[kTkMuon][kEta]->Fill(lsb_eta * muIt.hwEta());
    controlHistos_[kTkMuon][kIso]->Fill(muIt.hwIso());
    controlHistos_[kTkMuon][kQual]->Fill(muIt.hwQual());
    controlHistos_[kTkMuon][kZ0]->Fill(lsb_z0 * muIt.hwZ0());
    controlHistos_[kTkMuon][kD0]->Fill(lsb_d0 * muIt.hwD0());
  }
}

fillEfficiencyHistos()

void L1TPhase2MuonOffline::fillEfficiencyHistos ( )

protected

FOR TK MUONS

Definition at line 286 of file L1TPhase2MuonOffline.cc.

References DMR_cfg::cut, cuts_, efficiencyDen_, efficiencyNum_, effTypes_, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), gmtSAMuonPairs_, gmtTkMuonPairs_, kEffPt, kSAMuon, kTkMuon, submitPVResolutionJobs::q, and ALCARECOEcalPhiSym_cff::var.

Referenced by analyze().

                                                {
  for (const auto& muIt : gmtSAMuonPairs_) {
    auto eta = muIt.etaRegion();
    for (const auto var : effTypes_) {
      auto varToFill = muIt.getVar(var);
      for (const auto& cut : cuts_) {
        const auto q = cut.second;
        efficiencyDen_[kSAMuon][eta][q][var]->Fill(varToFill);
        if (muIt.gmtPt() < 0)
          continue;  // there is not an assciated gmt muon
        if (muIt.gmtQual() < q * 4)
          continue;  //quality requirements
        const auto gmtPtCut = cut.first;
        if (var != kEffPt && muIt.gmtPt() < gmtPtCut)
          continue;  // pt requirement
        efficiencyNum_[kSAMuon][eta][q][var]->Fill(varToFill);
      }
    }
  }

  for (const auto& muIt : gmtTkMuonPairs_) {
    auto eta = muIt.etaRegion();
    for (const auto var : effTypes_) {
      auto varToFill = muIt.getVar(var);
      for (const auto& cut : cuts_) {
        const auto q = cut.second;
        efficiencyDen_[kTkMuon][eta][q][var]->Fill(varToFill);
        if (muIt.gmtPt() < 0)
          continue;  // there is not an assciated gmt muon
        if (muIt.gmtQual() < q * 4)
          continue;  //quality requirements
        const auto gmtPtCut = cut.first;
        if (var != kEffPt && muIt.gmtPt() < gmtPtCut)
          continue;  // pt requirement
        efficiencyNum_[kTkMuon][eta][q][var]->Fill(varToFill);
      }
    }
  }
}

fillResolutionHistos()

void L1TPhase2MuonOffline::fillResolutionHistos ( )

protected

Definition at line 327 of file L1TPhase2MuonOffline.cc.

References PVValHelper::eta, dqm::impl::MonitorElement::Fill(), gmtSAMuonPairs_, gmtTkMuonPairs_, kSAMuon, kTkMuon, submitPVResolutionJobs::q, qualLevels_, resolutionHistos_, resTypes_, and ALCARECOEcalPhiSym_cff::var.

Referenced by analyze().

                                                {
  for (const auto& muIt : gmtSAMuonPairs_) {
    if (muIt.gmtPt() < 0)
      continue;

    auto eta = muIt.etaRegion();
    for (const auto q : qualLevels_) {
      if (muIt.gmtQual() < q * 4)
        continue;
      for (const auto var : resTypes_) {
        auto varToFill = muIt.getDeltaVar(var);

        resolutionHistos_[kSAMuon][eta][q][var]->Fill(varToFill);
      }
    }
  }

  for (const auto& muIt : gmtTkMuonPairs_) {
    if (muIt.gmtPt() < 0)
      continue;

    auto eta = muIt.etaRegion();
    for (const auto q : qualLevels_) {
      if (muIt.gmtQual() < q * 4)
        continue;
      for (const auto var : resTypes_) {
        auto varToFill = muIt.getDeltaVar(var);

        resolutionHistos_[kTkMuon][eta][q][var]->Fill(varToFill);
      }
    }
  }
}

getHistBinsEff()

std::vector< float > L1TPhase2MuonOffline::getHistBinsEff ( EffType  eff )

private

Definition at line 401 of file L1TPhase2MuonOffline.cc.

References L1TMuonDQMOffline_cfi::effVsEtaBins, effVsEtaBins_, L1TMuonDQMOffline_cfi::effVsPhiBins, effVsPhiBins_, L1TMuonDQMOffline_cfi::effVsPtBins, effVsPtBins_, kEffEta, kEffPhi, and kEffPt.

Referenced by bookEfficiencyHistos().

                                                                 {
  if (eff == kEffPt) {
    std::vector<float> effVsPtBins(effVsPtBins_.begin(), effVsPtBins_.end());
    return effVsPtBins;
  }
  if (eff == kEffPhi) {
    std::vector<float> effVsPhiBins(effVsPhiBins_.begin(), effVsPhiBins_.end());
    return effVsPhiBins;
  }
  if (eff == kEffEta) {
    std::vector<float> effVsEtaBins(effVsEtaBins_.begin(), effVsEtaBins_.end());
    return effVsEtaBins;
  }
  return {0., 1.};
}

getHistBinsRes()

std::tuple< int, double, double > L1TPhase2MuonOffline::getHistBinsRes ( ResType  res )

private

Definition at line 417 of file L1TPhase2MuonOffline.cc.

References kRes1OverPt, kResCh, kResEta, kResPhi, kResPt, and kResQOverPt.

Referenced by bookResolutionHistos().

                                                                              {
  if (res == kResPt)
    return {50, -2., 2.};
  if (res == kRes1OverPt)
    return {50, -2., 2.};
  if (res == kResQOverPt)
    return {50, -2., 2.};
  if (res == kResPhi)
    return {96, -0.2, 0.2};
  if (res == kResEta)
    return {100, -0.1, 0.1};
  if (res == kResCh)
    return {5, -2, 3};
  return {1, 0, 1};
}

getMuonGmtPairs()

void L1TPhase2MuonOffline::getMuonGmtPairs ( edm::Handle< l1t::MuonBxCollection > &  gmtCands )

private

matchMuonsToGen()

void L1TPhase2MuonOffline::matchMuonsToGen ( std::vector< const reco::GenParticle *>  genmus )

private

Definition at line 362 of file L1TPhase2MuonOffline.cc.

References GenMuonGMTPair::dR2(), gmtSAMuon_, gmtSAMuonPairs_, gmtTkMuon_, gmtTkMuonPairs_, and maxGmtMuonDR_.

Referenced by analyze().

                                                                                   {
  gmtSAMuonPairs_.clear();
  gmtTkMuonPairs_.clear();

  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::matchMuonsToGen() " << endl;

  for (const reco::GenParticle* gen : genmus) {
    edm::LogInfo("L1TPhase2MuonOffline") << "Looping on genmus: " << gen << endl;
    GenMuonGMTPair pairBestCand(&(*gen), nullptr);
    float dr2Best = maxGmtMuonDR_ * maxGmtMuonDR_;
    for (auto& muIt : *gmtSAMuon_) {
      GenMuonGMTPair pairTmpCand(&(*gen), &(muIt));
      float dr2Tmp = pairTmpCand.dR2();
      if (dr2Tmp < dr2Best) {
        dr2Best = dr2Tmp;
        pairBestCand = pairTmpCand;
      }
    }
    gmtSAMuonPairs_.emplace_back(pairBestCand);

    GenMuonGMTPair pairBestCand2(&(*gen), nullptr);
    dr2Best = maxGmtMuonDR_ * maxGmtMuonDR_;
    for (auto& tkmuIt : *gmtTkMuon_) {
      GenMuonGMTPair pairTmpCand(&(*gen), &(tkmuIt));
      float dr2Tmp = pairTmpCand.dR2();
      if (dr2Tmp < dr2Best) {
        dr2Best = dr2Tmp;
        pairBestCand2 = pairTmpCand;
      }
    }
    gmtTkMuonPairs_.emplace_back(pairBestCand2);
  }
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::matchMuonsToGen() gmtSAMuons: "
                                       << gmtSAMuonPairs_.size() << endl;
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::matchMuonsToGen() gmtTkMuons: "
                                       << gmtTkMuonPairs_.size() << endl;
  edm::LogInfo("L1TPhase2MuonOffline") << "L1TPhase2MuonOffline::matchMuonsToGen() END " << endl;
}

Member Data Documentation

controlHistos_

MonitorElement* L1TPhase2MuonOffline::controlHistos_[kNMuTypes][kNVarTypes]

private

Definition at line 119 of file L1TPhase2MuonOffline.h.

Referenced by bookControlHistos(), and fillControlHistos().

cuts_

std::vector<std::pair<int, QualLevel> > L1TPhase2MuonOffline::cuts_

private

Definition at line 124 of file L1TPhase2MuonOffline.h.

Referenced by fillEfficiencyHistos().

cutsVPSet_

const std::vector<edm::ParameterSet> L1TPhase2MuonOffline::cutsVPSet_

private

Definition at line 103 of file L1TPhase2MuonOffline.h.

efficiencyDen_

MonitorElement* L1TPhase2MuonOffline::efficiencyDen_[kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes]

private

Definition at line 117 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), and fillEfficiencyHistos().

efficiencyNum_

MonitorElement* L1TPhase2MuonOffline::efficiencyNum_[kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes]

private

Definition at line 116 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), and fillEfficiencyHistos().

effTypes_

const std::vector<EffType> L1TPhase2MuonOffline::effTypes_

private

Definition at line 86 of file L1TPhase2MuonOffline.h.

Referenced by fillEfficiencyHistos().

effVsEtaBins_

const std::vector<double> L1TPhase2MuonOffline::effVsEtaBins_

private

Definition at line 107 of file L1TPhase2MuonOffline.h.

Referenced by getHistBinsEff().

effVsPhiBins_

const std::vector<double> L1TPhase2MuonOffline::effVsPhiBins_

private

Definition at line 106 of file L1TPhase2MuonOffline.h.

Referenced by getHistBinsEff().

effVsPtBins_

const std::vector<double> L1TPhase2MuonOffline::effVsPtBins_

private

Definition at line 105 of file L1TPhase2MuonOffline.h.

Referenced by getHistBinsEff().

etaNames_

std::map<EtaRegion, std::string> L1TPhase2MuonOffline::etaNames_

private

Definition at line 97 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), and bookResolutionHistos().

etaRegions_

const std::vector<EtaRegion> L1TPhase2MuonOffline::etaRegions_

private

Definition at line 89 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), and bookResolutionHistos().

genparticles_

edm::Handle<std::vector<reco::GenParticle> > L1TPhase2MuonOffline::genparticles_

private

Definition at line 80 of file L1TPhase2MuonOffline.h.

Referenced by analyze().

genParticleToken_

edm::EDGetTokenT<std::vector<reco::GenParticle> > L1TPhase2MuonOffline::genParticleToken_

private

Definition at line 76 of file L1TPhase2MuonOffline.h.

Referenced by analyze().

gmtMuonToken_

edm::EDGetTokenT<l1t::SAMuonCollection> L1TPhase2MuonOffline::gmtMuonToken_

private

Definition at line 74 of file L1TPhase2MuonOffline.h.

Referenced by analyze().

gmtSAMuon_

edm::Handle<l1t::SAMuonCollection> L1TPhase2MuonOffline::gmtSAMuon_

private

Definition at line 78 of file L1TPhase2MuonOffline.h.

Referenced by analyze(), fillControlHistos(), and matchMuonsToGen().

gmtSAMuonPairs_

std::vector<GenMuonGMTPair> L1TPhase2MuonOffline::gmtSAMuonPairs_

private

Definition at line 122 of file L1TPhase2MuonOffline.h.

Referenced by fillEfficiencyHistos(), fillResolutionHistos(), and matchMuonsToGen().

gmtTkMuon_

edm::Handle<l1t::TrackerMuonCollection> L1TPhase2MuonOffline::gmtTkMuon_

private

Definition at line 79 of file L1TPhase2MuonOffline.h.

Referenced by analyze(), fillControlHistos(), and matchMuonsToGen().

gmtTkMuonPairs_

std::vector<GenMuonGMTPair> L1TPhase2MuonOffline::gmtTkMuonPairs_

private

Definition at line 123 of file L1TPhase2MuonOffline.h.

Referenced by fillEfficiencyHistos(), fillResolutionHistos(), and matchMuonsToGen().

gmtTkMuonToken_

edm::EDGetTokenT<l1t::TrackerMuonCollection> L1TPhase2MuonOffline::gmtTkMuonToken_

private

Definition at line 75 of file L1TPhase2MuonOffline.h.

Referenced by analyze().

histFolder_

const std::string L1TPhase2MuonOffline::histFolder_

private

Definition at line 102 of file L1TPhase2MuonOffline.h.

Referenced by bookControlHistos(), bookEfficiencyHistos(), and bookResolutionHistos().

lsb_d0

const float L1TPhase2MuonOffline::lsb_d0 = Phase2L1GMT::LSBSAd0

private

Definition at line 130 of file L1TPhase2MuonOffline.h.

Referenced by fillControlHistos().

lsb_eta

const float L1TPhase2MuonOffline::lsb_eta = Phase2L1GMT::LSBeta

private

Definition at line 128 of file L1TPhase2MuonOffline.h.

Referenced by fillControlHistos().

lsb_phi

const float L1TPhase2MuonOffline::lsb_phi = Phase2L1GMT::LSBphi

private

Definition at line 127 of file L1TPhase2MuonOffline.h.

Referenced by fillControlHistos().

lsb_pt

const float L1TPhase2MuonOffline::lsb_pt = Phase2L1GMT::LSBpt

private

Definition at line 126 of file L1TPhase2MuonOffline.h.

Referenced by fillControlHistos().

lsb_z0

const float L1TPhase2MuonOffline::lsb_z0 = Phase2L1GMT::LSBSAz0

private

Definition at line 129 of file L1TPhase2MuonOffline.h.

Referenced by fillControlHistos().

maxGmtMuonDR_

const double L1TPhase2MuonOffline::maxGmtMuonDR_

private

Definition at line 109 of file L1TPhase2MuonOffline.h.

Referenced by matchMuonsToGen().

muonNames_

std::map<MuType, std::string> L1TPhase2MuonOffline::muonNames_

private

Definition at line 99 of file L1TPhase2MuonOffline.h.

Referenced by bookControlHistos(), bookEfficiencyHistos(), and bookResolutionHistos().

muonTypes_

const std::vector<MuType> L1TPhase2MuonOffline::muonTypes_

private

Definition at line 85 of file L1TPhase2MuonOffline.h.

Referenced by bookHistograms().

qualLevels_

const std::vector<QualLevel> L1TPhase2MuonOffline::qualLevels_

private

Definition at line 90 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), bookResolutionHistos(), and fillResolutionHistos().

qualNames_

std::map<QualLevel, std::string> L1TPhase2MuonOffline::qualNames_

private

Definition at line 98 of file L1TPhase2MuonOffline.h.

Referenced by bookEfficiencyHistos(), and bookResolutionHistos().

resLabels_

std::map<ResType, std::string> L1TPhase2MuonOffline::resLabels_

private

Definition at line 96 of file L1TPhase2MuonOffline.h.

Referenced by bookResolutionHistos().

resNames_

std::map<ResType, std::string> L1TPhase2MuonOffline::resNames_

private

Definition at line 95 of file L1TPhase2MuonOffline.h.

Referenced by bookResolutionHistos().

resolutionHistos_

MonitorElement* L1TPhase2MuonOffline::resolutionHistos_[kNMuTypes][kNEtaRegions][kNQualLevels][kNResTypes]

private

Definition at line 118 of file L1TPhase2MuonOffline.h.

Referenced by bookResolutionHistos(), and fillResolutionHistos().

resTypes_

const std::vector<ResType> L1TPhase2MuonOffline::resTypes_

private

Definition at line 87 of file L1TPhase2MuonOffline.h.

Referenced by bookResolutionHistos(), and fillResolutionHistos().