CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes
L1TMuonDQMOffline Class Reference

#include <L1TMuonDQMOffline.h>

Inheritance diagram for L1TMuonDQMOffline:
one::DQMEDAnalyzer< T > one::dqmimplementation::DQMBaseClass< T... >

Public Types

enum  Control {
  kCtrlTagPt, kCtrlTagEta, kCtrlTagPhi, kCtrlProbePt,
  kCtrlProbeEta, kCtrlProbePhi, kCtrlTagProbeDr, kCtrlTagHltDr,
  kCtrlMuonGmtDeltaR, kCtrlNTightVsAll, kCtrlNProbesVsTight
}
 
enum  EffType { kEffPt, kEffPhi, kEffEta, kEffVtx }
 
enum  EtaRegion {
  kEtaRegionAll, kEtaRegionBmtf, kEtaRegionOmtf, kEtaRegionEmtf,
  kEtaRegionOut
}
 
enum  QualLevel { kQualAll, kQualOpen, kQualDouble, kQualSingle }
 
enum  ResType {
  kResPt, kRes1OverPt, kResQOverPt, kResPhi,
  kResEta, kResCh
}
 

Public Member Functions

 L1TMuonDQMOffline (const edm::ParameterSet &ps)
 
 ~L1TMuonDQMOffline () override
 
- Public Member Functions inherited from one::DQMEDAnalyzer< T >
 DQMEDAnalyzer ()=default
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
 
 ~DQMEDAnalyzer () override=default
 

Protected Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c) override
 
virtual void bookControlHistos (DQMStore::IBooker &)
 
virtual void bookEfficiencyHistos (DQMStore::IBooker &ibooker)
 
void bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &iSetup) override
 
virtual void bookResolutionHistos (DQMStore::IBooker &ibooker)
 
void dqmBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override
 

Private Types

typedef std::tuple< EffType, int, EtaRegionm_histoKeyEffDenVarType
 
typedef std::pair< int, QualLevelm_histoKeyEffNumEtaType
 
typedef std::tuple< EffType, int, EtaRegion, QualLevelm_histoKeyEffNumVarType
 
typedef std::tuple< ResType, EtaRegion, QualLevelm_histoKeyResType
 

Private Member Functions

std::vector< float > getHistBinsEff (EffType eff)
 
std::tuple< int, double, double > getHistBinsRes (ResType res)
 
void getMuonGmtPairs (edm::Handle< l1t::MuonBxCollection > &gmtCands)
 
const unsigned int getNVertices (edm::Handle< reco::VertexCollection > &vertex)
 
const reco::Vertex getPrimaryVertex (edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)
 
void getProbeMuons (edm::Handle< edm::TriggerResults > &trigResults, edm::Handle< trigger::TriggerEvent > &trigEvent)
 
void getTightMuons (edm::Handle< reco::MuonCollection > &muons, const reco::Vertex &vertex)
 
double matchHlt (edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon *mu)
 

Private Attributes

edm::EDGetTokenT< reco::BeamSpotm_BsInputTag
 
std::map< Control, MonitorElement * > m_ControlHistos
 
std::vector< std::pair< int, QualLevel > > m_cuts
 
std::vector< edm::ParameterSetm_cutsVPSet
 
std::map< int, MonitorElement * > m_EfficiencyDenEtaHistos
 
std::map< EtaRegion, MonitorElement * > m_EfficiencyDenPtHistos
 
std::map< std::tuple< EffType, int, EtaRegion >, MonitorElement * > m_EfficiencyDenVarHistos
 
std::map< std::pair< int, QualLevel >, MonitorElement * > m_EfficiencyNumEtaHistos
 
std::map< std::tuple< EffType, int, EtaRegion, QualLevel >, MonitorElement * > m_EfficiencyNumVarHistos
 
std::map< EffType, std::string > m_effLabelStrings
 
std::map< EffType, std::string > m_effStrings
 
const std::vector< EffTypem_effTypes
 
std::vector< double > m_effVsEtaBins
 
std::vector< double > m_effVsPhiBins
 
std::vector< double > m_effVsPtBins
 
std::vector< double > m_effVsVtxBins
 
const std::vector< EtaRegionm_etaRegions
 
std::map< EtaRegion, std::string > m_etaStrings
 
edm::EDGetTokenT< l1t::MuonBxCollectionm_GmtInputTag
 
std::string m_HistFolder
 
HLTConfigProvider m_hltConfig
 
BXVector< l1t::Muonm_L1tL1tMuons
 
std::vector< l1t::MuonBxCollectionm_L1tMuons
 
float m_maxGmtMuonDR
 
float m_maxHltMuonDR
 
float m_minTagProbeDR
 
std::vector< MuonGmtPairm_MuonGmtPairs
 
edm::EDGetTokenT< reco::MuonCollectionm_MuonInputTag
 
std::vector< const reco::Muon * > m_ProbeMuons
 
PropagateToMuon m_propagator
 
const std::vector< QualLevelm_qualLevelsRes
 
std::map< QualLevel, std::string > m_qualStrings
 
std::vector< reco::MuonCollectionm_RecoMuons
 
std::vector< reco::Muonm_RecoRecoMuons
 
double m_recoToL1PtCutFactor
 
std::map< ResType, std::string > m_resLabelStrings
 
std::map< std::tuple< ResType, EtaRegion, QualLevel >, MonitorElement * > m_ResolutionHistos
 
std::map< ResType, std::string > m_resStrings
 
const std::vector< ResTypem_resTypes
 
double m_TagPtCut
 
std::vector< const reco::Muon * > m_TightMuons
 
std::vector< int > m_trigIndices
 
edm::EDGetTokenT< trigger::TriggerEventm_trigInputTag
 
std::vector< std::string > m_trigNames
 
std::string m_trigProcess
 
edm::EDGetTokenT< edm::TriggerResultsm_trigProcess_token
 
bool m_useAtVtxCoord
 
bool m_verbose
 
edm::EDGetTokenT< reco::VertexCollectionm_VtxInputTag
 

Detailed Description

Definition at line 63 of file L1TMuonDQMOffline.h.

Member Typedef Documentation

Definition at line 104 of file L1TMuonDQMOffline.h.

typedef std::pair<int, QualLevel> L1TMuonDQMOffline::m_histoKeyEffNumEtaType
private

Definition at line 103 of file L1TMuonDQMOffline.h.

Definition at line 102 of file L1TMuonDQMOffline.h.

Definition at line 101 of file L1TMuonDQMOffline.h.

Member Enumeration Documentation

Constructor & Destructor Documentation

L1TMuonDQMOffline::L1TMuonDQMOffline ( const edm::ParameterSet ps)

