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

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

143 { }

Member Function Documentation

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

Definition at line 179 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.

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

301  {
302  if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Control Plot Histos" << endl;
303 
304  ibooker.setCurrentFolder(m_HistFolder+"/control_variables");
305 
306  m_ControlHistos[kCtrlMuonGmtDeltaR] = ibooker.book1D("MuonGmtDeltaR", "MuonGmtDeltaR; #DeltaR", 50, 0., 0.5);
307  m_ControlHistos[kCtrlNTightVsAll] = ibooker.book2D("NTightVsAll", "NTightVsAll; # muons; # tight muons", 20, -0.5, 19.5, 16, -0.5, 15.5);
308  m_ControlHistos[kCtrlNProbesVsTight] = ibooker.book2D("NProbesVsTight", "NProbesVsTight; # tight muons; # probe muons", 8, -0.5, 7.5, 8, -0.5, 7.5);
309 
310  m_ControlHistos[kCtrlTagPt] = ibooker.book1D("TagMuonPt", "TagMuonPt; p_{T}", 50, 0., 100.);
311  m_ControlHistos[kCtrlTagPhi] = ibooker.book1D("TagMuonPhi", "TagMuonPhi; #phi", 66, -3.3, 3.3);
312  m_ControlHistos[kCtrlTagEta] = ibooker.book1D("TagMuonEta", "TagMuonEta; #eta", 50, -2.5, 2.5);
313 
314  m_ControlHistos[kCtrlProbePt] = ibooker.book1D("ProbeMuonPt", "ProbeMuonPt; p_{T}", 50, 0., 100.);
315  m_ControlHistos[kCtrlProbePhi] = ibooker.book1D("ProbeMuonPhi", "ProbeMuonPhi; #phi", 66, -3.3, 3.3);
316  m_ControlHistos[kCtrlProbeEta] = ibooker.book1D("ProbeMuonEta", "ProbeMuonEta; #eta", 50, -2.5, 2.5);
317 
318  m_ControlHistos[kCtrlTagProbeDr] = ibooker.book1D("TagMuonProbeMuonDeltaR", "TagMuonProbeMuonDeltaR; #DeltaR", 50, 0.,5.0);
319  m_ControlHistos[kCtrlTagHltDr] = ibooker.book1D("TagMuonHltDeltaR", "TagMuonHltDeltaR;#DeltaR", 55, 0., 0.11);
320 }
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 323 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().

323  {
324  ibooker.setCurrentFolder(m_HistFolder+"/numerators_and_denominators");
325 
326  for(const auto var : m_effTypes) {
327  auto histBins = getHistBinsEff(var);
328  // histograms for eta variable get a special treatment
329  if (var == kEffEta) {
330  for (const auto cut : m_cuts) {
331  const auto gmtPtCut = cut.first;
332  const auto qualLevel = cut.second;
333  std::string name = "effDen_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut);
334  m_EfficiencyDenEtaHistos[gmtPtCut] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
335  name = "effNum_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_qualStrings[qualLevel];
336  m_histoKeyEffNumEtaType histoKeyEffNumEta = {gmtPtCut, qualLevel};
337  m_EfficiencyNumEtaHistos[histoKeyEffNumEta] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
338  }
339  } else {
340  for (const auto etaReg : m_etaRegions) {
341  // denominator histograms for pt variable get a special treatment
342  if (var == kEffPt) {
343  std::string name = "effDen_"+m_effStrings[var]+"_"+m_etaStrings[etaReg];
344  m_EfficiencyDenPtHistos[etaReg] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
345  } else {
346  for (const auto cut : m_cuts) {
347  const int gmtPtCut = cut.first;
348  std::string name = "effDen_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_etaStrings[etaReg];
349  m_histoKeyEffDenVarType histoKeyEffDenVar = {var, gmtPtCut, etaReg};
350  m_EfficiencyDenVarHistos[histoKeyEffDenVar] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
351  }
352  }
353  for (const auto cut : m_cuts) {
354  const auto gmtPtCut = cut.first;
355  const auto qualLevel = cut.second;
356  std::string name = "effNum_"+m_effStrings[var]+"_"+std::to_string(gmtPtCut)+"_"+m_etaStrings[etaReg]+"_"+m_qualStrings[qualLevel];
357  m_histoKeyEffNumVarType histoKeyEffNum = {var, gmtPtCut, etaReg, qualLevel};
358  m_EfficiencyNumVarHistos[histoKeyEffNum] = ibooker.book1D(name, name+";"+m_effLabelStrings[var], histBins.size()-1, &histBins[0]);
359  }
360  }
361  }
362  }
363 }
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 153 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().

