|
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
|
#ifndef JetTagPlotter_H
#define JetTagPlotter_H
#include "TH1F.h"
#include <string>
#include <vector>
#include "DQMOffline/RecoB/interface/FlavourHistorgrams.h"
#include "DQMOffline/RecoB/interface/BaseBTagPlotter.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/BTauReco/interface/JetTag.h"
#include "DQMServices/Core/interface/DQMStore.h"
class JetTagPlotter : public BaseBTagPlotter {
public:
JetTagPlotter(const std::string& tagName,
const EtaPtBin& etaPtBin,
const edm::ParameterSet& pSet,
unsigned int mc,
bool willFinalize,
DQMStore::IBooker& ibook,
bool doCTagPlots = false,
bool doDifferentialPlots = false,
double discrCut = -999.);
~JetTagPlotter() override;
void analyzeTag(); //added to fill the jet multiplicity on data
void analyzeTag(float w); //added to fill the jet multiplicity on mc
//void analyzeTag (const reco::JetTag & jetTag, const double & jec, const int & jetFlavour);
void analyzeTag(const reco::JetTag& jetTag, double jec, int jetFlavour, float w = 1);
//void analyzeTag (const reco::Jet & jet, const double & jec, const float& discriminator, const int& jetFlavour);
void analyzeTag(const reco::Jet& jet, double jec, float discriminator, int jetFlavour, float w = 1);
// final computation, plotting, printing .......
void finalize(DQMStore::IBooker& ibook_, DQMStore::IGetter& igetter_) override;
// get "2d" histograms for misid. vs. b-eff
EffPurFromHistos& getEffPurFromHistos() { return *effPurFromHistos_; }
void epsPlot(const std::string& name) override;
void psPlot(const std::string& name) override;
int nBinEffPur() const { return nBinEffPur_; }
double startEffPur() const { return startEffPur_; }
double endEffPur() const { return endEffPur_; }
protected:
// binning and bounds
// 1) for 'efficiency' versus discriminator cut histos
double discrStart_;
double discrEnd_;
int nBinEffPur_;
double startEffPur_;
double endEffPur_;
unsigned int mcPlots_;
bool willFinalize_;
bool doCTagPlots_;
// Differential plots: efficiency vs. variable for cut on discrimator > cutValue_
bool doDifferentialPlots_;
double cutValue_;
std::vector<int> nJets_;
// jet multiplicity
std::unique_ptr<FlavourHistograms<int>> jetMultiplicity_;
// for the misid vs. eff plots
std::unique_ptr<EffPurFromHistos> effPurFromHistos_;
std::unique_ptr<FlavourHistograms<int>> dJetFlav_;
// Discriminator: again with reasonable binning
std::unique_ptr<FlavourHistograms<double>> dDiscriminator_;
// reconstructed jet momentum
std::unique_ptr<FlavourHistograms<double>> dJetRecMomentum_;
// reconstructed jet transverse momentum
std::unique_ptr<FlavourHistograms<double>> dJetRecPt_;
// reconstructed jet eta
std::unique_ptr<FlavourHistograms<double>> dJetRecPseudoRapidity_;
// reconstructed jet phi
std::unique_ptr<FlavourHistograms<double>> dJetRecPhi_;
// jet Phi larger than requested discrimnator cut
std::unique_ptr<FlavourHistograms<double>> dJetPhiDiscrCut_;
// jet Eta larger than requested discrimnator cut
std::unique_ptr<FlavourHistograms<double>> dJetPseudoRapidityDiscrCut_;
};
#endif
|