Definition at line 95 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, kEffEta, kEffPhi, kEffPt, kEffVtx, kEtaRegionAll, kEtaRegionBmtf, kEtaRegionEmtf, kEtaRegionOmtf, kQualAll, kQualDouble, kQualOpen, kQualSingle, kRes1OverPt, kResCh, kResEta, kResPhi, kResPt, kResQOverPt, m_BsInputTag, m_cuts, m_cutsVPSet, m_effLabelStrings, m_effStrings, m_effVsEtaBins, m_effVsPhiBins, m_effVsPtBins, m_effVsVtxBins, m_etaRegions, m_etaStrings, m_GmtInputTag, m_HistFolder, m_maxGmtMuonDR, m_maxHltMuonDR, m_minTagProbeDR, m_MuonInputTag, m_qualLevelsRes, m_qualStrings, m_recoToL1PtCutFactor, m_resLabelStrings, m_resStrings, m_resTypes, m_TagPtCut, m_trigInputTag, m_trigNames, m_trigProcess, m_trigProcess_token, m_useAtVtxCoord, m_verbose, and m_VtxInputTag.

95  :
101  m_effStrings({ {kEffPt, "pt"}, {kEffPhi, "phi"}, {kEffEta, "eta"}, {kEffVtx, "vtx"} }),
102  m_effLabelStrings({ {kEffPt, "p_{T} (GeV)"}, {kEffPhi, "#phi"}, {kEffEta, "#eta"}, {kEffVtx, "# vertices"} }),
103  m_resStrings({ {kResPt, "pt"}, {kRes1OverPt, "1overpt"}, {kResQOverPt, "qoverpt"}, {kResPhi, "phi"}, {kResEta, "eta"}, {kResCh, "charge"} }),
104  m_resLabelStrings({ {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}"} }),
105  m_etaStrings({ {kEtaRegionAll, "etaMin0_etaMax2p4"}, {kEtaRegionBmtf, "etaMin0_etaMax0p83"}, {kEtaRegionOmtf, "etaMin0p83_etaMax1p24"}, {kEtaRegionEmtf, "etaMin1p24_etaMax2p4"} }),
106  m_qualStrings({ {kQualAll, "qualAll"}, {kQualOpen, "qualOpen"}, {kQualDouble, "qualDouble"}, {kQualSingle, "qualSingle"} }),
107  m_verbose(ps.getUntrackedParameter<bool>("verbose")),
108  m_HistFolder(ps.getUntrackedParameter<string>("histFolder")),
109  m_TagPtCut(ps.getUntrackedParameter<double>("tagPtCut")),
110  m_recoToL1PtCutFactor(ps.getUntrackedParameter<double>("recoToL1PtCutFactor")),
111  m_cutsVPSet(ps.getUntrackedParameter<std::vector<edm::ParameterSet>>("cuts")),
112  m_MuonInputTag(consumes<reco::MuonCollection>(ps.getUntrackedParameter<InputTag>("muonInputTag"))),
113  m_GmtInputTag(consumes<l1t::MuonBxCollection>(ps.getUntrackedParameter<InputTag>("gmtInputTag"))),
114  m_VtxInputTag(consumes<VertexCollection>(ps.getUntrackedParameter<InputTag>("vtxInputTag"))),
115  m_BsInputTag(consumes<BeamSpot>(ps.getUntrackedParameter<InputTag>("bsInputTag"))),
116  m_trigInputTag(consumes<trigger::TriggerEvent>(ps.getUntrackedParameter<InputTag>("trigInputTag"))),
117  m_trigProcess(ps.getUntrackedParameter<string>("trigProcess")),
118  m_trigProcess_token(consumes<edm::TriggerResults>(ps.getUntrackedParameter<InputTag>("trigProcess_token"))),
119  m_trigNames(ps.getUntrackedParameter<vector<string> >("triggerNames")),
120  m_effVsPtBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPtBins")),
121  m_effVsPhiBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPhiBins")),
122  m_effVsEtaBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsEtaBins")),
123  m_effVsVtxBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsVtxBins")),
124  m_useAtVtxCoord(ps.getUntrackedParameter<bool>("useL1AtVtxCoord")),
125  m_maxGmtMuonDR(0.3),
126  m_minTagProbeDR(0.5),
127  m_maxHltMuonDR(0.1)
128 {
129  if (m_verbose) cout << "[L1TMuonDQMOffline:] ____________ Storage initialization ____________ " << endl;
130 
131  for (const auto cutsPSet : m_cutsVPSet) {
132  const auto qCut = cutsPSet.getUntrackedParameter<int>("qualCut");
133  QualLevel qLevel = kQualAll;
134  if (qCut > 11) {
135  qLevel = kQualSingle;
136  } else if (qCut > 7) {
137  qLevel = kQualDouble;
138  } else if (qCut > 3) {
139  qLevel = kQualOpen;
140  }
141  m_cuts.emplace_back(std::make_pair(cutsPSet.getUntrackedParameter<int>("ptCut"), qLevel));
142  }
143 }
std::map< EffType, std::string > m_effStrings
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< reco::BeamSpot > m_BsInputTag
edm::EDGetTokenT< reco::MuonCollection > m_MuonInputTag
std::vector< double > m_effVsVtxBins
std::map< ResType, std::string > m_resLabelStrings
std::map< EtaRegion, std::string > m_etaStrings
edm::EDGetTokenT< edm::TriggerResults > m_trigProcess_token
const std::vector< EtaRegion > m_etaRegions
std::vector< double > m_effVsPhiBins
std::string m_HistFolder
std::vector< std::string > m_trigNames
std::map< QualLevel, std::string > m_qualStrings
std::string m_trigProcess
std::vector< edm::ParameterSet > m_cutsVPSet
std::map< EffType, std::string > m_effLabelStrings
edm::EDGetTokenT< l1t::MuonBxCollection > m_GmtInputTag
std::vector< double > m_effVsPtBins
std::vector< std::pair< int, QualLevel > > m_cuts
const std::vector< EffType > m_effTypes
const std::vector< QualLevel > m_qualLevelsRes
std::map< ResType, std::string > m_resStrings
std::vector< double > m_effVsEtaBins
const std::vector< ResType > m_resTypes
edm::EDGetTokenT< reco::VertexCollection > m_VtxInputTag
edm::EDGetTokenT< trigger::TriggerEvent > m_trigInputTag
PropagateToMuon m_propagator
L1TMuonDQMOffline::~L1TMuonDQMOffline ( )
override

Definition at line 146 of file L1TMuonDQMOffline.cc.

146 { }

Member Function Documentation

void L1TMuonDQMOffline::analyze ( const edm::Event e,
const edm::EventSetup c 
)
overrideprotected

Definition at line 182 of file L1TMuonDQMOffline.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, gather_cfg::cout, TkAlMuonSelectors_cfi::cut, edm::Event::getByToken(), getMuonGmtPairs(), getNVertices(), getPrimaryVertex(), getProbeMuons(), getTightMuons(), kEffEta, kEffPt, kEffVtx, kEtaRegionAll, kEtaRegionOut, kQualAll, kResPt, m_BsInputTag, m_cuts, m_EfficiencyDenEtaHistos, m_EfficiencyDenPtHistos, m_EfficiencyDenVarHistos, m_EfficiencyNumEtaHistos, m_EfficiencyNumVarHistos, m_effTypes, m_GmtInputTag, m_MuonGmtPairs, m_MuonInputTag, m_qualLevelsRes, m_recoToL1PtCutFactor, m_ResolutionHistos, m_resTypes, m_trigInputTag, m_trigProcess_token, m_verbose, m_VtxInputTag, extraflags_cff::muons, impactParameterTagInfos_cfi::primaryVertex, and JetChargeProducer_cfi::var.

182  {
183 
185  iEvent.getByToken(m_MuonInputTag, muons);
187  iEvent.getByToken(m_BsInputTag, beamSpot);
189  iEvent.getByToken(m_VtxInputTag, vertex);
191  iEvent.getByToken(m_GmtInputTag,gmtCands);
192  Handle<edm::TriggerResults> trigResults;
193  iEvent.getByToken(m_trigProcess_token,trigResults);
195  iEvent.getByToken(m_trigInputTag,trigEvent);
196 
197  const auto nVtx = getNVertices(vertex);
198  const Vertex primaryVertex = getPrimaryVertex(vertex,beamSpot);
199 
200  getTightMuons(muons,primaryVertex);
201  getProbeMuons(trigResults,trigEvent); // CB add flag to run on orthogonal datasets (no T&P)
202 
203  getMuonGmtPairs(gmtCands);
204 
205  if (m_verbose) cout << "[L1TMuonDQMOffline:] Computing efficiencies" << endl;
206 
207  vector<MuonGmtPair>::const_iterator muonGmtPairsIt = m_MuonGmtPairs.begin();
208  vector<MuonGmtPair>::const_iterator muonGmtPairsEnd = m_MuonGmtPairs.end();
209 
210  // To fill once for global eta and once for TF eta region of the L1T muon.
211  // The second entry is a placeholder and will be replaced by the TF eta region of the L1T muon.
212  std::array<EtaRegion, 2> regsToFill { {kEtaRegionAll, kEtaRegionAll} };
213 
214  for(; muonGmtPairsIt!=muonGmtPairsEnd; ++muonGmtPairsIt) {
215  // Fill the resolution histograms
216  if( (muonGmtPairsIt->etaRegion() != kEtaRegionOut) && (muonGmtPairsIt->gmtPt() > 0) ){
217  regsToFill[1] = muonGmtPairsIt->etaRegion();
219  for (const auto var : m_resTypes) {
220  const auto varToFill = muonGmtPairsIt->getDeltaVar(var);
221  std::get<0>(histoKeyRes) = var;
222  // Fill for the global eta and for TF eta region that the probe muon is in
223  for (const auto regToFill : regsToFill) {
224  std::get<1>(histoKeyRes) = regToFill;
225  for (const auto qualLevel : m_qualLevelsRes) {
226  // This assumes that the qualLevel enum has increasing qualities
227  // HW quality levels can be 0, 4, 8, or 12
228  int qualCut = qualLevel * 4;
229  if (muonGmtPairsIt->gmtQual() >= qualCut) {
230  std::get<2>(histoKeyRes) = qualLevel;
231  m_ResolutionHistos[histoKeyRes]->Fill(varToFill);
232  }
233  }
234  }
235  }
236  }
237 
238  // Fill the efficiency numerator and denominator histograms
239  if (muonGmtPairsIt->etaRegion() != kEtaRegionOut) {
240  unsigned int cutsCounter = 0;
241  for (const auto cut : m_cuts) {
242  const auto gmtPtCut = cut.first;
243  const auto qualLevel = cut.second;
244  const bool gmtAboveCut = (muonGmtPairsIt->gmtPt() > gmtPtCut);
245 
246  // default keys
247  m_histoKeyEffDenVarType histoKeyEffDenVar = {kEffPt, gmtPtCut, kEtaRegionAll};
248  m_histoKeyEffNumVarType histoKeyEffNumVar = {kEffPt, gmtPtCut, kEtaRegionAll, qualLevel};
249 
250  regsToFill[1] = muonGmtPairsIt->etaRegion();
251  for(const auto var : m_effTypes) {
252  if(var != kEffPt){
253  if (muonGmtPairsIt->pt() < m_recoToL1PtCutFactor * gmtPtCut) break; // efficiency at plateau
254  }
255  double varToFill;
256  if (var == kEffVtx) {
257  varToFill = static_cast<double>(nVtx);
258  } else {
259  varToFill = muonGmtPairsIt->getVar(var);
260  }
261  // Fill denominators
262  if (var == kEffEta) {
263  m_EfficiencyDenEtaHistos[gmtPtCut]->Fill(varToFill);
264  } else {
265  std::get<0>(histoKeyEffDenVar) = var;
266  // Fill for the global eta and for TF eta region that the probe muon is in
267  for (const auto regToFill : regsToFill) {
268  if (var == kEffPt) {
269  if (cutsCounter == 0) {
270  m_EfficiencyDenPtHistos[regToFill]->Fill(varToFill);
271  }
272  } else {
273  std::get<2>(histoKeyEffDenVar) = regToFill;
274  m_EfficiencyDenVarHistos[histoKeyEffDenVar]->Fill(varToFill);
275  }
276  }
277  }
278  // Fill numerators
279  std::get<0>(histoKeyEffNumVar) = var;
280  // This assumes that the qualLevel enum has increasing qualities
281  if (gmtAboveCut && muonGmtPairsIt->gmtQual() >= qualLevel * 4) {
282  if (var == kEffEta) {
283  m_histoKeyEffNumEtaType histoKeyEffNumEta = {gmtPtCut, qualLevel};
284  m_EfficiencyNumEtaHistos[histoKeyEffNumEta]->Fill(varToFill);
285  } else {
286  std::get<3>(histoKeyEffNumVar) = qualLevel;
287  // Fill for the global eta and for TF eta region that the probe muon is in
288  for (const auto regToFill : regsToFill) {
289  std::get<2>(histoKeyEffNumVar) = regToFill;
290  m_EfficiencyNumVarHistos[histoKeyEffNumVar]->Fill(varToFill);
291  }
292  }
293  }
294  }
295  ++cutsCounter;
296  }
297  }
298  }
299 
300  if (m_verbose) cout << "[L1TMuonDQMOffline:] Computation finished" << endl;
301 }
edm::EDGetTokenT< reco::BeamSpot > m_BsInputTag
edm::EDGetTokenT< reco::MuonCollection > m_MuonInputTag
edm::EDGetTokenT< edm::TriggerResults > m_trigProcess_token
std::map< std::tuple< ResType, EtaRegion, QualLevel >, MonitorElement * > m_ResolutionHistos
std::map< EtaRegion, MonitorElement * > m_EfficiencyDenPtHistos
std::tuple< EffType, int, EtaRegion, QualLevel > m_histoKeyEffNumVarType
std::pair< int, QualLevel > m_histoKeyEffNumEtaType
std::tuple< ResType, EtaRegion, QualLevel > m_histoKeyResType
void getProbeMuons(edm::Handle< edm::TriggerResults > &trigResults, edm::Handle< trigger::TriggerEvent > &trigEvent)
const unsigned int getNVertices(edm::Handle< reco::VertexCollection > &vertex)
int iEvent
Definition: GenABIO.cc:230
edm::EDGetTokenT< l1t::MuonBxCollection > m_GmtInputTag
std::tuple< EffType, int, EtaRegion > m_histoKeyEffDenVarType
std::map< std::pair< int, QualLevel >, MonitorElement * > m_EfficiencyNumEtaHistos
std::map< int, MonitorElement * > m_EfficiencyDenEtaHistos
std::vector< std::pair< int, QualLevel > > m_cuts
const std::vector< EffType > m_effTypes
const std::vector< QualLevel > m_qualLevelsRes
std::vector< MuonGmtPair > m_MuonGmtPairs
std::map< std::tuple< EffType, int, EtaRegion >, MonitorElement * > m_EfficiencyDenVarHistos
void getMuonGmtPairs(edm::Handle< l1t::MuonBxCollection > &gmtCands)
const reco::Vertex getPrimaryVertex(edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)
const std::vector< ResType > m_resTypes
edm::EDGetTokenT< reco::VertexCollection > m_VtxInputTag
edm::EDGetTokenT< trigger::TriggerEvent > m_trigInputTag
std::map< std::tuple< EffType, int, EtaRegion, QualLevel >, MonitorElement * > m_EfficiencyNumVarHistos
void getTightMuons(edm::Handle< reco::MuonCollection > &muons, const reco::Vertex &vertex)
void L1TMuonDQMOffline::bookControlHistos ( DQMStore::IBooker ibooker)
protectedvirtual

Definition at line 304 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, kCtrlMuonGmtDeltaR, kCtrlNProbesVsTight, kCtrlNTightVsAll, kCtrlProbeEta, kCtrlProbePhi, kCtrlProbePt, kCtrlTagEta, kCtrlTagHltDr, kCtrlTagPhi, kCtrlTagProbeDr, kCtrlTagPt, m_ControlHistos, m_HistFolder, m_verbose, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

304  {
305  if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Control Plot Histos" << endl;
306 
307  ibooker.setCurrentFolder(m_HistFolder+"/control_variables");
308 
309  m_ControlHistos[kCtrlMuonGmtDeltaR] = ibooker.book1D("MuonGmtDeltaR", "MuonGmtDeltaR; #DeltaR", 50, 0., 0.5);
310  m_ControlHistos[kCtrlNTightVsAll] = ibooker.book2D("NTightVsAll", "NTightVsAll; # muons; # tight muons", 20, -0.5, 19.5, 16, -0.5, 15.5);
311  m_ControlHistos[kCtrlNProbesVsTight] = ibooker.book2D("NProbesVsTight", "NProbesVsTight; # tight muons; # probe muons", 8, -0.5, 7.5, 8, -0.5, 7.5);
312 
313  m_ControlHistos[kCtrlTagPt] = ibooker.book1D("TagMuonPt", "TagMuonPt; p_{T}", 50, 0., 100.);
314  m_ControlHistos[kCtrlTagPhi] = ibooker.book1D("TagMuonPhi", "TagMuonPhi; #phi", 66, -3.3, 3.3);
315  m_ControlHistos[kCtrlTagEta] = ibooker.book1D("TagMuonEta", "TagMuonEta; #eta", 50, -2.5, 2.5);
316 
317  m_ControlHistos[kCtrlProbePt] = ibooker.book1D("ProbeMuonPt", "ProbeMuonPt; p_{T}", 50, 0., 100.);
318  m_ControlHistos[kCtrlProbePhi] = ibooker.book1D("ProbeMuonPhi", "ProbeMuonPhi; #phi", 66, -3.3, 3.3);
319  m_ControlHistos[kCtrlProbeEta] = ibooker.book1D("ProbeMuonEta", "ProbeMuonEta; #eta", 50, -2.5, 2.5);
320 
321  m_ControlHistos[kCtrlTagProbeDr] = ibooker.book1D("TagMuonProbeMuonDeltaR", "TagMuonProbeMuonDeltaR; #DeltaR", 50, 0.,5.0);
322  m_ControlHistos[kCtrlTagHltDr] = ibooker.book1D("TagMuonHltDeltaR", "TagMuonHltDeltaR;#DeltaR", 55, 0., 0.11);
323 }
std::string m_HistFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:109
std::map< Control, MonitorElement * > m_ControlHistos
void L1TMuonDQMOffline::bookEfficiencyHistos ( DQMStore::IBooker ibooker)
protectedvirtual

Definition at line 326 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), TkAlMuonSelectors_cfi::cut, getHistBinsEff(), kEffEta, kEffPt, m_cuts, m_EfficiencyDenEtaHistos, m_EfficiencyDenPtHistos, m_EfficiencyDenVarHistos, m_EfficiencyNumEtaHistos, m_EfficiencyNumVarHistos, m_effLabelStrings, m_effStrings, m_effTypes, m_etaRegions, m_etaStrings, m_HistFolder, m_qualStrings, dataset::name, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and JetChargeProducer_cfi::var.

Referenced by bookHistograms().

326  {
327  ibooker.setCurrentFolder(m_HistFolder+"/numerators_and_denominators");
328 
329  for(const auto var : m_effTypes) {
330  auto histBins = getHistBinsEff(var);
331  // histograms for eta variable get a special treatment
332  if (var == kEffEta) {
333  for (const auto cut : m_cuts) {
334  const auto gmtPtCut = cut.first;
335  const auto qualLevel = cut.second;
336  std::string name = "effDen_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut);
337  m_EfficiencyDenEtaHistos[gmtPtCut] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
338  name = "effNum_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_qualStrings[qualLevel];
339  m_histoKeyEffNumEtaType histoKeyEffNumEta = {gmtPtCut, qualLevel};
340  m_EfficiencyNumEtaHistos[histoKeyEffNumEta] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
341  }
342  } else {
343  for (const auto etaReg : m_etaRegions) {
344  // denominator histograms for pt variable get a special treatment
345  if (var == kEffPt) {
346  std::string name = "effDen_"+m_effStrings[var]+"_"+m_etaStrings[etaReg];
347  m_EfficiencyDenPtHistos[etaReg] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
348  } else {
349  for (const auto cut : m_cuts) {
350  const int gmtPtCut = cut.first;
351  std::string name = "effDen_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_etaStrings[etaReg];
352  m_histoKeyEffDenVarType histoKeyEffDenVar = {var, gmtPtCut, etaReg};
353  m_EfficiencyDenVarHistos[histoKeyEffDenVar] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
354  }
355  }
356  for (const auto cut : m_cuts) {
357  const auto gmtPtCut = cut.first;
358  const auto qualLevel = cut.second;
359  std::string name = "effNum_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_etaStrings[etaReg]+"_"+m_qualStrings[qualLevel];
360  m_histoKeyEffNumVarType histoKeyEffNum = {var, gmtPtCut, etaReg, qualLevel};
361  m_EfficiencyNumVarHistos[histoKeyEffNum] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
362  }
363  }
364  }
365  }
366 }
std::map< EffType, std::string > m_effStrings
std::map< EtaRegion, std::string > m_etaStrings
const std::vector< EtaRegion > m_etaRegions
std::map< EtaRegion, MonitorElement * > m_EfficiencyDenPtHistos
std::tuple< EffType, int, EtaRegion, QualLevel > m_histoKeyEffNumVarType
std::string m_HistFolder
std::pair< int, QualLevel > m_histoKeyEffNumEtaType
std::vector< float > getHistBinsEff(EffType eff)
std::map< QualLevel, std::string > m_qualStrings
std::map< EffType, std::string > m_effLabelStrings
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
std::tuple< EffType, int, EtaRegion > m_histoKeyEffDenVarType
std::map< std::pair< int, QualLevel >, MonitorElement * > m_EfficiencyNumEtaHistos
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
std::map< int, MonitorElement * > m_EfficiencyDenEtaHistos
std::vector< std::pair< int, QualLevel > > m_cuts
const std::vector< EffType > m_effTypes
std::map< std::tuple< EffType, int, EtaRegion >, MonitorElement * > m_EfficiencyDenVarHistos
std::map< std::tuple< EffType, int, EtaRegion, QualLevel >, MonitorElement * > m_EfficiencyNumVarHistos
void L1TMuonDQMOffline::bookHistograms ( DQMStore::IBooker ibooker,
const edm::Run run,
const edm::EventSetup iSetup 
)
overrideprotected

Definition at line 156 of file L1TMuonDQMOffline.cc.

References bookControlHistos(), bookEfficiencyHistos(), bookResolutionHistos(), gather_cfg::cout, createfilelist::int, m_hltConfig, m_trigIndices, m_trigNames, m_verbose, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().

156  {
157  //book histos
158  bookControlHistos(ibooker);
159  bookEfficiencyHistos(ibooker);
160  bookResolutionHistos(ibooker);
161 
162  vector<string>::const_iterator trigNamesIt = m_trigNames.begin();
163  vector<string>::const_iterator trigNamesEnd = m_trigNames.end();
164 
165  for (; trigNamesIt!=trigNamesEnd; ++trigNamesIt) {
166  TString tNameTmp = TString(*trigNamesIt); // use TString as it handles regex
167  TRegexp tNamePattern = TRegexp(tNameTmp,true);
168  int tIndex = -1;
169 
170  for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
171  TString tmpName = TString(m_hltConfig.triggerName(ipath));
172  if (tmpName.Contains(tNamePattern)) {
173  tIndex = int(ipath);
174  m_trigIndices.push_back(tIndex);
175  }
176  }
177  if (tIndex < 0 && m_verbose) cout << "[L1TMuonDQMOffline:] Warning: Could not find trigger " << (*trigNamesIt) << endl;
178  }
179 }
unsigned int size() const
number of trigger paths in trigger table
const std::string & triggerName(unsigned int triggerIndex) const
std::vector< std::string > m_trigNames
virtual void bookControlHistos(DQMStore::IBooker &)
virtual void bookResolutionHistos(DQMStore::IBooker &ibooker)
HLTConfigProvider m_hltConfig
std::vector< int > m_trigIndices
virtual void bookEfficiencyHistos(DQMStore::IBooker &ibooker)
void L1TMuonDQMOffline::bookResolutionHistos ( DQMStore::IBooker ibooker)
protectedvirtual

