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