153  {
154  //book histos
155  bookControlHistos(ibooker);
156  bookEfficiencyHistos(ibooker);
157  bookResolutionHistos(ibooker);
158 
159  vector<string>::const_iterator trigNamesIt = m_trigNames.begin();
160  vector<string>::const_iterator trigNamesEnd = m_trigNames.end();
161 
162  for (; trigNamesIt!=trigNamesEnd; ++trigNamesIt) {
163  TString tNameTmp = TString(*trigNamesIt); // use TString as it handles regex
164  TRegexp tNamePattern = TRegexp(tNameTmp,true);
165  int tIndex = -1;
166 
167  for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
168  TString tmpName = TString(m_hltConfig.triggerName(ipath));
169  if (tmpName.Contains(tNamePattern)) {
170  tIndex = int(ipath);
171  m_trigIndices.push_back(tIndex);
172  }
173  }
174  if (tIndex < 0 && m_verbose) cout << "[L1TMuonDQMOffline:] Warning: Could not find trigger " << (*trigNamesIt) << endl;
175  }
176 }
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 365 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().

365  {
366  if(m_verbose) cout << "[L1TMuonOffline:] Booking Resolution Plot Histos" << endl;
367  ibooker.setCurrentFolder(m_HistFolder+"/resolution");
368 
369  for (const auto var : m_resTypes) {
370  auto nbins = std::get<0>(getHistBinsRes(var));
371  auto xmin = std::get<1>(getHistBinsRes(var));
372  auto xmax = std::get<2>(getHistBinsRes(var));
373  for (const auto etaReg : m_etaRegions) {
374  for (const auto qualLevel : m_qualLevelsRes) {
375  m_histoKeyResType histoKeyRes = {var, etaReg, qualLevel};
376  std::string name = "resolution_"+m_resStrings[var]+"_"+m_etaStrings[etaReg]+"_"+m_qualStrings[qualLevel];
377  m_ResolutionHistos[histoKeyRes] = ibooker.book1D(name, name+";"+m_resLabelStrings[var], nbins, xmin, xmax);
378  }
379  }
380  }
381 }
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 145 of file L1TMuonDQMOffline.cc.

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

145  {
146  if (m_verbose) cout << "[L1TMuonDQMOffline:] Called beginRun." << endl;
147  bool changed = true;
148  m_hltConfig.init(run,iSetup,m_trigProcess,changed);
149  m_propagator.init(iSetup);
150 }
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 565 of file L1TMuonDQMOffline.cc.

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

Referenced by bookEfficiencyHistos().

565  {
566  if (eff == kEffPt) {
567  std::vector<float> effVsPtBins(m_effVsPtBins.begin(), m_effVsPtBins.end());
568  return effVsPtBins;
569  }
570  if (eff == kEffPhi) {
571  std::vector<float> effVsPhiBins(m_effVsPhiBins.begin(), m_effVsPhiBins.end());
572  return effVsPhiBins;
573  }
574  if (eff == kEffEta) {
575  std::vector<float> effVsEtaBins(m_effVsEtaBins.begin(), m_effVsEtaBins.end());
576  return effVsEtaBins;
577  }
578  if (eff == kEffVtx) {
579  std::vector<float> effVsVtxBins(m_effVsVtxBins.begin(), m_effVsVtxBins.end());
580  return effVsVtxBins;
581  }
582  return {0., 1.};
583 }
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 585 of file L1TMuonDQMOffline.cc.

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