Definition at line 368 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), gather_cfg::cout, getHistBinsRes(), m_etaRegions, m_etaStrings, m_HistFolder, m_qualLevelsRes, m_qualStrings, m_resLabelStrings, m_ResolutionHistos, m_resStrings, m_resTypes, m_verbose, dataset::name, pileupCalc::nbins, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, JetChargeProducer_cfi::var, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

Referenced by bookHistograms().

368  {
369  if(m_verbose) cout << "[L1TMuonOffline:] Booking Resolution Plot Histos" << endl;
370  ibooker.setCurrentFolder(m_HistFolder+"/resolution");
371 
372  for (const auto var : m_resTypes) {
373  auto nbins = std::get<0>(getHistBinsRes(var));
374  auto xmin = std::get<1>(getHistBinsRes(var));
375  auto xmax = std::get<2>(getHistBinsRes(var));
376  for (const auto etaReg : m_etaRegions) {
377  for (const auto qualLevel : m_qualLevelsRes) {
378  m_histoKeyResType histoKeyRes = {var, etaReg, qualLevel};
379  std::string name = "resolution_"+m_resStrings[var]+"_"+m_etaStrings[etaReg]+"_"+m_qualStrings[qualLevel];
380  m_ResolutionHistos[histoKeyRes] = ibooker.book1D(name, name+";"+m_resLabelStrings[var], nbins, xmin, xmax);
381  }
382  }
383  }
384 }
std::map< ResType, std::string > m_resLabelStrings
std::map< EtaRegion, std::string > m_etaStrings
std::map< std::tuple< ResType, EtaRegion, QualLevel >, MonitorElement * > m_ResolutionHistos
const std::vector< EtaRegion > m_etaRegions
std::string m_HistFolder
std::tuple< ResType, EtaRegion, QualLevel > m_histoKeyResType
std::map< QualLevel, std::string > m_qualStrings
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
std::tuple< int, double, double > getHistBinsRes(ResType res)
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
const std::vector< QualLevel > m_qualLevelsRes
std::map< ResType, std::string > m_resStrings
const std::vector< ResType > m_resTypes
void L1TMuonDQMOffline::dqmBeginRun ( const edm::Run run,
const edm::EventSetup iSetup 
)
overrideprotected

