CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
IPTagPlotter< Container, Base > Class Template Reference

#include <IPTagPlotter.h>

Inheritance diagram for IPTagPlotter< Container, Base >:
BaseTagInfoPlotter BaseBTagPlotter

Public Member Functions

void analyzeTag (const reco::BaseTagInfo *baseTagInfo, const double &jec, const int &jetFlavour)
 
void analyzeTag (const reco::BaseTagInfo *baseTagInfo, const double &jec, const int &jetFlavour, const float &w)
 
void epsPlot (const std::string &name)
 
virtual void finalize (DQMStore::IBooker &ibook_, DQMStore::IGetter &igetter_)
 
reco::TrackBase::TrackQuality highestTrackQual (const reco::Track *track) const
 
 IPTagPlotter (const std::string &tagName, const EtaPtBin &etaPtBin, const edm::ParameterSet &pSet, const unsigned int &mc, const bool &wf, DQMStore::IBooker &ibook)
 
void psPlot (const std::string &name)
 
 ~IPTagPlotter ()
 
- Public Member Functions inherited from BaseTagInfoPlotter
virtual void analyzeTag (const std::vector< const reco::BaseTagInfo * > &tagInfos, const double &jec, const int &jetFlavour)
 
virtual void analyzeTag (const std::vector< const reco::BaseTagInfo * > &tagInfos, const double &jec, const int &jetFlavour, const float &w)
 
 BaseTagInfoPlotter (const std::string &tagName, const EtaPtBin &etaPtBin)
 
virtual void setEventSetup (const edm::EventSetup &setup)
 
virtual std::vector< std::string > tagInfoRequirements () const
 
virtual ~BaseTagInfoPlotter ()
 
- Public Member Functions inherited from BaseBTagPlotter
 BaseBTagPlotter (const std::string &tagName, const EtaPtBin &etaPtBin)
 
const EtaPtBin & etaPtBin ()
 
virtual ~BaseBTagPlotter ()
 

Private Attributes

EffPurFromHistos * effPurFromHistos [4]
 
double endEffPur_
 
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
 
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
 
TrackIPHistograms< double > * ghostTrackWeightHisto
 
double lowerIPBound
 
double lowerIPEBound
 
double lowerIPSBound
 
bool makeQualityPlots_
 
double maxDecayLength
 
double maxJetDistance
 
unsigned int mcPlots_
 
double minDecayLength
 
double minJetDistance
 
int nBinEffPur_
 
int nBinsIP
 
int nBinsIPE
 
int nBinsIPS
 
FlavourHistograms2D< double,
int > * 		selectedTrackMultVsJetPtHisto
 
FlavourHistograms< int > * selectedTrackQualHisto
 
double startEffPur_
 
TrackIPHistograms< double > * tkcntHistosDecayLengthVal2D [5]
 
TrackIPHistograms< double > * tkcntHistosDecayLengthVal3D [5]
 
TrackIPHistograms< double > * tkcntHistosErr2D [5]
 
TrackIPHistograms< double > * tkcntHistosErr3D [5]
 
TrackIPHistograms< double > * tkcntHistosJetDistSign2D [5]
 
TrackIPHistograms< double > * tkcntHistosJetDistSign3D [5]
 
TrackIPHistograms< double > * tkcntHistosJetDistVal2D [5]
 
TrackIPHistograms< double > * tkcntHistosJetDistVal3D [5]
 
TrackIPHistograms< float > * tkcntHistosProb2D [5]
 
TrackIPHistograms< float > * tkcntHistosProb3D [5]
 
TrackIPHistograms< double > * tkcntHistosSig2D [5]
 
TrackIPHistograms< double > * tkcntHistosSig3D [5]
 
TrackIPHistograms< double > * tkcntHistosTkNChiSqr2D [5]
 
TrackIPHistograms< double > * tkcntHistosTkNChiSqr3D [5]
 
TrackIPHistograms< int > * tkcntHistosTkNHits2D [5]
 
TrackIPHistograms< int > * tkcntHistosTkNHits3D [5]
 
TrackIPHistograms< int > * tkcntHistosTkNPixelHits2D [5]
 
TrackIPHistograms< int > * tkcntHistosTkNPixelHits3D [5]
 
TrackIPHistograms< float > * tkcntHistosTkProbIPneg2D
 
TrackIPHistograms< float > * tkcntHistosTkProbIPneg3D
 
TrackIPHistograms< float > * tkcntHistosTkProbIPpos2D
 
TrackIPHistograms< float > * tkcntHistosTkProbIPpos3D
 
TrackIPHistograms< double > * tkcntHistosTkPt2D [5]
 
TrackIPHistograms< double > * tkcntHistosTkPt3D [5]
 
TrackIPHistograms< double > * tkcntHistosVal2D [5]
 
TrackIPHistograms< double > * tkcntHistosVal3D [5]
 
FlavourHistograms2D< double,
int > * 		trackMultVsJetPtHisto
 
FlavourHistograms< int > * trackQualHisto
 
FlavourHistograms< int > * trkNbr2D
 
FlavourHistograms< int > * trkNbr3D
 
double upperIPBound
 
double upperIPEBound
 
double upperIPSBound
 
bool willFinalize_
 

Additional Inherited Members

- Protected Attributes inherited from BaseBTagPlotter
const EtaPtBin etaPtBin_
 
const std::string tagName_
 
const std::string theExtensionString
 

Detailed Description

template<class Container, class Base>
class IPTagPlotter< Container, Base >

Definition at line 11 of file IPTagPlotter.h.

Constructor & Destructor Documentation

template<class Container , class Base >
IPTagPlotter< Container, Base >::IPTagPlotter ( const std::string &  tagName,
const EtaPtBin &  etaPtBin,
const edm::ParameterSet &  pSet,
const unsigned int &  mc,
const bool &  wf,
DQMStore::IBooker &  ibook 
)

Definition at line 8 of file IPTagPlotter_cc.h.

References IPTagPlotter< Container, Base >::ghostTrackDistanceSignHisto, IPTagPlotter< Container, Base >::ghostTrackDistanceValuHisto, IPTagPlotter< Container, Base >::ghostTrackWeightHisto, IPTagPlotter< Container, Base >::lowerIPBound, IPTagPlotter< Container, Base >::lowerIPEBound, IPTagPlotter< Container, Base >::lowerIPSBound, IPTagPlotter< Container, Base >::makeQualityPlots_, IPTagPlotter< Container, Base >::nBinsIP, IPTagPlotter< Container, Base >::nBinsIPE, IPTagPlotter< Container, Base >::nBinsIPS, IPTagPlotter< Container, Base >::selectedTrackMultVsJetPtHisto, IPTagPlotter< Container, Base >::selectedTrackQualHisto, AlCaHLTBitMon_QueryRunRegistry::string, BaseBTagPlotter::theExtensionString, IPTagPlotter< Container, Base >::tkcntHistosDecayLengthVal2D, IPTagPlotter< Container, Base >::tkcntHistosDecayLengthVal3D, IPTagPlotter< Container, Base >::tkcntHistosErr2D, IPTagPlotter< Container, Base >::tkcntHistosErr3D, IPTagPlotter< Container, Base >::tkcntHistosJetDistVal2D, IPTagPlotter< Container, Base >::tkcntHistosJetDistVal3D, IPTagPlotter< Container, Base >::tkcntHistosProb2D, IPTagPlotter< Container, Base >::tkcntHistosProb3D, IPTagPlotter< Container, Base >::tkcntHistosSig2D, IPTagPlotter< Container, Base >::tkcntHistosSig3D, IPTagPlotter< Container, Base >::tkcntHistosTkNChiSqr2D, IPTagPlotter< Container, Base >::tkcntHistosTkNChiSqr3D, IPTagPlotter< Container, Base >::tkcntHistosTkNHits2D, IPTagPlotter< Container, Base >::tkcntHistosTkNHits3D, IPTagPlotter< Container, Base >::tkcntHistosTkNPixelHits2D, IPTagPlotter< Container, Base >::tkcntHistosTkNPixelHits3D, IPTagPlotter< Container, Base >::tkcntHistosTkProbIPneg2D, IPTagPlotter< Container, Base >::tkcntHistosTkProbIPneg3D, IPTagPlotter< Container, Base >::tkcntHistosTkProbIPpos2D, IPTagPlotter< Container, Base >::tkcntHistosTkProbIPpos3D, IPTagPlotter< Container, Base >::tkcntHistosTkPt2D, IPTagPlotter< Container, Base >::tkcntHistosTkPt3D, IPTagPlotter< Container, Base >::tkcntHistosVal2D, IPTagPlotter< Container, Base >::tkcntHistosVal3D, IPTagPlotter< Container, Base >::trackMultVsJetPtHisto, IPTagPlotter< Container, Base >::trackQualHisto, IPTagPlotter< Container, Base >::trkNbr2D, IPTagPlotter< Container, Base >::trkNbr3D, IPTagPlotter< Container, Base >::upperIPBound, IPTagPlotter< Container, Base >::upperIPEBound, IPTagPlotter< Container, Base >::upperIPSBound, and IPTagPlotter< Container, Base >::willFinalize_.

