CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
FourVectorHLTOffline.h
Go to the documentation of this file.
1 #ifndef FOURVECTORHLTOFFLINE_H
2 #define FOURVECTORHLTOFFLINE_H
3 // -*- C++ -*-
4 //
5 // Package: FourVectorHLTOffline
6 // Class: FourVectorHLTOffline
7 //
16 //
17 // Original Author: Jeffrey Berryhill
18 // Created: June 2008
19 // Rewritten by: Vladimir Rekovic
20 // Date: May 2009
21 //
22 // $Id: FourVectorHLTOffline.h,v 1.64 2010/10/26 04:55:28 wmtan Exp $
23 //
24 //
25 // system include files
26 #include <memory>
27 #include <unistd.h>
28 
29 
30 // user include files
31 #include "FWCore/Framework/interface/Frameworkfwd.h"
32 #include "FWCore/Framework/interface/EDAnalyzer.h"
33 
34 #include "FWCore/Framework/interface/Event.h"
35 #include "FWCore/Framework/interface/MakerMacros.h"
36 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
37 #include "DataFormats/Common/interface/TriggerResults.h"
38 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
39 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
40 
41 #include "DQMServices/Core/interface/DQMStore.h"
42 #include "DQMServices/Core/interface/MonitorElement.h"
43 #include "FWCore/ServiceRegistry/interface/Service.h"
44 #include "FWCore/MessageLogger/interface/MessageLogger.h"
45 #include "FWCore/ParameterSet/interface/ParameterSet.h"
46 
47 // added VR
48 #include "FWCore/Framework/interface/EDAnalyzer.h"
49 #include "DataFormats/Common/interface/Handle.h"
50 #include "FWCore/Framework/interface/Run.h"
51 #include "FWCore/Framework/interface/MakerMacros.h"
52 #include "FWCore/Framework/interface/ESHandle.h"
53 
54 #include "FWCore/ParameterSet/interface/ParameterSet.h"
55 
56 #include "DataFormats/HLTReco/interface/TriggerObject.h"
57 #include "DataFormats/Common/interface/TriggerResults.h"
58 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
59 #include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
60 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
61 
62 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
63 #include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"
64 #include "DataFormats/EgammaCandidates/interface/Photon.h"
65 #include "DataFormats/TrackReco/interface/TrackFwd.h"
66 #include "DataFormats/TrackReco/interface/Track.h"
67 #include "DataFormats/TrackReco/interface/TrackExtra.h"
68 #include "DataFormats/MuonReco/interface/MuonFwd.h"
69 #include "DataFormats/MuonReco/interface/Muon.h"
70 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
71 #include "DataFormats/JetReco/interface/CaloJet.h"
72 #include "DataFormats/TauReco/interface/CaloTauFwd.h"
73 #include "DataFormats/TauReco/interface/CaloTau.h"
74 #include "DataFormats/METReco/interface/CaloMETCollection.h"
75 #include "DataFormats/METReco/interface/CaloMET.h"
76 #include "DataFormats/BTauReco/interface/JetTag.h"
77 #include "DataFormats/MuonReco/interface/MuonSelectors.h"
78 #include "RecoJets/JetProducers/interface/JetIDHelper.h"
79 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
80 
81 #include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
82 #include "DataFormats/JetReco/interface/JetFloatAssociation.h"
83 
84 #include "DataFormats/TauReco/interface/PFTau.h"
85 #include "DataFormats/TauReco/interface/PFTauFwd.h"
86 #include "DataFormats/TauReco/interface/PFTauDiscriminator.h"
87 
88 /* MC
89 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
90 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
91 #include "DataFormats/JetReco/interface/GenJetCollection.h"
92 #include "DataFormats/METReco/interface/GenMETCollection.h"
93 #include "DataFormats/METReco/interface/GenMET.h"
94 */
95 #include "DataFormats/Math/interface/deltaR.h"
96 #include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
97 
98 #include "DQMServices/Core/interface/MonitorElement.h"
99 
100 
101 
102 #include <iostream>
103 #include <fstream>
104 #include <vector>
105 
106 namespace edm {
107  class TriggerNames;
108 }
109 
110 typedef std::multimap<float,int> fimmap ;
111 typedef std::set<fimmap , std::less<fimmap> > mmset;
112 
113 class FourVectorHLTOffline : public edm::EDAnalyzer {
114 
115  public:
116  explicit FourVectorHLTOffline(const edm::ParameterSet&);
117  ~FourVectorHLTOffline();
118 
119  void cleanDRMatchSet(mmset& tempSet);
120 
121  edm::Handle<trigger::TriggerEvent> fTriggerObj;
122  edm::Handle<edm::TriggerResults> fTriggerResults;
123  edm::Handle<reco::BeamSpot> fBeamSpotHandle;
124 
125  private:
126  virtual void beginJob() ;
127  virtual void analyze(const edm::Event&, const edm::EventSetup&);
128  virtual void endJob() ;
129 
130  // BeginRun
131  void beginRun(const edm::Run& run, const edm::EventSetup& c);
132 
133  // EndRun
134  void endRun(const edm::Run& run, const edm::EventSetup& c);
135  void fillHltMatrix(const edm::TriggerNames & triggerNames);
136  void setupHltMatrix(const std::string& label, std::vector<std::string> & paths);
137 
138  void setupHltLsPlots();
139  void setupHltBxPlots();
140  void countHLTPathHitsEndLumiBlock(const int & lumi);
141  void countHLTGroupHitsEndLumiBlock(const int & lumi);
142  void countHLTGroupL1HitsEndLumiBlock(const int & lumi);
143  void countHLTGroupBXHitsEndLumiBlock(const int & lumi);
144  int getTriggerTypeParsePathName(const std::string & pathname);
145  const std::string getL1ConditionModuleName(const std::string & pathname);
146  bool hasL1Passed(const std::string & pathname, const edm::TriggerNames & triggerNames);
147  bool hasHLTPassed(const std::string & pathname, const edm::TriggerNames& triggerNames);
148  int getHltThresholdFromName(const std::string & pathname);
149 
150  void selectMuons(const edm::Handle<reco::MuonCollection> & muonHandle);
151  bool isVBTFMuon(const reco::Muon& muon);
152  void selectElectrons(const edm::Event& iEvent, const edm::EventSetup& iSetup, const edm::Handle<reco::GsfElectronCollection> & eleHandle);
153  void selectPhotons(const edm::Handle<reco::PhotonCollection> & phoHandle);
154  void selectJets(const edm::Event& iEvent,const edm::Handle<reco::CaloJetCollection> & jetHandle);
155  void selectMet(const edm::Handle<reco::CaloMETCollection> & metHandle);
156  void selectTaus(const edm::Event& iEvent);
157  void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c);
158  void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c);
159 
160  // ----------member data ---------------------------
161  int nev_;
162  DQMStore * dbe_;
163  bool fLumiFlag;
164  bool fIsSetup;
165 
166  // JetID helper
167  reco::helper::JetIDHelper *jetID;
168 
169  // Tau discriminators
170  std::string tauDscrmtrLabel1_;
171  std::string tauDscrmtrLabel2_;
172  std::string tauDscrmtrLabel3_;
173 
174  MonitorElement* ME_HLTAll_LS;
175  MonitorElement* ME_HLT_BX;
176  MonitorElement* ME_HLT_CUSTOM_BX;
177  std::vector<MonitorElement*> v_ME_HLTAll_LS;
178  std::vector<MonitorElement*> v_ME_Total_BX;
179  std::vector<MonitorElement*> v_ME_Total_BX_Norm;
180 
181  std::string pathsSummaryFolder_;
182  std::string pathsSummaryHLTCorrelationsFolder_;
183  std::string pathsSummaryFilterEfficiencyFolder_;
184  std::string pathsSummaryFilterCountsFolder_;
185  std::string pathsSummaryHLTPathsPerLSFolder_;
186  std::string pathsIndividualHLTPathsPerLSFolder_;
187  std::string pathsSummaryHLTPathsPerBXFolder_;
188  std::string fCustomBXPath;
189 
190  std::vector<std::string> fGroupName;
191 
192  reco::MuonCollection * fSelectedMuons;
193  edm::Handle<reco::MuonCollection> fSelMuonsHandle;
194 
195  reco::GsfElectronCollection * fSelectedElectrons;
196  edm::Handle<reco::GsfElectronCollection> fSelElectronsHandle;
197 
198  reco::PhotonCollection * fSelectedPhotons;
199  edm::Handle<reco::PhotonCollection> fSelPhotonsHandle;
200 
201  reco::CaloJetCollection * fSelectedJets;
202  edm::Handle<reco::CaloJetCollection> fSelJetsHandle;
203 
204  reco::CaloMETCollection * fSelectedMet;
205  edm::Handle<reco::CaloMETCollection> fSelMetHandle;
206 
207  //reco::CaloTauCollection * fSelectedTaus;
208  //edm::Handle<reco::CaloTauCollection> fSelTausHandle;
209  reco::PFTauCollection * fSelectedTaus;
210  edm::Handle<reco::PFTauCollection> fSelTausHandle;
211 
212  unsigned int nLS_;
213  double LSsize_ ;
214  double thresholdFactor_ ;
215  unsigned int referenceBX_;
216  unsigned int Nbx_;
217 
218  bool plotAll_;
219  bool doCombineRuns_;
220  int currentRun_;
221 
222  unsigned int nBins_;
223  unsigned int nBinsDR_;
224  unsigned int nBins2D_;
225  unsigned int nBinsOneOverEt_;
226  double ptMin_ ;
227  double ptMax_ ;
228  double dRMax_ ;
229 
230  double electronEtaMax_;
231  double electronEtMin_;
232  double electronDRMatch_;
233  double electronL1DRMatch_;
234  double muonEtaMax_;
235  double muonEtMin_;
236  double muonDRMatch_;
237  double muonL1DRMatch_;
238  double tauEtaMax_;
239  double tauEtMin_;
240  double tauDRMatch_;
241  double tauL1DRMatch_;
242  double jetEtaMax_;
243  double jetEtMin_;
244  double jetDRMatch_;
245  double jetL1DRMatch_;
246  double bjetEtaMax_;
247  double bjetEtMin_;
248  double bjetDRMatch_;
249  double bjetL1DRMatch_;
250  double photonEtaMax_;
251  double photonEtMin_;
252  double photonDRMatch_;
253  double photonL1DRMatch_;
254  double trackEtaMax_;
255  double trackEtMin_;
256  double trackDRMatch_;
257  double trackL1DRMatch_;
258  double metEtaMax_;
259  double metMin_;
260  double metDRMatch_;
261  double metL1DRMatch_;
262  double htEtaMax_;
263  double htMin_;
264  double htDRMatch_;
265  double htL1DRMatch_;
266  double sumEtMin_;
267 
268  // Muon quality cuts
270  double dxyCut_;
271  double normalizedChi2Cut_;
272  int trackerHitsCut_;
273  int pixelHitsCut_;
274  int muonHitsCut_;
275  bool isAlsoTrackerMuon_;
276  int nMatchesCut_;
277 
278  // Electron quality cuts
280  float eleMaxOver3x3_;
281  // Ecal Barrel
282  float dr03TkSumPtEB_;
283  float dr04EcalRecHitSumEtEB_;
284  float dr04HcalTowerSumEtEB_;
285  float hadronicOverEmEB_;
286  float deltaPhiSuperClusterTrackAtVtxEB_;
287  float deltaEtaSuperClusterTrackAtVtxEB_;
288  float sigmaIetaIetaEB_;
289  //spikes
290  float sigmaIetaIetaSpikesEB_;
291  // Ecal Endcap
292  float dr03TkSumPtEC_;
293  float dr04EcalRecHitSumEtEC_;
294  float dr04HcalTowerSumEtEC_;
295  float hadronicOverEmEC_;
296  float deltaPhiSuperClusterTrackAtVtxEC_;
297  float deltaEtaSuperClusterTrackAtVtxEC_;
298  float sigmaIetaIetaEC_;
299  //spikes
300  float sigmaIetaIetaSpikesEC_;
301 
302  // Jet quality cuts
304  float emEnergyFractionJet_;
305  float fHPDJet_;
306  int n90Jet_;
307 
308  std::vector<std::pair<std::string, std::string> > custompathnamepairs_;
309 
310  std::vector <std::vector <std::string> > triggerFilters_;
311  std::vector <std::vector <uint> > triggerFilterIndices_;
312  std::vector <std::pair<std::string, float> > fPathTempCountPair;
313  std::vector <std::pair<std::string, std::vector<int> > > fPathBxTempCountPair;
314  std::vector <std::pair<std::string, float> > fGroupTempCountPair;
315  std::vector <std::pair<std::string, float> > fGroupL1TempCountPair;
316  std::vector <std::pair<std::string, std::vector<std::string> > > fGroupNamePathsPair;
317 
318  std::vector<std::string> specialPaths_;
319 
320  std::string dirname_;
321  std::string processname_;
322  std::string muonRecoCollectionName_;
323  bool monitorDaemon_;
324  int theHLTOutputType;
325  edm::InputTag triggerSummaryLabel_;
326  edm::InputTag triggerResultsLabel_;
327  edm::InputTag recHitsEBTag_, recHitsEETag_;
328  HLTConfigProvider hltConfig_;
329  // data across paths
330  MonitorElement* scalersSelect;
331  // helper class to store the data path
332 
333  class PathInfo {
334 
335  PathInfo():
336  pathIndex_(-1), denomPathName_("unset"), pathName_("unset"), l1pathName_("unset"), filterName_("unset"), processName_("unset"), objectType_(-1) {};
337 
338  public:
339 
340  void setFilterHistos(MonitorElement* const filters)
341  {
342  filters_ = filters;
343  }
344 
345  void setHistos(
346 
347  MonitorElement* const NOn,
348  MonitorElement* const onEtOn,
349  MonitorElement* const onOneOverEtOn,
350  MonitorElement* const onEtavsonPhiOn,
351  MonitorElement* const NOff,
352  MonitorElement* const offEtOff,
353  MonitorElement* const offEtavsoffPhiOff,
354  MonitorElement* const NL1,
355  MonitorElement* const l1EtL1,
356  MonitorElement* const l1Etavsl1PhiL1,
357  MonitorElement* const NL1On,
358  MonitorElement* const l1EtL1On,
359  MonitorElement* const l1Etavsl1PhiL1On,
360  MonitorElement* const NL1Off,
361  MonitorElement* const offEtL1Off,
362  MonitorElement* const offEtavsoffPhiL1Off,
363  MonitorElement* const NOnOff,
364  MonitorElement* const offEtOnOff,
365  MonitorElement* const offEtavsoffPhiOnOff,
366  MonitorElement* const NL1OnUM,
367  MonitorElement* const l1EtL1OnUM,
368  MonitorElement* const l1Etavsl1PhiL1OnUM,
369  MonitorElement* const NL1OffUM,
370  MonitorElement* const offEtL1OffUM,
371  MonitorElement* const offEtavsoffPhiL1OffUM,
372  MonitorElement* const NOnOffUM,
373  MonitorElement* const offEtOnOffUM,
374  MonitorElement* const offEtavsoffPhiOnOffUM,
375  MonitorElement* const offDRL1Off,
376  MonitorElement* const offDROnOff,
377  MonitorElement* const l1DRL1On)
378  {
379 
380  NOn_ = NOn;
381  onEtOn_ = onEtOn;
382  onOneOverEtOn_ = onOneOverEtOn;
383  onEtavsonPhiOn_ = onEtavsonPhiOn;
384  NOff_ = NOff;
385  offEtOff_ = offEtOff;
386  offEtavsoffPhiOff_ = offEtavsoffPhiOff;
387  NL1_ = NL1;
388  l1EtL1_ = l1EtL1;
389  l1Etavsl1PhiL1_ = l1Etavsl1PhiL1;
390  NL1On_ = NL1On;
391  l1EtL1On_ = l1EtL1On;
392  l1Etavsl1PhiL1On_ = l1Etavsl1PhiL1On;
393  NL1Off_ = NL1Off;
394  offEtL1Off_ = offEtL1Off;
395  offEtavsoffPhiL1Off_ = offEtavsoffPhiL1Off;
396  NOnOff_ = NOnOff;
397  offEtOnOff_ = offEtOnOff;
398  offEtavsoffPhiOnOff_ = offEtavsoffPhiOnOff;
399  NL1OnUM_ = NL1OnUM;
400  l1EtL1OnUM_ = l1EtL1OnUM;
401  l1Etavsl1PhiL1OnUM_ = l1Etavsl1PhiL1OnUM;
402  NL1OffUM_ = NL1OffUM;
403  offEtL1OffUM_ = offEtL1OffUM;
404  offEtavsoffPhiL1OffUM_ = offEtavsoffPhiL1OffUM;
405  NOnOffUM_ = NOnOffUM;
406  offEtOnOffUM_ = offEtOnOffUM;
407  offEtavsoffPhiOnOffUM_ = offEtavsoffPhiOnOffUM;
408  offDRL1Off_ = offDRL1Off;
409  offDROnOff_ = offDROnOff;
410  l1DRL1On_ = l1DRL1On;
411 
412  }
413 
414  MonitorElement * getNOnHisto() {
415  return NOn_;
416  }
417  MonitorElement * getOnEtOnHisto() {
418  return onEtOn_;
419  }
420  MonitorElement * getOnOneOverEtOnHisto() {
421  return onOneOverEtOn_;
422  }
423  MonitorElement * getOnEtaVsOnPhiOnHisto() {
424  return onEtavsonPhiOn_;
425  }
426  MonitorElement * getNOffHisto() {
427  return NOff_;
428  }
429  MonitorElement * getOffEtOffHisto() {
430  return offEtOff_;
431  }
432  MonitorElement * getOffEtaVsOffPhiOffHisto() {
433  return offEtavsoffPhiOff_;
434  }
435  MonitorElement * getNL1Histo() {
436  return NL1_;
437  }
438  MonitorElement * getL1EtL1Histo() {
439  return l1EtL1_;
440  }
441  MonitorElement * getL1EtaVsL1PhiL1Histo() {
442  return l1Etavsl1PhiL1_;
443  }
444  MonitorElement * getNL1OnHisto() {
445  return NL1On_;
446  }
447  MonitorElement * getL1EtL1OnHisto() {
448  return l1EtL1On_;
449  }
450  MonitorElement * getL1EtaVsL1PhiL1OnHisto() {
451  return l1Etavsl1PhiL1On_;
452  }
453  MonitorElement * getNL1OffHisto() {
454  return NL1Off_;
455  }
456  MonitorElement * getOffEtL1OffHisto() {
457  return offEtL1Off_;
458  }
459  MonitorElement * getOffEtaVsOffPhiL1OffHisto() {
460  return offEtavsoffPhiL1Off_;
461  }
462  MonitorElement * getNOnOffHisto() {
463  return NOnOff_;
464  }
465  MonitorElement * getOffEtOnOffHisto() {
466  return offEtOnOff_;
467  }
468  MonitorElement * getOffEtaVsOffPhiOnOffHisto() {
469  return offEtavsoffPhiOnOff_;
470  }
471  MonitorElement * getNL1OnUMHisto() {
472  return NL1OnUM_;
473  }
474  MonitorElement * getL1EtL1OnUMHisto() {
475  return l1EtL1OnUM_;
476  }
477  MonitorElement * getL1EtaVsL1PhiL1OnUMHisto() {
478  return l1Etavsl1PhiL1OnUM_;
479  }
480  MonitorElement * getNL1OffUMHisto() {
481  return NL1OffUM_;
482  }
483  MonitorElement * getOffEtL1OffUMHisto() {
484  return offEtL1OffUM_;
485  }
486  MonitorElement * getOffEtaVsOffPhiL1OffUMHisto() {
487  return offEtavsoffPhiL1OffUM_;
488  }
489  MonitorElement * getNOnOffUMHisto() {
490  return NOnOffUM_;
491  }
492  MonitorElement * getOffEtOnOffUMHisto() {
493  return offEtOnOffUM_;
494  }
495  MonitorElement * getOffEtaVsOffPhiOnOffUMHisto() {
496  return offEtavsoffPhiOnOffUM_;
497  }
498  MonitorElement * getOffDRL1OffHisto() {
499  return offDRL1Off_;
500  }
501  MonitorElement * getOffDROnOffHisto() {
502  return offDROnOff_;
503  }
504  MonitorElement * getL1DROnL1Histo() {
505  return l1DRL1On_;
506  }
507  MonitorElement * getFiltersHisto() {
508  return filters_;
509  }
510  const std::string getLabel(void ) const {
511  return filterName_;
512  }
513  void setLabel(std::string labelName){
514  filterName_ = labelName;
515  return;
516  }
517  const std::string & getPath(void ) const {
518  return pathName_;
519  }
520  const std::string & getl1Path(void ) const {
521  return l1pathName_;
522  }
523  const int getL1ModuleIndex(void ) const {
524  return l1ModuleIndex_;
525  }
526  const std::string & getDenomPath(void ) const {
527  return denomPathName_;
528  }
529  const std::string & getProcess(void ) const {
530  return processName_;
531  }
532  const int getObjectType(void ) const {
533  return objectType_;
534  }
535 
536  const edm::InputTag getTag(void) const{
537  edm::InputTag tagName(filterName_,"",processName_);
538  return tagName;
539  }
540 
541  ~PathInfo() {};
542 
543  PathInfo(std::string denomPathName, std::string pathName, std::string l1pathName, int l1ModuleIndex, std::string filterName, std::string processName, size_t type, float ptmin, float ptmax, float hltThreshold, float l1Threshold):
544 
545  denomPathName_(denomPathName),
546  pathName_(pathName),
547  l1pathName_(l1pathName),
548  l1ModuleIndex_(l1ModuleIndex),
549  filterName_(filterName),
550  processName_(processName),
551  objectType_(type),
552  NOn_(0), onEtOn_(0), onOneOverEtOn_(0), onEtavsonPhiOn_(0),
553  NOff_(0), offEtOff_(0), offEtavsoffPhiOff_(0),
554  NL1_(0), l1EtL1_(0), l1Etavsl1PhiL1_(0),
555  NL1On_(0), l1EtL1On_(0), l1Etavsl1PhiL1On_(0),
556  NL1Off_(0), offEtL1Off_(0), offEtavsoffPhiL1Off_(0),
557  NOnOff_(0), offEtOnOff_(0), offEtavsoffPhiOnOff_(0),
558  NL1OnUM_(0), l1EtL1OnUM_(0), l1Etavsl1PhiL1OnUM_(0),
559  NL1OffUM_(0), offEtL1OffUM_(0), offEtavsoffPhiL1OffUM_(0),
560  NOnOffUM_(0), offEtOnOffUM_(0), offEtavsoffPhiOnOffUM_(0),
561  offDRL1Off_(0), offDROnOff_(0), l1DRL1On_(0), filters_(0),
562  ptmin_(ptmin), ptmax_(ptmax),
563  hltThreshold_(hltThreshold), l1Threshold_(l1Threshold)
564 
565  {
566  };
567 
568  PathInfo(std::string denomPathName, std::string pathName, std::string l1pathName, std::string filterName, std::string processName, size_t type,
569  MonitorElement *NOn,
570  MonitorElement *onEtOn,
571  MonitorElement *onOneOverEtOn,
572  MonitorElement *onEtavsonPhiOn,
573  MonitorElement *NOff,
574  MonitorElement *offEtOff,
575  MonitorElement *offEtavsoffPhiOff,
576  MonitorElement *NL1,
577  MonitorElement *l1EtL1,
578  MonitorElement *l1Etavsl1PhiL1,
579  MonitorElement *NL1On,
580  MonitorElement *l1EtL1On,
581  MonitorElement *l1Etavsl1PhiL1On,
582  MonitorElement *NL1Off,
583  MonitorElement *offEtL1Off,
584  MonitorElement *offEtavsoffPhiL1Off,
585  MonitorElement *NOnOff,
586  MonitorElement *offEtOnOff,
587  MonitorElement *offEtavsoffPhiOnOff,
588  MonitorElement *NL1OnUM,
589  MonitorElement *l1EtL1OnUM,
590  MonitorElement *l1Etavsl1PhiL1OnUM,
591  MonitorElement *NL1OffUM,
592  MonitorElement *offEtL1OffUM,
593  MonitorElement *offEtavsoffPhiL1OffUM,
594  MonitorElement *NOnOffUM,
595  MonitorElement *offEtOnOffUM,
596  MonitorElement *offEtavsoffPhiOnOffUM,
597  MonitorElement *offDRL1Off,
598  MonitorElement *offDROnOff,
599  MonitorElement *l1DRL1On,
600  MonitorElement *filters,
601  float ptmin, float ptmax
602  ):
603 
604  denomPathName_(denomPathName),
605  pathName_(pathName), l1pathName_(l1pathName),
606  filterName_(filterName), processName_(processName), objectType_(type),
607  NOn_(NOn), onEtOn_(onEtOn), onOneOverEtOn_(onOneOverEtOn), onEtavsonPhiOn_(onEtavsonPhiOn),
608  NOff_(NOff), offEtOff_(offEtOff), offEtavsoffPhiOff_(offEtavsoffPhiOff),
609  NL1_(NL1), l1EtL1_(l1EtL1), l1Etavsl1PhiL1_(l1Etavsl1PhiL1),
610  NL1On_(NL1On), l1EtL1On_(l1EtL1On), l1Etavsl1PhiL1On_(l1Etavsl1PhiL1On),
611  NL1Off_(NL1Off), offEtL1Off_(offEtL1Off), offEtavsoffPhiL1Off_(offEtavsoffPhiL1Off),
612  NOnOff_(NOnOff), offEtOnOff_(offEtOnOff), offEtavsoffPhiOnOff_(offEtavsoffPhiOnOff),
613  NL1OnUM_(NL1OnUM), l1EtL1OnUM_(l1EtL1OnUM), l1Etavsl1PhiL1OnUM_(l1Etavsl1PhiL1OnUM),
614  NL1OffUM_(NL1OffUM), offEtL1OffUM_(offEtL1OffUM), offEtavsoffPhiL1OffUM_(offEtavsoffPhiL1OffUM),
615  NOnOffUM_(NOnOffUM), offEtOnOffUM_(offEtOnOffUM), offEtavsoffPhiOnOffUM_(offEtavsoffPhiOnOffUM),
616  offDRL1Off_(offDRL1Off),
617  offDROnOff_(offDROnOff),
618  l1DRL1On_(l1DRL1On),
619  filters_(filters),
620  ptmin_(ptmin), ptmax_(ptmax)
621  {
622  };
623 
624  bool operator==(const std::string& v)
625  {
626  return v==filterName_;
627  }
628 
629  bool operator!=(const std::string& v)
630  {
631  return v!=filterName_;
632  }
633 
634  float getPtMin() const { return ptmin_; }
635  float getPtMax() const { return ptmax_; }
636  float getHltThreshold() const { return hltThreshold_; }
637  float getL1Threshold() const { return l1Threshold_; }
638 
639  std::vector< std::pair<std::string,unsigned int> > filtersAndIndices;
640 
641 
642  private:
643 
644  int pathIndex_;
645  std::string denomPathName_;
646  std::string pathName_;
647  std::string l1pathName_;
648  int l1ModuleIndex_;
649  std::string filterName_;
650  std::string processName_;
651  int objectType_;
652 
653  // we don't own this data
654  MonitorElement *NOn_, *onEtOn_, *onOneOverEtOn_, *onEtavsonPhiOn_;
655  MonitorElement *NOff_, *offEtOff_, *offEtavsoffPhiOff_;
656  MonitorElement *NL1_, *l1EtL1_, *l1Etavsl1PhiL1_;
657  MonitorElement *NL1On_, *l1EtL1On_, *l1Etavsl1PhiL1On_;
658  MonitorElement *NL1Off_, *offEtL1Off_, *offEtavsoffPhiL1Off_;
659  MonitorElement *NOnOff_, *offEtOnOff_, *offEtavsoffPhiOnOff_;
660  MonitorElement *NL1OnUM_, *l1EtL1OnUM_, *l1Etavsl1PhiL1OnUM_;
661  MonitorElement *NL1OffUM_, *offEtL1OffUM_, *offEtavsoffPhiL1OffUM_;
662  MonitorElement *NOnOffUM_, *offEtOnOffUM_, *offEtavsoffPhiOnOffUM_;
663  MonitorElement *offDRL1Off_, *offDROnOff_, *l1DRL1On_;
664  MonitorElement *filters_;
665 
666  float ptmin_, ptmax_;
667  float hltThreshold_, l1Threshold_;
668 
669  const int index() {
670  return pathIndex_;
671  }
672  const int type() {
673  return objectType_;
674  }
675 
676 
677  };
678 
679 
680  public:
681 
682  // simple collection - just
683  class PathInfoCollection: public std::vector<PathInfo> {
684  public:
685 
686  PathInfoCollection(): std::vector<PathInfo>()
687  {};
688  std::vector<PathInfo>::iterator find(std::string pathName) {
689  return std::find(begin(), end(), pathName);
690  }
691  };
692 
693  PathInfoCollection hltPaths_;
694 
695  PathInfoCollection hltPathsDiagonal_;
696 
697 };
698 
699 
700 class BaseMonitor
701 {
702  public:
703  virtual void clearSets( void ) = 0;
704  virtual void setPath(FourVectorHLTOffline::PathInfoCollection::iterator v) = 0;
705  virtual void monitorOffline( void ) = 0;
706  virtual void monitorL1( const int l1Index, FourVectorHLTOffline* fv) = 0;
707  virtual void monitorOnline(const int hltIndex, const int l1Index, FourVectorHLTOffline* fv) = 0;
708 
709  virtual void matchL1Offline(const trigger::TriggerObject& l1FV, FourVectorHLTOffline* fv, const int& NL1, unsigned int& NL1OffUM) = 0;
710  virtual void matchOnlineL1(const trigger::TriggerObject& onlineFV, const int& l1Index, FourVectorHLTOffline* fv, const int& NOn) = 0;
711  virtual void matchOnlineOffline(const trigger::TriggerObject& onlineFV, FourVectorHLTOffline* fv, const int& NOn) = 0;
712 
713  virtual void fillL1Match(FourVectorHLTOffline* fv) = 0;
714  virtual void fillOnlineMatch(const int l1Index, FourVectorHLTOffline* fv) = 0;
715 
716  virtual bool isTriggerType(int t) = 0;
717  virtual ~BaseMonitor(){}
718 
719 };
720 
721 template <class T>
722 class objMonData:public BaseMonitor {
723 public:
724  objMonData() { EtaMax_= 2.5; EtMin_=3.0; GenJetsFlag_ = false; BJetsFlag_ = false; fL2MuFlag = false; }
725  void setLimits(float etaMax, float etMin, float drMatch, float l1drMatch, float dRRange, float thresholdFactor)
726  {
727  EtaMax_= etaMax;
728  EtMin_= etMin;
729  DRMatch_= drMatch;
730  L1DRMatch_= l1drMatch;
731  DRRange_ = dRRange;
732  thresholdFactor_ = thresholdFactor;
733  }
734  void setTriggerType(std::vector<int> trigType) { triggerType_ = trigType; }
735  void pushTriggerType(int trigType) { triggerType_.push_back(trigType); }
736  void setL1TriggerType(std::vector<int> trigType) { l1triggerType_ = trigType; }
737  void pushL1TriggerType(int trigType) { l1triggerType_.push_back(trigType); }
738  void setPath(FourVectorHLTOffline::PathInfoCollection::iterator v) { v_ = v; }
739  void setReco(edm::Handle<T> offColl) { offColl_ = offColl; }
740  void setRecoB(edm::Handle<reco::JetTagCollection> offCollB) { offCollB_ = offCollB; }
741  void setRecoMu(edm::Handle<reco::MuonCollection> offCollMu) { offCollMu_ = offCollMu; }
742  void setRecoEle(edm::Handle<reco::GsfElectronCollection> offCollEle) { offCollEle_ = offCollEle; }
743 
744 
745 
746  // Monitor methods
747 
748  void monitorOffline();
749  void monitorL1(const int l1Index, FourVectorHLTOffline* fv);
750  void monitorOnline(const int hltIndex, const int l1Index, FourVectorHLTOffline* fv);
751  void matchL1Offline(const trigger::TriggerObject& l1FV, FourVectorHLTOffline* fv, const int& NL1, unsigned int& NL1OffUM);
752  void matchOnlineL1(const trigger::TriggerObject& onlineFV, const int& l1Index, FourVectorHLTOffline* fv, const int& NOn);
753  void matchOnlineOffline(const trigger::TriggerObject& onlineFV, FourVectorHLTOffline* fv, const int& NOn);
754  void fillOnlineMatch(const int l1Index, FourVectorHLTOffline* fv);
755  void fillOnOffMatch(FourVectorHLTOffline* fv);
756  void fillOnL1Match(const int l1Index, FourVectorHLTOffline* fv);
757  void fillL1Match(FourVectorHLTOffline* fv);
758  void fillL1OffMatch(FourVectorHLTOffline* fv);
759 
760  void clearSets();
761 
762  bool isTriggerType(int t);
763  bool isL1TriggerType(int t);
764 
765  mmset L1OffDRMatchSet;
766  mmset L1MCDRMatchSet;
767  mmset OnOffDRMatchSet;
768  mmset OnMCDRMatchSet;
769  mmset OnL1DRMatchSet;
770  mmset OffMCDRMatchSet;
771 
772 
773  void setBJetsFlag(bool flag)
774  {
775  BJetsFlag_ = flag;
776  }
777  void setL2MuFlag(bool flag)
778  {
779  fL2MuFlag = flag;
780  }
781 
782 
783 private:
784 
785  int pdgId_;
786  int pdgStatus_;
787 
788  float EtaMax_;
789  float EtMin_;
790 
791  float DRMatch_;
792  float L1DRMatch_;
793  float DRRange_;
794  float thresholdFactor_;
795 
796  bool GenJetsFlag_;
797  bool BJetsFlag_;
798  bool fL2MuFlag;
799 
800  std::vector<int> triggerType_;
801  std::vector<int> l1triggerType_;
802 
803  edm::Handle<T> offColl_;
804  edm::Handle<reco::JetTagCollection> offCollB_;
805  edm::Handle<reco::MuonCollection> offCollMu_;
806  edm::Handle<reco::GsfElectronCollection> offCollEle_;
807 
808  FourVectorHLTOffline::PathInfoCollection::iterator v_;
809 
810 };
811 
812 
813 template <class T>
814 bool objMonData<T>::isTriggerType(int t)
815 {
816  bool rc = false;
817 
818  for(std::vector<int>::const_iterator it = triggerType_.begin(); it != triggerType_.end(); ++it)
819  {
820 
821  if(t == *it) { rc = true; break; }
822 
823  } // end for
824 
825  if (t==0) rc = true;
826 
827  return rc;
828 
829 }
830 
831 
832 template <class T>
833 bool objMonData<T>::isL1TriggerType(int t)
834 {
835  bool rc = false;
836 
837  for(std::vector<int>::const_iterator it = l1triggerType_.begin(); it != l1triggerType_.end(); ++it)
838  {
839 
840  if(fabs(t) == fabs(*it)) { rc = true; break; }
841 
842  } // end for
843 
844  return rc;
845 
846 }
847 
848 
849 
850 
851 template <class T>
852 void objMonData<T>::monitorOffline()
853 {
854 
855  if(! isTriggerType(v_->getObjectType()) ) return;
856 
857  unsigned int NOff = 0;
858 
859  if( offCollB_.isValid()) {
860  typedef typename reco::JetTagCollection::const_iterator const_iterator;
861  for( const_iterator iter = offCollB_->begin(), iend = offCollB_->end(); iter != iend; ++iter )
862  {
863 
864  float recoEta = (*iter).first->eta();
865  float recoPhi = (*iter).first->phi();
866  float recoPt = (*iter).first->pt();
867 
868 
869  if (fabs(recoEta) <= EtaMax_ && recoPt >= EtMin_ )
870  {
871 
872  NOff++;
873  v_->getOffEtOffHisto()->Fill(recoPt);
874  if(recoPt >= thresholdFactor_*v_->getHltThreshold())
875  v_->getOffEtaVsOffPhiOffHisto()->Fill(recoEta, recoPhi);
876 
877  }
878  /*
879  else {
880 
881  continue;
882 
883  }
884  */
885 
886  }
887 
888  }
889  else if(offCollEle_.isValid()) {
890 
891  typedef typename reco::GsfElectronCollection::const_iterator const_iterator;
892  for( const_iterator iter = offCollEle_->begin(), iend = offCollEle_->end(); iter != iend; ++iter )
893  {
894 
895  if (fabs(iter->eta()) <= EtaMax_ && iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()) >= EtMin_ )
896  {
897 
898  NOff++;
899  v_->getOffEtOffHisto()->Fill(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()));
900 
901  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
902  v_->getOffEtaVsOffPhiOffHisto()->Fill(iter->eta(), iter->phi());
903 
904  }
905  /*
906  else {
907 
908  continue;
909 
910  }
911  */
912 
913  }
914 
915  } // end else if
916  else if(offColl_.isValid()) {
917 
918  typedef typename T::const_iterator const_iterator;
919  for( const_iterator iter = offColl_->begin(), iend = offColl_->end(); iter != iend; ++iter )
920  {
921 
922  if (fabs(iter->eta()) <= EtaMax_ && iter->pt() >= EtMin_ )
923  {
924 
925  NOff++;
926  v_->getOffEtOffHisto()->Fill(iter->pt());
927 
928  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
929  v_->getOffEtaVsOffPhiOffHisto()->Fill(iter->eta(), iter->phi());
930 
931  }
932  /*
933  else {
934 
935  continue;
936 
937  }
938  */
939 
940  }
941 
942  } // end else if
943 
944  if(NOff>0)v_->getNOffHisto()->Fill(NOff);
945 
946 }
947 
948 
949 template <class T>
950 void objMonData<T>::monitorL1(const int l1Index, FourVectorHLTOffline* fv)
951 {
952 
953  if ( l1Index >= fv->fTriggerObj->sizeFilters() ) return;
954 
955  unsigned int NL1=0;
956  unsigned int NL1OffUM=0;
957 
958  const trigger::TriggerObjectCollection & toc(fv->fTriggerObj->getObjects());
959  const trigger::Vids & idtype = fv->fTriggerObj->filterIds(l1Index);
960  const trigger::Keys & l1k = fv->fTriggerObj->filterKeys(l1Index);
961  bool l1accept = l1k.size() > 0;
962 
963  if(!l1accept) return;
964 
965  trigger::Vids::const_iterator idtypeiter = idtype.begin();
966 
967  for (trigger::Keys::const_iterator l1ki = l1k.begin(); l1ki !=l1k.end(); ++l1ki ) {
968 
969  trigger::TriggerObject l1FV = toc[*l1ki];
970 
971  if(isL1TriggerType(*idtypeiter))
972  {
973 
974  NL1++;
975 
976  v_->getL1EtL1Histo()->Fill(l1FV.pt());
977  v_->getL1EtaVsL1PhiL1Histo()->Fill(l1FV.eta(), l1FV.phi());
978 
979  matchL1Offline(l1FV, fv, NL1, NL1OffUM);
980 
981  } // end if isL1TriggerType
982 
983  ++idtypeiter;
984 
985  } // end for l1ki
986 
987  if(NL1 > 0) v_->getNL1Histo()->Fill(NL1);
988  if(NL1OffUM > 0) v_->getNL1OffUMHisto()->Fill(NL1OffUM);
989 
990 }
991 
992 
993 template <class T>
994 void objMonData<T>::matchL1Offline(const trigger::TriggerObject& l1FV, FourVectorHLTOffline* fv, const int& NL1, unsigned int& NL1OffUM)
995 {
996 
997  fimmap L1OffDRMatchMap;
998 
999  if (offCollB_.isValid()) {
1000 
1001  int j=0;
1002  typedef typename reco::JetTagCollection::const_iterator const_iterator;
1003  for( const_iterator iter = offCollB_->begin(), iend = offCollB_->end(); iter != iend; ++iter )
1004  {
1005 
1006  float recoEta = (*iter).first->eta();
1007  float recoPhi = (*iter).first->phi();
1008  float recoPt = (*iter).first->pt();
1009 
1010  if (fabs(recoEta) <= EtaMax_ && recoPt >= EtMin_ )
1011  {
1012 
1013  // fill UM histos (no matching required)
1014  if(NL1 == 1) {
1015 
1016  NL1OffUM++;
1017  v_->getOffEtL1OffUMHisto()->Fill(recoPt);
1018 
1019  if(recoPt >= thresholdFactor_*v_->getHltThreshold())
1020  v_->getOffEtaVsOffPhiL1OffUMHisto()->Fill(recoEta,recoPhi);
1021 
1022  }
1023 
1024  // make maps of matched objects
1025  float dR = reco::deltaR(recoEta,recoPhi,l1FV.eta(),l1FV.phi());
1026  if ( dR < DRRange_)
1027  {
1028 
1029  L1OffDRMatchMap.insert(std::pair<float,int>(dR,j));
1030 
1031  }
1032 
1033  }
1034 
1035  j++;
1036 
1037  }
1038 
1039  }
1040  else if (offCollMu_.isValid()) {
1041 
1042  int j=0;
1043  typedef typename reco::MuonCollection::const_iterator const_iterator;
1044  for( const_iterator iter = offCollMu_->begin(), iend = offCollMu_->end(); iter != iend; ++iter )
1045  {
1046 
1047  // get Eta, Phi of the MuonDetectorTrack,
1048  // looking at the detector most inner Position
1049  // This should be close to what L1 sees
1050  float recoEta = iter->outerTrack()->innerPosition().eta();
1051  float recoPhi = iter->outerTrack()->innerPosition().phi();
1052  float recoPt = iter->pt();
1053 
1054  if (fabs(recoEta) <= EtaMax_ && recoPt >= EtMin_ )
1055  {
1056 
1057  // fill UM histos (no matching required)
1058  if(NL1 == 1) {
1059 
1060  NL1OffUM++;
1061  v_->getOffEtL1OffUMHisto()->Fill(recoPt);
1062 
1063  if(recoPt >= thresholdFactor_*v_->getHltThreshold())
1064  v_->getOffEtaVsOffPhiL1OffUMHisto()->Fill(recoEta,recoPhi);
1065 
1066  }
1067 
1068  // make maps of matched objects
1069  float dR = reco::deltaR(recoEta,recoPhi,l1FV.eta(),l1FV.phi());
1070  if ( dR < DRRange_)
1071  {
1072 
1073  L1OffDRMatchMap.insert(std::pair<float,int>(dR,j));
1074 
1075  }
1076 
1077  }
1078 
1079  j++;
1080 
1081  }
1082 
1083  }
1084  else if (offCollEle_.isValid()) {
1085 
1086  int j=0;
1087  typedef typename reco::GsfElectronCollection::const_iterator const_iterator;
1088  for( const_iterator iter = offCollEle_->begin(), iend = offCollEle_->end(); iter != iend; ++iter )
1089  {
1090 
1091  if (fabs(iter->eta()) <= EtaMax_ && iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()) >= EtMin_ )
1092  {
1093 
1094  // fill UM histos (no matching required)
1095  if(NL1 == 1) {
1096 
1097  NL1OffUM++;
1098  v_->getOffEtL1OffUMHisto()->Fill(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()));
1099 
1100  if(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()) >= thresholdFactor_*v_->getHltThreshold())
1101  v_->getOffEtaVsOffPhiL1OffUMHisto()->Fill(iter->eta(),iter->phi());
1102 
1103  }
1104 
1105  // make maps of matched objects
1106  float dR = reco::deltaR(iter->eta(),iter->phi(),l1FV.eta(),l1FV.phi());
1107  if ( dR < DRRange_)
1108  {
1109 
1110  L1OffDRMatchMap.insert(std::pair<float,int>(dR,j));
1111 
1112  }
1113 
1114  }
1115 
1116  j++;
1117 
1118  }
1119 
1120  }
1121  else if (offColl_.isValid()) {
1122 
1123  int j=0;
1124  typedef typename T::const_iterator const_iterator;
1125  for( const_iterator iter = offColl_->begin(), iend = offColl_->end(); iter != iend; ++iter )
1126  {
1127 
1128  if (fabs(iter->eta()) <= EtaMax_ && iter->pt() >= EtMin_ )
1129  {
1130 
1131  // fill UM histos (no matching required)
1132  if(NL1 == 1) {
1133 
1134  NL1OffUM++;
1135  v_->getOffEtL1OffUMHisto()->Fill(iter->pt());
1136 
1137  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1138  v_->getOffEtaVsOffPhiL1OffUMHisto()->Fill(iter->eta(),iter->phi());
1139 
1140  }
1141 
1142  // make maps of matched objects
1143  float dR = reco::deltaR(iter->eta(),iter->phi(),l1FV.eta(),l1FV.phi());
1144  if ( dR < DRRange_)
1145  {
1146 
1147  L1OffDRMatchMap.insert(std::pair<float,int>(dR,j));
1148 
1149  }
1150 
1151  }
1152 
1153  j++;
1154 
1155  }
1156 
1157  }
1158  if(! L1OffDRMatchMap.empty()) L1OffDRMatchSet.insert(L1OffDRMatchMap);
1159 
1160 }
1161 
1162 
1163 template <class T>
1164 void objMonData<T>::monitorOnline(const int hltIndex, const int l1Index, FourVectorHLTOffline* fv)
1165 {
1166 
1167  if(! isTriggerType(v_->getObjectType()) ) return;
1168 
1169  // Get keys of objects passed by the last filter
1170  const trigger::Keys & k = fv->fTriggerObj->filterKeys(hltIndex);
1171 
1172  const trigger::TriggerObjectCollection & toc(fv->fTriggerObj->getObjects());
1173 
1174  unsigned int NOn=0;
1175 
1176  // Loop over HLT objects
1177  for (trigger::Keys::const_iterator ki = k.begin(); ki !=k.end(); ++ki ) {
1178 
1179  trigger::TriggerObject onlineFV = toc[*ki];
1180 
1181  NOn++;
1182 
1183  v_->getOnEtOnHisto()->Fill(onlineFV.pt());
1184  v_->getOnOneOverEtOnHisto()->Fill(1./onlineFV.pt());
1185  v_->getOnEtaVsOnPhiOnHisto()->Fill(onlineFV.eta(), onlineFV.phi());
1186 
1187  matchOnlineL1(onlineFV,l1Index, fv, NOn);
1188  matchOnlineOffline(onlineFV,fv, NOn);
1189 
1190  } // end loop over HLT objects
1191 
1192  if(NOn>0) v_->getNOnHisto()->Fill(NOn);
1193 
1194 }
1195 
1196 template <class T>
1197 void objMonData<T>::matchOnlineL1(const trigger::TriggerObject& onlineFV, const int& l1Index, FourVectorHLTOffline* fv, const int& NOn)
1198 {
1199 
1200  if ( l1Index >= fv->fTriggerObj->sizeFilters() ) return;
1201 
1202  unsigned int NOnL1UM=0;
1203 
1204  const trigger::TriggerObjectCollection & toc(fv->fTriggerObj->getObjects());
1205  const trigger::Vids & idtype = fv->fTriggerObj->filterIds(l1Index);
1206  const trigger::Keys & l1k = fv->fTriggerObj->filterKeys(l1Index);
1207 
1208  fimmap OnL1DRMatchMap;
1209  int j=0;
1210  trigger::Vids::const_iterator idtypeiter = idtype.begin();
1211  for (trigger::Keys::const_iterator l1ki = l1k.begin(); l1ki !=l1k.end(); ++l1ki )
1212  {
1213 
1214  if(isL1TriggerType(*idtypeiter))
1215  {
1216 
1217  trigger::TriggerObject l1FV = toc[*l1ki];
1218 
1219  // fill UM histos (no matching required)
1220  if(NOn == 1) {
1221 
1222  NOnL1UM++;
1223  v_->getL1EtL1OnUMHisto()->Fill(l1FV.pt());
1224  v_->getL1EtaVsL1PhiL1OnUMHisto()->Fill(l1FV.eta(),l1FV.phi());
1225 
1226  }
1227 
1228 
1229  float dR = reco::deltaR(l1FV.eta(),l1FV.phi(),onlineFV.eta(),onlineFV.phi());
1230 
1231  if ( dR < DRRange_)
1232  {
1233 
1234  OnL1DRMatchMap.insert(std::pair<float,int>(dR,j));
1235 
1236  }
1237 
1238  } // end if isL1TriggerType
1239 
1240  ++idtypeiter;
1241  j++;
1242 
1243  } // end for
1244 
1245  if(! OnL1DRMatchMap.empty()) OnL1DRMatchSet.insert(OnL1DRMatchMap);
1246 
1247 }
1248 
1249 
1250 template <class T>
1251 void objMonData<T>::matchOnlineOffline(const trigger::TriggerObject& onlineFV, FourVectorHLTOffline* fv, const int& NOn)
1252 {
1253 
1254  unsigned int NOnOffUM=0;
1255 
1256  fimmap OnOffDRMatchMap;
1257 
1258  if (offCollB_.isValid()) {
1259 
1260  int j=0;
1261  typedef typename reco::JetTagCollection::const_iterator const_iterator;
1262  for( const_iterator iter = offCollB_->begin(), iend = offCollB_->end(); iter != iend; ++iter )
1263  {
1264 
1265  float recoEta = (*iter).first->eta();
1266  float recoPhi = (*iter).first->phi();
1267  float recoPt = (*iter).first->pt();
1268 
1269  if (fabs(recoEta) <= EtaMax_ && recoPt >= EtMin_ )
1270  {
1271 
1272  // fill UM histos (no matching required)
1273  if(NOn == 1) {
1274 
1275  NOnOffUM++;
1276  v_->getOffEtOnOffUMHisto()->Fill(recoPt);
1277 
1278  if(recoPt >= thresholdFactor_*v_->getHltThreshold())
1279  v_->getOffEtaVsOffPhiOnOffUMHisto()->Fill(recoEta,recoPhi);
1280 
1281  }
1282 
1283  // make maps of matched objects
1284  float dR = reco::deltaR(recoEta,recoPhi,onlineFV.eta(),onlineFV.phi());
1285  if ( dR < DRRange_)
1286  {
1287 
1288  OnOffDRMatchMap.insert(std::pair<float,int>(dR,j));
1289 
1290  }
1291 
1292  }
1293 
1294  j++;
1295 
1296  }
1297 
1298  }
1299  else if (offCollMu_.isValid() && fL2MuFlag) {
1300 
1301  int j=0;
1302  typedef typename reco::MuonCollection::const_iterator const_iterator;
1303  for( const_iterator iter = offCollMu_->begin(), iend = offCollMu_->end(); iter != iend; ++iter )
1304  {
1305 
1306  // get Eta, Phi of the MuonDetectorTrack,
1307  // looking at the detector most inner Position
1308  // This should be close to what L1 sees
1309  float recoEta = iter->outerTrack()->innerPosition().eta();
1310  float recoPhi = iter->outerTrack()->innerPosition().phi();
1311  float recoPt = iter->pt();
1312 
1313  if (fabs(recoEta) <= EtaMax_ && recoPt >= EtMin_ )
1314  {
1315  // fill UM histos (no matching required)
1316  if(NOn == 1) {
1317 
1318  NOnOffUM++;
1319  v_->getOffEtOnOffUMHisto()->Fill(iter->pt());
1320 
1321  if(recoPt >= thresholdFactor_*v_->getHltThreshold())
1322  v_->getOffEtaVsOffPhiOnOffUMHisto()->Fill(iter->eta(),iter->phi());
1323 
1324  }
1325 
1326  // make maps of matched objects
1327  float dR = reco::deltaR(recoEta,recoPhi,onlineFV.eta(),onlineFV.phi());
1328  if ( dR < DRRange_)
1329  {
1330 
1331  OnOffDRMatchMap.insert(std::pair<float,int>(dR,j));
1332 
1333  }
1334 
1335  }
1336 
1337  j++;
1338 
1339 
1340  }
1341 
1342  }
1343  else if (offCollEle_.isValid()) {
1344 
1345  int j=0;
1346 
1347  typedef typename reco::GsfElectronCollection::const_iterator const_iterator;
1348  for( const_iterator iter = offCollEle_->begin(), iend = offCollEle_->end(); iter != iend; ++iter )
1349  {
1350 
1351  if (fabs(iter->eta()) <= EtaMax_ && iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()) >= EtMin_ )
1352  {
1353 
1354  // fill UM histos (no matching required)
1355  if(NOn == 1) {
1356 
1357  NOnOffUM++;
1358  v_->getOffEtOnOffUMHisto()->Fill(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()));
1359 
1360  if(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()) >= thresholdFactor_*v_->getHltThreshold())
1361  v_->getOffEtaVsOffPhiOnOffUMHisto()->Fill(iter->eta(),iter->phi());
1362 
1363  }
1364 
1365  // make maps of matched objects
1366  float dR = reco::deltaR(iter->eta(),iter->phi(),onlineFV.eta(),onlineFV.phi());
1367  if ( dR < DRRange_)
1368  {
1369 
1370  OnOffDRMatchMap.insert(std::pair<float,int>(dR,j));
1371 
1372  }
1373 
1374  }
1375 
1376  j++;
1377 
1378 
1379  }
1380 
1381  }
1382  else if (offColl_.isValid()) {
1383 
1384  int j=0;
1385 
1386  typedef typename T::const_iterator const_iterator;
1387  for( const_iterator iter = offColl_->begin(), iend = offColl_->end(); iter != iend; ++iter )
1388  {
1389 
1390  if (fabs(iter->eta()) <= EtaMax_ && iter->pt() >= EtMin_ )
1391  {
1392 
1393  // fill UM histos (no matching required)
1394  if(NOn == 1) {
1395 
1396  NOnOffUM++;
1397  v_->getOffEtOnOffUMHisto()->Fill(iter->pt());
1398 
1399  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1400  v_->getOffEtaVsOffPhiOnOffUMHisto()->Fill(iter->eta(),iter->phi());
1401 
1402  }
1403 
1404  // make maps of matched objects
1405  float dR = reco::deltaR(iter->eta(),iter->phi(),onlineFV.eta(),onlineFV.phi());
1406  if ( dR < DRRange_)
1407  {
1408 
1409  OnOffDRMatchMap.insert(std::pair<float,int>(dR,j));
1410 
1411  }
1412 
1413  }
1414 
1415  j++;
1416 
1417 
1418  }
1419 
1420  }
1421 
1422  if(! OnOffDRMatchMap.empty()) OnOffDRMatchSet.insert(OnOffDRMatchMap);
1423 
1424 }
1425 
1426 template <class T>
1427 void objMonData<T>::fillL1OffMatch(FourVectorHLTOffline* fv)
1428 {
1429 
1430  float NL1Off=0;
1431 
1432  if(L1OffDRMatchSet.size() > 1) {
1433 
1434  LogDebug("FourVectorHLTOffline") << " Cleaning L1Off mmset" << std::endl;
1435  fv->cleanDRMatchSet(L1OffDRMatchSet);
1436 
1437  }
1438  // clean the set L1-Off
1439  // now fill histos
1440  for ( mmset::iterator setIter = L1OffDRMatchSet.begin( ); setIter != L1OffDRMatchSet.end( ); setIter++ )
1441  {
1442 
1443  fimmap tempMap = *setIter;
1444 
1445  fimmap::iterator it = tempMap.begin();
1446  int i = (*it).second ;
1447  float dR = (*it).first;
1448  v_->getOffDRL1OffHisto()->Fill(dR);
1449 
1450  if (dR > L1DRMatch_) continue;
1451  if( offCollB_.isValid()) {
1452 
1453  typedef typename reco::JetTagCollection::const_iterator const_iterator;
1454  const_iterator iter = offCollB_->begin();
1455  for (int count = 0; count < i; count++) iter++;
1456 
1457 
1458  NL1Off++;
1459  v_->getOffEtL1OffHisto()->Fill((*iter).first->pt());
1460  if((*iter).first->pt() >= thresholdFactor_*v_->getHltThreshold())
1461  v_->getOffEtaVsOffPhiL1OffHisto()->Fill((*iter).first->eta(),(*iter).first->phi());
1462 
1463 
1464  }
1465  else if( offCollMu_.isValid()) {
1466 
1467  typedef typename reco::MuonCollection::const_iterator const_iterator;
1468  const_iterator iter = offCollMu_->begin();
1469  for (int count = 0; count < i; count++) iter++;
1470 
1471 
1472  NL1Off++;
1473  v_->getOffEtL1OffHisto()->Fill(iter->pt());
1474  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1475  v_->getOffEtaVsOffPhiL1OffHisto()->Fill(iter->outerTrack()->innerPosition().eta(),iter->outerTrack()->innerPosition().phi());
1476 
1477  }
1478  else if( offCollEle_.isValid()) {
1479 
1480  typedef typename reco::GsfElectronCollection::const_iterator const_iterator;
1481  const_iterator iter = offCollEle_->begin();
1482  for (int count = 0; count < i; count++) iter++;
1483 
1484 
1485  NL1Off++;
1486  v_->getOffEtL1OffHisto()->Fill(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()));
1487  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1488  v_->getOffEtaVsOffPhiL1OffHisto()->Fill(iter->eta(),iter->phi());
1489 
1490  }
1491  else if( offColl_.isValid()) {
1492 
1493  typedef typename T::const_iterator const_iterator;
1494  const_iterator iter = offColl_->begin();
1495  for (int count = 0; count < i; count++) iter++;
1496 
1497 
1498  NL1Off++;
1499  v_->getOffEtL1OffHisto()->Fill(iter->pt());
1500  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1501  v_->getOffEtaVsOffPhiL1OffHisto()->Fill(iter->eta(),iter->phi());
1502 
1503  }
1504 
1505  }
1506 
1507  if(NL1Off > 0) v_->getNL1OffHisto()->Fill(NL1Off);
1508 
1509 }
1510 
1511 
1512 template <class T>
1513 void objMonData<T>::fillOnOffMatch(FourVectorHLTOffline* fv)
1514 {
1515 
1516  unsigned int NOnOff=0;
1517 
1518  // clean the set L1-Off
1519  if(OnOffDRMatchSet.size() > 1){
1520 
1521  LogDebug("FourVectorHLTOffline") << " Cleaning OnOff mmset" << std::endl;
1522  fv->cleanDRMatchSet(OnOffDRMatchSet);
1523 
1524  }
1525  // now fill histos
1526  for ( mmset::iterator setIter = OnOffDRMatchSet.begin( ); setIter != OnOffDRMatchSet.end( ); setIter++ )
1527  {
1528 
1529 
1530  fimmap tempMap = *setIter;
1531 
1532  fimmap::iterator it = tempMap.begin();
1533  int i = (*it).second ;
1534  float dR = (*it).first;
1535  v_->getOffDROnOffHisto()->Fill(dR);
1536 
1537 
1538  if (dR > DRMatch_) continue;
1539 
1540  if( offCollB_.isValid()) {
1541 
1542 
1543  typedef typename reco::JetTagCollection::const_iterator const_iterator;
1544  const_iterator iter = offCollB_->begin();
1545  for (int count = 0; count < i; count++) iter++;
1546 
1547 
1548  NOnOff++;
1549  v_->getOffEtOnOffHisto()->Fill((*iter).first->pt());
1550  if((*iter).first->pt() >= thresholdFactor_*v_->getHltThreshold())
1551  v_->getOffEtaVsOffPhiOnOffHisto()->Fill((*iter).first->eta(),(*iter).first->phi());
1552 
1553  }
1554  else if( offCollMu_.isValid() && fL2MuFlag) {
1555 
1556  typedef typename reco::MuonCollection::const_iterator const_iterator;
1557  const_iterator iter = offCollMu_->begin();
1558  for (int count = 0; count < i; count++) iter++;
1559 
1560 
1561  NOnOff++;
1562  v_->getOffEtOnOffHisto()->Fill(iter->pt());
1563  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1564  v_->getOffEtaVsOffPhiOnOffHisto()->Fill(iter->outerTrack()->innerPosition().eta(),iter->outerTrack()->innerPosition().phi());
1565 
1566  }
1567  else if( offCollEle_.isValid()) {
1568 
1569  typedef typename reco::GsfElectronCollection::const_iterator const_iterator;
1570  const_iterator iter = offCollEle_->begin();
1571  for (int count = 0; count < i; count++) iter++;
1572 
1573  NOnOff++;
1574  v_->getOffEtOnOffHisto()->Fill(iter->superCluster()->energy()*sin(iter->superCluster()->position().Theta()));
1575  if(iter->superCluster()->energy()*fabs(sin(iter->superCluster()->position().Theta())) >= thresholdFactor_*v_->getHltThreshold())
1576  v_->getOffEtaVsOffPhiOnOffHisto()->Fill(iter->eta(),iter->phi());
1577 
1578  }
1579  else if( offColl_.isValid()) {
1580 
1581  typedef typename T::const_iterator const_iterator;
1582  const_iterator iter = offColl_->begin();
1583  for (int count = 0; count < i; count++) iter++;
1584 
1585  NOnOff++;
1586  v_->getOffEtOnOffHisto()->Fill(iter->pt());
1587  if(iter->pt() >= thresholdFactor_*v_->getHltThreshold())
1588  v_->getOffEtaVsOffPhiOnOffHisto()->Fill(iter->eta(),iter->phi());
1589 
1590  }
1591 
1592  }
1593 
1594  v_->getNOnOffHisto()->Fill(NOnOff);
1595 
1596 }
1597 
1598 
1599 template <class T>
1600 void objMonData<T>::fillOnL1Match(const int l1Index, FourVectorHLTOffline* fv)
1601 {
1602 
1603  const trigger::TriggerObjectCollection & toc(fv->fTriggerObj->getObjects());
1604  const trigger::Keys & l1k = fv->fTriggerObj->filterKeys(l1Index);
1605 
1606  unsigned int NOnL1=0;
1607 
1608  // clean the set On-L1
1609  if(OnL1DRMatchSet.size() > 1) {
1610 
1611  LogDebug("FourVectorHLTOffline") << " Cleaning L1On mmset" << std::endl;
1612  fv->cleanDRMatchSet(OnL1DRMatchSet);
1613 
1614  }
1615  // now fill histos
1616  for ( mmset::iterator setIter = OnL1DRMatchSet.begin( ); setIter != OnL1DRMatchSet.end( ); setIter++ )
1617  {
1618 
1619  fimmap tempMap = *setIter;
1620 
1621  fimmap::iterator it = tempMap.begin();
1622  int i = (*it).second ;
1623  float dR = (*it).first;
1624  v_->getL1DROnL1Histo()->Fill(dR);
1625 
1626  if (dR > L1DRMatch_) continue;
1627 
1628  trigger::Keys::const_iterator l1ki = l1k.begin();
1629  for (int count = 0; count < i; count++) l1ki++;
1630 
1631  NOnL1++;
1632  v_->getL1EtL1OnHisto()->Fill(toc[*l1ki].pt());
1633  v_->getL1EtaVsL1PhiL1OnHisto()->Fill(toc[*l1ki].eta(),toc[*l1ki].phi());
1634 
1635  }
1636 
1637  v_->getNL1OnHisto()->Fill(NOnL1);
1638 
1639 }
1640 
1641 template <class T>
1642 void objMonData<T>::fillOnlineMatch(const int l1Index, FourVectorHLTOffline* fv)
1643 {
1644 
1645  if(! isTriggerType(v_->getObjectType()) ) return;
1646  fillOnOffMatch(fv);
1647 
1648  if ( l1Index >= fv->fTriggerObj->sizeFilters() ) return;
1649  fillOnL1Match(l1Index, fv);
1650 
1651 }
1652 
1653 template <class T>
1654 void objMonData<T>::fillL1Match(FourVectorHLTOffline* fv)
1655 {
1656 
1657  fillL1OffMatch(fv);
1658 
1659 }
1660 
1661 template <class T>
1662 void objMonData<T>::clearSets()
1663 {
1664 
1665  L1OffDRMatchSet.clear();
1666  L1MCDRMatchSet.clear();
1667  OnOffDRMatchSet.clear();
1668  OnMCDRMatchSet.clear();
1669  OnL1DRMatchSet.clear();
1670  OffMCDRMatchSet.clear();
1671 
1672 }
1673 
1674 
1675 
1676 #endif
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:510
FourVectorHLTOffline::PathInfo::getL1EtaVsL1PhiL1OnUMHisto
MonitorElement * getL1EtaVsL1PhiL1OnUMHisto()
Definition: FourVectorHLTOffline.h:477
FourVectorHLTOffline::setupHltMatrix
void setupHltMatrix(const std::string &label, std::vector< std::string > &paths)
Definition: FourVectorHLTOffline.cc:1477
FourVectorHLTOffline::jetID
reco::helper::JetIDHelper * jetID
Definition: FourVectorHLTOffline.h:167
ExpressReco_HICollisions_FallBack.etMin
tuple etMin
Definition: ExpressReco_HICollisions_FallBack.py:25165
type
type
Definition: HCALResponse.h:22
i
int i
Definition: DBlmapReader.cc:9
objMonData::offColl_
edm::Handle< T > offColl_
Definition: FourVectorHLTOnline.h:651
FourVectorHLTOffline::hasL1Passed
bool hasL1Passed(const std::string &pathname, const edm::TriggerNames &triggerNames)
Definition: FourVectorHLTOffline.cc:2269
BaseMonitor::fillL1Match
virtual void fillL1Match(FourVectorHLTOffline *fv)=0
flag
long int flag
Definition: mlp_lapack.h:47
FourVectorHLTOffline::selectPhotons
void selectPhotons(const edm::Handle< reco::PhotonCollection > &phoHandle)
Definition: FourVectorHLTOffline.cc:2421
reco::PFTauCollection
std::vector< PFTau > PFTauCollection
collection of PFTau objects
Definition: PFTauFwd.h:9
BaseMonitor::~BaseMonitor
virtual ~BaseMonitor()
Definition: FourVectorHLTOffline.h:717
FourVectorHLTOffline::countHLTPathHitsEndLumiBlock
void countHLTPathHitsEndLumiBlock(const int &lumi)
Definition: FourVectorHLTOffline.cc:2097
BaseMonitor::monitorOnline
virtual void monitorOnline(const int hltIndex, const int l1Index, FourVectorHLTOnline *fv)=0
FourVectorHLTOffline::countHLTGroupBXHitsEndLumiBlock
void countHLTGroupBXHitsEndLumiBlock(const int &lumi)
Definition: FourVectorHLTOffline.cc:1799
FourVectorHLTOffline::fSelElectronsHandle
edm::Handle< reco::GsfElectronCollection > fSelElectronsHandle
Definition: FourVectorHLTOffline.h:196
FourVectorHLTOffline::fSelMetHandle
edm::Handle< reco::CaloMETCollection > fSelMetHandle
Definition: FourVectorHLTOffline.h:205
corrections_cfg.tagName
tuple tagName
Definition: corrections_cfg.py:139
objMonData::v_
FourVectorHLTOffline::PathInfoCollection::iterator v_
Definition: FourVectorHLTOffline.h:808
label
const std::string & label
Definition: MVAComputer.cc:186
BaseMonitor::fillOnlineMatch
virtual void fillOnlineMatch(const int l1Index, FourVectorHLTOnline *fv)=0
edm::AssociationVector::const_iterator
transient_vector_type::const_iterator const_iterator
Definition: AssociationVector.h:95
FourVectorHLTOffline::triggerFilterIndices_
std::vector< std::vector< uint > > triggerFilterIndices_
Definition: FourVectorHLTOffline.h:311
objMonData::isL1TriggerType
bool isL1TriggerType(int t)
Definition: FourVectorHLTOnline.h:678
FourVectorHLTOffline::PathInfo::getL1EtaVsL1PhiL1OnHisto
MonitorElement * getL1EtaVsL1PhiL1OnHisto()
Definition: FourVectorHLTOffline.h:450
FourVectorHLTOffline::fSelectedMet
reco::CaloMETCollection * fSelectedMet
Definition: FourVectorHLTOffline.h:204
FourVectorHLTOffline::hltPaths_
PathInfoCollection hltPaths_
Definition: FourVectorHLTOffline.h:693
FourVectorHLTOffline::PathInfo::getOffEtaVsOffPhiL1OffUMHisto
MonitorElement * getOffEtaVsOffPhiL1OffUMHisto()
Definition: FourVectorHLTOffline.h:486
objMonData::isTriggerType
bool isTriggerType(int t)
Definition: FourVectorHLTOnline.h:659
trigger::TriggerObject::phi
float phi() const
Definition: TriggerObject.h:60
FourVectorHLTOffline::fSelectedTaus
reco::PFTauCollection * fSelectedTaus
Definition: FourVectorHLTOffline.h:209
FourVectorHLTOffline::PathInfoCollection::find
std::vector< PathInfo >::iterator find(std::string pathName)
Definition: FourVectorHLTOffline.h:688
FourVectorHLTOffline::endRun
void endRun(const edm::Run &run, const edm::EventSetup &c)
EndRun.
Definition: FourVectorHLTOffline.cc:1323
FourVectorHLTOffline::pathsSummaryFilterEfficiencyFolder_
std::string pathsSummaryFilterEfficiencyFolder_
Definition: FourVectorHLTOffline.h:183
objMonData::setRecoMu
void setRecoMu(edm::Handle< reco::MuonCollection > offCollMu)
Definition: FourVectorHLTOffline.h:741
fjr2json.lumi
tuple lumi
Definition: fjr2json.py:41
FourVectorHLTOffline::pathsSummaryHLTPathsPerLSFolder_
std::string pathsSummaryHLTPathsPerLSFolder_
Definition: FourVectorHLTOffline.h:185
FourVectorHLTOffline::getTriggerTypeParsePathName
int getTriggerTypeParsePathName(const std::string &pathname)
Definition: FourVectorHLTOffline.cc:2211
FourVectorHLTOffline::fillHltMatrix
void fillHltMatrix(const edm::TriggerNames &triggerNames)
Definition: FourVectorHLTOffline.cc:1551
FourVectorHLTOffline::fPathBxTempCountPair
std::vector< std::pair< std::string, std::vector< int > > > fPathBxTempCountPair
Definition: FourVectorHLTOffline.h:313
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
FourVectorHLTOffline::specialPaths_
std::vector< std::string > specialPaths_
Definition: FourVectorHLTOffline.h:318
objMonData::matchOnlineL1
void matchOnlineL1(const trigger::TriggerObject &onlineFV, const int &l1Index, FourVectorHLTOnline *fv, const int &NOn)
Definition: FourVectorHLTOnline.h:789
FourVectorHLTOffline::FourVectorHLTOffline
FourVectorHLTOffline(const edm::ParameterSet &)
Definition: FourVectorHLTOffline.cc:15
objMonData::l1triggerType_
std::vector< int > l1triggerType_
Definition: FourVectorHLTOnline.h:649
FourVectorHLTOffline::cleanDRMatchSet
void cleanDRMatchSet(mmset &tempSet)
Clean DR Match Set.
Definition: FourVectorHLTOffline.cc:1331
fimmap
std::multimap< float, int > fimmap
Definition: FourVectorHLTOnline.h:73
reco_skim_cfg_mod.filterName
tuple filterName
Definition: reco_skim_cfg_mod.py:159
FourVectorHLTOffline::PathInfo::PathInfo
PathInfo(std::string denomPathName, std::string pathName, std::string l1pathName, std::string filterName, std::string processName, size_t type, MonitorElement *NOn, MonitorElement *onEtOn, MonitorElement *onOneOverEtOn, MonitorElement *onEtavsonPhiOn, MonitorElement *NOff, MonitorElement *offEtOff, MonitorElement *offEtavsoffPhiOff, MonitorElement *NL1, MonitorElement *l1EtL1, MonitorElement *l1Etavsl1PhiL1, MonitorElement *NL1On, MonitorElement *l1EtL1On, MonitorElement *l1Etavsl1PhiL1On, MonitorElement *NL1Off, MonitorElement *offEtL1Off, MonitorElement *offEtavsoffPhiL1Off, MonitorElement *NOnOff, MonitorElement *offEtOnOff, MonitorElement *offEtavsoffPhiOnOff, MonitorElement *NL1OnUM, MonitorElement *l1EtL1OnUM, MonitorElement *l1Etavsl1PhiL1OnUM, MonitorElement *NL1OffUM, MonitorElement *offEtL1OffUM, MonitorElement *offEtavsoffPhiL1OffUM, MonitorElement *NOnOffUM, MonitorElement *offEtOnOffUM, MonitorElement *offEtavsoffPhiOnOffUM, MonitorElement *offDRL1Off, MonitorElement *offDROnOff, MonitorElement *l1DRL1On, MonitorElement *filters, float ptmin, float ptmax)
Definition: FourVectorHLTOffline.h:568
objMonData::triggerType_
std::vector< int > triggerType_
Definition: FourVectorHLTOnline.h:648
FourVectorHLTOffline::fSelMuonsHandle
edm::Handle< reco::MuonCollection > fSelMuonsHandle
Definition: FourVectorHLTOffline.h:193
objMonData::offCollEle_
edm::Handle< reco::GsfElectronCollection > offCollEle_
Definition: FourVectorHLTOffline.h:806
trigger::TriggerObject::eta
float eta() const
Definition: TriggerObject.h:59
FourVectorHLTOffline::fBeamSpotHandle
edm::Handle< reco::BeamSpot > fBeamSpotHandle
Definition: FourVectorHLTOffline.h:123
FourVectorHLTOffline::PathInfo::operator!=
bool operator!=(const std::string &v)
Definition: FourVectorHLTOffline.h:629
BaseMonitor::monitorOffline
virtual void monitorOffline(void)=0
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
FourVectorHLTOffline::fSelTausHandle
edm::Handle< reco::PFTauCollection > fSelTausHandle
Definition: FourVectorHLTOffline.h:210
objMonData::fillOnlineMatch
void fillOnlineMatch(const int l1Index, FourVectorHLTOnline *fv)
Definition: FourVectorHLTOnline.h:891
objMonData::offCollB_
edm::Handle< reco::JetTagCollection > offCollB_
Definition: FourVectorHLTOffline.h:804
BaseMonitor::matchOnlineOffline
virtual void matchOnlineOffline(const trigger::TriggerObject &onlineFV, FourVectorHLTOffline *fv, const int &NOn)=0
FourVectorHLTOffline::PathInfo::getPath
const std::string & getPath(void) const
Definition: FourVectorHLTOffline.h:517
reco::GsfElectronCollection
std::vector< GsfElectron > GsfElectronCollection
collection of GsfElectron objects
Definition: GsfElectronFwd.h:14
eta
T eta() const
Definition: newBasic3DVector.h:122
FourVectorHLTOffline::PathInfo::getTag
const edm::InputTag getTag(void) const
Definition: FourVectorHLTOffline.h:536
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
FourVectorHLTOffline::triggerSummaryLabel_
edm::InputTag triggerSummaryLabel_
Definition: FourVectorHLTOffline.h:325
FourVectorHLTOffline::countHLTGroupL1HitsEndLumiBlock
void countHLTGroupL1HitsEndLumiBlock(const int &lumi)
Definition: FourVectorHLTOffline.cc:1938
objMonData::setRecoB
void setRecoB(edm::Handle< reco::JetTagCollection > offCollB)
Definition: FourVectorHLTOffline.h:740
BaseMonitor::isTriggerType
virtual bool isTriggerType(int t)=0
FourVectorHLTOffline::hltPathsDiagonal_
PathInfoCollection hltPathsDiagonal_
Definition: FourVectorHLTOffline.h:695
FourVectorHLTOffline::custompathnamepairs_
std::vector< std::pair< std::string, std::string > > custompathnamepairs_
Definition: FourVectorHLTOffline.h:308
objMonData::setPath
void setPath(FourVectorHLTOffline::PathInfoCollection::iterator v)
Definition: FourVectorHLTOffline.h:738
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
mmset
std::set< fimmap, std::less< fimmap > > mmset
Definition: FourVectorHLTOnline.h:74
FourVectorHLTOffline::fSelPhotonsHandle
edm::Handle< reco::PhotonCollection > fSelPhotonsHandle
Definition: FourVectorHLTOffline.h:199
reco::deltaR
double deltaR(double eta1, double phi1, double eta2, double phi2)
Definition: deltaR.h:19
FourVectorHLTOffline::fSelectedElectrons
reco::GsfElectronCollection * fSelectedElectrons
Definition: FourVectorHLTOffline.h:195
trigger::TriggerObject
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:24
pathName
const std::string * pathName() const
Definition: HLTadd.h:31
FourVectorHLTOffline::PathInfo::getL1ModuleIndex
const int getL1ModuleIndex(void) const
Definition: FourVectorHLTOffline.h:523
FourVectorHLTOffline::PathInfo::getOffEtaVsOffPhiOffHisto
MonitorElement * getOffEtaVsOffPhiOffHisto()
Definition: FourVectorHLTOffline.h:432
objMonData::v_
FourVectorHLTOnline::PathInfoCollection::iterator v_
Definition: FourVectorHLTOnline.h:653
objMonData::matchOnlineOffline
void matchOnlineOffline(const trigger::TriggerObject &onlineFV, FourVectorHLTOffline *fv, const int &NOn)
Definition: FourVectorHLTOffline.h:1251
FourVectorHLTOffline::getL1ConditionModuleName
const std::string getL1ConditionModuleName(const std::string &pathname)
Definition: FourVectorHLTOffline.cc:2242
iEvent
int iEvent
Definition: GenABIO.cc:243
objMonData::pushTriggerType
void pushTriggerType(int trigType)
Definition: FourVectorHLTOffline.h:735
objMonData::monitorOnline
void monitorOnline(const int hltIndex, const int l1Index, FourVectorHLTOnline *fv)
Definition: FourVectorHLTOnline.h:747
FourVectorHLTOffline::fTriggerResults
edm::Handle< edm::TriggerResults > fTriggerResults
Definition: FourVectorHLTOffline.h:122
FourVectorHLTOffline::selectJets
void selectJets(const edm::Event &iEvent, const edm::Handle< reco::CaloJetCollection > &jetHandle)
Definition: FourVectorHLTOffline.cc:2453
FourVectorHLTOffline::v_ME_Total_BX
std::vector< MonitorElement * > v_ME_Total_BX
Definition: FourVectorHLTOffline.h:178
objMonData::matchL1Offline
void matchL1Offline(const trigger::TriggerObject &l1FV, FourVectorHLTOffline *fv, const int &NL1, unsigned int &NL1OffUM)
Definition: FourVectorHLTOffline.h:994
FourVectorHLTOffline::pathsSummaryHLTPathsPerBXFolder_
std::string pathsSummaryHLTPathsPerBXFolder_
Definition: FourVectorHLTOffline.h:187
BaseMonitor::matchL1Offline
virtual void matchL1Offline(const trigger::TriggerObject &l1FV, FourVectorHLTOffline *fv, const int &NL1, unsigned int &NL1OffUM)=0
FourVectorHLTOffline::fPathTempCountPair
std::vector< std::pair< std::string, float > > fPathTempCountPair
Definition: FourVectorHLTOffline.h:312
j
int j
Definition: DBlmapReader.cc:9
FourVectorHLTOffline::PathInfo::getProcess
const std::string & getProcess(void) const
Definition: FourVectorHLTOffline.h:529
FourVectorHLTOffline::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
Definition: FourVectorHLTOffline.cc:1785
FourVectorHLTOffline::hltConfig_
HLTConfigProvider hltConfig_
Definition: FourVectorHLTOffline.h:328
BaseMonitor::matchOnlineL1
virtual void matchOnlineL1(const trigger::TriggerObject &onlineFV, const int &l1Index, FourVectorHLTOnline *fv, const int &NOn)=0
FourVectorHLTOffline::PathInfo::setLabel
void setLabel(std::string labelName)
Definition: FourVectorHLTOffline.h:513
end
#define end
Definition: vmac.h:38
FourVectorHLTOffline::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: FourVectorHLTOffline.cc:236
FourVectorHLTOffline::PathInfo::setHistos
void setHistos(MonitorElement *const NOn, MonitorElement *const onEtOn, MonitorElement *const onOneOverEtOn, MonitorElement *const onEtavsonPhiOn, MonitorElement *const NOff, MonitorElement *const offEtOff, MonitorElement *const offEtavsoffPhiOff, MonitorElement *const NL1, MonitorElement *const l1EtL1, MonitorElement *const l1Etavsl1PhiL1, MonitorElement *const NL1On, MonitorElement *const l1EtL1On, MonitorElement *const l1Etavsl1PhiL1On, MonitorElement *const NL1Off, MonitorElement *const offEtL1Off, MonitorElement *const offEtavsoffPhiL1Off, MonitorElement *const NOnOff, MonitorElement *const offEtOnOff, MonitorElement *const offEtavsoffPhiOnOff, MonitorElement *const NL1OnUM, MonitorElement *const l1EtL1OnUM, MonitorElement *const l1Etavsl1PhiL1OnUM, MonitorElement *const NL1OffUM, MonitorElement *const offEtL1OffUM, MonitorElement *const offEtavsoffPhiL1OffUM, MonitorElement *const NOnOffUM, MonitorElement *const offEtOnOffUM, MonitorElement *const offEtavsoffPhiOnOffUM, MonitorElement *const offDRL1Off, MonitorElement *const offDROnOff, MonitorElement *const l1DRL1On)
Definition: FourVectorHLTOffline.h:345
FourVectorHLTOffline::fSelectedPhotons
reco::PhotonCollection * fSelectedPhotons
Definition: FourVectorHLTOffline.h:198
FourVectorHLTOffline::v_ME_Total_BX_Norm
std::vector< MonitorElement * > v_ME_Total_BX_Norm
Definition: FourVectorHLTOffline.h:179
objMonData::fillL1OffMatch
void fillL1OffMatch(FourVectorHLTOffline *fv)
Definition: FourVectorHLTOffline.h:1427
FourVectorHLTOffline::selectTaus
void selectTaus(const edm::Event &iEvent)
Definition: FourVectorHLTOffline.cc:2506
FourVectorHLTOffline::PathInfo::filtersAndIndices
std::vector< std::pair< std::string, unsigned int > > filtersAndIndices
Definition: FourVectorHLTOffline.h:639
FourVectorHLTOffline::PathInfo::operator==
bool operator==(const std::string &v)
Definition: FourVectorHLTOffline.h:624
FourVectorHLTOffline::PathInfo::getObjectType
const int getObjectType(void) const
Definition: FourVectorHLTOffline.h:532
FourVectorHLTOffline::isVBTFMuon
bool isVBTFMuon(const reco::Muon &muon)
Definition: FourVectorHLTOffline.cc:2589
FourVectorHLTOffline::ME_HLT_BX
MonitorElement * ME_HLT_BX
Definition: FourVectorHLTOffline.h:175
ExpressReco_HICollisions_FallBack.pt
tuple pt
Definition: ExpressReco_HICollisions_FallBack.py:1070
FourVectorHLTOffline::selectElectrons
void selectElectrons(const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::Handle< reco::GsfElectronCollection > &eleHandle)
Definition: FourVectorHLTOffline.cc:2334
gen::k
int k[5][pyjets_maxn]
Definition: Pythia6Hadronizer.cc:86
FourVectorHLTOffline::fSelectedMuons
reco::MuonCollection * fSelectedMuons
Definition: FourVectorHLTOffline.h:192
FourVectorHLTOffline::triggerResultsLabel_
edm::InputTag triggerResultsLabel_
Definition: FourVectorHLTOffline.h:326
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:83
BaseMonitor::clearSets
virtual void clearSets(void)=0
objMonData::monitorL1
void monitorL1(const int l1Index, FourVectorHLTOnline *fv)
Definition: FourVectorHLTOnline.h:695
objMonData::setLimits
void setLimits(float etaMax, float etMin, float drMatch, float l1drMatch, float dRRange, float thresholdFactor)
Definition: FourVectorHLTOffline.h:725
FourVectorHLTOffline::beginRun
void beginRun(const edm::Run &run, const edm::EventSetup &c)
Definition: FourVectorHLTOffline.cc:777
reco::CaloMETCollection
std::vector< reco::CaloMET > CaloMETCollection
collection of CaloMET objects
Definition: CaloMETCollection.h:21
FourVectorHLTOffline::ME_HLT_CUSTOM_BX
MonitorElement * ME_HLT_CUSTOM_BX
Definition: FourVectorHLTOffline.h:176
FourVectorHLTOffline::triggerFilters_
std::vector< std::vector< std::string > > triggerFilters_
Definition: FourVectorHLTOffline.h:310
objMonData::fillOnOffMatch
void fillOnOffMatch(FourVectorHLTOffline *fv)
Definition: FourVectorHLTOffline.h:1513
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:22
FourVectorHLTOffline::PathInfo::getOffEtaVsOffPhiOnOffUMHisto
MonitorElement * getOffEtaVsOffPhiOnOffUMHisto()
Definition: FourVectorHLTOffline.h:495
FourVectorHLTOffline::fSelJetsHandle
edm::Handle< reco::CaloJetCollection > fSelJetsHandle
Definition: FourVectorHLTOffline.h:202
FourVectorHLTOffline::PathInfo::getLabel
const std::string getLabel(void) const
Definition: FourVectorHLTOffline.h:510
objMonData::setL1TriggerType
void setL1TriggerType(std::vector< int > trigType)
Definition: FourVectorHLTOffline.h:736
reco::PhotonCollection
std::vector< Photon > PhotonCollection
collectin of Photon objects
Definition: PhotonFwd.h:9
FourVectorHLTOffline::fGroupTempCountPair
std::vector< std::pair< std::string, float > > fGroupTempCountPair
Definition: FourVectorHLTOffline.h:314
FourVectorHLTOffline::fGroupNamePathsPair
std::vector< std::pair< std::string, std::vector< std::string > > > fGroupNamePathsPair
Definition: FourVectorHLTOffline.h:316
objMonData::setRecoEle
void setRecoEle(edm::Handle< reco::GsfElectronCollection > offCollEle)
Definition: FourVectorHLTOffline.h:742
objMonData::fillL1Match
void fillL1Match(FourVectorHLTOffline *fv)
Definition: FourVectorHLTOffline.h:1654
objMonData::setL2MuFlag
void setL2MuFlag(bool flag)
Definition: FourVectorHLTOffline.h:777
BaseMonitor::setPath
virtual void setPath(FourVectorHLTOnline::PathInfoCollection::iterator v)=0
ptmin
double ptmin
Definition: HydjetWrapper.h:86
FourVectorHLTOffline::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
Definition: FourVectorHLTOffline.cc:1777
FourVectorHLTOffline::pathsSummaryHLTCorrelationsFolder_
std::string pathsSummaryHLTCorrelationsFolder_
Definition: FourVectorHLTOffline.h:182
FourVectorHLTOffline::pathsSummaryFilterCountsFolder_
std::string pathsSummaryFilterCountsFolder_
Definition: FourVectorHLTOffline.h:184
FourVectorHLTOffline::PathInfo::getl1Path
const std::string & getl1Path(void) const
Definition: FourVectorHLTOffline.h:520
FourVectorHLTOffline::pathsIndividualHLTPathsPerLSFolder_
std::string pathsIndividualHLTPathsPerLSFolder_
Definition: FourVectorHLTOffline.h:186
begin
#define begin
Definition: vmac.h:31
FourVectorHLTOffline::PathInfo::getOffEtaVsOffPhiOnOffHisto
MonitorElement * getOffEtaVsOffPhiOnOffHisto()
Definition: FourVectorHLTOffline.h:468
FourVectorHLTOffline::selectMet
void selectMet(const edm::Handle< reco::CaloMETCollection > &metHandle)
Definition: FourVectorHLTOffline.cc:2483
FourVectorHLTOffline::PathInfo::getOffEtaVsOffPhiL1OffHisto
MonitorElement * getOffEtaVsOffPhiL1OffHisto()
Definition: FourVectorHLTOffline.h:459
FourVectorHLTOffline::ME_HLTAll_LS
MonitorElement * ME_HLTAll_LS
Definition: FourVectorHLTOffline.h:174
FourVectorHLTOffline::fGroupName
std::vector< std::string > fGroupName
Definition: FourVectorHLTOffline.h:190
objMonData::fillOnL1Match
void fillOnL1Match(const int l1Index, FourVectorHLTOnline *fv)
Definition: FourVectorHLTOnline.h:849
FourVectorHLTOffline::countHLTGroupHitsEndLumiBlock
void countHLTGroupHitsEndLumiBlock(const int &lumi)
Definition: FourVectorHLTOffline.cc:2018
BaseMonitor::monitorL1
virtual void monitorL1(const int l1Index, FourVectorHLTOnline *fv)=0
FourVectorHLTOffline::fGroupL1TempCountPair
std::vector< std::pair< std::string, float > > fGroupL1TempCountPair
Definition: FourVectorHLTOffline.h:315
ExpressReco_HICollisions_FallBack.etaMax
tuple etaMax
Definition: ExpressReco_HICollisions_FallBack.py:9492
FourVectorHLTOffline::fSelectedJets
reco::CaloJetCollection * fSelectedJets
Definition: FourVectorHLTOffline.h:201
objMonData::pushL1TriggerType
void pushL1TriggerType(int trigType)
Definition: FourVectorHLTOffline.h:737
FourVectorHLTOffline::PathInfo::PathInfo
PathInfo(std::string denomPathName, std::string pathName, std::string l1pathName, int l1ModuleIndex, std::string filterName, std::string processName, size_t type, float ptmin, float ptmax, float hltThreshold, float l1Threshold)
Definition: FourVectorHLTOffline.h:543
FourVectorHLTOffline::getHltThresholdFromName
int getHltThresholdFromName(const std::string &pathname)
Definition: FourVectorHLTOffline.cc:2551
FourVectorHLTOffline::fTriggerObj
edm::Handle< trigger::TriggerEvent > fTriggerObj
Definition: FourVectorHLTOffline.h:121
objMonData::setReco
void setReco(edm::Handle< T > offColl)
Definition: FourVectorHLTOffline.h:739
FourVectorHLTOffline::scalersSelect
MonitorElement * scalersSelect
Definition: FourVectorHLTOffline.h:330
objMonData::setTriggerType
void setTriggerType(std::vector< int > trigType)
Definition: FourVectorHLTOffline.h:734
FourVectorHLTOffline::v_ME_HLTAll_LS
std::vector< MonitorElement * > v_ME_HLTAll_LS
Definition: FourVectorHLTOffline.h:177
FourVectorHLTOffline::PathInfo::setFilterHistos
void setFilterHistos(MonitorElement *const filters)
Definition: FourVectorHLTOffline.h:340
v
mathSSE::Vec4< T > v
Definition: newBasic3DVector.h:210
trigger::Vids
std::vector< int > Vids
Definition: TriggerTypeDefs.h:24
edm::Run
Definition: Run.h:31
objMonData::offCollMu_
edm::Handle< reco::MuonCollection > offCollMu_
Definition: FourVectorHLTOffline.h:805
reco::CaloJetCollection
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
Definition: CaloJetCollection.h:16
FourVectorHLTOffline::selectMuons
void selectMuons(const edm::Handle< reco::MuonCollection > &muonHandle)
Definition: FourVectorHLTOffline.cc:2309
FourVectorHLTOffline::hasHLTPassed
bool hasHLTPassed(const std::string &pathname, const edm::TriggerNames &triggerNames)
Definition: FourVectorHLTOffline.cc:2295
FourVectorHLTOffline::PathInfo::getDenomPath
const std::string & getDenomPath(void) const
Definition: FourVectorHLTOffline.h:526
phi
Definition: DDAxes.h:10
objMonData::setBJetsFlag
void setBJetsFlag(bool flag)
Definition: FourVectorHLTOffline.h:773
Generated for CMSSW Reference Manual by doxygen 1.8.5