Definition at line 148 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, PropagateToMuon::init(), HLTConfigProvider::init(), m_hltConfig, m_propagator, m_trigProcess, and m_verbose.

148  {
149  if (m_verbose) cout << "[L1TMuonDQMOffline:] Called beginRun." << endl;
150  bool changed = true;
151  m_hltConfig.init(run,iSetup,m_trigProcess,changed);
152  m_propagator.init(iSetup);
153 }
std::string m_trigProcess
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
HLTConfigProvider m_hltConfig
void init(const edm::EventSetup &iSetup)
Call this method at the beginning of each run, to initialize geometry, magnetic field and propagators...
PropagateToMuon m_propagator
std::vector< float > L1TMuonDQMOffline::getHistBinsEff ( EffType  eff)
private

Definition at line 568 of file L1TMuonDQMOffline.cc.

References kEffEta, kEffPhi, kEffPt, kEffVtx, m_effVsEtaBins, m_effVsPhiBins, m_effVsPtBins, and m_effVsVtxBins.

Referenced by bookEfficiencyHistos().

568  {
569  if (eff == kEffPt) {
570  std::vector<float> effVsPtBins(m_effVsPtBins.begin(), m_effVsPtBins.end());
571  return effVsPtBins;
572  }
573  if (eff == kEffPhi) {
574  std::vector<float> effVsPhiBins(m_effVsPhiBins.begin(), m_effVsPhiBins.end());
575  return effVsPhiBins;
576  }
577  if (eff == kEffEta) {
578  std::vector<float> effVsEtaBins(m_effVsEtaBins.begin(), m_effVsEtaBins.end());
579  return effVsEtaBins;
580  }
581  if (eff == kEffVtx) {
582  std::vector<float> effVsVtxBins(m_effVsVtxBins.begin(), m_effVsVtxBins.end());
583  return effVsVtxBins;
584  }
585  return {0., 1.};
586 }
std::vector< double > m_effVsVtxBins
std::vector< double > m_effVsPhiBins
std::vector< double > m_effVsPtBins
std::vector< double > m_effVsEtaBins
std::tuple< int, double, double > L1TMuonDQMOffline::getHistBinsRes ( ResType  res)
private

