|
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
|
#ifndef IPTagPlotter_H
#define IPTagPlotter_H
#include "DQMOffline/RecoB/interface/BaseTagInfoPlotter.h"
#include "DQMOffline/RecoB/interface/TrackIPHistograms.h"
#include "DQMOffline/RecoB/interface/FlavourHistorgrams2D.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/BTauReco/interface/IPTagInfo.h"
template <class Container, class Base>
class IPTagPlotter : public BaseTagInfoPlotter {
public:
IPTagPlotter(const std::string& tagName,
const EtaPtBin& etaPtBin,
const edm::ParameterSet& pSet,
unsigned int mc,
bool wf,
DQMStore::IBooker& ibook);
~IPTagPlotter() override;
void analyzeTag(const reco::BaseTagInfo* baseTagInfo, double jec, int jetFlavour, float w = 1) override;
void finalize(DQMStore::IBooker& ibook_, DQMStore::IGetter& igetter_) override;
void epsPlot(const std::string& name) override;
void psPlot(const std::string& name) override;
reco::TrackBase::TrackQuality highestTrackQual(const reco::Track* track) const;
private:
int nBinEffPur_;
double startEffPur_;
double endEffPur_;
unsigned int mcPlots_;
bool willFinalize_;
bool makeQualityPlots_;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosSig3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosSig2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosErr3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosErr2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosVal3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosVal2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosDecayLengthVal2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosDecayLengthVal3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosJetDistVal2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosJetDistVal3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosJetDistSign2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosJetDistSign3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosTkNChiSqr2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosTkNChiSqr3D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosTkPt2D;
std::vector<std::unique_ptr<TrackIPHistograms<double>>> tkcntHistosTkPt3D;
std::vector<std::unique_ptr<TrackIPHistograms<int>>> tkcntHistosTkNHits2D;
std::vector<std::unique_ptr<TrackIPHistograms<int>>> tkcntHistosTkNHits3D;
std::vector<std::unique_ptr<TrackIPHistograms<int>>> tkcntHistosTkNPixelHits2D;
std::vector<std::unique_ptr<TrackIPHistograms<int>>> tkcntHistosTkNPixelHits3D;
std::unique_ptr<FlavourHistograms<int>> trkNbr3D, trkNbr2D;
double lowerIPSBound, upperIPSBound, lowerIPBound, upperIPBound, lowerIPEBound, upperIPEBound;
int nBinsIPS, nBinsIP, nBinsIPE;
double minDecayLength, maxDecayLength, minJetDistance, maxJetDistance;
std::vector<std::unique_ptr<EffPurFromHistos>> effPurFromHistos;
std::vector<std::unique_ptr<TrackIPHistograms<float>>> tkcntHistosProb3D;
std::vector<std::unique_ptr<TrackIPHistograms<float>>> tkcntHistosProb2D;
std::unique_ptr<TrackIPHistograms<float>> tkcntHistosTkProbIPneg2D, tkcntHistosTkProbIPpos2D;
std::unique_ptr<TrackIPHistograms<float>> tkcntHistosTkProbIPneg3D, tkcntHistosTkProbIPpos3D;
std::unique_ptr<TrackIPHistograms<double>> ghostTrackWeightHisto;
std::unique_ptr<TrackIPHistograms<double>> ghostTrackDistanceValuHisto, ghostTrackDistanceSignHisto;
std::unique_ptr<FlavourHistograms<int>> trackQualHisto;
std::unique_ptr<FlavourHistograms<int>> selectedTrackQualHisto;
std::unique_ptr<FlavourHistograms2D<double, int>> trackMultVsJetPtHisto;
std::unique_ptr<FlavourHistograms2D<double, int>> selectedTrackMultVsJetPtHisto;
};
#include "DQMOffline/RecoB/interface/IPTagPlotter_cc.h"
#endif
|