CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
BDHadronTrackMonitoringHarvester Class Reference

#include <BDHadronTrackMonitoringHarvester.h>

Inheritance diagram for BDHadronTrackMonitoringHarvester:

Public Member Functions

 BDHadronTrackMonitoringHarvester (const edm::ParameterSet &pSet)
 
 ~BDHadronTrackMonitoringHarvester () override
 

Private Member Functions

void beginJob () override
 
void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 

Private Attributes

MonitorElement * nTrk_absolute_bjet
 
MonitorElement * nTrk_absolute_cjet
 
MonitorElement * nTrk_absolute_dusgjet
 
MonitorElement * nTrk_relative_bjet
 
MonitorElement * nTrk_relative_cjet
 
MonitorElement * nTrk_relative_dusgjet
 
MonitorElement * nTrk_std_bjet
 
MonitorElement * nTrk_std_cjet
 
MonitorElement * nTrk_std_dusgjet
 

Detailed Description

Top level steering routine for b tag performance harvesting.

Definition at line 24 of file BDHadronTrackMonitoringHarvester.h.

Constructor & Destructor Documentation

BDHadronTrackMonitoringHarvester::BDHadronTrackMonitoringHarvester ( const edm::ParameterSet &  pSet)
explicit

Definition at line 11 of file BDHadronTrackMonitoringHarvester.cc.

11 {}
BDHadronTrackMonitoringHarvester::~BDHadronTrackMonitoringHarvester ( )
override

Definition at line 13 of file BDHadronTrackMonitoringHarvester.cc.

13 {}

Member Function Documentation

void BDHadronTrackMonitoringHarvester::beginJob ( void  )
overrideprivate

Definition at line 15 of file BDHadronTrackMonitoringHarvester.cc.

15 {}
void BDHadronTrackMonitoringHarvester::dqmEndJob ( DQMStore::IBooker &  ibook,
DQMStore::IGetter &  iget 
)
overrideprivate

Definition at line 17 of file BDHadronTrackMonitoringHarvester.cc.

References DQMStore::IBooker::book1D(), DEFINE_FWK_MODULE, DQMStore::IGetter::get(), MonitorElement::getMean(), MonitorElement::getRMS(), mps_fire::i, SiStripPI::max, nTrk_absolute_bjet, nTrk_absolute_cjet, nTrk_absolute_dusgjet, nTrk_relative_bjet, nTrk_relative_cjet, nTrk_relative_dusgjet, nTrk_std_bjet, nTrk_std_cjet, nTrk_std_dusgjet, DQMStore::IGetter::removeElement(), MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinContent(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), RecoBTag::setTDRStyle(), and BDHadronTrackMonitoringAnalyzer::TrkHistCat.