Definition at line 588 of file L1TMuonDQMOffline.cc.

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

Referenced by bookResolutionHistos().

588  {
589  if (res == kResPt) return {50, -2., 2.};
590  if (res == kRes1OverPt) return {50, -2., 2.};
591  if (res == kResQOverPt) return {50, -2., 2.};
592  if (res == kResPhi) return {96, -0.2, 0.2};
593  if (res == kResEta) return {100, -0.1, 0.1};
594  if (res == kResCh) return {5, -2, 3};
595  return {1, 0, 1};
596 }
Definition: Electron.h:6
void L1TMuonDQMOffline::getMuonGmtPairs ( edm::Handle< l1t::MuonBxCollection > &  gmtCands)
private

Definition at line 505 of file L1TMuonDQMOffline.cc.

References BXVector< T >::begin(), gather_cfg::cout, MuonGmtPair::dR(), BXVector< T >::end(), MuonGmtPair::getVar(), kCtrlMuonGmtDeltaR, kCtrlProbeEta, kCtrlProbePhi, kCtrlProbePt, kEffEta, kEffPhi, kEffPt, m_ControlHistos, m_maxGmtMuonDR, m_MuonGmtPairs, m_ProbeMuons, m_propagator, m_useAtVtxCoord, and m_verbose.

