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