|
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
|
#ifndef RecoParticleFlow_Benchmark_MatchCandidateBenchmark_h
#define RecoParticleFlow_Benchmark_MatchCandidateBenchmark_h
#include "DQMOffline/PFTau/interface/Benchmark.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <vector>
/// To plot Candidate quantities
class MatchCandidateBenchmark : public Benchmark {
public:
MatchCandidateBenchmark(Mode mode);
~MatchCandidateBenchmark() override;
/// book histograms
void setup(DQMStore::IBooker &b);
void setup(DQMStore::IBooker &b, const edm::ParameterSet ¶meterSet);
/// fill histograms with a given particle
void fillOne(const reco::Candidate &candidate, const reco::Candidate &matchedCandidate);
void fillOne(const reco::Candidate &candidate,
const reco::Candidate &matchedCandidate,
const edm::ParameterSet ¶meterSet);
protected:
TH2F *delta_et_Over_et_VS_et_;
TH2F *delta_et_VS_et_;
TH2F *delta_eta_VS_et_;
TH2F *delta_phi_VS_et_;
TH2F *BRdelta_et_Over_et_VS_et_;
TH2F *ERdelta_et_Over_et_VS_et_;
std::vector<TH1F *> pTRes_;
std::vector<TH1F *> BRpTRes_;
std::vector<TH1F *> ERpTRes_;
std::vector<float> ptBins_;
bool histogramBooked_;
double eta_min_barrel_;
double eta_max_barrel_;
double eta_min_endcap_;
double eta_max_endcap_;
private:
void computePtBins(const edm::ParameterSet &, const edm::ParameterSet &);
bool inEtaRange(double, bool);
inline bool inBarrelRange(double value) { return inEtaRange(value, true); }
inline bool inEndcapRange(double value) { return inEtaRange(value, false); }
};
#endif
|