Referenced by analyze().

505  {
506 
507  m_MuonGmtPairs.clear();
508  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting muon GMT pairs" << endl;
509 
510  vector<const reco::Muon*>::const_iterator probeMuIt = m_ProbeMuons.begin();
511  vector<const reco::Muon*>::const_iterator probeMuEnd = m_ProbeMuons.end();
512 
514  l1t::MuonBxCollection::const_iterator gmtEnd = gmtCands->end(0);
515 
516  for (; probeMuIt!=probeMuEnd; ++probeMuIt) {
517  MuonGmtPair pairBestCand((*probeMuIt), nullptr, m_propagator, m_useAtVtxCoord);
518 
519  // Fill the control histograms with the probe muon kinematic variables used
520  m_ControlHistos[kCtrlProbeEta]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffEta));
521  m_ControlHistos[kCtrlProbePhi]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffPhi));
522  m_ControlHistos[kCtrlProbePt]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffPt));
523 
524  gmtIt = gmtCands->begin(0); // use only on L1T muons from BX 0
525 
526  for(; gmtIt!=gmtEnd; ++gmtIt) {
527  MuonGmtPair pairTmpCand((*probeMuIt), &(*gmtIt), m_propagator, m_useAtVtxCoord);
528 
529  if ( (pairTmpCand.dR() < m_maxGmtMuonDR) && (pairTmpCand.dR() < pairBestCand.dR() ) ) {
530  pairBestCand = pairTmpCand;
531  }
532 
533  }
534  m_MuonGmtPairs.push_back(pairBestCand);
535  m_ControlHistos[kCtrlMuonGmtDeltaR]->Fill(pairBestCand.dR());
536  }
537 }
const_iterator end(int bx) const
std::vector< const reco::Muon * > m_ProbeMuons
std::vector< MuonGmtPair > m_MuonGmtPairs
PropagateToMuon m_propagator
const_iterator begin(int bx) const
std::map< Control, MonitorElement * > m_ControlHistos
std::vector< Muon >::const_iterator const_iterator
Definition: BXVector.h:20
const unsigned int L1TMuonDQMOffline::getNVertices ( edm::Handle< reco::VertexCollection > &  vertex)
private

Definition at line 387 of file L1TMuonDQMOffline.cc.

References edm::HandleBase::isValid().

Referenced by analyze().

387  {
388  unsigned int nVtx = 0;
389 
390  if (vertex.isValid()) {
391  for (const auto vertexIt : *vertex) {
392  if (vertexIt.isValid() && !vertexIt.isFake()) {
393  ++nVtx;
394  }
395  }
396  }
397  return nVtx;
398 }
bool isValid() const
Definition: HandleBase.h:74
const reco::Vertex L1TMuonDQMOffline::getPrimaryVertex ( edm::Handle< reco::VertexCollection > &  vertex,
edm::Handle< reco::BeamSpot > &  beamSpot 
)
private

Definition at line 401 of file L1TMuonDQMOffline.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), edm::HandleBase::isValid(), reco::BeamSpot::position(), impactParameterTagInfos_cfi::primaryVertex, and reco::BeamSpot::sigmaZ().

Referenced by analyze().

402  {
403  Vertex::Point posVtx;
404  Vertex::Error errVtx;
405 
406  bool hasPrimaryVertex = false;
407 
408  if (vertex.isValid()) {
409  vector<Vertex>::const_iterator vertexIt = vertex->begin();
410  vector<Vertex>::const_iterator vertexEnd = vertex->end();
411 
412  for (;vertexIt!=vertexEnd;++vertexIt) {
413  if (vertexIt->isValid() && !vertexIt->isFake()) {
414  posVtx = vertexIt->position();
415  errVtx = vertexIt->error();
416  hasPrimaryVertex = true;
417  break;
418  }
419  }
420  }
421 
422  if ( !hasPrimaryVertex ) {
423  posVtx = beamSpot->position();
424  errVtx(0,0) = beamSpot->BeamWidthX();
425  errVtx(1,1) = beamSpot->BeamWidthY();
426  errVtx(2,2) = beamSpot->sigmaZ();
427  }
428  const Vertex primaryVertex(posVtx,errVtx);
429  return primaryVertex;
430 }
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:86
bool isValid() const
Definition: HandleBase.h:74
double sigmaZ() const
sigma z
Definition: BeamSpot.h:80
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:88
const Point & position() const
position
Definition: BeamSpot.h:62
void L1TMuonDQMOffline::getProbeMuons ( edm::Handle< edm::TriggerResults > &  trigResults,
edm::Handle< trigger::TriggerEvent > &  trigEvent 
)
private