Referenced by bookResolutionHistos().

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

Definition at line 502 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().

502  {
503 
504  m_MuonGmtPairs.clear();
505  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting muon GMT pairs" << endl;
506 
507  vector<const reco::Muon*>::const_iterator probeMuIt = m_ProbeMuons.begin();
508  vector<const reco::Muon*>::const_iterator probeMuEnd = m_ProbeMuons.end();
509 
511  l1t::MuonBxCollection::const_iterator gmtEnd = gmtCands->end(0);
512 
513  for (; probeMuIt!=probeMuEnd; ++probeMuIt) {
514  MuonGmtPair pairBestCand((*probeMuIt), nullptr, m_propagator, m_useAtVtxCoord);
515 
516  // Fill the control histograms with the probe muon kinematic variables used
517  m_ControlHistos[kCtrlProbeEta]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffEta));
518  m_ControlHistos[kCtrlProbePhi]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffPhi));
519  m_ControlHistos[kCtrlProbePt]->Fill(pairBestCand.getVar(L1TMuonDQMOffline::kEffPt));
520 
521  gmtIt = gmtCands->begin(0); // use only on L1T muons from BX 0
522 
523  for(; gmtIt!=gmtEnd; ++gmtIt) {
524  MuonGmtPair pairTmpCand((*probeMuIt), &(*gmtIt), m_propagator, m_useAtVtxCoord);
525 
526  if ( (pairTmpCand.dR() < m_maxGmtMuonDR) && (pairTmpCand.dR() < pairBestCand.dR() ) ) {
527  pairBestCand = pairTmpCand;
528  }
529 
530  }
531  m_MuonGmtPairs.push_back(pairBestCand);
532  m_ControlHistos[kCtrlMuonGmtDeltaR]->Fill(pairBestCand.dR());
533  }
534 }
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 384 of file L1TMuonDQMOffline.cc.

References edm::HandleBase::isValid().

Referenced by analyze().

384  {
385  unsigned int nVtx = 0;
386 
387  if (vertex.isValid()) {
388  for (const auto vertexIt : *vertex) {
389  if (vertexIt.isValid() && !vertexIt.isFake()) {
390  ++nVtx;
391  }
392  }
393  }
394  return nVtx;
395 }
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 398 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().

399  {
400  Vertex::Point posVtx;
401  Vertex::Error errVtx;
402 
403  bool hasPrimaryVertex = false;
404 
405  if (vertex.isValid()) {
406  vector<Vertex>::const_iterator vertexIt = vertex->begin();
407  vector<Vertex>::const_iterator vertexEnd = vertex->end();
408 
409  for (;vertexIt!=vertexEnd;++vertexIt) {
410  if (vertexIt->isValid() && !vertexIt->isFake()) {
411  posVtx = vertexIt->position();
412  errVtx = vertexIt->error();
413  hasPrimaryVertex = true;
414  break;
415  }
416  }
417  }
418 
419  if ( !hasPrimaryVertex ) {
420  posVtx = beamSpot->position();
421  errVtx(0,0) = beamSpot->BeamWidthX();
422  errVtx(1,1) = beamSpot->BeamWidthY();
423  errVtx(2,2) = beamSpot->sigmaZ();
424  }
425  const Vertex primaryVertex(posVtx,errVtx);
426  return primaryVertex;
427 }
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 446 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().

