CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
TrackingQualityChecker.cc
Go to the documentation of this file.
1 #include "DQM/TrackingMonitorClient/interface/TrackingQualityChecker.h"
2 #include "FWCore/ServiceRegistry/interface/Service.h"
3 #include "DQMServices/Core/interface/DQMStore.h"
4 #include "DQMServices/Core/interface/MonitorElement.h"
5 #include "DQMServices/Core/interface/QReport.h"
6 
7 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
8 
9 #include "CalibTracker/SiStripCommon/interface/TkDetMap.h"
10 
11 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
12 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
13 #include "Geometry/Records/interface/IdealGeometryRecord.h"
14 
15 #include "DQM/TrackingMonitorClient/interface/TrackingUtility.h"
16 
17 #include "FWCore/MessageLogger/interface/MessageLogger.h"
18 #include "FWCore/ParameterSet/interface/ParameterSet.h"
19 
20 
21 #include <iomanip>
22 //
23 // -- Constructor
24 //
25 TrackingQualityChecker::TrackingQualityChecker(edm::ParameterSet const& ps) :
26  pSet_(ps)
27 {
28  edm::LogInfo("TrackingQualityChecker") << " Creating TrackingQualityChecker " << "\n" ;
29 
30  bookedTrackingGlobalStatus_ = false;
31  bookedTrackingLSStatus_ = false;
32 
33  if(!edm::Service<TkDetMap>().isAvailable()){
34  edm::LogError("TkHistoMap") <<
35  "\n------------------------------------------"
36  "\nUnAvailable Service TkHistoMap: please insert in the configuration file an instance like"
37  "\n\tprocess.TkDetMap = cms.Service(\"TkDetMap\")"
38  "\n------------------------------------------";
39  }
40 
41  TopFolderName_ = pSet_.getUntrackedParameter<std::string>("TopFolderName","Tracking");
42 
43  TrackingMEs tracking_mes;
44  std::vector<edm::ParameterSet> TrackingGlobalQualityMEs = pSet_.getParameter< std::vector<edm::ParameterSet> >("TrackingGlobalQualityPSets" );
45  for ( auto meQTset : TrackingGlobalQualityMEs ) {
46 
47  std::string QTname = meQTset.getParameter<std::string>("QT");
48  tracking_mes.HistoDir = meQTset.getParameter<std::string>("dir");
49  tracking_mes.HistoName = meQTset.getParameter<std::string>("name");
50  // std::cout << "[TrackingQualityChecker::TrackingQualityChecker] inserting " << QTname << " in TrackingMEsMap" << std::endl;
51  TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>(QTname, tracking_mes));
52  }
53  // std::cout << "[TrackingQualityChecker::TrackingQualityChecker] created TrackingMEsMap" << std::endl;
54 
55  TrackingLSMEs tracking_ls_mes;
56  std::vector<edm::ParameterSet> TrackingLSQualityMEs = pSet_.getParameter< std::vector<edm::ParameterSet> >("TrackingLSQualityPSets" );
57  for ( auto meQTset : TrackingLSQualityMEs ) {
58 
59  std::string QTname = meQTset.getParameter<std::string>("QT");
60  tracking_ls_mes.HistoLSDir = meQTset.exists("LSdir") ? meQTset.getParameter<std::string>("LSdir") : "";
61  tracking_ls_mes.HistoLSName = meQTset.exists("LSname") ? meQTset.getParameter<std::string>("LSname") : "";
62  tracking_ls_mes.HistoLSLowerCut = meQTset.exists("LSlowerCut") ? meQTset.getParameter<double>("LSlowerCut") : -1.;
63  tracking_ls_mes.HistoLSUpperCut = meQTset.exists("LSupperCut") ? meQTset.getParameter<double>("LSupperCut") : -1.;
64  tracking_ls_mes.TrackingFlag = 0;
65 
66  // std::cout << "[TrackingQualityChecker::TrackingQualityChecker] inserting " << QTname << " in TrackingMEsMap" << std::endl;
67  TrackingLSMEsMap.insert(std::pair<std::string, TrackingLSMEs>(QTname, tracking_ls_mes));
68  }
69  // std::cout << "[TrackingQualityChecker::TrackingQualityChecker] created TrackingLSMEsMap" << std::endl;
70 
71 }
72 //
73 // -- Destructor
74 //
75 TrackingQualityChecker::~TrackingQualityChecker() {
76  edm::LogInfo("TrackingQualityChecker") << " Deleting TrackingQualityChecker " << "\n" ;
77 }
78 //
79 // -- create reportSummary MEs
80 //
81 void TrackingQualityChecker::bookGlobalStatus(DQMStore* dqm_store) {
82 
83  // std::cout << "[TrackingQualityChecker::bookGlobalStatus] already booked ? " << (bookedTrackingGlobalStatus_ ? "yes" : "nope") << std::endl;
84 
85  if (!bookedTrackingGlobalStatus_) {
86  dqm_store->cd();
87  edm::LogInfo("TrackingQualityChecker") << " booking TrackingQualityStatus" << "\n";
88 
89  std::string tracking_dir = "";
90  TrackingUtility::getTopFolderPath(dqm_store, TopFolderName_, tracking_dir);
91  dqm_store->setCurrentFolder(TopFolderName_+"/EventInfo");
92 
93  TrackGlobalSummaryReportGlobal = dqm_store->bookFloat("reportSummary");
94 
95  std::string hname, htitle;
96  hname = "reportSummaryMap";
97  htitle = "Tracking Report Summary Map";
98 
99  size_t nQT = TrackingMEsMap.size();
100  // std::cout << "[TrackingQualityChecker::bookGlobalStatus] nQT: " << nQT << std::endl;
101  TrackGlobalSummaryReportMap = dqm_store->book2D(hname, htitle, nQT,0.5,float(nQT)+0.5,1,0.5,1.5);
102  TrackGlobalSummaryReportMap->setAxisTitle("Track Quality Type", 1);
103  TrackGlobalSummaryReportMap->setAxisTitle("QTest Flag", 2);
104  size_t ibin =0;
105  for ( auto meQTset : TrackingMEsMap ) {
106  TrackGlobalSummaryReportMap->setBinLabel(ibin+1,meQTset.first);
107  ibin++;
108  }
109 
110  dqm_store->setCurrentFolder(TopFolderName_+"/EventInfo/reportSummaryContents");
111 
112  for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
113  it != TrackingMEsMap.end(); it++) {
114  std::string meQTname = it->first;
115  it->second.TrackingFlag = dqm_store->bookFloat("Track"+meQTname);
116  // std::cout << "[TrackingQualityChecker::bookGlobalStatus] " << it->first << " exists ? " << it->second.TrackingFlag << std::endl;
117  // std::cout << "[TrackingQualityChecker::bookGlobalStatus] DONE w/ TrackingMEsMap" << std::endl;
118  }
119 
120  bookedTrackingGlobalStatus_ = true;
121  dqm_store->cd();
122  }
123 }
124 
125 void TrackingQualityChecker::bookLSStatus(DQMStore* dqm_store) {
126 
127  // std::cout << "[TrackingQualityChecker::bookLSStatus] already booked ? " << (bookedTrackingLSStatus_ ? "yes" : "nope") << std::endl;
128 
129  if (!bookedTrackingLSStatus_) {
130  dqm_store->cd();
131  edm::LogInfo("TrackingQualityChecker") << " booking TrackingQualityStatus" << "\n";
132 
133  std::string tracking_dir = "";
134  TrackingUtility::getTopFolderPath(dqm_store, TopFolderName_, tracking_dir);
135  dqm_store->setCurrentFolder(TopFolderName_+"/EventInfo");
136 
137  TrackLSSummaryReportGlobal = dqm_store->bookFloat("reportSummary");
138 
139  std::string hname, htitle;
140  hname = "reportSummaryMap";
141  htitle = "Tracking Report Summary Map";
142 
143  // size_t nQT = TrackingLSMEsMap.size();
144  // std::cout << "[TrackingQualityChecker::bookLSStatus] nQT: " << nQT << std::endl;
145 
146  dqm_store->setCurrentFolder(TopFolderName_+"/EventInfo/reportSummaryContents");
147  for (std::map<std::string, TrackingLSMEs>::iterator it = TrackingLSMEsMap.begin();
148  it != TrackingLSMEsMap.end(); it++) {
149  std::string meQTname = it->first;
150  it->second.TrackingFlag = dqm_store->bookFloat("Track"+meQTname);
151  // std::cout << "[TrackingQualityChecker::bookLSStatus] " << it->first << " exists ? " << it->second.TrackingFlag << std::endl;
152  // std::cout << "[TrackingQualityChecker::bookLSStatus] DONE w/ TrackingLSMEsMap" << std::endl;
153  }
154 
155  bookedTrackingLSStatus_ = true;
156  dqm_store->cd();
157  }
158 }
159 
160 //
161 // -- Fill Dummy Status
162 //
163 void TrackingQualityChecker::fillDummyGlobalStatus(){
164  // std::cout << "[TrackingQualityChecker::fillDummyGlobalStatus] starting ..." << std::endl;
165 
166  resetGlobalStatus();
167  // std::cout << "[TrackingQualityChecker::fillDummyGlobalStatus] already booked ? " << (bookedTrackingGlobalStatus_ ? "yes" : "nope") << std::endl;
168  if (bookedTrackingGlobalStatus_) {
169 
170  TrackGlobalSummaryReportGlobal->Fill(-1.0);
171 
172  for (int ibin = 1; ibin < TrackGlobalSummaryReportMap->getNbinsX()+1; ibin++) {
173  fillStatusHistogram(TrackGlobalSummaryReportMap, ibin, 1, -1.0);
174  }
175 
176  for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
177  it != TrackingMEsMap.end(); it++)
178  it->second.TrackingFlag->Fill(-1.0);
179  // std::cout << "[TrackingQualityChecker::fillDummyGlobalStatus] DONE w/ TrackingMEsMap" << std::endl;
180 
181  }
182 }
183 void TrackingQualityChecker::fillDummyLSStatus(){
184  // std::cout << "[TrackingQualityChecker::fillDummyLSStatus] starting ..." << std::endl;
185 
186  resetLSStatus();
187  // std::cout << "[TrackingQualityChecker::fillDummyLSStatus] already booked ? " << (bookedTrackingLSStatus_ ? "yes" : "nope") << std::endl;
188  if (bookedTrackingLSStatus_) {
189 
190  TrackLSSummaryReportGlobal->Fill(-1.0);
191  for (std::map<std::string, TrackingLSMEs>::iterator it = TrackingLSMEsMap.begin();
192  it != TrackingLSMEsMap.end(); it++)
193  it->second.TrackingFlag->Fill(-1.0);
194  // std::cout << "[TrackingQualityChecker::fillDummyLSStatus] DONE w/ TrackingLSMEsMap" << std::endl;
195 
196  }
197 }
198 
199 //
200 // -- Reset Status
201 //
202 void TrackingQualityChecker::resetGlobalStatus() {
203 
204  // std::cout << "[TrackingQualityChecker::resetGlobalStatus] already booked ? " << (bookedTrackingGlobalStatus_ ? "yes" : "nope") << std::endl;
205  if (bookedTrackingGlobalStatus_) {
206 
207  TrackGlobalSummaryReportGlobal -> Reset();
208  TrackGlobalSummaryReportMap -> Reset();
209 
210  for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
211  it != TrackingMEsMap.end(); it++) {
212  MonitorElement* me = it->second.TrackingFlag;
213  // std::cout << "[TrackingQualityChecker::resetGlobalStatus] " << it->second.HistoName << " exist ? " << ( it->second.TrackingFlag == NULL ? "nope" : "yes" ) << " ---> " << me << std::endl;
214  me->Reset();
215  }
216  // std::cout << "[TrackingQualityChecker::resetGlobalStatus] DONE w/ TrackingMEsMap" << std::endl;
217 
218  }
219 }
220 void TrackingQualityChecker::resetLSStatus() {
221 
222  // std::cout << "[TrackingQualityChecker::resetLSStatus] already booked ? " << (bookedTrackingLSStatus_ ? "yes" : "nope") << std::endl;
223  if (bookedTrackingLSStatus_) {
224 
225  TrackLSSummaryReportGlobal -> Reset();
226  for (std::map<std::string, TrackingLSMEs>::iterator it = TrackingLSMEsMap.begin();
227  it != TrackingLSMEsMap.end(); it++) {
228  MonitorElement* me = it->second.TrackingFlag;
229  // std::cout << "[TrackingQualityChecker::resetLSStatus] " << it->second.HistoLSName << " exist ? " << ( it->second.TrackingFlag == NULL ? "nope" : "yes" ) << " ---> " << me << std::endl;
230  me->Reset();
231  }
232  // std::cout << "[TrackingQualityChecker::resetLSStatus] DONE w/ TrackingLSMEsMap" << std::endl;
233 
234  }
235 }
236 
237 //
238 // -- Fill Status
239 //
240 void TrackingQualityChecker::fillGlobalStatus(DQMStore* dqm_store) {
241 
242  // std::cout << "[TrackingQualityChecker::fillGlobalStatus] already booked ? " << (bookedTrackingGlobalStatus_ ? "yes" : "nope") << std::endl;
243  if (!bookedTrackingGlobalStatus_) bookGlobalStatus(dqm_store);
244 
245  fillDummyGlobalStatus();
246  fillTrackingStatus(dqm_store);
247  // std::cout << "[TrackingQualityChecker::fillGlobalStatus] DONE" << std::endl;
248  dqm_store->cd();
249 }
250 
251 void TrackingQualityChecker::fillLSStatus(DQMStore* dqm_store) {
252 
253  // std::cout << "[TrackingQualityChecker::fillLSStatus] already booked ? " << (bookedTrackingLSStatus_ ? "yes" : "nope") << std::endl;
254  if (!bookedTrackingLSStatus_) bookLSStatus(dqm_store);
255 
256  fillDummyLSStatus();
257  fillTrackingStatusAtLumi(dqm_store);
258  // std::cout << "[TrackingQualityChecker::fillLSStatus] DONE" << std::endl;
259  dqm_store->cd();
260 }
261 
262 //
263 // -- Fill Tracking Status
264 //
265 void TrackingQualityChecker::fillTrackingStatus(DQMStore* dqm_store) {
266 
267  float gstatus = 0.0;
268 
269  dqm_store->cd();
270  if (!TrackingUtility::goToDir(dqm_store, TopFolderName_)) return;
271 
272 
273  int ibin = 0;
274  for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
275  it != TrackingMEsMap.end(); it++) {
276 
277  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] ME: " << it->first << " [" << it->second.TrackingFlag->getFullname() << "] flag: " << it->second.TrackingFlag->getFloatValue() << std::endl;
278 
279  ibin++;
280 
281  std::string localMEdirpath = it->second.HistoDir;
282  std::string MEname = it->second.HistoName;
283 
284  std::vector<MonitorElement*> tmpMEvec = dqm_store->getContents(dqm_store->pwd()+"/"+localMEdirpath);
285  MonitorElement* me = NULL;
286 
287  size_t nMEs = 0;
288  for ( auto ime : tmpMEvec ) {
289  std::string name = ime->getName();
290  if ( name.find(MEname) != std::string::npos) {
291  me = ime;
292  nMEs++;
293  }
294  }
295  // only one ME found
296  if (nMEs == 1) {
297  float status = 0.;
298  for ( auto ime : tmpMEvec ) {
299  std::string name = ime->getName();
300  if ( name.find(MEname) != std::string::npos) {
301  me = ime;
302  }
303  }
304  if (!me) continue;
305  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] status: " << status << std::endl;
306  std::vector<QReport *> qt_reports = me->getQReports();
307  size_t nQTme = qt_reports.size();
308  if (nQTme != 0) {
309  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] qt_reports: " << qt_reports.size() << std::endl;
310  // loop on possible QTs
311  for ( auto iQT : qt_reports ) {
312  status += iQT->getQTresult();
313  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] iQT: " << iQT->getQRName() << std::endl;
314  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] MEname: " << MEname << " status: " << iQT->getQTresult() << " exists ? " << (it->second.TrackingFlag ? "yes " : "no ") << it->second.TrackingFlag << std::endl;
315  }
316  status = status/float(nQTme);
317  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] MEname: " << MEname << " status: " << status << std::endl;
318  it->second.TrackingFlag->Fill(status);
319  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] TrackGlobalSummaryReportMap: " << TrackGlobalSummaryReportMap << std::endl;
320  fillStatusHistogram(TrackGlobalSummaryReportMap, ibin, 1, status);
321  }
322 
323  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] gstatus: " << gstatus << " x status: " << status << std::endl;
324  if ( status < 0. ) gstatus = -1.;
325  else gstatus += status;
326  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] ===> gstatus: " << gstatus << std::endl;
327  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] ME: " << it->first << " [" << it->second.TrackingFlag->getFullname() << "] flag: " << it->second.TrackingFlag->getFloatValue() << std::endl;
328 
329  } else { // more than 1 ME w/ the same root => they need to be considered together
330  float status = 1.;
331  for ( auto ime : tmpMEvec ) {
332  float tmp_status = 1.;
333  std::string name = ime->getName();
334  if ( name.find(MEname) != std::string::npos) {
335  me = ime;
336 
337  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] status: " << status << std::endl;
338  std::vector<QReport *> qt_reports = me->getQReports();
339  size_t nQTme = qt_reports.size();
340  if (nQTme != 0) {
341  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] qt_reports: " << qt_reports.size() << std::endl;
342  // loop on possible QTs
343  for ( auto iQT : qt_reports ) {
344  tmp_status += iQT->getQTresult();
345  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] iQT: " << iQT->getQRName() << std::endl;
346  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] MEname: " << MEname << " status: " << iQT->getQTresult() << " exists ? " << (it->second.TrackingFlag ? "yes " : "no ") << it->second.TrackingFlag << std::endl;
347  }
348  tmp_status = tmp_status/float(nQTme);
349  }
350  }
351  status = fminf(tmp_status,status);
352  }
353  if ( status < 0. ) gstatus = -1.;
354  else gstatus += status;
355  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] MEname: " << MEname << " status: " << status << std::endl;
356  it->second.TrackingFlag->Fill(status);
357  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] TrackGlobalSummaryReportMap: " << TrackGlobalSummaryReportMap << std::endl;
358 
359  fillStatusHistogram(TrackGlobalSummaryReportMap, ibin, 1, status);
360  }
361  }
362 
363  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] gstatus: " << gstatus << std::endl;
364  size_t nQT = TrackingMEsMap.size();
365  if (gstatus < 1.) gstatus = -1.;
366  else gstatus = gstatus/float(nQT);
367 
368  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] ===> gstatus: " << gstatus << std::endl;
369  TrackGlobalSummaryReportGlobal->Fill(gstatus);
370  dqm_store->cd();
371 
372  // std::cout << "[TrackingQualityChecker::fillTrackingStatus] DONE" << std::endl;
373 
374 }
375 
376 //
377 // -- Fill Report Summary Map
378 //
379  void TrackingQualityChecker::fillStatusHistogram(MonitorElement* me, int xbin, int ybin, float val){
380  if (me && me->kind() == MonitorElement::DQM_KIND_TH2F) {
381  TH2F* th2d = me->getTH2F();
382  th2d->SetBinContent(xbin, ybin, val);
383  }
384  }
385 
386 // Fill Tracking Status MEs at the Lumi block
387 //
388 void TrackingQualityChecker::fillTrackingStatusAtLumi(DQMStore* dqm_store){
389 
390  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] starting .. " << std::endl;
391  float gstatus = 1.0;
392 
393  dqm_store->cd();
394  if (!TrackingUtility::goToDir(dqm_store, TopFolderName_)) return;
395 
396 
397  int ibin = 0;
398  for (std::map<std::string, TrackingLSMEs>::iterator it = TrackingLSMEsMap.begin();
399  it != TrackingLSMEsMap.end(); it++) {
400 
401  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] ME: " << it->first << " [" << it->second.TrackingFlag->getFullname() << "] flag: " << it->second.TrackingFlag->getFloatValue() << std::endl;
402 
403  ibin++;
404 
405  std::string localMEdirpath = it->second.HistoLSDir;
406  std::string MEname = it->second.HistoLSName;
407  float lower_cut = it->second.HistoLSLowerCut;
408  float upper_cut = it->second.HistoLSUpperCut;
409 
410  float status = 1.0;
411 
412  std::vector<MonitorElement*> tmpMEvec = dqm_store->getContents(dqm_store->pwd()+"/"+localMEdirpath);
413  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] tmpMEvec: " << tmpMEvec.size() << std::endl;
414 
415  MonitorElement* me = NULL;
416 
417  size_t nMEs = 0;
418  for ( auto ime : tmpMEvec ) {
419  std::string name = ime->getName();
420  if ( name.find(MEname) != std::string::npos) {
421  me = ime;
422  nMEs++;
423  }
424  }
425  // only one ME found
426  if (nMEs == 1) {
427  for ( auto ime : tmpMEvec ) {
428  std::string name = ime->getName();
429  if ( name.find(MEname) != std::string::npos) {
430  me = ime;
431  }
432  }
433  if (!me) continue;
434 
435  if (me->kind() == MonitorElement::DQM_KIND_TH1F) {
436  float x_mean = me->getMean();
437  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] MEname: " << MEname << " x_mean: " << x_mean << std::endl;
438  if (x_mean <= lower_cut || x_mean > upper_cut) status = 0.0;
439  else status = 1.0;
440  }
441  } else { // more than 1 ME w/ the same root => they need to be considered together
442  for ( auto ime : tmpMEvec ) {
443  float tmp_status = 1.;
444  std::string name = ime->getName();
445  if ( name.find(MEname) != std::string::npos) {
446  me = ime;
447  if (!me) continue;
448 
449  if (me->kind() == MonitorElement::DQM_KIND_TH1F) {
450  float x_mean = me->getMean();
451  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] MEname: " << MEname << "[" << me->getName() << "] x_mean: " << x_mean << std::endl;
452  if (x_mean <= lower_cut || x_mean > upper_cut) tmp_status = 0.0;
453  else tmp_status = 1.0;
454  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] tmp_status: " << tmp_status << std::endl;
455  }
456  }
457  status = fminf(tmp_status,status);
458  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] ==> status: " << status << std::endl;
459  } // loop on tmpMEvec
460  }
461  it->second.TrackingFlag->Fill(status);
462  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] ===> status: " << status << " [" << gstatus << "]" << std::endl;
463  if (status == 0.0) gstatus = -1.0;
464  else gstatus = gstatus * status;
465  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] ===> gstatus: " << gstatus << std::endl;
466  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] ME: " << it->first << " [" << it->second.TrackingFlag->getFullname() << "] flag: " << it->second.TrackingFlag->getFloatValue() << std::endl;
467  }
468  TrackLSSummaryReportGlobal->Fill(gstatus);
469  dqm_store->cd();
470 
471  // std::cout << "[TrackingQualityChecker::fillTrackingStatusAtLumi] DONE" << std::endl;
472 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
TrackingQualityChecker::fillTrackingStatus
void fillTrackingStatus(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:265
TrackingQualityChecker::TrackGlobalSummaryReportMap
MonitorElement * TrackGlobalSummaryReportMap
Definition: TrackingQualityChecker.h:64
TrackingUtility::goToDir
static bool goToDir(DQMStore *dqm_store, std::string name)
Definition: TrackingUtility.cc:182
TrackingQualityChecker::fillGlobalStatus
void fillGlobalStatus(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:240
DQMStore::cd
void cd(void)
go to top directory (ie. root)
Definition: DQMStore.cc:561
TrackingQualityChecker::TrackingMEsMap
std::map< std::string, TrackingMEs > TrackingMEsMap
Definition: TrackingQualityChecker.h:61
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:901
TrackingQualityChecker::bookLSStatus
void bookLSStatus(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:125
NULL
#define NULL
Definition: scimark2.h:8
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:681
TrackingUtility::getTopFolderPath
static void getTopFolderPath(DQMStore *dqm_store, std::string top_dir, std::string &path)
Definition: TrackingUtility.cc:231
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DQMStore::bookFloat
MonitorElement * bookFloat(const char *name)
Book float.
Definition: DQMStore.cc:809
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
TrackingQualityChecker::TrackingQualityChecker
TrackingQualityChecker(edm::ParameterSet const &ps)
Definition: TrackingQualityChecker.cc:25
TrackingQualityChecker::TrackLSSummaryReportGlobal
MonitorElement * TrackLSSummaryReportGlobal
Definition: TrackingQualityChecker.h:67
TrackingQualityChecker::fillStatusHistogram
void fillStatusHistogram(MonitorElement *, int xbin, int ybin, float val)
Definition: TrackingQualityChecker.cc:379
MonitorElement::kind
Kind kind(void) const
Get the type of the monitor element.
Definition: MonitorElement.h:100
DQMStore::getContents
std::vector< MonitorElement * > getContents(const std::string &path) const
Definition: DQMStore.cc:1652
TrackingQualityChecker::bookGlobalStatus
void bookGlobalStatus(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:81
MonitorElement::getQReports
std::vector< QReport * > getQReports(void) const
get map of QReports
Definition: MonitorElement.cc:566
TrackingQualityChecker::fillTrackingStatusAtLumi
void fillTrackingStatusAtLumi(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:388
TrackingQualityChecker::fillLSStatus
void fillLSStatus(DQMStore *dqm_store)
Definition: TrackingQualityChecker.cc:251
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:706
TrackingQualityChecker::TrackGlobalSummaryReportGlobal
MonitorElement * TrackGlobalSummaryReportGlobal
Definition: TrackingQualityChecker.h:65
HcalObjRepresent::Reset
void Reset(std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:620
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1522
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1000
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:928
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:434
TrackingQualityChecker::TrackingLSMEsMap
std::map< std::string, TrackingLSMEs > TrackingLSMEsMap
Definition: TrackingQualityChecker.h:62
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
DQMStore::pwd
const std::string & pwd(void) const
Definition: DQMStore.cc:556
Generated for CMSSW Reference Manual by doxygen 1.8.5