Definition at line 449 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, particleFlow_cfi::dEta, SoftLeptonByDistance_cfi::distance, particleFlow_cfi::dPhi, PVValHelper::eta, kCtrlNProbesVsTight, kCtrlTagEta, kCtrlTagHltDr, kCtrlTagPhi, kCtrlTagProbeDr, kCtrlTagPt, m_ControlHistos, m_minTagProbeDR, m_ProbeMuons, m_TagPtCut, m_TightMuons, m_verbose, matchHlt(), phi, EnergyCorrector::pt, and mathSSE::sqrt().

Referenced by analyze().

450  {
451 
452  if (m_verbose) cout << "[L1TMuonDQMOffline:] getting probe muons" << endl;
453  m_ProbeMuons.clear();
454  std::vector<const reco::Muon*> tagMuonsInHist;
455 
456  tagMuonsInHist.clear();
457 
458  vector<const reco::Muon*>::const_iterator probeCandIt = m_TightMuons.begin();
459  vector<const reco::Muon*>::const_iterator tightMuonsEnd = m_TightMuons.end();
460 
461  for (; probeCandIt!=tightMuonsEnd; ++probeCandIt) {
462  bool isProbe = false;
463  vector<const reco::Muon*>::const_iterator tagCandIt = m_TightMuons.begin();
464  float deltar = 0.;
465 
466  for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
467  bool tagMuonAlreadyInHist = false;
468  bool tagHasTrig = false;
469  float eta = (*tagCandIt)->eta();
470  float phi = (*tagCandIt)->phi();
471  float pt = (*tagCandIt)->pt();
472  float dEta = eta - (*probeCandIt)->eta();
473  float dPhi = phi - (*probeCandIt)->phi();
474  deltar = sqrt(dEta*dEta + dPhi*dPhi);
475 
476  if ( (*tagCandIt) == (*probeCandIt) || deltar<m_minTagProbeDR ) continue; // CB has a little bias for closed-by muons
477  auto matchHltDeltaR = matchHlt(trigEvent,(*tagCandIt));
478  tagHasTrig = (matchHltDeltaR < m_maxHltMuonDR) && (pt > m_TagPtCut);
479  isProbe |= tagHasTrig;
480  if (tagHasTrig) {
481  if (std::distance(m_TightMuons.begin(), m_TightMuons.end()) > 2 ) {
482  for (vector<const reco::Muon*>::const_iterator tagMuonsInHistIt = tagMuonsInHist.begin(); tagMuonsInHistIt!=tagMuonsInHist.end(); ++tagMuonsInHistIt) {
483  if ( (*tagCandIt) == (*tagMuonsInHistIt) ) {
484  tagMuonAlreadyInHist = true;
485  break;
486  }
487  }
488  if (tagMuonAlreadyInHist == false) tagMuonsInHist.push_back((*tagCandIt));
489  }
490  if (tagMuonAlreadyInHist == false) {
491  m_ControlHistos[kCtrlTagEta]->Fill(eta);
492  m_ControlHistos[kCtrlTagPhi]->Fill(phi);
493  m_ControlHistos[kCtrlTagPt]->Fill(pt);
494  m_ControlHistos[kCtrlTagProbeDr]->Fill(deltar);
495  m_ControlHistos[kCtrlTagHltDr]->Fill(matchHltDeltaR);
496  }
497  }
498  }
499  if (isProbe) m_ProbeMuons.push_back((*probeCandIt));
500  }
502 }
std::vector< const reco::Muon * > m_ProbeMuons
double matchHlt(edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon *mu)
T sqrt(T t)
Definition: SSEVec.h:18
Definition: deltar.py:1
std::vector< const reco::Muon * > m_TightMuons
std::map< Control, MonitorElement * > m_ControlHistos
void L1TMuonDQMOffline::getTightMuons ( edm::Handle< reco::MuonCollection > &  muons,
const reco::Vertex vertex 
)
private

Definition at line 433 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, muon::isTightMuon(), kCtrlNTightVsAll, m_ControlHistos, m_TightMuons, and m_verbose.

Referenced by analyze().

433  {
434 
435  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting tight muons" << endl;
436  m_TightMuons.clear();
437  MuonCollection::const_iterator muonIt = muons->begin();
438  MuonCollection::const_iterator muonEnd = muons->end();
439 
440  for(; muonIt!=muonEnd; ++muonIt) {
441  if (muon::isTightMuon((*muonIt), vertex)) {
442  m_TightMuons.push_back(&(*muonIt));
443  }
444  }
445  m_ControlHistos[kCtrlNTightVsAll]->Fill(muons->size(), m_TightMuons.size());
446 }
std::vector< const reco::Muon * > m_TightMuons
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
std::map< Control, MonitorElement * > m_ControlHistos
double L1TMuonDQMOffline::matchHlt ( edm::Handle< trigger::TriggerEvent > &  triggerEvent,
const reco::Muon mu 
)
private

Definition at line 540 of file L1TMuonDQMOffline.cc.

References deltaR(), trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), m_hltConfig, m_trigIndices, m_trigProcess, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), and TriggerAnalyzer::trigObjs.

Referenced by getProbeMuons().

540  {
541 
542  double matchDeltaR = 9999;
543 
544  TriggerObjectCollection trigObjs = triggerEvent->getObjects();
545 
546  vector<int>::const_iterator trigIndexIt = m_trigIndices.begin();
547  vector<int>::const_iterator trigIndexEnd = m_trigIndices.end();
548 
549  for(; trigIndexIt!=trigIndexEnd; ++trigIndexIt) {
550  const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
551  const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt))-2;
552  const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex],"",m_trigProcess));
553 
554  if (hltFilterIndex < triggerEvent->sizeFilters()) {
555  const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
556  const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
557  const unsigned nTriggers = triggerVids.size();
558  for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
559  const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
560  double dRtmp = deltaR((*mu),trigObject);
561  if (dRtmp < matchDeltaR) matchDeltaR = dRtmp;
562  }
563  }
564  }
565  return matchDeltaR;
566 }
unsigned int size() const
number of trigger paths in trigger table
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:111
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
Definition: TriggerEvent.h:123
std::string m_trigProcess
const Vids & filterIds(trigger::size_type index) const
Definition: TriggerEvent.h:110
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
std::vector< size_type > Keys
HLTConfigProvider m_hltConfig
std::vector< int > Vids
std::vector< int > m_trigIndices

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> L1TMuonDQMOffline::m_BsInputTag
private

Definition at line 150 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::map<Control, MonitorElement*> L1TMuonDQMOffline::m_ControlHistos
private
std::vector<std::pair<int, QualLevel> > L1TMuonDQMOffline::m_cuts
private

Definition at line 125 of file L1TMuonDQMOffline.h.

Referenced by analyze(), bookEfficiencyHistos(), and L1TMuonDQMOffline().

std::vector<edm::ParameterSet> L1TMuonDQMOffline::m_cutsVPSet
private

Definition at line 146 of file L1TMuonDQMOffline.h.

Referenced by L1TMuonDQMOffline().

std::map<int, MonitorElement*> L1TMuonDQMOffline::m_EfficiencyDenEtaHistos
private