17  {
18  // ***********************
19  //
20  // Book all new histograms.
21  //
22  // ***********************
23  RecoBTag::setTDRStyle();
24  ibook.setCurrentFolder("BDHadronTracks/JetContent");
25 
26  // b jets
27  // absolute average number of tracks
28  nTrk_absolute_bjet = ibook.book1D("nTrk_absolute_bjet", "absolute average number of tracks in b jets", 6, -0.5, 5.5);
29  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
30  nTrk_absolute_bjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
31  }
32  nTrk_absolute_bjet->setAxisRange(0, 5, 2);
33  nTrk_absolute_bjet->setAxisTitle("average number of tracks", 2);
34 
35  // relative (in percent) average number of tracks
36  nTrk_relative_bjet = ibook.book1D("nTrk_relative_bjet", "relative average number of tracks in b jets", 6, -0.5, 5.5);
37  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
38  nTrk_relative_bjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
39  }
40  nTrk_relative_bjet->setAxisRange(0, 1, 2);
41  nTrk_relative_bjet->setAxisTitle("average fraction of tracks", 2);
42 
43  // standard deviation of number of tracks
44  nTrk_std_bjet = ibook.book1D("nTrk_std_bjet", "RMS of number of tracks in b jets", 6, -0.5, 5.5);
45  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
46  nTrk_std_bjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
47  }
48  nTrk_std_bjet->setAxisRange(0, 3, 2);
49  nTrk_std_bjet->setAxisTitle("RMS of number of tracks", 2);
50 
51  // c jets
52  nTrk_absolute_cjet = ibook.book1D("nTrk_absolute_cjet", "absolute average number of tracks in c jets", 6, -0.5, 5.5);
53  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
54  nTrk_absolute_cjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
55  }
56  nTrk_absolute_cjet->setAxisRange(0, 5, 2);
57  nTrk_absolute_cjet->setAxisTitle("average number of tracks", 2);
58 
59  nTrk_relative_cjet = ibook.book1D("nTrk_relative_cjet", "relative average number of tracks in c jets", 6, -0.5, 5.5);
60  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
61  nTrk_relative_cjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
62  }
63  nTrk_relative_cjet->setAxisRange(0, 1, 2);
64  nTrk_relative_cjet->setAxisTitle("average fraction of tracks", 2);
65 
66  nTrk_std_cjet = ibook.book1D("nTrk_std_cjet", "RMS of number of tracks in c jets", 6, -0.5, 5.5);
67  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
68  nTrk_std_cjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
69  }
70  nTrk_std_cjet->setAxisRange(0, 3, 2);
71  nTrk_std_cjet->setAxisTitle("RMS of number of tracks", 2);
72 
73  // udsg jets
74  nTrk_absolute_dusgjet =
75  ibook.book1D("nTrk_absolute_dusgjet", "absolute average number of tracks in dusg jets", 6, -0.5, 5.5);
76  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
77  nTrk_absolute_dusgjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
78  }
79  nTrk_absolute_dusgjet->setAxisRange(0, 5, 2);
80  nTrk_absolute_dusgjet->setAxisTitle("average number of tracks", 2);
81 
82  nTrk_relative_dusgjet =
83  ibook.book1D("nTrk_relative_dusgjet", "relative average number of tracks in dusg jets", 6, -0.5, 5.5);
84  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
85  nTrk_relative_dusgjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
86  }
87  nTrk_relative_dusgjet->setAxisRange(0, 1, 2);
88  nTrk_relative_dusgjet->setAxisTitle("average fraction of tracks", 2);
89 
90  nTrk_std_dusgjet = ibook.book1D("nTrk_std_dusgjet", "RMS of number of tracks in dusg jets", 6, -0.5, 5.5);
91  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
92  nTrk_std_dusgjet->setBinLabel(i + 1, BDHadronTrackMonitoringAnalyzer::TrkHistCat[i], 1);
93  }
94  nTrk_std_dusgjet->setAxisRange(0, 3, 2);
95  nTrk_std_dusgjet->setAxisTitle("RMS of number of tracks", 2);
96 
97  // ***********************
98  //
99  // get all the old histograms
100  //
101  // ***********************
102 
103  // b jets
104  MonitorElement *nTrk_bjet[6];
105  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
106  nTrk_bjet[i] = iget.get("BDHadronTracks/JetContent/nTrk_bjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
107  }
108  MonitorElement *nTrkAll_bjet = iget.get("BDHadronTracks/JetContent/nTrkAll_bjet");
109 
110  // c jets
111  MonitorElement *nTrk_cjet[6];
112  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
113  nTrk_cjet[i] = iget.get("BDHadronTracks/JetContent/nTrk_cjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
114  }
115  MonitorElement *nTrkAll_cjet = iget.get("BDHadronTracks/JetContent/nTrkAll_cjet");
116 
117  // dusg jets
118  MonitorElement *nTrk_dusgjet[6];
119  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
120  nTrk_dusgjet[i] =
121  iget.get("BDHadronTracks/JetContent/nTrk_dusgjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
122  }
123  MonitorElement *nTrkAll_dusgjet = iget.get("BDHadronTracks/JetContent/nTrkAll_dusgjet");
124 
125  // ***********************
126  //
127  // Calculate contents of new histograms
128  //
129  // ***********************
130 
131  // b jets
132  float mean_bjets[6];
133  float std_bjets[6];
134  float meanAll_bjets;
135  meanAll_bjets = std::max(0.01, nTrkAll_bjet->getMean(1));
136  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
137  mean_bjets[i] = nTrk_bjet[i]->getMean(1); // mean number of tracks per category
138  std_bjets[i] = nTrk_bjet[i]->getRMS(1);
139  nTrk_absolute_bjet->setBinContent(i + 1, mean_bjets[i]);
140  nTrk_relative_bjet->setBinContent(i + 1, mean_bjets[i] / meanAll_bjets);
141  nTrk_std_bjet->setBinContent(i + 1, std_bjets[i]);
142  }
143 
144  // c jets
145  float mean_cjets[6];
146  float std_cjets[6];
147  float meanAll_cjets;
148  meanAll_cjets = std::max(0.01, nTrkAll_cjet->getMean(1));
149  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
150  mean_cjets[i] = nTrk_cjet[i]->getMean(1); // mean number of tracks per category
151  std_cjets[i] = nTrk_cjet[i]->getRMS(1);
152  nTrk_absolute_cjet->setBinContent(i + 1, mean_cjets[i]);
153  nTrk_relative_cjet->setBinContent(i + 1, mean_cjets[i] / meanAll_cjets);
154  nTrk_std_cjet->setBinContent(i + 1, std_cjets[i]);
155  }
156 
157  // dusg jets
158  float mean_dusgjets[6];
159  float std_dusgjets[6];
160  float meanAll_dusgjets;
161  meanAll_dusgjets = std::max(0.01, nTrkAll_dusgjet->getMean(1));
162  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
163  mean_dusgjets[i] = nTrk_dusgjet[i]->getMean(1); // mean number of tracks per category
164  std_dusgjets[i] = nTrk_dusgjet[i]->getRMS(1);
165  nTrk_absolute_dusgjet->setBinContent(i + 1, mean_dusgjets[i]);
166  nTrk_relative_dusgjet->setBinContent(i + 1, mean_dusgjets[i] / meanAll_dusgjets);
167  nTrk_std_dusgjet->setBinContent(i + 1, std_dusgjets[i]);
168  }
169 
170  // ***********************
171  //
172  // Remove histograms with nTracks that we do not want to show
173  //
174  // ***********************
175 
176  // b jets
177  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
178  iget.removeElement("nTrk_bjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
179  }
180  iget.removeElement("nTrkAll_bjet");
181 
182  // c jets
183  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
184  iget.removeElement("nTrk_cjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
185  }
186  iget.removeElement("nTrkAll_cjet");
187 
188  // dusg jets
189  for (unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
190  iget.removeElement("nTrk_dusgjet_" + BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
191  }
192  iget.removeElement("nTrkAll_dusgjet");
193 }
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:866
mps_fire.i
i
Definition: mps_fire.py:338
MonitorElement::setAxisRange
void setAxisRange(double xmin, double xmax, int axis=1)
set x-, y- or z-axis range (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:964
RecoBTag::setTDRStyle
TStyle * setTDRStyle()
Definition: Tools.cc:346
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:942
MonitorElement::getMean
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:722
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
BDHadronTrackMonitoringAnalyzer::TrkHistCat
static const std::vector< std::string > TrkHistCat
Definition: BDHadronTrackMonitoringAnalyzer.h:61
DQMStore::IGetter::get
MonitorElement * get(std::string const &path)
Definition: DQMStore.cc:303
MonitorElement::getRMS
double getRMS(int axis=1) const
get RMS of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:735
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:973
DQMStore::IGetter::removeElement
void removeElement(Args &&...args)
Definition: DQMStore.h:199

Member Data Documentation

MonitorElement* BDHadronTrackMonitoringHarvester::nTrk_absolute_bjet
private

Definition at line 35 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement* BDHadronTrackMonitoringHarvester::nTrk_absolute_cjet
private

Definition at line 40 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement* BDHadronTrackMonitoringHarvester::nTrk_absolute_dusgjet
private

Definition at line 45 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_relative_bjet
private

Definition at line 35 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_relative_cjet
private

Definition at line 40 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_relative_dusgjet
private

Definition at line 45 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_std_bjet
private

Definition at line 35 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_std_cjet
private

Definition at line 40 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

MonitorElement * BDHadronTrackMonitoringHarvester::nTrk_std_dusgjet
private

Definition at line 45 of file BDHadronTrackMonitoringHarvester.h.

Referenced by dqmEndJob().

Generated for CMSSW Reference Manual by doxygen 1.8.11