|
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
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
|
// -*- C++ -*-
// MuonIsolationDQM.h
// Package: Muon Isolation DQM
// Class: MuonIsolationDQM
//
/*
Description: Muon Isolation DQM class
NOTE: The static member variable declarations *should* include the key word "static", but
I haven't found an elegant way to initalize the vectors. Static primatives (e.g. int,
float, ...) and simple static objects are easy to initialze. Outside of the class
decleration, you would write
int MuonIsolationDQM::CONST_INT = 5;
FooType MuonIsolationDQM::CONST_FOOT = Foo(constructor_argument);
but you can't do this if you want to, say, initalize a std::vector with a bunch of
different values. So, you can't make them static and you have to initialize them using
a member method. To keep it consistent, I've just initialized them all in the same
method, even the simple types.
*/
//
// Original Author: "C. Jess Riedel", UC Santa Barbara
// Created: Tue Jul 17 15:58:24 CDT 2007
//
//Base class
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
//Member types
#include "FWCore/Utilities/interface/InputTag.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
//Other include files
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
//----------------------------------------
//Forward declarations
class TH1;
class TH1I;
class TH1D;
class TH2;
class TProfile;
//------------------------------------------
// Class Declaration: MuonIsolationDQM
//--------------------------------------
class MuonIsolationDQM : public DQMEDAnalyzer {
//---------namespace and typedefs--------------
typedef edm::View<reco::Muon>::const_iterator MuonIterator;
typedef edm::RefToBase<reco::Muon> MuonBaseRef;
typedef edm::Handle<reco::IsoDepositMap> MuIsoDepHandle;
typedef const reco::IsoDeposit MuIsoDepRef;
public:
//---------methods----------------------------
explicit MuonIsolationDQM(const edm::ParameterSet&);
~MuonIsolationDQM() override;
void analyze(const edm::Event&, const edm::EventSetup&) override;
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
private:
//---------methods----------------------------
void InitStatics();
void RecordData(const reco::Muon& muon); //Fills Histograms with info from single muo
// void doPFIsoPlots(MuonIterator muon); //Fills Histograms with PF info from single muo (only for GLB)
void InitHistos(); //adds title, bin information to member histograms
void FillHistos(int); //Fills histograms with data
void FillNVtxHistos(int);
void NormalizeHistos(); //Normalize to number of muons
//----- helper methods
int GetNVtxBin(int);
TH1* GetTH1FromMonitorElement(MonitorElement* me);
//----------Static Variables---------------
int vtxBin_;
double vtxMin_;
double vtxMax_;
//Collection labels
edm::EDGetTokenT<reco::VertexCollection> theVertexCollectionLabel_;
edm::EDGetTokenT<edm::View<reco::Muon> > theMuonCollectionLabel_;
//root file name
std::string rootfilename;
// Directories within the rootfile
std::string dirName;
//Histogram parameters
static const int NUM_VARS = 48; // looking at R03 and R05. Total of 54 histos.
static const int NUM_VARS_2D = 10; // looking only at R03. Total of 8 TH2F.
static const int NUM_VARS_NVTX = 6;
double L_BIN_WIDTH; //large bins
double S_BIN_WIDTH; //small bins
int LOG_BINNING_ENABLED; //pseudo log binning for profile plots
int NUM_LOG_BINS;
double LOG_BINNING_RATIO;
bool requireGLBMuon;
bool requireSTAMuon;
bool requireTRKMuon;
std::string title_sam;
std::string title_cone;
// std::string title_cd;
std::vector<std::string> main_titles; //[NUM_VARS]
std::vector<std::string> axis_titles; //[NUM_VARS]
std::vector<std::string> names; //[NUM_VARS]
std::vector<std::vector<double> > param; //[NUM_VARS][3]
std::vector<int> isContinuous; //[NUM_VARS]
std::vector<std::string> titles_2D; //[NUM_VARS]
std::vector<std::string> names_2D; //[NUM_VARS]
std::vector<std::string> main_titles_NVtxs;
std::vector<std::string> names_NVtxs;
std::vector<std::string> axis_titles_NVtxs;
//---------------Dynamic Variables---------------------
//The Data
double theData[NUM_VARS];
double theData2D[NUM_VARS_2D];
double theDataNVtx[NUM_VARS_NVTX];
//Histograms
MonitorElement* h_nMuons;
std::vector<MonitorElement*> h_1D; //[NUM_VARS]
std::vector<MonitorElement*> h_2D; //[NUM_VARS_2D]
std::vector<MonitorElement*> h_1D_NVTX; //[NUM_VARS_NVTX]
// std::vector<MonitorElement*> cd_plots;//[NUM_VARS]
//Counters
int nEvents;
int nSTAMuons;
int nGLBMuons;
int nTRKMuons;
//enums for monitorElement
enum { NOAXIS, XAXIS, YAXIS, ZAXIS };
};
|