|
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
|
#ifndef RecoParticleFlow_Benchmark_Benchmark_h
#define RecoParticleFlow_Benchmark_Benchmark_h
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include <string>
class TH1;
class TH1F;
class TH2;
class TProfile;
class TH2F;
class TDirectory;
/// abstract base class
class Benchmark {
public:
typedef dqm::legacy::DQMStore DQMStore;
class PhaseSpace {
public:
int n;
float m;
float M;
PhaseSpace() : n(1), m(0), M(1) {}
PhaseSpace(int n, float m, float M) : n(n), m(m), M(M) {}
};
enum Mode { DEFAULT, DQMOFFLINE, VALIDATION };
Benchmark(Mode mode = DEFAULT)
: dir_(nullptr), mode_(mode), ptMin_(0), ptMax_(10e10), etaMin_(-10), etaMax_(10), phiMin_(-10), phiMax_(10) {}
virtual ~Benchmark() noexcept(false);
void setParameters(Mode mode) { mode_ = mode; }
void setRange(float ptMin, float ptMax, float etaMin, float etaMax, float phiMin, float phiMax) {
ptMin_ = ptMin;
ptMax_ = ptMax;
etaMin_ = etaMin;
etaMax_ = etaMax;
phiMin_ = phiMin;
phiMax_ = phiMax;
}
bool isInRange(float pt, float eta, float phi) const {
return (pt > ptMin_ && pt < ptMax_ && eta > etaMin_ && eta < etaMax_ && phi > phiMin_ && phi < phiMax_);
}
virtual void setDirectory(TDirectory *dir);
/// write to the TFile, in plain ROOT mode. No need to call this function in DQM mode
void write();
protected:
/// book a 1D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not.
//TH1F *book1D(const char *histname, const char *title,
/// book a 1D histogram, either through IBooker or plain root
TH1F *book1D(DQMStore::IBooker &b, const char *histname, const char *title, int nbins, float xmin, float xmax);
/// book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not.
//TH2F *book2D(const char *histname, const char *title,
/// book a 2D histogram, either through IBooker or plain root
TH2F *book2D(DQMStore::IBooker &b,
const char *histname,
const char *title,
int nbinsx,
float xmin,
float xmax,
int nbinsy,
float ymin,
float ymax);
/// book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not.
//TH2F *book2D(const char *histname, const char *title,
/// book a 2D histogram, either through IBooker or plain root
TH2F *book2D(DQMStore::IBooker &b,
const char *histname,
const char *title,
int nbinsx,
float *xbins,
int nbinsy,
float ymin,
float ymax);
/// book a TProfile histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not.
//TProfile *bookProfile(const char *histname, const char *title,
/// book a TProfile, either through IBooker or plain root
TProfile *bookProfile(DQMStore::IBooker &b,
const char *histname,
const char *title,
int nbinsx,
float xmin,
float xmax,
float ymin,
float ymax,
const char *option);
/// book a TProfile histogram, either with DQM or plain root depending if DQM_
/// has been initialized in a child analyzer or not.
// TProfile *bookProfile(const char *histname, const char *title,
/// book a TProfile, either through IBooker or plain root
TProfile *bookProfile(DQMStore::IBooker &b,
const char *histname,
const char *title,
int nbinsx,
float *xbins,
float ymin,
float ymax,
const char *option);
TDirectory *dir_;
Mode mode_;
float ptMin_;
float ptMax_;
float etaMin_;
float etaMax_;
float phiMin_;
float phiMax_;
};
#endif
|