447  {
448 
449  if (m_verbose) cout << "[L1TMuonDQMOffline:] getting probe muons" << endl;
450  m_ProbeMuons.clear();
451  std::vector<const reco::Muon*> tagMuonsInHist;
452 
453  tagMuonsInHist.clear();
454 
455  vector<const reco::Muon*>::const_iterator probeCandIt = m_TightMuons.begin();
456  vector<const reco::Muon*>::const_iterator tightMuonsEnd = m_TightMuons.end();
457 
458  for (; probeCandIt!=tightMuonsEnd; ++probeCandIt) {
459  bool isProbe = false;
460  vector<const reco::Muon*>::const_iterator tagCandIt = m_TightMuons.begin();
461  float deltar = 0.;
462 
463  for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
464  bool tagMuonAlreadyInHist = false;
465  bool tagHasTrig = false;
466  float eta = (*tagCandIt)->eta();
467  float phi = (*tagCandIt)->phi();
468  float pt = (*tagCandIt)->pt();
469  float dEta = eta - (*probeCandIt)->eta();
470  float dPhi = phi - (*probeCandIt)->phi();
471  deltar = sqrt(dEta*dEta + dPhi*dPhi);
472 
473  if ( (*tagCandIt) == (*probeCandIt) || deltar<m_minTagProbeDR ) continue; // CB has a little bias for closed-by muons
474  auto matchHltDeltaR = matchHlt(trigEvent,(*tagCandIt));
475  tagHasTrig = (matchHltDeltaR < m_maxHltMuonDR) && (pt > m_TagPtCut);
476  isProbe |= tagHasTrig;
477  if (tagHasTrig) {
478  if (std::distance(m_TightMuons.begin(), m_TightMuons.end()) > 2 ) {
479  for (vector<const reco::Muon*>::const_iterator tagMuonsInHistIt = tagMuonsInHist.begin(); tagMuonsInHistIt!=tagMuonsInHist.end(); ++tagMuonsInHistIt) {
480  if ( (*tagCandIt) == (*tagMuonsInHistIt) ) {
481  tagMuonAlreadyInHist = true;
482  break;
483  }
484  }
485  if (tagMuonAlreadyInHist == false) tagMuonsInHist.push_back((*tagCandIt));
486  }
487  if (tagMuonAlreadyInHist == false) {
488  m_ControlHistos[kCtrlTagEta]->Fill(eta);
489  m_ControlHistos[kCtrlTagPhi]->Fill(phi);
490  m_ControlHistos[kCtrlTagPt]->Fill(pt);
491  m_ControlHistos[kCtrlTagProbeDr]->Fill(deltar);
492  m_ControlHistos[kCtrlTagHltDr]->Fill(matchHltDeltaR);
493  }
494  }
495  }
496  if (isProbe) m_ProbeMuons.push_back((*probeCandIt));
497  }
499 }
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 430 of file L1TMuonDQMOffline.cc.

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

Referenced by analyze().

430  {
431 
432  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting tight muons" << endl;
433  m_TightMuons.clear();
434  MuonCollection::const_iterator muonIt = muons->begin();
435  MuonCollection::const_iterator muonEnd = muons->end();
436 
437  for(; muonIt!=muonEnd; ++muonIt) {
438  if (muon::isTightMuon((*muonIt), vertex)) {
439  m_TightMuons.push_back(&(*muonIt));
440  }
441  }
442  m_ControlHistos[kCtrlNTightVsAll]->Fill(muons->size(), m_TightMuons.size());
443 }
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 537 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().

537  {
538 
539  double matchDeltaR = 9999;
540 
541  TriggerObjectCollection trigObjs = triggerEvent->getObjects();
542 
543  vector<int>::const_iterator trigIndexIt = m_trigIndices.begin();
544  vector<int>::const_iterator trigIndexEnd = m_trigIndices.end();
545 
546  for(; trigIndexIt!=trigIndexEnd; ++trigIndexIt) {
547  const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
548  const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt))-2;
549  const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex],"",m_trigProcess));
550 
551  if (hltFilterIndex < triggerEvent->sizeFilters()) {
552  const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
553  const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
554  const unsigned nTriggers = triggerVids.size();
555  for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
556  const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
557  double dRtmp = deltaR((*mu),trigObject);
558  if (dRtmp < matchDeltaR) matchDeltaR = dRtmp;
559  }
560  }
561  }
562  return matchDeltaR;
563 }
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().