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