Definition at line 111 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::map<EtaRegion, MonitorElement*> L1TMuonDQMOffline::m_EfficiencyDenPtHistos
private

Definition at line 110 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::map<std::tuple<EffType, int, EtaRegion>, MonitorElement*> L1TMuonDQMOffline::m_EfficiencyDenVarHistos
private

Definition at line 109 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::map<std::pair<int, QualLevel>, MonitorElement*> L1TMuonDQMOffline::m_EfficiencyNumEtaHistos
private

Definition at line 108 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::map<std::tuple<EffType, int, EtaRegion, QualLevel>, MonitorElement*> L1TMuonDQMOffline::m_EfficiencyNumVarHistos
private

Definition at line 107 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::map<EffType, std::string> L1TMuonDQMOffline::m_effLabelStrings
private

Definition at line 135 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos(), and L1TMuonDQMOffline().

std::map<EffType, std::string> L1TMuonDQMOffline::m_effStrings
private

Definition at line 134 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos(), and L1TMuonDQMOffline().

const std::vector<EffType> L1TMuonDQMOffline::m_effTypes
private

Definition at line 128 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookEfficiencyHistos().

std::vector<double> L1TMuonDQMOffline::m_effVsEtaBins
private

Definition at line 158 of file L1TMuonDQMOffline.h.

Referenced by getHistBinsEff(), and L1TMuonDQMOffline().

std::vector<double> L1TMuonDQMOffline::m_effVsPhiBins
private

Definition at line 157 of file L1TMuonDQMOffline.h.

Referenced by getHistBinsEff(), and L1TMuonDQMOffline().

std::vector<double> L1TMuonDQMOffline::m_effVsPtBins
private

Definition at line 156 of file L1TMuonDQMOffline.h.

Referenced by getHistBinsEff(), and L1TMuonDQMOffline().

std::vector<double> L1TMuonDQMOffline::m_effVsVtxBins
private

Definition at line 159 of file L1TMuonDQMOffline.h.

Referenced by getHistBinsEff(), and L1TMuonDQMOffline().

const std::vector<EtaRegion> L1TMuonDQMOffline::m_etaRegions
private
std::map<EtaRegion, std::string> L1TMuonDQMOffline::m_etaStrings
private
edm::EDGetTokenT<l1t::MuonBxCollection> L1TMuonDQMOffline::m_GmtInputTag
private

Definition at line 148 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::string L1TMuonDQMOffline::m_HistFolder
private
HLTConfigProvider L1TMuonDQMOffline::m_hltConfig
private

Definition at line 93 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), dqmBeginRun(), and matchHlt().

BXVector<l1t::Muon> L1TMuonDQMOffline::m_L1tL1tMuons
private

Definition at line 123 of file L1TMuonDQMOffline.h.

std::vector<l1t::MuonBxCollection> L1TMuonDQMOffline::m_L1tMuons
private

Definition at line 121 of file L1TMuonDQMOffline.h.

float L1TMuonDQMOffline::m_maxGmtMuonDR
private

Definition at line 164 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and L1TMuonDQMOffline().

float L1TMuonDQMOffline::m_maxHltMuonDR
private

Definition at line 166 of file L1TMuonDQMOffline.h.

Referenced by L1TMuonDQMOffline().

float L1TMuonDQMOffline::m_minTagProbeDR
private

Definition at line 165 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons(), and L1TMuonDQMOffline().

std::vector<MuonGmtPair> L1TMuonDQMOffline::m_MuonGmtPairs
private

Definition at line 118 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and getMuonGmtPairs().

edm::EDGetTokenT<reco::MuonCollection> L1TMuonDQMOffline::m_MuonInputTag
private

Definition at line 147 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::vector<const reco::Muon*> L1TMuonDQMOffline::m_ProbeMuons
private

Definition at line 117 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and getProbeMuons().

PropagateToMuon L1TMuonDQMOffline::m_propagator
private

Definition at line 95 of file L1TMuonDQMOffline.h.

Referenced by dqmBeginRun(), and getMuonGmtPairs().

const std::vector<QualLevel> L1TMuonDQMOffline::m_qualLevelsRes
private

Definition at line 131 of file L1TMuonDQMOffline.h.

Referenced by analyze(), bookResolutionHistos(), and L1TMuonDQMOffline().

std::map<QualLevel, std::string> L1TMuonDQMOffline::m_qualStrings
private
std::vector<reco::MuonCollection> L1TMuonDQMOffline::m_RecoMuons
private

Definition at line 120 of file L1TMuonDQMOffline.h.

std::vector<reco::Muon> L1TMuonDQMOffline::m_RecoRecoMuons
private

Definition at line 122 of file L1TMuonDQMOffline.h.

double L1TMuonDQMOffline::m_recoToL1PtCutFactor
private

Definition at line 145 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::map<ResType, std::string> L1TMuonDQMOffline::m_resLabelStrings
private

Definition at line 137 of file L1TMuonDQMOffline.h.

Referenced by bookResolutionHistos(), and L1TMuonDQMOffline().

std::map<std::tuple<ResType, EtaRegion, QualLevel>, MonitorElement*> L1TMuonDQMOffline::m_ResolutionHistos
private

Definition at line 112 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookResolutionHistos().

std::map<ResType, std::string> L1TMuonDQMOffline::m_resStrings
private

Definition at line 136 of file L1TMuonDQMOffline.h.

Referenced by bookResolutionHistos(), and L1TMuonDQMOffline().

const std::vector<ResType> L1TMuonDQMOffline::m_resTypes
private

Definition at line 129 of file L1TMuonDQMOffline.h.

Referenced by analyze(), bookResolutionHistos(), and L1TMuonDQMOffline().

double L1TMuonDQMOffline::m_TagPtCut
private

Definition at line 144 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons(), and L1TMuonDQMOffline().

std::vector<const reco::Muon*> L1TMuonDQMOffline::m_TightMuons
private

Definition at line 116 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons(), and getTightMuons().

std::vector<int> L1TMuonDQMOffline::m_trigIndices
private

Definition at line 161 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), and matchHlt().

edm::EDGetTokenT<trigger::TriggerEvent> L1TMuonDQMOffline::m_trigInputTag
private

Definition at line 151 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::vector<std::string> L1TMuonDQMOffline::m_trigNames
private

Definition at line 155 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), and L1TMuonDQMOffline().

std::string L1TMuonDQMOffline::m_trigProcess
private

Definition at line 152 of file L1TMuonDQMOffline.h.

Referenced by dqmBeginRun(), L1TMuonDQMOffline(), and matchHlt().

edm::EDGetTokenT<edm::TriggerResults> L1TMuonDQMOffline::m_trigProcess_token
private

Definition at line 153 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

bool L1TMuonDQMOffline::m_useAtVtxCoord
private

Definition at line 163 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and L1TMuonDQMOffline().

bool L1TMuonDQMOffline::m_verbose
private
edm::EDGetTokenT<reco::VertexCollection> L1TMuonDQMOffline::m_VtxInputTag
private

Definition at line 149 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().