|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
#ifndef DQMOffline_Trigger_HLTMuonMatchAndPlot_H
#define DQMOffline_Trigger_HLTMuonMatchAndPlot_H
/** \class HLTMuonMatchAndPlot
* Match reconstructed muons to HLT objects and plot efficiencies.
*
* Note that this is not a true EDAnalyzer;
*
* Documentation available on the CMS TWiki:
* https://twiki.cern.ch/twiki/bin/view/CMS/MuonHLTOfflinePerformance
*
*
* \author J. Slaunwhite, Jeff Klukas
*/
// Base Class Headers
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/Math/interface/deltaR.h"
#include <vector>
#include "TFile.h"
#include "TNtuple.h"
#include "TString.h"
#include "TPRegexp.h"
//////////////////////////////////////////////////////////////////////////////
//////// Typedefs and Constants //////////////////////////////////////////////
using LorentzVector = math::XYZTLorentzVector;
const double NOMATCH = 999.;
const std::string EFFICIENCY_SUFFIXES[2] = {"denom", "numer"};
//////////////////////////////////////////////////////////////////////////////
//////// HLTMuonMatchAndPlot Class Definition ////////////////////////////////
class HLTMuonMatchAndPlot {
public:
typedef dqm::legacy::DQMStore DQMStore;
typedef dqm::legacy::MonitorElement MonitorElement;
/// Constructor
HLTMuonMatchAndPlot(const edm::ParameterSet &, std::string, std::string, bool);
// Analyzer Methods
void beginRun(DQMStore::IBooker &, const edm::Run &, const edm::EventSetup &);
void analyze(edm::Handle<reco::MuonCollection> &,
edm::Handle<reco::BeamSpot> &,
edm::Handle<reco::VertexCollection> &,
edm::Handle<trigger::TriggerEvent> &,
edm::Handle<edm::TriggerResults> &,
const edm::TriggerNames &);
void endRun(const edm::Run &, const edm::EventSetup &);
// Helper Methods
bool fillEdges(size_t &nBins, float *&edges, const std::vector<double> &binning);
template <class T>
void fillMapFromPSet(std::map<std::string, T> &, const edm::ParameterSet &, const std::string &);
template <class T1, class T2>
std::vector<size_t> matchByDeltaR(const std::vector<T1> &, const std::vector<T2> &, const double maxDeltaR = NOMATCH);
private:
// Internal Methods
void book1D(DQMStore::IBooker &, std::string, const std::string &, std::string);
void book2D(DQMStore::IBooker &, const std::string &, const std::string &, const std::string &, const std::string &);
reco::MuonCollection selectedMuons(const reco::MuonCollection &, const reco::BeamSpot &, bool);
trigger::TriggerObjectCollection selectedTriggerObjects(const trigger::TriggerObjectCollection &,
const trigger::TriggerEvent &);
// Input from Configuration File
std::string hltProcessName_;
std::string destination_;
std::vector<std::string> requiredTriggers_;
std::map<std::string, std::vector<double> > binParams_;
std::map<std::string, double> plotCuts_;
edm::ParameterSet targetParams_;
edm::ParameterSet probeParams_;
// Member Variables
std::string triggerLevel_;
unsigned int cutMinPt_;
std::string hltPath_;
std::string moduleLabel_;
bool isLastFilter_;
std::map<std::string, MonitorElement *> hists_;
double targetMuonEtaMax_;
double targetMuonEtaMin_;
bool targetIsMuonGlb_;
double targetZ0Cut_;
double targetD0Cut_;
double targetptCutZ_;
double targetptCutJpsi_;
double probeMuonEtaMax_;
double probeMuonEtaMin_;
bool probeIsMuonGlb_;
double probeZ0Cut_;
double probeD0Cut_;
double triggerEtaMaxCut_;
double triggerEtaMinCut_;
};
#endif
|