10  :
11  BaseTagInfoPlotter(tagName, etaPtBin),
12  nBinEffPur_(pSet.getParameter<int>("nBinEffPur")),
13  startEffPur_(pSet.getParameter<double>("startEffPur")),
14  endEffPur_(pSet.getParameter<double>("endEffPur")),
15  mcPlots_(mc), willFinalize_(wf),
16  makeQualityPlots_(pSet.getParameter<bool>("QualityPlots")),
17  lowerIPSBound(pSet.getParameter<double>("LowerIPSBound")),
18  upperIPSBound(pSet.getParameter<double>("UpperIPSBound")),
19  lowerIPBound(pSet.getParameter<double>("LowerIPBound")),
20  upperIPBound(pSet.getParameter<double>("UpperIPBound")),
21  lowerIPEBound(pSet.getParameter<double>("LowerIPEBound")),
22  upperIPEBound(pSet.getParameter<double>("UpperIPEBound")),
23  nBinsIPS(pSet.getParameter<int>("NBinsIPS")),
24  nBinsIP(pSet.getParameter<int>("NBinsIP")),
25  nBinsIPE(pSet.getParameter<int>("NBinsIPE")),
26  minDecayLength(pSet.getParameter<double>("MinDecayLength")),
27  maxDecayLength(pSet.getParameter<double>("MaxDecayLength")),
28  minJetDistance(pSet.getParameter<double>("MinJetDistance")),
29  maxJetDistance(pSet.getParameter<double>("MaxJetDistance"))
30 {
31  const std::string trackIPDir(theExtensionString.substr(1));
32 
33  if (willFinalize_) return;
34 
35  trkNbr3D = new TrackIPHistograms<int>
36  ("selTrksNbr_3D" + theExtensionString, "Number of selected tracks for 3D IPS", 31, -0.5, 30.5,
37  false, true, true, "b", trackIPDir ,mc, makeQualityPlots_, ibook_);
38 
39  trkNbr2D = new TrackIPHistograms<int>
40  ("selTrksNbr_2D" + theExtensionString, "Number of selected tracks for 2D IPS", 31, -0.5, 30.5,
41  false, true, true, "b", trackIPDir ,mc, makeQualityPlots_, ibook_);
42 
43  // IP significance
44  // 3D
45  tkcntHistosSig3D[4] = new TrackIPHistograms<double>
46  ("ips_3D" + theExtensionString, "3D IP significance",
47  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
48 
49  tkcntHistosSig3D[0] = new TrackIPHistograms<double>
50  ("ips1_3D" + theExtensionString, "3D IP significance 1.trk",
51  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir,mc, makeQualityPlots_, ibook_) ;
52 
53  tkcntHistosSig3D[1] = new TrackIPHistograms<double>
54  ("ips2_3D" + theExtensionString, "3D IP significance 2.trk",
55  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
56 
57  tkcntHistosSig3D[2] = new TrackIPHistograms<double>
58  ("ips3_3D" + theExtensionString, "3D IP significance 3.trk",
59  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
60 
61  tkcntHistosSig3D[3] = new TrackIPHistograms<double>
62  ("ips4_3D" + theExtensionString, "3D IP significance 4.trk",
63  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
64 
65  //2D
66  tkcntHistosSig2D[4] = new TrackIPHistograms<double>
67  ("ips_2D" + theExtensionString, "2D IP significance",
68  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
69 
70  tkcntHistosSig2D[0] = new TrackIPHistograms<double>
71  ("ips1_2D" + theExtensionString, "2D IP significance 1.trk",
72  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
73 
74  tkcntHistosSig2D[1] = new TrackIPHistograms<double>
75  ("ips2_2D" + theExtensionString, "2D IP significance 2.trk",
76  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
77 
78  tkcntHistosSig2D[2] = new TrackIPHistograms<double>
79  ("ips3_2D" + theExtensionString, "2D IP significance 3.trk",
80  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
81 
82  tkcntHistosSig2D[3] = new TrackIPHistograms<double>
83  ("ips4_2D" + theExtensionString, "2D IP significance 4.trk",
84  nBinsIPS, lowerIPSBound, upperIPSBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
85 
86  // IP value
87  //3D
88  tkcntHistosVal3D[4] = new TrackIPHistograms<double>
89  ("ip_3D" + theExtensionString, "3D IP value",
90  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
91 
92  tkcntHistosVal3D[0] = new TrackIPHistograms<double>
93  ("ip1_3D" + theExtensionString, "3D IP value 1.trk",
94  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
95 
96  tkcntHistosVal3D[1] = new TrackIPHistograms<double>
97  ("ip2_3D" + theExtensionString, "3D IP value 2.trk",
98  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
99 
100  tkcntHistosVal3D[2] = new TrackIPHistograms<double>
101  ("ip3_3D" + theExtensionString, "3D IP value 3.trk",
102  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
103 
104  tkcntHistosVal3D[3] = new TrackIPHistograms<double>
105  ("ip4_3D" + theExtensionString, "3D IP value 4.trk",
106  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
107 
108  //2D
109  tkcntHistosVal2D[4] = new TrackIPHistograms<double>
110  ("ip_2D" + theExtensionString, "2D IP value",
111  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
112 
113  tkcntHistosVal2D[0] = new TrackIPHistograms<double>
114  ("ip1_2D" + theExtensionString, "2D IP value 1.trk",
115  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
116 
117  tkcntHistosVal2D[1] = new TrackIPHistograms<double>
118  ("ip2_2D" + theExtensionString, "2D IP value 2.trk",
119  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
120 
121  tkcntHistosVal2D[2] = new TrackIPHistograms<double>
122  ("ip3_2D" + theExtensionString, "2D IP value 3.trk",
123  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
124 
125  tkcntHistosVal2D[3] = new TrackIPHistograms<double>
126  ("ip4_2D" + theExtensionString, "2D IP value 4.trk",
127  nBinsIP, lowerIPBound, upperIPBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
128 
129 
130  // IP error
131  // 3D
132  tkcntHistosErr3D[4] = new TrackIPHistograms<double>
133  ("ipe_3D" + theExtensionString, "3D IP error",
134  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
135 
136  tkcntHistosErr3D[0] = new TrackIPHistograms<double>
137  ("ipe1_3D" + theExtensionString, "3D IP error 1.trk",
138  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
139 
140  tkcntHistosErr3D[1] = new TrackIPHistograms<double>
141  ("ipe2_3D" + theExtensionString, "3D IP error 2.trk",
142  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
143 
144  tkcntHistosErr3D[2] = new TrackIPHistograms<double>
145  ("ipe3_3D" + theExtensionString, "3D IP error 3.trk",
146  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
147 
148  tkcntHistosErr3D[3] = new TrackIPHistograms<double>
149  ("ipe4_3D" + theExtensionString, "3D IP error 4.trk",
150  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
151 
152  //2D
153  tkcntHistosErr2D[4] = new TrackIPHistograms<double>
154  ("ipe_2D" + theExtensionString, "2D IP error",
155  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
156 
157  tkcntHistosErr2D[0] = new TrackIPHistograms<double>
158  ("ipe1_2D" + theExtensionString, "2D IP error 1.trk",
159  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
160 
161  tkcntHistosErr2D[1] = new TrackIPHistograms<double>
162  ("ipe2_2D" + theExtensionString, "2D IP error 2.trk",
163  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
164 
165  tkcntHistosErr2D[2] = new TrackIPHistograms<double>
166  ("ipe3_2D" + theExtensionString, "2D IP error 3.trk",
167  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
168 
169  tkcntHistosErr2D[3] = new TrackIPHistograms<double>
170  ("ipe4_2D" + theExtensionString, "2D IP error 4.trk",
171  nBinsIPE, lowerIPEBound, upperIPEBound, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
172 
173  // decay length
174  tkcntHistosDecayLengthVal2D[4] = new TrackIPHistograms<double>
175  ("decLen_2D" + theExtensionString, "Decay Length 2D",
176  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
177 
178  tkcntHistosDecayLengthVal2D[0] = new TrackIPHistograms<double>
179  ("decLen1_2D" + theExtensionString, "2D Decay Length 1.trk",
180  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
181 
182  tkcntHistosDecayLengthVal2D[1] = new TrackIPHistograms<double>
183  ("decLen2_2D" + theExtensionString, "2D Decay Length 2.trk",
184  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
185 
186  tkcntHistosDecayLengthVal2D[2] = new TrackIPHistograms<double>
187  ("decLen3_2D" + theExtensionString, "2D Decay Length 3.trk",
188  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
189 
190  tkcntHistosDecayLengthVal2D[3] = new TrackIPHistograms<double>
191  ("decLen4_2D" + theExtensionString, "2D Decay Length 4.trk",
192  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
193 
194  tkcntHistosDecayLengthVal3D[4] = new TrackIPHistograms<double>
195  ("decLen_3D" + theExtensionString, "3D Decay Length",
196  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
197 
198  tkcntHistosDecayLengthVal3D[0] = new TrackIPHistograms<double>
199  ("decLen1_3D" + theExtensionString, "3D Decay Length 1.trk",
200  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
201 
202  tkcntHistosDecayLengthVal3D[1] = new TrackIPHistograms<double>
203  ("decLen2_3D" + theExtensionString, "3D Decay Length 2.trk",
204  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
205 
206  tkcntHistosDecayLengthVal3D[2] = new TrackIPHistograms<double>
207  ("decLen3_3D" + theExtensionString, "3D Decay Length 3.trk",
208  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
209 
210  tkcntHistosDecayLengthVal3D[3] = new TrackIPHistograms<double>
211  ("decLen4_3D" + theExtensionString, "3D Decay Length 4.trk",
212  50, 0.0, 5.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
213 
214  // jet distance
215  tkcntHistosJetDistVal2D[4] = new TrackIPHistograms<double>
216  ("jetDist_2D" + theExtensionString, "JetDistance 2D",
217  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
218 
219  tkcntHistosJetDistVal2D[0] = new TrackIPHistograms<double>
220  ("jetDist1_2D" + theExtensionString, "JetDistance 2D 1.trk",
221  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
222 
223  tkcntHistosJetDistVal2D[1] = new TrackIPHistograms<double>
224  ("jetDist2_2D" + theExtensionString, "JetDistance 2D 2.trk",
225  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
226 
227  tkcntHistosJetDistVal2D[2] = new TrackIPHistograms<double>
228  ("jetDist3_2D" + theExtensionString, "JetDistance 2D 3.trk",
229  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
230 
231  tkcntHistosJetDistVal2D[3] = new TrackIPHistograms<double>
232  ("jetDist4_2D" + theExtensionString, "JetDistance 2D 4.trk",
233  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
234 
235  tkcntHistosJetDistVal3D[4] = new TrackIPHistograms<double>
236  ("jetDist_3D" + theExtensionString, "JetDistance 3D",
237  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
238 
239  tkcntHistosJetDistVal3D[0] = new TrackIPHistograms<double>
240  ("jetDist1_3D" + theExtensionString, "JetDistance 3D 1.trk",
241  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
242 
243  tkcntHistosJetDistVal3D[1] = new TrackIPHistograms<double>
244  ("jetDist2_3D" + theExtensionString, "JetDistance 3D 2.trk",
245  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
246 
247  tkcntHistosJetDistVal3D[2] = new TrackIPHistograms<double>
248  ("jetDist3_3D" + theExtensionString, "JetDistance 3D 3.trk",
249  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
250 
251  tkcntHistosJetDistVal3D[3] = new TrackIPHistograms<double>
252  ("jetDist4_3D" + theExtensionString, "JetDistance 3D 4.trk",
253  50, -0.1, 0.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
254  /* // the jet distance significance return always 0 => disabled
255  tkcntHistosJetDistSign2D[4] = new TrackIPHistograms<double>
256  ("jetDistSig_2D" + theExtensionString, "JetDistance Sign 2D",
257  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
258 
259  tkcntHistosJetDistSign2D[0] = new TrackIPHistograms<double>
260  ("jetDistSig1_2D" + theExtensionString, "JetDistance Sign 2D 1.trk",
261  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
262 
263  tkcntHistosJetDistSign2D[1] = new TrackIPHistograms<double>
264  ("jetDistSig2_2D" + theExtensionString, "JetDistance Sign 2D 2.trk",
265  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
266 
267  tkcntHistosJetDistSign2D[2] = new TrackIPHistograms<double>
268  ("jetDistSig3_2D" + theExtensionString, "JetDistance Sign 2D 3.trk",
269  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
270 
271  tkcntHistosJetDistSign2D[3] = new TrackIPHistograms<double>
272  ("jetDistSig4_2D" + theExtensionString, "JetDistance Sign 2D 4.trk",
273  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
274 
275  tkcntHistosJetDistSign3D[4] = new TrackIPHistograms<double>
276  ("jetDistSig_3D" + theExtensionString, "JetDistance Sign 3D",
277  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
278 
279  tkcntHistosJetDistSign3D[0] = new TrackIPHistograms<double>
280  ("jetDistSig1_3D" + theExtensionString, "JetDistance Sign 3D 1.trk",
281  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
282 
283  tkcntHistosJetDistSign3D[1] = new TrackIPHistograms<double>
284  ("jetDistSig2_3D" + theExtensionString, "JetDistance Sign 3D 2.trk",
285  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
286 
287  tkcntHistosJetDistSign3D[2] = new TrackIPHistograms<double>
288  ("jetDistSig3_3D" + theExtensionString, "JetDistance Sign 3D 3.trk",
289  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
290 
291  tkcntHistosJetDistSign3D[3] = new TrackIPHistograms<double>
292  ("jetDistSig4_3D" + theExtensionString, "JetDistance Sign 3D 4.trk",
293  100, -0.001, 0.001, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
294  */
295  // track chi-squared
296  tkcntHistosTkNChiSqr2D[4] = new TrackIPHistograms<double>
297  ("tkNChiSqr_2D" + theExtensionString, "Normalized Chi Squared 2D",
298  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
299 
300  tkcntHistosTkNChiSqr2D[0] = new TrackIPHistograms<double>
301  ("tkNChiSqr1_2D" + theExtensionString, "Normalized Chi Squared 2D 1.trk",
302  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
303 
304  tkcntHistosTkNChiSqr2D[1] = new TrackIPHistograms<double>
305  ("tkNChiSqr2_2D" + theExtensionString, "Normalized Chi Squared 2D 2.trk",
306  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
307 
308  tkcntHistosTkNChiSqr2D[2] = new TrackIPHistograms<double>
309  ("tkNChiSqr3_2D" + theExtensionString, "Normalized Chi Squared 2D 3.trk",
310  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
311 
312  tkcntHistosTkNChiSqr2D[3] = new TrackIPHistograms<double>
313  ("tkNChiSqr4_2D" + theExtensionString, "Normalized Chi Squared 2D 4.trk",
314  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
315 
316  tkcntHistosTkNChiSqr3D[4] = new TrackIPHistograms<double>
317  ("tkNChiSqr_3D" + theExtensionString, "Normalized Chi Squared 3D",
318  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
319 
320  tkcntHistosTkNChiSqr3D[0] = new TrackIPHistograms<double>
321  ("tkNChiSqr1_3D" + theExtensionString, "Normalized Chi Squared 3D 1.trk",
322  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
323 
324  tkcntHistosTkNChiSqr3D[1] = new TrackIPHistograms<double>
325  ("tkNChiSqr2_3D" + theExtensionString, "Normalized Chi Squared 3D 2.trk",
326  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
327 
328  tkcntHistosTkNChiSqr3D[2] = new TrackIPHistograms<double>
329  ("tkNChiSqr3_3D" + theExtensionString, "Normalized Chi Squared 3D 3.trk",
330  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
331 
332  tkcntHistosTkNChiSqr3D[3] = new TrackIPHistograms<double>
333  ("tkNChiSqr4_3D" + theExtensionString, "Normalized Chi Squared 3D 4.trk",
334  50, -0.1, 10.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
335 
336  // track pT
337  tkcntHistosTkPt2D[4] = new TrackIPHistograms<double>
338  ("tkPt_2D" + theExtensionString, "Track Pt 2D",
339  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
340 
341  tkcntHistosTkPt2D[0] = new TrackIPHistograms<double>
342  ("tkPt1_2D" + theExtensionString, "Track Pt 2D 1.trk",
343  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
344 
345  tkcntHistosTkPt2D[1] = new TrackIPHistograms<double>
346  ("tkPt2_2D" + theExtensionString, "Track Pt 2D 2.trk",
347  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
348 
349  tkcntHistosTkPt2D[2] = new TrackIPHistograms<double>
350  ("tkPt3_2D" + theExtensionString, "Track Pt 2D 3.trk",
351  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
352 
353  tkcntHistosTkPt2D[3] = new TrackIPHistograms<double>
354  ("tkPt4_2D" + theExtensionString, "Track Pt 2D 4.trk",
355  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
356 
357  tkcntHistosTkPt3D[4] = new TrackIPHistograms<double>
358  ("tkPt_3D" + theExtensionString, "Track Pt 3D",
359  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
360 
361  tkcntHistosTkPt3D[0] = new TrackIPHistograms<double>
362  ("tkPt1_3D" + theExtensionString, "Track Pt 3D 1.trk",
363  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
364 
365  tkcntHistosTkPt3D[1] = new TrackIPHistograms<double>
366  ("tkPt2_3D" + theExtensionString, "Track Pt 3D 2.trk",
367  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
368 
369  tkcntHistosTkPt3D[2] = new TrackIPHistograms<double>
370  ("tkPt3_3D" + theExtensionString, "Track Pt 3D 3.trk",
371  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
372 
373  tkcntHistosTkPt3D[3] = new TrackIPHistograms<double>
374  ("tkPt4_3D" + theExtensionString, "Track Pt 3D 4.trk",
375  50, -0.1, 50.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
376 
377  // track nHits
378  tkcntHistosTkNHits2D[4] = new TrackIPHistograms<int>
379  ("tkNHits_2D" + theExtensionString, "Track NHits 2D",
380  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
381 
382  tkcntHistosTkNHits2D[0] = new TrackIPHistograms<int>
383  ("tkNHits1_2D" + theExtensionString, "Track NHits 2D 1.trk",
384  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
385 
386  tkcntHistosTkNHits2D[1] = new TrackIPHistograms<int>
387  ("tkNHits2_2D" + theExtensionString, "Track NHits 2D 2.trk",
388  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
389 
390  tkcntHistosTkNHits2D[2] = new TrackIPHistograms<int>
391  ("tkNHits3_2D" + theExtensionString, "Track NHits 2D 3.trk",
392  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
393 
394  tkcntHistosTkNHits2D[3] = new TrackIPHistograms<int>
395  ("tkNHits4_2D" + theExtensionString, "Track NHits 2D 4.trk",
396  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
397 
398  tkcntHistosTkNHits3D[4] = new TrackIPHistograms<int>
399  ("tkNHits_3D" + theExtensionString, "Track NHits 3D",
400  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
401 
402  tkcntHistosTkNHits3D[0] = new TrackIPHistograms<int>
403  ("tkNHits1_3D" + theExtensionString, "Track NHits 3D 1.trk",
404  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
405 
406  tkcntHistosTkNHits3D[1] = new TrackIPHistograms<int>
407  ("tkNHits2_3D" + theExtensionString, "Track NHits 3D 2.trk",
408  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
409 
410  tkcntHistosTkNHits3D[2] = new TrackIPHistograms<int>
411  ("tkNHits3_3D" + theExtensionString, "Track NHits 3D 3.trk",
412  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
413 
414  tkcntHistosTkNHits3D[3] = new TrackIPHistograms<int>
415  ("tkNHits4_3D" + theExtensionString, "Track NHits 3D 4.trk",
416  31, -0.5, 30.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
417 
418  //Pixel hits
419  tkcntHistosTkNPixelHits2D[4] = new TrackIPHistograms<int>
420  ("tkNPixelHits_2D" + theExtensionString, "Track NPixelHits 2D",
421  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
422 
423  tkcntHistosTkNPixelHits2D[0] = new TrackIPHistograms<int>
424  ("tkNPixelHits1_2D" + theExtensionString, "Track NPixelHits 2D 1.trk",
425  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
426 
427  tkcntHistosTkNPixelHits2D[1] = new TrackIPHistograms<int>
428  ("tkNPixelHits2_2D" + theExtensionString, "Track NPixelHits 2D 2.trk",
429  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
430 
431  tkcntHistosTkNPixelHits2D[2] = new TrackIPHistograms<int>
432  ("tkNPixelHits3_2D" + theExtensionString, "Track NPixelHits 2D 3.trk",
433  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
434 
435  tkcntHistosTkNPixelHits2D[3] = new TrackIPHistograms<int>
436  ("tkNPixelHits4_2D" + theExtensionString, "Track NPixelHits 2D 4.trk",
437  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
438 
439  tkcntHistosTkNPixelHits3D[4] = new TrackIPHistograms<int>
440  ("tkNPixelHits_3D" + theExtensionString, "Track NPixelHits 3D",
441  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
442 
443  tkcntHistosTkNPixelHits3D[0] = new TrackIPHistograms<int>
444  ("tkNPixelHits1_3D" + theExtensionString, "Track NPixelHits 3D 1.trk",
445  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
446 
447  tkcntHistosTkNPixelHits3D[1] = new TrackIPHistograms<int>
448  ("tkNPixelHits2_3D" + theExtensionString, "Track NPixelHits 3D 2.trk",
449  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
450 
451  tkcntHistosTkNPixelHits3D[2] = new TrackIPHistograms<int>
452  ("tkNPixelHits3_3D" + theExtensionString, "Track NPixelHits 3D 3.trk",
453  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
454 
455  tkcntHistosTkNPixelHits3D[3] = new TrackIPHistograms<int>
456  ("tkNPixelHits4_3D" + theExtensionString, "Track NPixelHits 3D 4.trk",
457  11, -0.5, 10.5, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
458 
459  // probability
460  tkcntHistosProb3D[4] = new TrackIPHistograms<float>
461  ("prob_3D" + theExtensionString, "3D IP probability",
462  50, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
463 
464  tkcntHistosProb3D[0] = new TrackIPHistograms<float>
465  ("prob1_3D" + theExtensionString, "3D IP probability 1.trk",
466  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
467 
468  tkcntHistosProb3D[1] = new TrackIPHistograms<float>
469  ("prob2_3D" + theExtensionString, "3D IP probability 2.trk",
470  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
471 
472  tkcntHistosProb3D[2] = new TrackIPHistograms<float>
473  ("prob3_3D" + theExtensionString, "3D IP probability 3.trk",
474  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
475 
476  tkcntHistosProb3D[3] = new TrackIPHistograms<float>
477  ("prob4_3D" + theExtensionString, "3D IP probability 4.trk",
478  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
479 
480  tkcntHistosProb2D[4] = new TrackIPHistograms<float>
481  ("prob_2D" + theExtensionString, "2D IP probability",
482  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
483 
484  tkcntHistosProb2D[0] = new TrackIPHistograms<float>
485  ("prob1_2D" + theExtensionString, "2D IP probability 1.trk",
486  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
487 
488  tkcntHistosProb2D[1] = new TrackIPHistograms<float>
489  ("prob2_2D" + theExtensionString, "2D IP probability 2.trk",
490  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
491 
492  tkcntHistosProb2D[2] = new TrackIPHistograms<float>
493  ("prob3_2D" + theExtensionString, "2D IP probability 3.trk",
494  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
495 
496  tkcntHistosProb2D[3] = new TrackIPHistograms<float>
497  ("prob4_2D" + theExtensionString, "2D IP probability 4.trk",
498  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
499 
500  //probability for tracks with IP value < 0 or IP value > 0
501  tkcntHistosTkProbIPneg2D = new TrackIPHistograms<float>
502  ("probIPneg_2D" + theExtensionString, "2D negative IP probability",
503  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
504  tkcntHistosTkProbIPpos2D = new TrackIPHistograms<float>
505  ("probIPpos_2D" + theExtensionString, "2D positive IP probability",
506  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
507  tkcntHistosTkProbIPneg3D = new TrackIPHistograms<float>
508  ("probIPneg_3D" + theExtensionString, "3D negative IP probability",
509  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
510  tkcntHistosTkProbIPpos3D = new TrackIPHistograms<float>
511  ("probIPpos_3D" + theExtensionString, "3D positive IP probability",
512  52, -1.04, 1.04, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
513 
514  //ghost Tracks and others
515  ghostTrackDistanceValuHisto = new TrackIPHistograms<double>
516  ("ghostTrackDist" + theExtensionString, "GhostTrackDistance",
517  50, 0.0, 0.1, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
518  ghostTrackDistanceSignHisto = new TrackIPHistograms<double>
519  ("ghostTrackDistSign" + theExtensionString, "GhostTrackDistance significance",
520  50, -5.0, 15.0, false, true, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
521  ghostTrackWeightHisto = new TrackIPHistograms<double>
522  ("ghostTrackWeight" + theExtensionString, "GhostTrack fit participation weight",
523  50, 0.0, 1.0, false, false, true, "b", trackIPDir, mc, makeQualityPlots_, ibook_) ;
524 
525  trackQualHisto = new FlavourHistograms<int>
526  ("trackQual" + theExtensionString, "Track Quality of Tracks Associated to Jets",
527  4, -1.5, 2.5, false, true, true, "b", trackIPDir, mc, ibook_);
528 
529  selectedTrackQualHisto = new FlavourHistograms<int>
530  ("selectedTrackQual" + theExtensionString, "Track Quality of Selected Tracks Associated to Jets",
531  4, -1.5, 2.5, false, true, true, "b", trackIPDir, mc, ibook_);
532 
533  trackMultVsJetPtHisto = new FlavourHistograms2D<double, int>
534  ("trackMultVsJetPt" + theExtensionString, "Track Multiplicity vs Jet Pt for Tracks Associated to Jets",
535  50, 0.0, 250.0, 21, -0.5, 30.5, false, trackIPDir, mc, true, ibook_);
536 
537  selectedTrackMultVsJetPtHisto = new FlavourHistograms2D<double, int>
538  ("selectedTrackMultVsJetPt" + theExtensionString, "Track Multiplicity vs Jet Pt for Selected Tracks Associated to Jets",
539  50, 0.0, 250.0, 21, -0.5, 20.5, false, trackIPDir, mc, true, ibook_);
540 
541 }
IPTagPlotter::tkcntHistosTkNChiSqr2D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr2D[5]
Definition: IPTagPlotter.h:53
IPTagPlotter::lowerIPEBound
double lowerIPEBound
Definition: IPTagPlotter.h:62
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
IPTagPlotter::tkcntHistosDecayLengthVal3D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal3D[5]
Definition: IPTagPlotter.h:48
IPTagPlotter::ghostTrackDistanceValuHisto
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosVal3D
TrackIPHistograms< double > * tkcntHistosVal3D[5]
Definition: IPTagPlotter.h:45
IPTagPlotter::lowerIPBound
double lowerIPBound
Definition: IPTagPlotter.h:62
IPTagPlotter::tkcntHistosTkNHits3D
TrackIPHistograms< int > * tkcntHistosTkNHits3D[5]
Definition: IPTagPlotter.h:58
IPTagPlotter::lowerIPSBound
double lowerIPSBound
Definition: IPTagPlotter.h:62
IPTagPlotter::tkcntHistosTkProbIPneg2D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosProb3D
TrackIPHistograms< float > * tkcntHistosProb3D[5]
Definition: IPTagPlotter.h:68
IPTagPlotter::ghostTrackWeightHisto
TrackIPHistograms< double > * ghostTrackWeightHisto
Definition: IPTagPlotter.h:72
IPTagPlotter::tkcntHistosTkProbIPneg3D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg3D
Definition: IPTagPlotter.h:71
IPTagPlotter::tkcntHistosTkPt2D
TrackIPHistograms< double > * tkcntHistosTkPt2D[5]
Definition: IPTagPlotter.h:55
BaseBTagPlotter::theExtensionString
const std::string theExtensionString
Definition: BaseBTagPlotter.h:32
IPTagPlotter::tkcntHistosDecayLengthVal2D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal2D[5]
Definition: IPTagPlotter.h:47
IPTagPlotter::mcPlots_
unsigned int mcPlots_
Definition: IPTagPlotter.h:37
IPTagPlotter::startEffPur_
double startEffPur_
Definition: IPTagPlotter.h:35
IPTagPlotter::tkcntHistosTkNHits2D
TrackIPHistograms< int > * tkcntHistosTkNHits2D[5]
Definition: IPTagPlotter.h:57
IPTagPlotter::trackQualHisto
FlavourHistograms< int > * trackQualHisto
Definition: IPTagPlotter.h:75
IPTagPlotter::endEffPur_
double endEffPur_
Definition: IPTagPlotter.h:36
IPTagPlotter::tkcntHistosVal2D
TrackIPHistograms< double > * tkcntHistosVal2D[5]
Definition: IPTagPlotter.h:46
IPTagPlotter::tkcntHistosJetDistVal3D
TrackIPHistograms< double > * tkcntHistosJetDistVal3D[5]
Definition: IPTagPlotter.h:50
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
IPTagPlotter::upperIPSBound
double upperIPSBound
Definition: IPTagPlotter.h:62
IPTagPlotter::trkNbr2D
FlavourHistograms< int > * trkNbr2D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosTkProbIPpos3D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos3D
Definition: IPTagPlotter.h:71
IPTagPlotter::selectedTrackMultVsJetPtHisto
FlavourHistograms2D< double, int > * selectedTrackMultVsJetPtHisto
Definition: IPTagPlotter.h:78
IPTagPlotter::maxDecayLength
double maxDecayLength
Definition: IPTagPlotter.h:64
IPTagPlotter::ghostTrackDistanceSignHisto
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::selectedTrackQualHisto
FlavourHistograms< int > * selectedTrackQualHisto
Definition: IPTagPlotter.h:76
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::trackMultVsJetPtHisto
FlavourHistograms2D< double, int > * trackMultVsJetPtHisto
Definition: IPTagPlotter.h:77
IPTagPlotter::trkNbr3D
FlavourHistograms< int > * trkNbr3D
Definition: IPTagPlotter.h:61
IPTagPlotter::maxJetDistance
double maxJetDistance
Definition: IPTagPlotter.h:64
IPTagPlotter::minDecayLength
double minDecayLength
Definition: IPTagPlotter.h:64
IPTagPlotter::tkcntHistosErr3D
TrackIPHistograms< double > * tkcntHistosErr3D[5]
Definition: IPTagPlotter.h:43
IPTagPlotter::upperIPEBound
double upperIPEBound
Definition: IPTagPlotter.h:62
BaseTagInfoPlotter::BaseTagInfoPlotter
BaseTagInfoPlotter(const std::string &tagName, const EtaPtBin &etaPtBin)
Definition: BaseTagInfoPlotter.h:15
IPTagPlotter::minJetDistance
double minJetDistance
Definition: IPTagPlotter.h:64
IPTagPlotter::tkcntHistosTkPt3D
TrackIPHistograms< double > * tkcntHistosTkPt3D[5]
Definition: IPTagPlotter.h:56
IPTagPlotter::upperIPBound
double upperIPBound
Definition: IPTagPlotter.h:62
IPTagPlotter::tkcntHistosProb2D
TrackIPHistograms< float > * tkcntHistosProb2D[5]
Definition: IPTagPlotter.h:69
IPTagPlotter::tkcntHistosTkNPixelHits3D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits3D[5]
Definition: IPTagPlotter.h:60
IPTagPlotter::tkcntHistosTkProbIPpos2D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::tkcntHistosTkNChiSqr3D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr3D[5]
Definition: IPTagPlotter.h:54
IPTagPlotter::willFinalize_
bool willFinalize_
Definition: IPTagPlotter.h:38
IPTagPlotter::tkcntHistosTkNPixelHits2D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits2D[5]
Definition: IPTagPlotter.h:59
IPTagPlotter::tkcntHistosErr2D
TrackIPHistograms< double > * tkcntHistosErr2D[5]
Definition: IPTagPlotter.h:44
IPTagPlotter::makeQualityPlots_
bool makeQualityPlots_
Definition: IPTagPlotter.h:39
IPTagPlotter::tkcntHistosJetDistVal2D
TrackIPHistograms< double > * tkcntHistosJetDistVal2D[5]
Definition: IPTagPlotter.h:49
template<class Container , class Base >
IPTagPlotter< Container, Base >::~IPTagPlotter ( )

Definition at line 544 of file IPTagPlotter_cc.h.

References gen::n.

545 {
546  if (willFinalize_) {
547  for(int n=1; n != 3; ++n) {
548  delete tkcntHistosSig3D[n];
549  delete effPurFromHistos[n-1];
550  }
551  for(int n=1; n != 3; ++n) {
552  delete tkcntHistosSig2D[n];
553  delete effPurFromHistos[n+1];
554  }
555  return;
556  }
557 
558  delete trkNbr3D;
559  delete trkNbr2D;
560  delete tkcntHistosTkProbIPneg2D;
561  delete tkcntHistosTkProbIPpos2D;
562  delete tkcntHistosTkProbIPneg3D;
563  delete tkcntHistosTkProbIPpos3D;
564  delete ghostTrackDistanceValuHisto;
565  delete ghostTrackDistanceSignHisto;
566  delete ghostTrackWeightHisto;
567  delete trackQualHisto;
568  delete selectedTrackQualHisto;
569  delete trackMultVsJetPtHisto;
570  delete selectedTrackMultVsJetPtHisto;
571 
572  for(int n=0; n != 5; ++n) {
573  delete tkcntHistosSig2D[n];
574  delete tkcntHistosSig3D[n];
575  delete tkcntHistosVal2D[n];
576  delete tkcntHistosVal3D[n];
577  delete tkcntHistosErr2D[n];
578  delete tkcntHistosErr3D[n];
579  delete tkcntHistosDecayLengthVal2D[n];
580  delete tkcntHistosDecayLengthVal3D[n];
581  delete tkcntHistosJetDistVal2D[n];
582  delete tkcntHistosJetDistVal3D[n];
583  //delete tkcntHistosJetDistSign2D[n];
584  //delete tkcntHistosJetDistSign3D[n];
585  delete tkcntHistosTkNChiSqr2D[n];
586  delete tkcntHistosTkNChiSqr3D[n];
587  delete tkcntHistosTkPt2D[n];
588  delete tkcntHistosTkPt3D[n];
589  delete tkcntHistosTkNHits2D[n];
590  delete tkcntHistosTkNHits3D[n];
591  delete tkcntHistosTkNPixelHits2D[n];
592  delete tkcntHistosTkNPixelHits3D[n];
593  delete tkcntHistosProb2D[n];
594  delete tkcntHistosProb3D[n];
595  }
596 }
IPTagPlotter::tkcntHistosTkNChiSqr2D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr2D[5]
Definition: IPTagPlotter.h:53
IPTagPlotter::tkcntHistosDecayLengthVal3D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal3D[5]
Definition: IPTagPlotter.h:48
IPTagPlotter::ghostTrackDistanceValuHisto
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosVal3D
TrackIPHistograms< double > * tkcntHistosVal3D[5]
Definition: IPTagPlotter.h:45
IPTagPlotter::tkcntHistosTkNHits3D
TrackIPHistograms< int > * tkcntHistosTkNHits3D[5]
Definition: IPTagPlotter.h:58
IPTagPlotter::tkcntHistosTkProbIPneg2D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosProb3D
TrackIPHistograms< float > * tkcntHistosProb3D[5]
Definition: IPTagPlotter.h:68
IPTagPlotter::ghostTrackWeightHisto
TrackIPHistograms< double > * ghostTrackWeightHisto
Definition: IPTagPlotter.h:72
IPTagPlotter::tkcntHistosTkProbIPneg3D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg3D
Definition: IPTagPlotter.h:71
IPTagPlotter::tkcntHistosTkPt2D
TrackIPHistograms< double > * tkcntHistosTkPt2D[5]
Definition: IPTagPlotter.h:55
IPTagPlotter::tkcntHistosDecayLengthVal2D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal2D[5]
Definition: IPTagPlotter.h:47
IPTagPlotter::effPurFromHistos
EffPurFromHistos * effPurFromHistos[4]
Definition: IPTagPlotter.h:66
IPTagPlotter::tkcntHistosTkNHits2D
TrackIPHistograms< int > * tkcntHistosTkNHits2D[5]
Definition: IPTagPlotter.h:57
IPTagPlotter::trackQualHisto
FlavourHistograms< int > * trackQualHisto
Definition: IPTagPlotter.h:75
IPTagPlotter::tkcntHistosVal2D
TrackIPHistograms< double > * tkcntHistosVal2D[5]
Definition: IPTagPlotter.h:46
IPTagPlotter::tkcntHistosJetDistVal3D
TrackIPHistograms< double > * tkcntHistosJetDistVal3D[5]
Definition: IPTagPlotter.h:50
IPTagPlotter::trkNbr2D
FlavourHistograms< int > * trkNbr2D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosTkProbIPpos3D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos3D
Definition: IPTagPlotter.h:71
IPTagPlotter::selectedTrackMultVsJetPtHisto
FlavourHistograms2D< double, int > * selectedTrackMultVsJetPtHisto
Definition: IPTagPlotter.h:78
IPTagPlotter::ghostTrackDistanceSignHisto
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::selectedTrackQualHisto
FlavourHistograms< int > * selectedTrackQualHisto
Definition: IPTagPlotter.h:76
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::trackMultVsJetPtHisto
FlavourHistograms2D< double, int > * trackMultVsJetPtHisto
Definition: IPTagPlotter.h:77
IPTagPlotter::trkNbr3D
FlavourHistograms< int > * trkNbr3D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosErr3D
TrackIPHistograms< double > * tkcntHistosErr3D[5]
Definition: IPTagPlotter.h:43
IPTagPlotter::tkcntHistosTkPt3D
TrackIPHistograms< double > * tkcntHistosTkPt3D[5]
Definition: IPTagPlotter.h:56
IPTagPlotter::tkcntHistosProb2D
TrackIPHistograms< float > * tkcntHistosProb2D[5]
Definition: IPTagPlotter.h:69
IPTagPlotter::tkcntHistosTkNPixelHits3D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits3D[5]
Definition: IPTagPlotter.h:60
IPTagPlotter::tkcntHistosTkProbIPpos2D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::tkcntHistosTkNChiSqr3D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr3D[5]
Definition: IPTagPlotter.h:54
IPTagPlotter::willFinalize_
bool willFinalize_
Definition: IPTagPlotter.h:38
IPTagPlotter::tkcntHistosTkNPixelHits2D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits2D[5]
Definition: IPTagPlotter.h:59
IPTagPlotter::tkcntHistosErr2D
TrackIPHistograms< double > * tkcntHistosErr2D[5]
Definition: IPTagPlotter.h:44
IPTagPlotter::tkcntHistosJetDistVal2D
TrackIPHistograms< double > * tkcntHistosJetDistVal2D[5]
Definition: IPTagPlotter.h:49

Member Function Documentation

template<class Container , class Base >
void IPTagPlotter< Container, Base >::analyzeTag ( const reco::BaseTagInfo *  baseTagInfo,
const double &  jec,
const int &  jetFlavour 
)
virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 599 of file IPTagPlotter_cc.h.

602 {
603  analyzeTag(baseTagInfo, jec, jetFlavour, 1.);
604 }
IPTagPlotter::analyzeTag
void analyzeTag(const reco::BaseTagInfo *baseTagInfo, const double &jec, const int &jetFlavour)
Definition: IPTagPlotter_cc.h:599
patTestJEC_cfi.jec
tuple jec
Definition: patTestJEC_cfi.py:5
template<class Container , class Base >
void IPTagPlotter< Container, Base >::analyzeTag ( const reco::BaseTagInfo *  baseTagInfo,
const double &  jec,
const int &  jetFlavour,
const float &  w 
)
virtual

Reimplemented from BaseTagInfoPlotter.

Definition at line 607 of file IPTagPlotter_cc.h.

References listHistos::decayLength, Exception, reco::Track::found(), reco::IPTagInfo< Container, Base >::hasProbabilities(), reco::TrackBase::hitPattern(), reco::IPTagInfo< Container, Base >::impactParameterData(), reco::btag::IP2DSig, reco::btag::IP3DSig, patTestJEC_cfi::jec, mag(), gen::n, reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidPixelHits(), reco::IPTagInfo< Container, Base >::primaryVertex(), reco::btag::Prob2D, reco::btag::Prob3D, reco::IPTagInfo< Container, Base >::probabilities(), reco::TrackBase::pt(), MetAnalyzer::pv(), TrackCollections2monitor_cff::selectedTracks, reco::IPTagInfo< Container, Base >::sortedIndexes(), reco::IPTagInfo< Container, Base >::sortedTracks(), dumpRecoGeometry_cfg::tagInfo, reco::btag::toTrack(), reco::TrackBase::undefQuality, and w.

610 {
611  // const reco::TrackIPTagInfo * tagInfo =
612  // dynamic_cast<const reco::TrackIPTagInfo *>(baseTagInfo);
613  const reco::IPTagInfo<Container, Base> * tagInfo =
614  dynamic_cast<const reco::IPTagInfo<Container, Base> *>(baseTagInfo);
615 
616  if (!tagInfo) {
617  throw cms::Exception("Configuration")
618  << "BTagPerformanceAnalyzer: Extended TagInfo not of type TrackIPTagInfo. " << std::endl;
619  }
620 
621  const GlobalPoint pv(tagInfo->primaryVertex()->position().x(),
622  tagInfo->primaryVertex()->position().y(),
623  tagInfo->primaryVertex()->position().z());
624 
625  const std::vector<reco::btag::TrackIPData>& ip = tagInfo->impactParameterData();
626 
627  std::vector<float> prob2d, prob3d;
628  if (tagInfo->hasProbabilities()) {
629  prob2d = tagInfo->probabilities(0);
630  prob3d = tagInfo->probabilities(1);
631  }
632 
633  std::vector<std::size_t> sortedIndices = tagInfo->sortedIndexes(reco::btag::IP2DSig);
634  std::vector<std::size_t> selectedIndices;
635  Container sortedTracks = tagInfo->sortedTracks(sortedIndices);
636  Container selectedTracks;
637  for(unsigned int n = 0; n != sortedIndices.size(); ++n) {
638  double decayLength = (ip[sortedIndices[n]].closestToJetAxis - pv).mag();
639  double jetDistance = ip[sortedIndices[n]].distanceToJetAxis.value();
640  if(decayLength > minDecayLength && decayLength < maxDecayLength &&
641  fabs(jetDistance) >= minJetDistance && fabs(jetDistance) < maxJetDistance ) {
642  selectedIndices.push_back(sortedIndices[n]);
643  selectedTracks.push_back(sortedTracks[n]);
644  }
645  }
646 
647  trkNbr2D->fill(jetFlavour, selectedIndices.size(),w);
648 
649  for(unsigned int n=0; n != selectedIndices.size(); ++n) {
650  const reco::Track * track = reco::btag::toTrack(selectedTracks[n]);
651  const reco::TrackBase::TrackQuality& trackQual = highestTrackQual(track);
652  tkcntHistosSig2D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.significance(), true,w);
653  tkcntHistosVal2D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.value(), true,w);
654  tkcntHistosErr2D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.error(), true,w);
655  const double& decayLen = (ip[selectedIndices[n]].closestToJetAxis - pv).mag();
656  tkcntHistosDecayLengthVal2D[4]->fill(jetFlavour, trackQual, decayLen, true,w);
657  tkcntHistosJetDistVal2D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.value(), true,w);
658  //tkcntHistosJetDistSign2D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.significance(), true,w);
659  tkcntHistosTkNChiSqr2D[4]->fill(jetFlavour, trackQual, track->normalizedChi2(), true,w);
660  tkcntHistosTkPt2D[4]->fill(jetFlavour, trackQual, track->pt(), true,w);
661  tkcntHistosTkNHits2D[4]->fill(jetFlavour, trackQual, track->found(), true,w);
662  tkcntHistosTkNPixelHits2D[4]->fill(jetFlavour, trackQual, track->hitPattern().numberOfValidPixelHits(), true,w);
663  if(n >= 4) continue;
664  tkcntHistosSig2D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.significance(), true,w);
665  tkcntHistosVal2D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.value(), true,w);
666  tkcntHistosErr2D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip2d.error(), true,w);
667  tkcntHistosDecayLengthVal2D[n]->fill(jetFlavour, trackQual, decayLen, true,w);
668  tkcntHistosJetDistVal2D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.value(), true,w);
669  //tkcntHistosJetDistSign2D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.significance(), true,w);
670  tkcntHistosTkNChiSqr2D[n]->fill(jetFlavour, trackQual, track->normalizedChi2(), true,w);
671  tkcntHistosTkPt2D[n]->fill(jetFlavour, trackQual, track->pt(), true,w);
672  tkcntHistosTkNHits2D[n]->fill(jetFlavour, trackQual, track->found(), true,w);
673  tkcntHistosTkNPixelHits2D[n]->fill(jetFlavour, trackQual, track->hitPattern().numberOfValidPixelHits(), true,w);
674  }
675  sortedIndices = tagInfo->sortedIndexes(reco::btag::Prob2D);
676  selectedIndices.clear();
677  sortedTracks = tagInfo->sortedTracks(sortedIndices);
678  selectedTracks.clear();
679  for(unsigned int n = 0; n != sortedIndices.size(); ++n) {
680  double decayLength = (ip[sortedIndices[n]].closestToJetAxis - pv).mag();
681  double jetDistance = ip[sortedIndices[n]].distanceToJetAxis.value();
682  if(decayLength > minDecayLength && decayLength < maxDecayLength &&
683  fabs(jetDistance) >= minJetDistance && fabs(jetDistance) < maxJetDistance ) {
684  selectedIndices.push_back(sortedIndices[n]);
685  selectedTracks.push_back(sortedTracks[n]);
686  }
687  }
688  for(unsigned int n=0; n != selectedIndices.size(); ++n) {
689  const reco::Track * track = reco::btag::toTrack(selectedTracks[n]);
690  const reco::TrackBase::TrackQuality& trackQual = highestTrackQual(track);
691  tkcntHistosProb2D[4]->fill(jetFlavour, trackQual, prob2d[selectedIndices[n]], true,w);
692  if(ip[selectedIndices[n]].ip2d.value() < 0) tkcntHistosTkProbIPneg2D->fill(jetFlavour, trackQual, prob2d[selectedIndices[n]], true,w);
693  else tkcntHistosTkProbIPpos2D->fill(jetFlavour, trackQual, prob2d[selectedIndices[n]], true,w);
694  if(n >= 4) continue;
695  tkcntHistosProb2D[n]->fill(jetFlavour, trackQual, prob2d[selectedIndices[n]], true,w);
696  }
697  for(unsigned int n=selectedIndices.size(); n < 4; ++n){
698  const reco::TrackBase::TrackQuality trackQual = reco::TrackBase::undefQuality;
699  tkcntHistosSig2D[n]->fill(jetFlavour, trackQual, lowerIPSBound-1.0, false,w);
700  tkcntHistosVal2D[n]->fill(jetFlavour, trackQual, lowerIPBound-1.0, false,w);
701  tkcntHistosErr2D[n]->fill(jetFlavour, trackQual, lowerIPEBound-1.0, false,w);
702  }
703  sortedIndices = tagInfo->sortedIndexes(reco::btag::IP3DSig);
704  selectedIndices.clear();
705  sortedTracks = tagInfo->sortedTracks(sortedIndices);
706  selectedTracks.clear();
707  for(unsigned int n = 0; n != sortedIndices.size(); ++n) {
708  double decayLength = (ip[sortedIndices[n]].closestToJetAxis - pv).mag();
709  double jetDistance = ip[sortedIndices[n]].distanceToJetAxis.value();
710  if(decayLength > minDecayLength && decayLength < maxDecayLength &&
711  fabs(jetDistance) >= minJetDistance && fabs(jetDistance) < maxJetDistance ) {
712  selectedIndices.push_back(sortedIndices[n]);
713  selectedTracks.push_back(sortedTracks[n]);
714  }
715  }
716 
717  trkNbr3D->fill(jetFlavour, selectedIndices.size(),w);
718  int nSelectedTracks = selectedIndices.size();
719 
720  for(unsigned int n=0; n != selectedIndices.size(); ++n) {
721  const reco::Track * track = reco::btag::toTrack(selectedTracks[n]);
722  const reco::TrackBase::TrackQuality& trackQual = highestTrackQual(track);
723  tkcntHistosSig3D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.significance(), true,w);
724  tkcntHistosVal3D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.value(), true,w);
725  tkcntHistosErr3D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.error(), true,w);
726  const double& decayLen = (ip[selectedIndices[n]].closestToJetAxis - pv).mag();
727  tkcntHistosDecayLengthVal3D[4]->fill(jetFlavour, trackQual, decayLen, true,w);
728  tkcntHistosJetDistVal3D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.value(), true,w);
729  //tkcntHistosJetDistSign3D[4]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.significance(), true,w);
730  tkcntHistosTkNChiSqr3D[4]->fill(jetFlavour, trackQual, track->normalizedChi2(), true,w);
731  tkcntHistosTkPt3D[4]->fill(jetFlavour, trackQual, track->pt(), true,w);
732  tkcntHistosTkNHits3D[4]->fill(jetFlavour, trackQual, track->found(), true,w);
733  tkcntHistosTkNPixelHits3D[4]->fill(jetFlavour, trackQual, track->hitPattern().numberOfValidPixelHits(), true,w);
734  //ghostTrack infos
735  ghostTrackDistanceValuHisto->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToGhostTrack.value(), true,w);
736  ghostTrackDistanceSignHisto->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToGhostTrack.significance(), true,w);
737  ghostTrackWeightHisto->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ghostTrackWeight, true,w);
738  selectedTrackQualHisto->fill(jetFlavour, trackQual,w);
739  if(n >= 4) continue;
740  tkcntHistosSig3D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.significance(), true,w);
741  tkcntHistosVal3D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.value(), true,w);
742  tkcntHistosErr3D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].ip3d.error(), true,w);
743  tkcntHistosDecayLengthVal3D[n]->fill(jetFlavour, trackQual, decayLen, true,w);
744  tkcntHistosJetDistVal3D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.value(), true,w);
745  //tkcntHistosJetDistSign3D[n]->fill(jetFlavour, trackQual, ip[selectedIndices[n]].distanceToJetAxis.significance(), true,w);
746  tkcntHistosTkNChiSqr3D[n]->fill(jetFlavour, trackQual, track->normalizedChi2(), true,w);
747  tkcntHistosTkPt3D[n]->fill(jetFlavour, trackQual, track->pt(), true,w);
748  tkcntHistosTkNHits3D[n]->fill(jetFlavour, trackQual, track->found(), true,w);
749  tkcntHistosTkNPixelHits3D[n]->fill(jetFlavour, trackQual, track->hitPattern().numberOfValidPixelHits(), true,w);
750  }
751  sortedIndices = tagInfo->sortedIndexes(reco::btag::Prob3D);
752  selectedIndices.clear();
753  sortedTracks = tagInfo->sortedTracks(sortedIndices);
754  selectedTracks.clear();
755  for(unsigned int n = 0; n != sortedIndices.size(); ++n) {
756  double decayLength = (ip[sortedIndices[n]].closestToJetAxis - pv).mag();
757  double jetDistance = ip[sortedIndices[n]].distanceToJetAxis.value();
758  if(decayLength > minDecayLength && decayLength < maxDecayLength &&
759  fabs(jetDistance) >= minJetDistance && fabs(jetDistance) < maxJetDistance ) {
760  selectedIndices.push_back(sortedIndices[n]);
761  selectedTracks.push_back(sortedTracks[n]);
762  }
763  }
764  for(unsigned int n=0; n != selectedIndices.size(); ++n) {
765  const reco::Track * track = reco::btag::toTrack(selectedTracks[n]);
766  const reco::TrackBase::TrackQuality& trackQual = highestTrackQual(track);
767  tkcntHistosProb3D[4]->fill(jetFlavour, trackQual, prob3d[selectedIndices[n]], true,w);
768  if(ip[selectedIndices[n]].ip3d.value() < 0) tkcntHistosTkProbIPneg3D->fill(jetFlavour, trackQual, prob3d[selectedIndices[n]], true,w);
769  else tkcntHistosTkProbIPpos3D->fill(jetFlavour, trackQual, prob3d[selectedIndices[n]], true,w);
770  if(n >= 4) continue;
771  tkcntHistosProb3D[n]->fill(jetFlavour, trackQual, prob3d[selectedIndices[n]], true,w);
772  }
773  for(unsigned int n=selectedIndices.size(); n < 4; ++n){
774  const reco::TrackBase::TrackQuality trackQual = reco::TrackBase::undefQuality;
775  tkcntHistosSig3D[n]->fill(jetFlavour, trackQual, lowerIPSBound-1.0, false,w);
776  tkcntHistosVal3D[n]->fill(jetFlavour, trackQual, lowerIPBound-1.0, false,w);
777  tkcntHistosErr3D[n]->fill(jetFlavour, trackQual, lowerIPEBound-1.0, false,w);
778  }
779  for(unsigned int n = 0; n != sortedTracks.size(); ++n) {
780  trackQualHisto->fill(jetFlavour, highestTrackQual(reco::btag::toTrack(sortedTracks[n])),w);
781  }
782 
783  //still need to implement weights in FlavourHistograms2D
784  trackMultVsJetPtHisto->fill(jetFlavour, tagInfo->jet()->pt()*jec, sortedTracks.size());
785  selectedTrackMultVsJetPtHisto->fill(jetFlavour, tagInfo->jet()->pt()*jec, nSelectedTracks); //tagInfo->selectedTracks().size());
786 }
IPTagPlotter::tkcntHistosTkNChiSqr2D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr2D[5]
Definition: IPTagPlotter.h:53
IPTagPlotter::lowerIPEBound
double lowerIPEBound
Definition: IPTagPlotter.h:62
TrackCollections2monitor_cff.selectedTracks
list selectedTracks
Definition: TrackCollections2monitor_cff.py:12
reco::IPTagInfo::probabilities
const std::vector< float > & probabilities(int ip) const
Definition: IPTagInfo.h:103
listHistos.decayLength
tuple decayLength
Definition: listHistos.py:101
IPTagPlotter::tkcntHistosDecayLengthVal3D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal3D[5]
Definition: IPTagPlotter.h:48
IPTagPlotter::ghostTrackDistanceValuHisto
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosVal3D
TrackIPHistograms< double > * tkcntHistosVal3D[5]
Definition: IPTagPlotter.h:45
IPTagPlotter::lowerIPBound
double lowerIPBound
Definition: IPTagPlotter.h:62
IPTagPlotter::tkcntHistosTkNHits3D
TrackIPHistograms< int > * tkcntHistosTkNHits3D[5]
Definition: IPTagPlotter.h:58
cond::ecalcond::Container
EcalChannelStatus Container
Definition: EcalChannelStatusPyWrapper.cc:31
w
const double w
Definition: UKUtility.cc:23
IPTagPlotter::lowerIPSBound
double lowerIPSBound
Definition: IPTagPlotter.h:62
reco::TrackBase::normalizedChi2
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:548
IPTagPlotter::tkcntHistosTkProbIPneg2D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg2D
Definition: IPTagPlotter.h:70
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:149
mag
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Definition: Basic3DVectorLD.h:147
IPTagPlotter::tkcntHistosProb3D
TrackIPHistograms< float > * tkcntHistosProb3D[5]
Definition: IPTagPlotter.h:68
IPTagPlotter::ghostTrackWeightHisto
TrackIPHistograms< double > * ghostTrackWeightHisto
Definition: IPTagPlotter.h:72
IPTagPlotter::tkcntHistosTkProbIPneg3D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg3D
Definition: IPTagPlotter.h:71
IPTagPlotter::tkcntHistosTkPt2D
TrackIPHistograms< double > * tkcntHistosTkPt2D[5]
Definition: IPTagPlotter.h:55
IPTagPlotter::tkcntHistosDecayLengthVal2D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal2D[5]
Definition: IPTagPlotter.h:47
FlavourHistograms::fill
void fill(const int &flavour, const T &variable) const
Definition: FlavourHistorgrams.h:306
IPTagPlotter::tkcntHistosTkNHits2D
TrackIPHistograms< int > * tkcntHistosTkNHits2D[5]
Definition: IPTagPlotter.h:57
dumpRecoGeometry_cfg.tagInfo
tuple tagInfo
Definition: dumpRecoGeometry_cfg.py:156
IPTagPlotter::trackQualHisto
FlavourHistograms< int > * trackQualHisto
Definition: IPTagPlotter.h:75
IPTagPlotter::tkcntHistosVal2D
TrackIPHistograms< double > * tkcntHistosVal2D[5]
Definition: IPTagPlotter.h:46
IPTagPlotter::tkcntHistosJetDistVal3D
TrackIPHistograms< double > * tkcntHistosJetDistVal3D[5]
Definition: IPTagPlotter.h:50
patTestJEC_cfi.jec
tuple jec
Definition: patTestJEC_cfi.py:5
reco::btag::toTrack
const reco::Track * toTrack(const reco::TrackBaseRef &t)
Definition: IPTagInfo.h:24
reco::IPTagInfo::primaryVertex
const edm::Ref< VertexCollection > & primaryVertex() const
Definition: IPTagInfo.h:133
IPTagPlotter::trkNbr2D
FlavourHistograms< int > * trkNbr2D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosTkProbIPpos3D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos3D
Definition: IPTagPlotter.h:71
FlavourHistograms2D::fill
void fill(const int &flavour, const T &variableX, const G &variableY) const
Definition: FlavourHistorgrams2D.h:433
TrackIPHistograms::fill
void fill(const int &flavour, const reco::TrackBase::TrackQuality &quality, const T &variable, const bool &hasTrack) const
Definition: TrackIPHistograms.h:92
IPTagPlotter::selectedTrackMultVsJetPtHisto
FlavourHistograms2D< double, int > * selectedTrackMultVsJetPtHisto
Definition: IPTagPlotter.h:78
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:608
IPTagPlotter::maxDecayLength
double maxDecayLength
Definition: IPTagPlotter.h:64
IPTagPlotter::ghostTrackDistanceSignHisto
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::selectedTrackQualHisto
FlavourHistograms< int > * selectedTrackQualHisto
Definition: IPTagPlotter.h:76
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::trackMultVsJetPtHisto
FlavourHistograms2D< double, int > * trackMultVsJetPtHisto
Definition: IPTagPlotter.h:77
reco::IPTagInfo::sortedIndexes
std::vector< size_t > sortedIndexes(btag::SortCriteria mode=reco::btag::IP3DSig) const
Definition: IPTagInfo.h:238
IPTagPlotter::trkNbr3D
FlavourHistograms< int > * trkNbr3D
Definition: IPTagPlotter.h:61
IPTagPlotter::maxJetDistance
double maxJetDistance
Definition: IPTagPlotter.h:64
IPTagPlotter::minDecayLength
double minDecayLength
Definition: IPTagPlotter.h:64
reco::IPTagInfo::impactParameterData
const std::vector< btag::TrackIPData > & impactParameterData() const
Definition: IPTagInfo.h:91
IPTagPlotter::tkcntHistosErr3D
TrackIPHistograms< double > * tkcntHistosErr3D[5]
Definition: IPTagPlotter.h:43
IPTagPlotter::minJetDistance
double minJetDistance
Definition: IPTagPlotter.h:64
IPTagPlotter::highestTrackQual
reco::TrackBase::TrackQuality highestTrackQual(const reco::Track *track) const
Definition: IPTagPlotter_cc.h:966
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:437
IPTagPlotter::tkcntHistosTkPt3D
TrackIPHistograms< double > * tkcntHistosTkPt3D[5]
Definition: IPTagPlotter.h:56
IPTagPlotter::tkcntHistosProb2D
TrackIPHistograms< float > * tkcntHistosProb2D[5]
Definition: IPTagPlotter.h:69
IPTagPlotter::tkcntHistosTkNPixelHits3D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits3D[5]
Definition: IPTagPlotter.h:60
IPTagPlotter::tkcntHistosTkProbIPpos2D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos2D
Definition: IPTagPlotter.h:70
reco::Track::found
unsigned short found() const
Number of valid hits on track.
Definition: Track.h:194
reco::HitPattern::numberOfValidPixelHits
int numberOfValidPixelHits() const
Definition: HitPattern.h:816
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::tkcntHistosTkNChiSqr3D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr3D[5]
Definition: IPTagPlotter.h:54
reco::IPTagInfo::sortedTracks
Container sortedTracks(const std::vector< size_t > &indexes) const
Definition: IPTagInfo.h:129
IPTagPlotter::tkcntHistosTkNPixelHits2D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits2D[5]
Definition: IPTagPlotter.h:59
IPTagPlotter::tkcntHistosErr2D
TrackIPHistograms< double > * tkcntHistosErr2D[5]
Definition: IPTagPlotter.h:44
reco::IPTagInfo::hasProbabilities
virtual bool hasProbabilities() const
Definition: IPTagInfo.h:85
IPTagPlotter::tkcntHistosJetDistVal2D
TrackIPHistograms< double > * tkcntHistosJetDistVal2D[5]
Definition: IPTagPlotter.h:49
template<class Container , class Base >
void IPTagPlotter< Container, Base >::epsPlot ( const std::string &  name)
virtual

Implements BaseBTagPlotter.

Definition at line 923 of file IPTagPlotter_cc.h.

References gen::n.

924 {
925  if (willFinalize_) {
926  for(int n=0; n != 4; ++n) effPurFromHistos[n]->epsPlot(name);
927  return;
928  }
929  trkNbr2D->epsPlot(name);
930  trkNbr3D->epsPlot(name);
931  tkcntHistosTkProbIPneg2D->epsPlot(name);
932  tkcntHistosTkProbIPpos2D->epsPlot(name);
933  tkcntHistosTkProbIPneg3D->epsPlot(name);
934  tkcntHistosTkProbIPpos3D->epsPlot(name);
935  ghostTrackDistanceValuHisto->epsPlot(name);
936  ghostTrackDistanceSignHisto->epsPlot(name);
937  ghostTrackWeightHisto->epsPlot(name);
938  for(int n=0; n != 5; ++n) {
939  tkcntHistosSig2D[n]->epsPlot(name);
940  tkcntHistosSig3D[n]->epsPlot(name);
941  tkcntHistosVal2D[n]->epsPlot(name);
942  tkcntHistosVal3D[n]->epsPlot(name);
943  tkcntHistosErr2D[n]->epsPlot(name);
944  tkcntHistosErr3D[n]->epsPlot(name);
945  tkcntHistosProb2D[n]->epsPlot(name);
946  tkcntHistosProb3D[n]->epsPlot(name);
947  tkcntHistosDecayLengthVal2D[n]->epsPlot(name);
948  tkcntHistosDecayLengthVal3D[n]->epsPlot(name);
949  tkcntHistosJetDistVal2D[n]->epsPlot(name);
950  tkcntHistosJetDistVal3D[n]->epsPlot(name);
951  //tkcntHistosJetDistSign2D[n]->epsPlot(name);
952  //tkcntHistosJetDistSign3D[n]->epsPlot(name);
953  tkcntHistosTkNChiSqr2D[n]->epsPlot(name);
954  tkcntHistosTkNChiSqr3D[n]->epsPlot(name);
955  tkcntHistosTkPt2D[n]->epsPlot(name);
956  tkcntHistosTkPt3D[n]->epsPlot(name);
957  tkcntHistosTkNHits2D[n]->epsPlot(name);
958  tkcntHistosTkNHits3D[n]->epsPlot(name);
959  tkcntHistosTkNPixelHits2D[n]->epsPlot(name);
960  tkcntHistosTkNPixelHits3D[n]->epsPlot(name);
961 
962  }
963 }
IPTagPlotter::tkcntHistosTkNChiSqr2D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr2D[5]
Definition: IPTagPlotter.h:53
IPTagPlotter::tkcntHistosDecayLengthVal3D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal3D[5]
Definition: IPTagPlotter.h:48
IPTagPlotter::ghostTrackDistanceValuHisto
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosVal3D
TrackIPHistograms< double > * tkcntHistosVal3D[5]
Definition: IPTagPlotter.h:45
IPTagPlotter::tkcntHistosTkNHits3D
TrackIPHistograms< int > * tkcntHistosTkNHits3D[5]
Definition: IPTagPlotter.h:58
IPTagPlotter::tkcntHistosTkProbIPneg2D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosProb3D
TrackIPHistograms< float > * tkcntHistosProb3D[5]
Definition: IPTagPlotter.h:68
IPTagPlotter::ghostTrackWeightHisto
TrackIPHistograms< double > * ghostTrackWeightHisto
Definition: IPTagPlotter.h:72
IPTagPlotter::tkcntHistosTkProbIPneg3D
TrackIPHistograms< float > * tkcntHistosTkProbIPneg3D
Definition: IPTagPlotter.h:71
IPTagPlotter::tkcntHistosTkPt2D
TrackIPHistograms< double > * tkcntHistosTkPt2D[5]
Definition: IPTagPlotter.h:55
IPTagPlotter::tkcntHistosDecayLengthVal2D
TrackIPHistograms< double > * tkcntHistosDecayLengthVal2D[5]
Definition: IPTagPlotter.h:47
IPTagPlotter::effPurFromHistos
EffPurFromHistos * effPurFromHistos[4]
Definition: IPTagPlotter.h:66
IPTagPlotter::tkcntHistosTkNHits2D
TrackIPHistograms< int > * tkcntHistosTkNHits2D[5]
Definition: IPTagPlotter.h:57
IPTagPlotter::tkcntHistosVal2D
TrackIPHistograms< double > * tkcntHistosVal2D[5]
Definition: IPTagPlotter.h:46
IPTagPlotter::tkcntHistosJetDistVal3D
TrackIPHistograms< double > * tkcntHistosJetDistVal3D[5]
Definition: IPTagPlotter.h:50
IPTagPlotter::trkNbr2D
FlavourHistograms< int > * trkNbr2D
Definition: IPTagPlotter.h:61
FlavourHistograms::epsPlot
void epsPlot(const std::string &name)
Definition: FlavourHistorgrams.h:503
IPTagPlotter::tkcntHistosTkProbIPpos3D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos3D
Definition: IPTagPlotter.h:71
IPTagPlotter::ghostTrackDistanceSignHisto
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::trkNbr3D
FlavourHistograms< int > * trkNbr3D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosErr3D
TrackIPHistograms< double > * tkcntHistosErr3D[5]
Definition: IPTagPlotter.h:43
IPTagPlotter::tkcntHistosTkPt3D
TrackIPHistograms< double > * tkcntHistosTkPt3D[5]
Definition: IPTagPlotter.h:56
IPTagPlotter::tkcntHistosProb2D
TrackIPHistograms< float > * tkcntHistosProb2D[5]
Definition: IPTagPlotter.h:69
IPTagPlotter::tkcntHistosTkNPixelHits3D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits3D[5]
Definition: IPTagPlotter.h:60
IPTagPlotter::epsPlot
void epsPlot(const std::string &name)
Definition: IPTagPlotter_cc.h:923
IPTagPlotter::tkcntHistosTkProbIPpos2D
TrackIPHistograms< float > * tkcntHistosTkProbIPpos2D
Definition: IPTagPlotter.h:70
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::tkcntHistosTkNChiSqr3D
TrackIPHistograms< double > * tkcntHistosTkNChiSqr3D[5]
Definition: IPTagPlotter.h:54
IPTagPlotter::willFinalize_
bool willFinalize_
Definition: IPTagPlotter.h:38
IPTagPlotter::tkcntHistosTkNPixelHits2D
TrackIPHistograms< int > * tkcntHistosTkNPixelHits2D[5]
Definition: IPTagPlotter.h:59
IPTagPlotter::tkcntHistosErr2D
TrackIPHistograms< double > * tkcntHistosErr2D[5]
Definition: IPTagPlotter.h:44
IPTagPlotter::tkcntHistosJetDistVal2D
TrackIPHistograms< double > * tkcntHistosJetDistVal2D[5]
Definition: IPTagPlotter.h:49
template<class Container , class Base >
void IPTagPlotter< Container, Base >::finalize ( DQMStore::IBooker &  ibook_,
DQMStore::IGetter &  igetter_ 
)
virtual

Implements BaseBTagPlotter.

Definition at line 789 of file IPTagPlotter_cc.h.

References bookConverter::compute(), gen::n, and AlCaHLTBitMon_QueryRunRegistry::string.

790 {
791  //
792  // final processing:
793  // produce the misid. vs. eff histograms
794  //
795  const std::string trackIPDir(theExtensionString.substr(1));
796 
797  tkcntHistosSig3D[1] = new TrackIPHistograms<double>
798  ("ips2_3D" + theExtensionString, "3D IP significance 2.trk",
799  nBinsIPS, lowerIPSBound, upperIPSBound, "b", trackIPDir, mcPlots_, makeQualityPlots_, igetter_) ;
800  effPurFromHistos[0] = new EffPurFromHistos (tkcntHistosSig3D[1],trackIPDir, mcPlots_, ibook_,
801  nBinEffPur_, startEffPur_,
802  endEffPur_);
803 
804  tkcntHistosSig3D[2] = new TrackIPHistograms<double>
805  ("ips3_3D" + theExtensionString, "3D IP significance 3.trk",
806  nBinsIPS, lowerIPSBound, upperIPSBound, "b", trackIPDir, mcPlots_, makeQualityPlots_, igetter_) ;
807  effPurFromHistos[1] = new EffPurFromHistos (tkcntHistosSig3D[2],trackIPDir, mcPlots_, ibook_,
808  nBinEffPur_, startEffPur_,
809  endEffPur_);
810 
811  tkcntHistosSig2D[1] = new TrackIPHistograms<double>
812  ("ips2_2D" + theExtensionString, "2D IP significance 2.trk",
813  nBinsIPS, lowerIPSBound, upperIPSBound, "b", trackIPDir, mcPlots_, makeQualityPlots_, igetter_) ;
814  effPurFromHistos[2] = new EffPurFromHistos (tkcntHistosSig2D[1],trackIPDir, mcPlots_, ibook_,
815  nBinEffPur_, startEffPur_,
816  endEffPur_);
817 
818  tkcntHistosSig2D[2] = new TrackIPHistograms<double>
819  ("ips3_2D" + theExtensionString, "2D IP significance 3.trk",
820  nBinsIPS, lowerIPSBound, upperIPSBound, "b", trackIPDir, mcPlots_, makeQualityPlots_, igetter_) ;
821  effPurFromHistos[3] = new EffPurFromHistos (tkcntHistosSig2D[2],trackIPDir, mcPlots_, ibook_,
822  nBinEffPur_, startEffPur_,
823  endEffPur_);
824 
825  for(int n=0; n != 4; ++n) effPurFromHistos[n]->compute(ibook_);
826 }
IPTagPlotter::lowerIPSBound
double lowerIPSBound
Definition: IPTagPlotter.h:62
BaseBTagPlotter::theExtensionString
const std::string theExtensionString
Definition: BaseBTagPlotter.h:32
IPTagPlotter::mcPlots_
unsigned int mcPlots_
Definition: IPTagPlotter.h:37
IPTagPlotter::effPurFromHistos
EffPurFromHistos * effPurFromHistos[4]
Definition: IPTagPlotter.h:66
IPTagPlotter::startEffPur_
double startEffPur_
Definition: IPTagPlotter.h:35
IPTagPlotter::endEffPur_
double endEffPur_
Definition: IPTagPlotter.h:36
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
IPTagPlotter::upperIPSBound
double upperIPSBound
Definition: IPTagPlotter.h:62
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::makeQualityPlots_
bool makeQualityPlots_
Definition: IPTagPlotter.h:39
template<class Container , class Base >
reco::TrackBase::TrackQuality IPTagPlotter< Container, Base >::highestTrackQual ( const reco::Track *  track) const

Definition at line 966 of file IPTagPlotter_cc.h.

References reco::TrackBase::highPurity, i, reco::TrackBase::quality(), and reco::TrackBase::undefQuality.

966  {
967  for(reco::TrackBase::TrackQuality i = reco::TrackBase::highPurity; i != reco::TrackBase::undefQuality; i = reco::TrackBase::TrackQuality(i - 1))
968  {
969  if(track->quality(i))
970  return i;
971  }
972 
973  return reco::TrackBase::undefQuality;
974 }
i
int i
Definition: DBlmapReader.cc:9
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:149
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:497
template<class Container , class Base >
void IPTagPlotter< Container, Base >::psPlot ( const std::string &  name)
virtual

Implements BaseBTagPlotter.

Definition at line 830 of file IPTagPlotter_cc.h.

References svgfig::canvas(), gen::n, RecoBTag::setTDRStyle(), and AlCaHLTBitMon_QueryRunRegistry::string.

831 {
832  const std::string cName("TrackIPPlots"+ theExtensionString);
833  RecoBTag::setTDRStyle()->cd();
834  TCanvas canvas(cName.c_str(), cName.c_str(), 600, 900);
835  canvas.UseCurrentStyle();
836  if (willFinalize_) {
837  for(int n=0; n != 2; ++n) {
838  canvas.Print((name + cName + ".ps").c_str());
839  canvas.Clear();
840  canvas.Divide(2,3);
841  canvas.cd(1);
842  effPurFromHistos[0+n]->discriminatorNoCutEffic()->plot();
843  canvas.cd(2);
844  effPurFromHistos[0+n]->discriminatorCutEfficScan()->plot();
845  canvas.cd(3);
846  effPurFromHistos[0+n]->plot();
847  canvas.cd(4);
848  effPurFromHistos[1+n]->discriminatorNoCutEffic()->plot();
849  canvas.cd(5);
850  effPurFromHistos[1+n]->discriminatorCutEfficScan()->plot();
851  canvas.cd(6);
852  effPurFromHistos[1+n]->plot();
853  }
854  return;
855  }
856 
857  canvas.Print((name + cName + ".ps[").c_str());
858  canvas.Clear();
859  canvas.Divide(2,3);
860 
861  canvas.cd(1);
862  trkNbr3D->plot();
863  canvas.cd(2);
864  tkcntHistosSig3D[4]->plot();
865  for(int n=0; n < 4; n++) {
866  canvas.cd(3+n);
867  tkcntHistosSig3D[n]->plot();
868  }
869 
870  canvas.Print((name + cName + ".ps").c_str());
871  canvas.Clear();
872  canvas.Divide(2,3);
873 
874  canvas.cd(1);
875  trkNbr3D->plot();
876  canvas.cd(2);
877  tkcntHistosProb3D[4]->plot();
878  for(int n=0; n < 4; n++) {
879  canvas.cd(3+n);
880  tkcntHistosProb3D[n]->plot();
881  }
882 
883  canvas.Print((name + cName + ".ps").c_str());
884  canvas.Clear();
885  canvas.Divide(2,3);
886  canvas.cd(1);
887  trkNbr2D->plot();
888  canvas.cd(2);
889  tkcntHistosSig2D[4]->plot();
890  for(int n=0; n != 4; ++n) {
891  canvas.cd(3+n);
892  tkcntHistosSig2D[n]->plot();
893  }
894 
895  canvas.Print((name + cName + ".ps").c_str());
896  canvas.Clear();
897  canvas.Divide(2,3);
898  canvas.cd(1);
899  trkNbr2D->plot();
900  canvas.cd(2);
901  tkcntHistosProb2D[4]->plot();
902  for(int n=0; n != 4; ++n) {
903  canvas.cd(3+n);
904  tkcntHistosProb2D[n]->plot();
905  }
906 
907  canvas.Print((name + cName + ".ps").c_str());
908  canvas.Clear();
909  canvas.Divide(1,3);
910  canvas.cd(1);
911  ghostTrackDistanceValuHisto->plot();
912  canvas.cd(2);
913  ghostTrackDistanceSignHisto->plot();
914  canvas.cd(3);
915  ghostTrackWeightHisto->plot();
916 
917  canvas.Print((name + cName + ".ps").c_str());
918  canvas.Print((name + cName + ".ps]").c_str());
919 }
FlavourHistograms::plot
void plot(TPad *theCanvas=0)
Definition: FlavourHistorgrams.h:367
IPTagPlotter::ghostTrackDistanceValuHisto
TrackIPHistograms< double > * ghostTrackDistanceValuHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosProb3D
TrackIPHistograms< float > * tkcntHistosProb3D[5]
Definition: IPTagPlotter.h:68
IPTagPlotter::ghostTrackWeightHisto
TrackIPHistograms< double > * ghostTrackWeightHisto
Definition: IPTagPlotter.h:72
BaseBTagPlotter::theExtensionString
const std::string theExtensionString
Definition: BaseBTagPlotter.h:32
RecoBTag::setTDRStyle
TStyle * setTDRStyle()
Definition: Tools.cc:287
IPTagPlotter::effPurFromHistos
EffPurFromHistos * effPurFromHistos[4]
Definition: IPTagPlotter.h:66
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
svgfig.canvas
def canvas
Definition: svgfig.py:481
IPTagPlotter::trkNbr2D
FlavourHistograms< int > * trkNbr2D
Definition: IPTagPlotter.h:61
EffPurFromHistos::discriminatorCutEfficScan
FlavourHistograms< double > * discriminatorCutEfficScan() const
Definition: EffPurFromHistos.h:61
EffPurFromHistos::plot
void plot(TPad *theCanvas=0)
Definition: EffPurFromHistos.cc:198
EffPurFromHistos::discriminatorNoCutEffic
FlavourHistograms< double > * discriminatorNoCutEffic() const
Definition: EffPurFromHistos.h:60
IPTagPlotter::ghostTrackDistanceSignHisto
TrackIPHistograms< double > * ghostTrackDistanceSignHisto
Definition: IPTagPlotter.h:73
IPTagPlotter::tkcntHistosSig3D
TrackIPHistograms< double > * tkcntHistosSig3D[5]
Definition: IPTagPlotter.h:41
IPTagPlotter::trkNbr3D
FlavourHistograms< int > * trkNbr3D
Definition: IPTagPlotter.h:61
IPTagPlotter::tkcntHistosProb2D
TrackIPHistograms< float > * tkcntHistosProb2D[5]
Definition: IPTagPlotter.h:69
IPTagPlotter::tkcntHistosSig2D
TrackIPHistograms< double > * tkcntHistosSig2D[5]
Definition: IPTagPlotter.h:42
IPTagPlotter::willFinalize_
bool willFinalize_
Definition: IPTagPlotter.h:38

Member Data Documentation

template<class Container, class Base>
EffPurFromHistos* IPTagPlotter< Container, Base >::effPurFromHistos[4]
private

Definition at line 66 of file IPTagPlotter.h.

template<class Container, class Base>
double IPTagPlotter< Container, Base >::endEffPur_
private

Definition at line 36 of file IPTagPlotter.h.

template<class Container, class Base>
TrackIPHistograms<double> * IPTagPlotter< Container, Base >::ghostTrackDistanceSignHisto
private

Definition at line 73 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::ghostTrackDistanceValuHisto
private

Definition at line 73 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::ghostTrackWeightHisto
private

Definition at line 72 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::lowerIPBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::lowerIPEBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::lowerIPSBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
bool IPTagPlotter< Container, Base >::makeQualityPlots_
private

Definition at line 39 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::maxDecayLength
private

Definition at line 64 of file IPTagPlotter.h.

template<class Container, class Base>
double IPTagPlotter< Container, Base >::maxJetDistance
private

Definition at line 64 of file IPTagPlotter.h.

template<class Container, class Base>
unsigned int IPTagPlotter< Container, Base >::mcPlots_
private

Definition at line 37 of file IPTagPlotter.h.

template<class Container, class Base>
double IPTagPlotter< Container, Base >::minDecayLength
private

Definition at line 64 of file IPTagPlotter.h.

template<class Container, class Base>
double IPTagPlotter< Container, Base >::minJetDistance
private

Definition at line 64 of file IPTagPlotter.h.

template<class Container, class Base>
int IPTagPlotter< Container, Base >::nBinEffPur_
private

Definition at line 34 of file IPTagPlotter.h.

template<class Container, class Base>
int IPTagPlotter< Container, Base >::nBinsIP
private

Definition at line 63 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
int IPTagPlotter< Container, Base >::nBinsIPE
private

Definition at line 63 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
int IPTagPlotter< Container, Base >::nBinsIPS
private

Definition at line 63 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms2D<double, int>* IPTagPlotter< Container, Base >::selectedTrackMultVsJetPtHisto
private

Definition at line 78 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms<int>* IPTagPlotter< Container, Base >::selectedTrackQualHisto
private

Definition at line 76 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::startEffPur_
private

Definition at line 35 of file IPTagPlotter.h.

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosDecayLengthVal2D[5]
private

Definition at line 47 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosDecayLengthVal3D[5]
private

Definition at line 48 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosErr2D[5]
private

Definition at line 44 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosErr3D[5]
private

Definition at line 43 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosJetDistSign2D[5]
private

Definition at line 51 of file IPTagPlotter.h.

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosJetDistSign3D[5]
private

Definition at line 52 of file IPTagPlotter.h.

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosJetDistVal2D[5]
private

Definition at line 49 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosJetDistVal3D[5]
private

Definition at line 50 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float>* IPTagPlotter< Container, Base >::tkcntHistosProb2D[5]
private

Definition at line 69 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float>* IPTagPlotter< Container, Base >::tkcntHistosProb3D[5]
private

Definition at line 68 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosSig2D[5]
private

Definition at line 42 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosSig3D[5]
private

Definition at line 41 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosTkNChiSqr2D[5]
private

Definition at line 53 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosTkNChiSqr3D[5]
private

Definition at line 54 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<int>* IPTagPlotter< Container, Base >::tkcntHistosTkNHits2D[5]
private

Definition at line 57 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<int>* IPTagPlotter< Container, Base >::tkcntHistosTkNHits3D[5]
private

Definition at line 58 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<int>* IPTagPlotter< Container, Base >::tkcntHistosTkNPixelHits2D[5]
private

Definition at line 59 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<int>* IPTagPlotter< Container, Base >::tkcntHistosTkNPixelHits3D[5]
private

Definition at line 60 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float>* IPTagPlotter< Container, Base >::tkcntHistosTkProbIPneg2D
private

Definition at line 70 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float>* IPTagPlotter< Container, Base >::tkcntHistosTkProbIPneg3D
private

Definition at line 71 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float> * IPTagPlotter< Container, Base >::tkcntHistosTkProbIPpos2D
private

Definition at line 70 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<float> * IPTagPlotter< Container, Base >::tkcntHistosTkProbIPpos3D
private

Definition at line 71 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosTkPt2D[5]
private

Definition at line 55 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosTkPt3D[5]
private

Definition at line 56 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosVal2D[5]
private

Definition at line 46 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
TrackIPHistograms<double>* IPTagPlotter< Container, Base >::tkcntHistosVal3D[5]
private

Definition at line 45 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms2D<double, int>* IPTagPlotter< Container, Base >::trackMultVsJetPtHisto
private

Definition at line 77 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms<int>* IPTagPlotter< Container, Base >::trackQualHisto
private

Definition at line 75 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms<int> * IPTagPlotter< Container, Base >::trkNbr2D
private

Definition at line 61 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
FlavourHistograms<int>* IPTagPlotter< Container, Base >::trkNbr3D
private

Definition at line 61 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::upperIPBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::upperIPEBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
double IPTagPlotter< Container, Base >::upperIPSBound
private

Definition at line 62 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

template<class Container, class Base>
bool IPTagPlotter< Container, Base >::willFinalize_
private

Definition at line 38 of file IPTagPlotter.h.

Referenced by IPTagPlotter< Container, Base >::IPTagPlotter().

Generated for CMSSW Reference Manual by doxygen 1.8.5