CMS 3D CMS Logo
SiPixelActionExecutor.cc
Go to the documentation of this file.
1 //#define printing false
2 //#define occupancyprinting false
3 
4 #include "DQM/SiPixelCommon/interface/SiPixelFolderOrganizer.h"
5 #include "DQM/SiPixelMonitorClient/interface/ANSIColors.h"
6 #include "DQM/SiPixelMonitorClient/interface/SiPixelActionExecutor.h"
7 #include "DQM/SiPixelMonitorClient/interface/SiPixelInformationExtractor.h"
8 #include "DQM/SiPixelMonitorClient/interface/SiPixelUtility.h"
9 #include "DQMServices/Core/interface/DQMStore.h"
10 #include "DataFormats/DetId/interface/DetId.h"
11 #include "DataFormats/TrackerCommon/interface/PixelBarrelName.h"
12 #include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
13 #include "DataFormats/TrackerCommon/interface/PixelEndcapName.h"
14 #include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
15 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
16 #include "FWCore/MessageLogger/interface/MessageLogger.h"
17 #include "FWCore/ParameterSet/interface/FileInPath.h"
18 
19 #include <cmath>
20 
21 #include <iostream>
22 using namespace std;
23 //=============================================================================================================
24 //
25 // -- Constructor
26 //
27 SiPixelActionExecutor::SiPixelActionExecutor(bool offlineXMLfile, bool Tier0Flag)
28  : offlineXMLfile_(offlineXMLfile), Tier0Flag_(Tier0Flag) {
29  edm::LogInfo("SiPixelActionExecutor") << " Creating SiPixelActionExecutor "
30  << "\n";
31  configParser_ = nullptr;
32  configWriter_ = nullptr;
33  ndet_ = 0;
34  // collationDone = false;
35 }
36 //=============================================================================================================
37 //
38 // -- Destructor
39 //
40 SiPixelActionExecutor::~SiPixelActionExecutor() {
41  edm::LogInfo("SiPixelActionExecutor") << " Deleting SiPixelActionExecutor "
42  << "\n";
43  if (configParser_)
44  delete configParser_;
45  if (configWriter_)
46  delete configWriter_;
47 }
48 //=============================================================================================================
49 //
50 // -- Read Configuration File
51 //
52 void SiPixelActionExecutor::readConfiguration() {
53  string localPath;
54  if (offlineXMLfile_)
55  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
56  else
57  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
58  if (configParser_ == nullptr) {
59  configParser_ = new SiPixelConfigParser();
60  configParser_->getDocument(edm::FileInPath(localPath).fullPath());
61  }
62 }
63 //=============================================================================================================
64 //
65 // -- Read Configuration File
66 //
67 bool SiPixelActionExecutor::readConfiguration(int &tkmap_freq,
68  int &sum_barrel_freq,
69  int &sum_endcap_freq,
70  int &sum_grandbarrel_freq,
71  int &sum_grandendcap_freq,
72  int &message_limit_,
73  int &source_type_,
74  int &calib_type_) {
75  // printing cout<<"Entering
76  // SiPixelActionExecutor::readConfiguration..."<<endl;
77  string localPath;
78  if (offlineXMLfile_)
79  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
80  else
81  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
82  if (configParser_ == nullptr) {
83  configParser_ = new SiPixelConfigParser();
84  configParser_->getDocument(edm::FileInPath(localPath).fullPath());
85  }
86 
87  if (!configParser_->getFrequencyForTrackerMap(tkmap_freq)) {
88  cout << "SiPixelActionExecutor::readConfiguration: Failed to read "
89  "TrackerMap configuration parameters!! ";
90  return false;
91  }
92  if (!configParser_->getFrequencyForBarrelSummary(sum_barrel_freq)) {
93  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Barrel Summary configuration parameters!! "
94  << "\n";
95  return false;
96  }
97  if (!configParser_->getFrequencyForEndcapSummary(sum_endcap_freq)) {
98  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Endcap Summary configuration parameters!! "
99  << "\n";
100  return false;
101  }
102  if (!configParser_->getFrequencyForGrandBarrelSummary(sum_grandbarrel_freq)) {
103  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Grand Barrel Summary configuration parameters!! "
104  << "\n";
105  return false;
106  }
107  if (!configParser_->getFrequencyForGrandEndcapSummary(sum_grandendcap_freq)) {
108  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Grand Endcap Summary configuration parameters!! "
109  << "\n";
110  return false;
111  }
112  if (!configParser_->getMessageLimitForQTests(message_limit_)) {
113  edm::LogInfo("SiPixelActionExecutor") << "Failed to read QTest Message Limit"
114  << "\n";
115  return false;
116  }
117  if (!configParser_->getSourceType(source_type_)) {
118  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Source Type"
119  << "\n";
120  return false;
121  }
122  if (!configParser_->getCalibType(calib_type_)) {
123  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Calib Type"
124  << "\n";
125  return false;
126  }
127  // printing cout<<"...leaving
128  // SiPixelActionExecutor::readConfiguration..."<<endl;
129  return true;
130 }
131 //=============================================================================================================
132 bool SiPixelActionExecutor::readConfiguration(int &tkmap_freq, int &summary_freq) {
133  // printing cout<<"Entering
134  // SiPixelActionExecutor::readConfiguration..."<<endl;
135  string localPath;
136  if (offlineXMLfile_)
137  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
138  else
139  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
140  if (configParser_ == nullptr) {
141  configParser_ = new SiPixelConfigParser();
142  configParser_->getDocument(edm::FileInPath(localPath).fullPath());
143  }
144 
145  if (!configParser_->getFrequencyForTrackerMap(tkmap_freq)) {
146  cout << "SiPixelActionExecutor::readConfiguration: Failed to read "
147  "TrackerMap configuration parameters!! ";
148  return false;
149  }
150  if (!configParser_->getFrequencyForBarrelSummary(summary_freq)) {
151  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Summary configuration parameters!! "
152  << "\n";
153  return false;
154  }
155  // printing cout<<"...leaving
156  // SiPixelActionExecutor::readConfiguration..."<<endl;
157  return true;
158 }
159 
160 //=============================================================================================================
161 void SiPixelActionExecutor::createSummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isUpgrade) {
162  // To be majorly overhauled and split into two, I guess.
163 
164  // cout<<"entering SiPixelActionExecutor::createSummary..."<<endl;
165  string barrel_structure_name;
166  vector<string> barrel_me_names;
167  string localPath;
168  if (offlineXMLfile_)
169  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
170  else
171  localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
172  // cout<<"*********************ATTENTION! LOCALPATH= "<<localPath<<endl;
173  if (configParser_ == nullptr) {
174  configParser_ = new SiPixelConfigParser();
175  configParser_->getDocument(edm::FileInPath(localPath).fullPath());
176  }
177  if (!configParser_->getMENamesForBarrelSummary(barrel_structure_name, barrel_me_names)) {
178  cout << "SiPixelActionExecutor::createSummary: Failed to read Barrel "
179  "Summary configuration parameters!! ";
180  return;
181  }
182  configParser_->getSourceType(source_type_);
183 
184  iBooker.setCurrentFolder("Pixel/");
185  iGetter.setCurrentFolder("Pixel/");
186  fillSummary(iBooker, iGetter, barrel_structure_name, barrel_me_names, true,
187  isUpgrade); // Barrel
188  iBooker.setCurrentFolder("Pixel/");
189  iGetter.setCurrentFolder("Pixel/");
190  string endcap_structure_name;
191  vector<string> endcap_me_names;
192  if (!configParser_->getMENamesForEndcapSummary(endcap_structure_name, endcap_me_names)) {
193  edm::LogInfo("SiPixelActionExecutor") << "Failed to read Endcap Summary configuration parameters!! "
194  << "\n";
195  return;
196  }
197 
198  // cout << "--- Processing endcap" << endl;
199 
200  iBooker.setCurrentFolder("Pixel/");
201  iGetter.setCurrentFolder("Pixel/");
202 
203  fillSummary(iBooker, iGetter, endcap_structure_name, endcap_me_names, false,
204  isUpgrade); // Endcap
205  iBooker.setCurrentFolder("Pixel/");
206  iGetter.setCurrentFolder("Pixel/");
207 
208  if (source_type_ == 0 || source_type_ == 5 || source_type_ == 20) { // do this only if RawData source is present
209  string federror_structure_name;
210  vector<string> federror_me_names;
211  if (!configParser_->getMENamesForFEDErrorSummary(federror_structure_name, federror_me_names)) {
212  cout << "SiPixelActionExecutor::createSummary: Failed to read FED Error "
213  "Summary configuration parameters!! ";
214  return;
215  }
216  iBooker.setCurrentFolder("Pixel/");
217  iGetter.setCurrentFolder("Pixel/");
218 
219  fillFEDErrorSummary(iBooker, iGetter, federror_structure_name, federror_me_names);
220  iBooker.setCurrentFolder("Pixel/");
221  iGetter.setCurrentFolder("Pixel/");
222  }
223  if (configWriter_)
224  delete configWriter_;
225  configWriter_ = nullptr;
226  // cout<<"leaving SiPixelActionExecutor::createSummary..."<<endl;
227 }
228 
229 //=============================================================================================================
230 void SiPixelActionExecutor::bookDeviations(DQMStore::IBooker &iBooker, bool isUpgrade) {
231  int nBPixModules;
232  if (isUpgrade) {
233  nBPixModules = 1184;
234  } else {
235  nBPixModules = 768;
236  }
237 
238  iBooker.cd();
239  iBooker.setCurrentFolder("Pixel/Barrel");
240  DEV_adc_Barrel = iBooker.book1D(
241  "DEV_adc_Barrel", "Deviation from reference;Module;<adc_ref>-<adc>", nBPixModules, 0., nBPixModules);
242  DEV_ndigis_Barrel = iBooker.book1D(
243  "DEV_ndigis_Barrel", "Deviation from reference;Module;<ndigis_ref>-<ndigis>", nBPixModules, 0., nBPixModules);
244  DEV_charge_Barrel = iBooker.book1D(
245  "DEV_charge_Barrel", "Deviation from reference;Module;<charge_ref>-<charge>", nBPixModules, 0., nBPixModules);
246  DEV_nclusters_Barrel = iBooker.book1D("DEV_nclusters_Barrel",
247  "Deviation from reference;Module;<nclusters_ref>-<nclusters>",
248  nBPixModules,
249  0.,
250  nBPixModules);
251  DEV_size_Barrel = iBooker.book1D(
252  "DEV_size_Barrel", "Deviation from reference;Module;<size_ref>-<size>", nBPixModules, 0., nBPixModules);
253  iBooker.cd();
254  iBooker.setCurrentFolder("Pixel/Endcap");
255  DEV_adc_Endcap = iBooker.book1D("DEV_adc_Endcap", "Deviation from reference;Module;<adc_ref>-<adc>", 672, 0., 672.);
256  DEV_ndigis_Endcap =
257  iBooker.book1D("DEV_ndigis_Endcap", "Deviation from reference;Module;<ndigis_ref>-<ndigis>", 672, 0., 672.);
258  DEV_charge_Endcap =
259  iBooker.book1D("DEV_charge_Endcap", "Deviation from reference;Module;<charge_ref>-<charge>", 672, 0., 672.);
260  DEV_nclusters_Endcap = iBooker.book1D(
261  "DEV_nclusters_Endcap", "Deviation from reference;Module;<nclusters_ref>-<nclusters>", 672, 0., 672.);
262  DEV_size_Endcap =
263  iBooker.book1D("DEV_size_Endcap", "Deviation from reference;Module;<size_ref>-<size>", 672, 0., 672.);
264  iBooker.cd();
265 }
266 
267 //=============================================================================================================
268 
269 void SiPixelActionExecutor::GetBladeSubdirs(DQMStore::IBooker &iBooker,
270  DQMStore::IGetter &iGetter,
271  vector<string> &blade_subdirs) {
272  blade_subdirs.clear();
273  vector<string> panels = iGetter.getSubdirs();
274  vector<string> modules;
275  for (vector<string>::const_iterator it = panels.begin(); it != panels.end(); it++) {
276  iGetter.cd(*it);
277  iBooker.cd(*it);
278  modules = iGetter.getSubdirs();
279  for (vector<string>::const_iterator m_it = modules.begin(); m_it != modules.end(); m_it++) {
280  blade_subdirs.push_back(*m_it);
281  }
282  }
283 }
284 
285 //=============================================================================================================
286 
287 void SiPixelActionExecutor::fillSummary(DQMStore::IBooker &iBooker,
288  DQMStore::IGetter &iGetter,
289  string dir_name,
290  vector<string> &me_names,
291  bool isbarrel,
292  bool isUpgrade) {
293  // cout<<"entering SiPixelActionExecutor::fillSummary..."<<endl;
294  string currDir = iBooker.pwd();
295  string prefix;
296  if (source_type_ == 0)
297  prefix = "SUMRAW";
298  else if (source_type_ == 1)
299  prefix = "SUMDIG";
300  else if (source_type_ == 2)
301  prefix = "SUMCLU";
302  else if (source_type_ == 3)
303  prefix = "SUMTRK";
304  else if (source_type_ == 4)
305  prefix = "SUMHIT";
306  else if (source_type_ >= 7 && source_type_ < 20)
307  prefix = "SUMCAL";
308  else if (source_type_ == 20)
309  prefix = "SUMOFF";
310  if (currDir.find(dir_name) != string::npos) {
311  vector<MonitorElement *> sum_mes;
312  for (vector<string>::const_iterator iv = me_names.begin(); iv != me_names.end(); iv++) {
313  if (source_type_ == 5 || source_type_ == 6) {
314  if ((*iv) == "errorType" || (*iv) == "NErrors" || (*iv) == "fullType" || (*iv) == "chanNmbr" ||
315  (*iv) == "TBMType" || (*iv) == "EvtNbr" || (*iv) == "evtSize" || (*iv) == "linkId" || (*iv) == "ROCId" ||
316  (*iv) == "DCOLId" || (*iv) == "PXId" || (*iv) == "ROCNmbr" || (*iv) == "TBMMessage" ||
317  (*iv) == "Type36Hitmap")
318  prefix = "SUMRAW";
319  else if ((*iv) == "ndigis" || (*iv) == "adc")
320  prefix = "SUMDIG";
321  else if ((*iv) == "nclusters" || (*iv) == "x" || (*iv) == "y" || (*iv) == "charge" || (*iv) == "size" ||
322  (*iv) == "sizeX" || (*iv) == "sizeY" || (*iv) == "minrow" || (*iv) == "maxrow" || (*iv) == "mincol" ||
323  (*iv) == "maxcol")
324  prefix = "SUMCLU";
325  if (currDir.find("Track") != string::npos)
326  prefix = "SUMTRK";
327  else if ((*iv) == "residualX" || (*iv) == "residualY")
328  prefix = "SUMTRK";
329  else if ((*iv) == "ClustX" || (*iv) == "ClustY" || (*iv) == "nRecHits" || (*iv) == "ErrorX" ||
330  (*iv) == "ErrorY")
331  prefix = "SUMHIT";
332  else if ((*iv) == "Gain1d" || (*iv) == "GainChi2NDF1d" || (*iv) == "GainChi2Prob1d" || (*iv) == "Pedestal1d" ||
333  (*iv) == "GainNPoints1d" || (*iv) == "GainHighPoint1d" || (*iv) == "GainLowPoint1d" ||
334  (*iv) == "GainEndPoint1d" || (*iv) == "GainFitResult2d" || (*iv) == "GainDynamicRange2d" ||
335  (*iv) == "GainSaturate2d" || (*iv) == "ScurveChi2ProbSummary" || (*iv) == "ScurveFitResultSummary" ||
336  (*iv) == "ScurveSigmasSummary" || (*iv) == "ScurveThresholdSummary" || (*iv) == "pixelAliveSummary" ||
337  (*iv) == "SiPixelErrorsCalibDigis")
338  prefix = "SUMCAL";
339  }
340  MonitorElement *temp;
341  string tag;
342  if ((*iv).find("residual") != string::npos) { // track residuals
343  tag = prefix + "_" + (*iv) + "_mean_" + currDir.substr(currDir.find(dir_name));
344  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
345  sum_mes.push_back(temp);
346  tag = prefix + "_" + (*iv) + "_RMS_" + currDir.substr(currDir.find(dir_name));
347  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
348  sum_mes.push_back(temp);
349  } else if (prefix == "SUMCAL") { // calibrations
350  if ((*iv) == "Gain1d" || (*iv) == "GainChi2NDF1d" || (*iv) == "GainChi2Prob1d" || (*iv) == "GainNPoints1d" ||
351  (*iv) == "GainHighPoint1d" || (*iv) == "GainLowPoint1d" || (*iv) == "GainEndPoint1d" ||
352  (*iv) == "GainDynamicRange2d" || (*iv) == "GainSaturate2d" || (*iv) == "Pedestal1d" ||
353  (*iv) == "ScurveChi2ProbSummary" || (*iv) == "ScurveFitResultSummary" || (*iv) == "ScurveSigmasSummary" ||
354  (*iv) == "ScurveThresholdSummary") {
355  tag = prefix + "_" + (*iv) + "_mean_" + currDir.substr(currDir.find(dir_name));
356  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
357  sum_mes.push_back(temp);
358  tag = prefix + "_" + (*iv) + "_RMS_" + currDir.substr(currDir.find(dir_name));
359  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
360  sum_mes.push_back(temp);
361  } else if ((*iv) == "SiPixelErrorsCalibDigis") {
362  tag = prefix + "_" + (*iv) + "_NCalibErrors_" + currDir.substr(currDir.find(dir_name));
363  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
364  sum_mes.push_back(temp);
365  } else if ((*iv) == "GainFitResult2d") {
366  tag = prefix + "_" + (*iv) + "_NNegativeFits_" + currDir.substr(currDir.find(dir_name));
367  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
368  sum_mes.push_back(temp);
369  } else if ((*iv) == "pixelAliveSummary") {
370  tag = prefix + "_" + (*iv) + "_FracOfPerfectPix_" + currDir.substr(currDir.find(dir_name));
371  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
372  sum_mes.push_back(temp);
373  tag = prefix + "_" + (*iv) + "_mean_" + currDir.substr(currDir.find(dir_name));
374  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
375  sum_mes.push_back(temp);
376  }
377  } else {
378  tag = prefix + "_" + (*iv) + "_" + currDir.substr(currDir.find(dir_name));
379  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
380  sum_mes.push_back(temp);
381  if ((*iv) == "ndigis") {
382  tag = prefix + "_" + (*iv) + "FREQ_" + currDir.substr(currDir.find(dir_name));
383  temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
384  sum_mes.push_back(temp);
385  }
386  if (prefix == "SUMDIG" && (*iv) == "adc") {
387  tag = "ALLMODS_" + (*iv) + "COMB_" + currDir.substr(currDir.find(dir_name));
388  temp = nullptr;
389  string fullpathname = iBooker.pwd() + "/" + tag;
390  temp = iGetter.get(fullpathname);
391  if (temp) {
392  temp->Reset();
393  } else {
394  temp = iBooker.book1D(tag.c_str(), tag.c_str(), 128, 0., 256.);
395  }
396  sum_mes.push_back(temp);
397  }
398  if (prefix == "SUMCLU" && (*iv) == "charge") {
399  tag = "ALLMODS_" + (*iv) + "COMB_" + currDir.substr(currDir.find(dir_name));
400  temp = nullptr;
401  string fullpathname = iBooker.pwd() + "/" + tag;
402  temp = iGetter.get(fullpathname);
403  if (temp) {
404  temp->Reset();
405  } else {
406  temp = iBooker.book1D(tag.c_str(), tag.c_str(), 100, 0.,
407  200.); // To look to get the size automatically
408  }
409  sum_mes.push_back(temp);
410  }
411  }
412  }
413  if (sum_mes.empty()) {
414  edm::LogInfo("SiPixelActionExecutor") << " Summary MEs can not be created"
415  << "\n";
416  return;
417  }
418  vector<string> subdirs = iGetter.getSubdirs();
419  // Blade
420  if (dir_name.find("Blade_") == 0)
421  GetBladeSubdirs(iBooker, iGetter, subdirs);
422 
423  int ndet = 0;
424  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
425  if (prefix != "SUMOFF" && (*it).find("Module_") == string::npos)
426  continue;
427  if (prefix == "SUMOFF" && (*it).find(isbarrel ? "Layer_" : "Disk_") == string::npos)
428  continue;
429  iBooker.cd(*it);
430  iGetter.cd(*it);
431  ndet++;
432 
433  vector<string> contents = iGetter.getMEs();
434 
435  for (vector<MonitorElement *>::const_iterator isum = sum_mes.begin(); isum != sum_mes.end(); isum++) {
436  for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++) {
437  string sname = ((*isum)->getName());
438  string tname = " ";
439  tname = sname.substr(7, (sname.find('_', 7) - 6));
440  if (sname.find("ALLMODS_adcCOMB_") != string::npos)
441  tname = "adc_";
442  if (sname.find("ALLMODS_chargeCOMB_") != string::npos)
443  tname = "charge_";
444  if (sname.find("_charge_") != string::npos && sname.find("Track_") == string::npos)
445  tname = "charge_";
446  if (sname.find("_nclusters_") != string::npos && sname.find("Track_") == string::npos)
447  tname = "nclusters_";
448  if (sname.find("_size_") != string::npos && sname.find("Track_") == string::npos)
449  tname = "size_";
450  if (sname.find("_charge_OffTrack_") != string::npos)
451  tname = "charge_OffTrack_";
452  if (sname.find("_nclusters_OffTrack_") != string::npos)
453  tname = "nclusters_OffTrack_";
454  if (sname.find("_size_OffTrack_") != string::npos)
455  tname = "size_OffTrack_";
456  if (sname.find("_sizeX_OffTrack_") != string::npos)
457  tname = "sizeX_OffTrack_";
458  if (sname.find("_sizeY_OffTrack_") != string::npos)
459  tname = "sizeY_OffTrack_";
460  if (sname.find("_charge_OnTrack_") != string::npos)
461  tname = "charge_OnTrack_";
462  if (sname.find("_nclusters_OnTrack_") != string::npos)
463  tname = "nclusters_OnTrack_";
464  if (sname.find("_size_OnTrack_") != string::npos)
465  tname = "size_OnTrack_";
466  if (sname.find("_sizeX_OnTrack_") != string::npos)
467  tname = "sizeX_OnTrack_";
468  if (sname.find("_sizeY_OnTrack_") != string::npos)
469  tname = "sizeY_OnTrack_";
470  if (tname.find("FREQ") != string::npos)
471  tname = "ndigis_";
472  if (((*im)).find(tname) == 0) {
473  string fullpathname = iBooker.pwd() + "/" + (*im);
474  MonitorElement *me = iGetter.get(fullpathname);
475 
476  if (me) {
477  if (sname.find("_charge") != string::npos && sname.find("Track_") == string::npos &&
478  me->getName().find("Track_") != string::npos)
479  continue;
480  if (sname.find("_nclusters_") != string::npos && sname.find("Track_") == string::npos &&
481  me->getName().find("Track_") != string::npos)
482  continue;
483  if (sname.find("_size") != string::npos && sname.find("Track_") == string::npos &&
484  me->getName().find("Track_") != string::npos)
485  continue;
486  // fill summary histos:
487  if (sname.find("_RMS_") != string::npos && sname.find("GainDynamicRange2d") == string::npos &&
488  sname.find("GainSaturate2d") == string::npos) {
489  (*isum)->Fill(ndet, me->getRMS());
490  } else if (sname.find("GainDynamicRange2d") != string::npos ||
491  sname.find("GainSaturate2d") != string::npos) {
492  float SumOfEntries = 0.;
493  float SumOfSquaredEntries = 0.;
494  int SumOfPixels = 0;
495  for (int cols = 1; cols != me->getNbinsX() + 1; cols++)
496  for (int rows = 1; rows != me->getNbinsY() + 1; rows++) {
497  SumOfEntries += me->getBinContent(cols, rows);
498  SumOfSquaredEntries += (me->getBinContent(cols, rows)) * (me->getBinContent(cols, rows));
499  SumOfPixels++;
500  }
501 
502  float MeanInZ = SumOfEntries / float(SumOfPixels);
503  float RMSInZ = sqrt(SumOfSquaredEntries / float(SumOfPixels));
504  if (sname.find("_mean_") != string::npos)
505  (*isum)->Fill(ndet, MeanInZ);
506  if (sname.find("_RMS_") != string::npos)
507  (*isum)->Fill(ndet, RMSInZ);
508  } else if (sname.find("_FracOfPerfectPix_") != string::npos) {
509  float nlast = me->getBinContent(me->getNbinsX());
510  float nall = (me->getTH1F())->Integral(1, 11);
511  (*isum)->Fill(ndet, nlast / nall);
512  } else if (sname.find("_NCalibErrors_") != string::npos || sname.find("FREQ_") != string::npos) {
513  float nall = me->getEntries();
514  (*isum)->Fill(ndet, nall);
515  } else if (sname.find("GainFitResult2d") != string::npos) {
516  int NegFitPixels = 0;
517  for (int cols = 1; cols != me->getNbinsX() + 1; cols++)
518  for (int rows = 1; rows != me->getNbinsY() + 1; rows++) {
519  if (me->getBinContent(cols, rows) < 0.)
520  NegFitPixels++;
521  }
522  (*isum)->Fill(ndet, float(NegFitPixels));
523  } else if (sname.find("ALLMODS_adcCOMB_") != string::npos ||
524  (sname.find("ALLMODS_chargeCOMB_") != string::npos &&
525  me->getName().find("Track_") == string::npos)) {
526  (*isum)->getTH1F()->Add(me->getTH1F());
527  } else if (sname.find("_NErrors_") != string::npos) {
528  string path1 = fullpathname;
529  path1 = path1.replace(path1.find("NErrors"), 7, "errorType");
530  MonitorElement *me1 = iGetter.get(path1);
531  bool notReset = true;
532  if (me1) {
533  for (int jj = 1; jj < 16; jj++) {
534  if (me1->getBinContent(jj) > 0.) {
535  if (jj == 6) { // errorType=30 (reset)
536  string path2 = path1;
537  path2 = path2.replace(path2.find("errorType"), 9, "TBMMessage");
538  MonitorElement *me2 = iGetter.get(path2);
539  if (me2)
540  if (me2->getBinContent(6) > 0. || me2->getBinContent(7) > 0.)
541  notReset = false;
542  }
543  }
544  }
545  }
546  if (notReset)
547  (*isum)->Fill(ndet, me1->getEntries());
548  } else if ((sname.find("_charge_") != string::npos && sname.find("Track_") == string::npos &&
549  me->getName().find("Track_") == string::npos) ||
550  (sname.find("_charge_") != string::npos && sname.find("_OnTrack_") != string::npos &&
551  me->getName().find("_OnTrack_") != string::npos) ||
552  (sname.find("_charge_") != string::npos && sname.find("_OffTrack_") != string::npos &&
553  me->getName().find("_OffTrack_") != string::npos) ||
554  (sname.find("_nclusters_") != string::npos && sname.find("Track_") == string::npos &&
555  me->getName().find("Track_") == string::npos) ||
556  (sname.find("_nclusters_") != string::npos && sname.find("_OnTrack_") != string::npos &&
557  me->getName().find("_OnTrack_") != string::npos) ||
558  (sname.find("_nclusters_") != string::npos && sname.find("_OffTrack_") != string::npos &&
559  me->getName().find("_OffTrack_") != string::npos) ||
560  (sname.find("_size") != string::npos && sname.find("Track_") == string::npos &&
561  me->getName().find("Track_") == string::npos) ||
562  (sname.find("_size") != string::npos && sname.find("_OnTrack_") != string::npos &&
563  me->getName().find("_OnTrack_") != string::npos) ||
564  (sname.find("_size") != string::npos && sname.find("_OffTrack_") != string::npos &&
565  me->getName().find("_OffTrack_") != string::npos)) {
566  (*isum)->Fill(ndet, me->getMean());
567  } else if (sname.find("_charge_") == string::npos && sname.find("_nclusters_") == string::npos &&
568  sname.find("_size") == string::npos) {
569  (*isum)->Fill(ndet, me->getMean());
570  }
571 
572  // set titles:
573  if (prefix == "SUMOFF") {
574  (*isum)->setAxisTitle(isbarrel ? "Ladders" : "Blades", 1);
575  } else if (sname.find("ALLMODS_adcCOMB_") != string::npos) {
576  (*isum)->setAxisTitle("Digi charge [ADC]", 1);
577  } else if (sname.find("ALLMODS_chargeCOMB_") != string::npos) {
578  (*isum)->setAxisTitle("Cluster charge [kilo electrons]", 1);
579  } else {
580  (*isum)->setAxisTitle("Modules", 1);
581  }
582  string title = " ";
583  if (sname.find("_RMS_") != string::npos) {
584  title = "RMS of " + sname.substr(7, (sname.find('_', 7) - 7)) + " per module";
585  } else if (sname.find("_FracOfPerfectPix_") != string::npos) {
586  title = "FracOfPerfectPix " + sname.substr(7, (sname.find('_', 7) - 7)) + " per module";
587  } else if (sname.find("_NCalibErrors_") != string::npos) {
588  title = "Number of CalibErrors " + sname.substr(7, (sname.find('_', 7) - 7)) + " per module";
589  } else if (sname.find("_NNegativeFits_") != string::npos) {
590  title = "Number of pixels with neg. fit result " + sname.substr(7, (sname.find('_', 7) - 7)) +
591  " per module";
592  } else if (sname.find("FREQ_") != string::npos) {
593  title = "NEvents with digis per module";
594  } else if (sname.find("ALLMODS_adcCOMB_") != string::npos) {
595  title = "NDigis";
596  } else if (sname.find("ALLMODS_chargeCOMB_") != string::npos) {
597  title = "NClusters";
598  } else if (sname.find("_NErrors_") != string::npos) {
599  if (prefix == "SUMOFF" && isbarrel)
600  title = "Total number of errors per Ladder";
601  else if (prefix == "SUMOFF" && !isbarrel)
602  title = "Total number of errors per Blade";
603  else
604  title = "Total number of errors per Module";
605  } else {
606  if (prefix == "SUMOFF")
607  title =
608  "Mean " + sname.substr(7, (sname.find('_', 7) - 7)) + (isbarrel ? " per Ladder" : " per Blade");
609  else
610  title = "Mean " + sname.substr(7, (sname.find('_', 7) - 7)) + " per Module";
611  }
612  (*isum)->setAxisTitle(title, 2);
613  }
614  break;
615  }
616  }
617  }
618  iBooker.goUp();
619  iGetter.setCurrentFolder(iBooker.pwd());
620  if (dir_name.find("Blade") == 0) {
621  iBooker.goUp(); // Going up a second time if we are processing the Blade
622  iGetter.setCurrentFolder(iBooker.pwd());
623  }
624  } // end for it (subdirs)
625  } else {
626  vector<string> subdirs = iGetter.getSubdirs();
627  // printing cout << "#\t" << iBooker.pwd() << endl;
628  if (isbarrel) {
629  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
630  // cout << "##\t" << iBooker.pwd() << "\t"
631  //<<
632  //(*it) << endl;
633  if ((iBooker.pwd()).find("Endcap") != string::npos ||
634  (iBooker.pwd()).find("AdditionalPixelErrors") != string::npos) {
635  iBooker.goUp();
636  iGetter.setCurrentFolder(iBooker.pwd());
637  }
638  iBooker.cd(*it);
639  iGetter.cd(*it);
640  if ((*it).find("Endcap") != string::npos || (*it).find("AdditionalPixelErrors") != string::npos)
641  continue;
642  fillSummary(iBooker, iGetter, dir_name, me_names, true,
643  isUpgrade); // Barrel
644  iBooker.goUp();
645  iGetter.setCurrentFolder(iBooker.pwd());
646  }
647  string grandbarrel_structure_name;
648  vector<string> grandbarrel_me_names;
649  if (!configParser_->getMENamesForGrandBarrelSummary(grandbarrel_structure_name, grandbarrel_me_names)) {
650  cout << "SiPixelActionExecutor::createSummary: Failed to read Grand "
651  "Barrel Summary configuration parameters!! ";
652  return;
653  }
654  fillGrandBarrelSummaryHistos(iBooker, iGetter, grandbarrel_me_names, isUpgrade);
655 
656  } else // Endcap
657  {
658  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
659  if ((iBooker.pwd()).find("Barrel") != string::npos ||
660  (iBooker.pwd()).find("AdditionalPixelErrors") != string::npos) {
661  iBooker.goUp();
662  iGetter.setCurrentFolder(iBooker.pwd());
663  }
664  iBooker.cd(*it);
665  iGetter.cd(*it);
666  if ((*it).find("Barrel") != string::npos || (*it).find("AdditionalPixelErrors") != string::npos)
667  continue;
668  fillSummary(iBooker, iGetter, dir_name, me_names, false,
669  isUpgrade); // Endcap
670  iBooker.goUp();
671  iGetter.setCurrentFolder(iBooker.pwd());
672  }
673  string grandendcap_structure_name;
674  vector<string> grandendcap_me_names;
675  if (!configParser_->getMENamesForGrandEndcapSummary(grandendcap_structure_name, grandendcap_me_names)) {
676  cout << "SiPixelActionExecutor::createSummary: Failed to read Grand "
677  "Endcap Summary configuration parameters!! ";
678  return;
679  }
680  fillGrandEndcapSummaryHistos(iBooker, iGetter, grandendcap_me_names, isUpgrade);
681  }
682  }
683  // cout<<"...leaving SiPixelActionExecutor::fillSummary!"<<endl;
684 }
685 
686 //=============================================================================================================
687 void SiPixelActionExecutor::fillFEDErrorSummary(DQMStore::IBooker &iBooker,
688  DQMStore::IGetter &iGetter,
689  string dir_name,
690  vector<string> &me_names) {
691  // printing cout<<"entering
692  // SiPixelActionExecutor::fillFEDErrorSummary..."<<endl;
693  string currDir = iBooker.pwd();
694  string prefix;
695  if (source_type_ == 0)
696  prefix = "SUMRAW";
697  else if (source_type_ == 20)
698  prefix = "SUMOFF";
699 
700  if (currDir.find(dir_name) != string::npos) {
701  vector<MonitorElement *> sum_mes;
702  for (vector<string>::const_iterator iv = me_names.begin(); iv != me_names.end(); iv++) {
703  bool isBooked = false;
704  vector<string> contents = iGetter.getMEs();
705  for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++)
706  if ((*im).find(*iv) != string::npos)
707  isBooked = true;
708  if (source_type_ == 5 || source_type_ == 6) {
709  if ((*iv) == "errorType" || (*iv) == "NErrors" || (*iv) == "fullType" || (*iv) == "chanNmbr" ||
710  (*iv) == "TBMType" || (*iv) == "EvtNbr" || (*iv) == "evtSize" || (*iv) == "linkId" || (*iv) == "ROCId" ||
711  (*iv) == "DCOLId" || (*iv) == "PXId" || (*iv) == "ROCNmbr" || (*iv) == "TBMMessage" ||
712  (*iv) == "Type36Hitmap" || (*iv) == "FedChLErr" || (*iv) == "FedChNErr" || (*iv) == "FedETypeNErr")
713  prefix = "SUMRAW";
714  }
715  if ((*iv) == "errorType" || (*iv) == "NErrors" || (*iv) == "fullType" || (*iv) == "chanNmbr" ||
716  (*iv) == "TBMType" || (*iv) == "EvtNbr" || (*iv) == "evtSize" || (*iv) == "linkId" || (*iv) == "ROCId" ||
717  (*iv) == "DCOLId" || (*iv) == "PXId" || (*iv) == "ROCNmbr" || (*iv) == "TBMMessage" ||
718  (*iv) == "Type36Hitmap") {
719  string tag = prefix + "_" + (*iv) + "_FEDErrors";
720  MonitorElement *temp = getFEDSummaryME(iBooker, iGetter, tag);
721  sum_mes.push_back(temp);
722  } else if ((*iv) == "FedChLErr" || (*iv) == "FedChNErr" || (*iv) == "FedETypeNErr") {
723  string tag = prefix + "_" + (*iv);
724  MonitorElement *temp;
725  if ((*iv) == "FedChLErr") {
726  if (!isBooked)
727  temp = iBooker.book2D("FedChLErr", "Type of last error", 40, -0.5, 39.5, 37, 0., 37.);
728  else {
729  string fullpathname = iBooker.pwd() + "/" + (*iv);
730  temp = iGetter.get(fullpathname);
731  temp->Reset();
732  }
733  } // If I don't reset this one, then I instead start adding error
734  // codes..
735  if ((*iv) == "FedChNErr") {
736  if (!isBooked)
737  temp = iBooker.book2D("FedChNErr", "Total number of errors", 40, -0.5, 39.5, 37, 0., 37.);
738  else {
739  string fullpathname = iBooker.pwd() + "/" + (*iv);
740  temp = iGetter.get(fullpathname);
741  temp->Reset();
742  }
743  } // If I don't reset this one, then I instead start adding error
744  // codes..
745  if ((*iv) == "FedETypeNErr") {
746  if (!isBooked) {
747  temp = iBooker.book2D("FedETypeNErr", "Number of each error type", 40, -0.5, 39.5, 21, 0., 21.);
748  temp->setBinLabel(1, "ROC of 25", 2);
749  temp->setBinLabel(2, "Gap word", 2);
750  temp->setBinLabel(3, "Dummy word", 2);
751  temp->setBinLabel(4, "FIFO full", 2);
752  temp->setBinLabel(5, "Timeout", 2);
753  temp->setBinLabel(6, "Stack full", 2);
754  temp->setBinLabel(7, "Pre-cal issued", 2);
755  temp->setBinLabel(8, "Trigger clear or sync", 2);
756  temp->setBinLabel(9, "No token bit", 2);
757  temp->setBinLabel(10, "Overflow", 2);
758  temp->setBinLabel(11, "FSM error", 2);
759  temp->setBinLabel(12, "Invalid #ROCs", 2);
760  temp->setBinLabel(13, "Event number", 2);
761  temp->setBinLabel(14, "Slink header", 2);
762  temp->setBinLabel(15, "Slink trailer", 2);
763  temp->setBinLabel(16, "Event size", 2);
764  temp->setBinLabel(17, "Invalid channel#", 2);
765  temp->setBinLabel(18, "ROC value", 2);
766  temp->setBinLabel(19, "Dcol or pixel value", 2);
767  temp->setBinLabel(20, "Readout order", 2);
768  temp->setBinLabel(21, "CRC error", 2);
769  } else {
770  string fullpathname = iBooker.pwd() + "/" + (*iv);
771  temp = iGetter.get(fullpathname);
772  temp->Reset();
773  } // If I don't reset this one, then I instead start adding error
774  // codes..
775  }
776  sum_mes.push_back(temp);
777  }
778  }
779  if (sum_mes.empty()) {
780  edm::LogInfo("SiPixelActionExecutor") << " Summary MEs can not be created"
781  << "\n";
782  return;
783  }
784  vector<string> subdirs = iGetter.getSubdirs();
785  int ndet = 0;
786  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
787  if ((*it).find("FED_") == string::npos)
788  continue;
789  iBooker.cd(*it);
790  iGetter.cd(*it);
791  string fedid = (*it).substr((*it).find("_") + 1);
792  std::istringstream isst;
793  isst.str(fedid);
794  isst >> ndet;
795  ndet++;
796  vector<string> contents = iGetter.getMEs();
797 
798  for (vector<MonitorElement *>::const_iterator isum = sum_mes.begin(); isum != sum_mes.end(); isum++) {
799  for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++) {
800  if (((*im).find("FedChNErr") != std::string::npos &&
801  (*isum)->getName().find("FedChNErr") != std::string::npos) ||
802  ((*im).find("FedChLErr") != std::string::npos &&
803  (*isum)->getName().find("FedChLErr") != std::string::npos) ||
804  ((*im).find("FedETypeNErr") != std::string::npos &&
805  (*isum)->getName().find("FedETypeNErr") != std::string::npos)) {
806  string fullpathname = iBooker.pwd() + "/" + (*im);
807  MonitorElement *me = iGetter.get(fullpathname);
808  if (me) {
809  for (int i = 0; i != 37; i++) {
810  if ((*im).find("FedETypeNErr") != std::string::npos && i < 21)
811  (*isum)->Fill(ndet - 1, i, me->getBinContent(i + 1));
812  else
813  (*isum)->Fill(ndet - 1, i, me->getBinContent(i + 1));
814  }
815  }
816  }
817  string sname = ((*isum)->getName());
818  string tname = " ";
819  tname = sname.substr(7, (sname.find('_', 7) - 6));
820  if (((*im)).find(tname) == 0) {
821  string fullpathname = iBooker.pwd() + "/" + (*im);
822  MonitorElement *me = iGetter.get(fullpathname);
823 
824  if (me) {
825  if (me->getMean() > 0.) {
826  if (sname.find("_NErrors_") != string::npos) {
827  string path1 = fullpathname;
828  path1 = path1.replace(path1.find("NErrors"), 7, "errorType");
829  MonitorElement *me1 = iGetter.get(path1);
830  bool notReset = true;
831  if (me1) {
832  for (int jj = 1; jj < 16; jj++) {
833  if (me1->getBinContent(jj) > 0.) {
834  if (jj == 6) { // errorType=30 (reset)
835  string path2 = path1;
836  path2 = path2.replace(path2.find("errorType"), 9, "TBMMessage");
837  MonitorElement *me2 = iGetter.get(path2);
838  if (me2)
839  if (me2->getBinContent(6) > 0. || me2->getBinContent(7) > 0.)
840  notReset = false;
841  }
842  }
843  }
844  }
845  if (notReset)
846  (*isum)->setBinContent(ndet, (*isum)->getBinContent(ndet) + me1->getEntries());
847  } else
848  (*isum)->setBinContent(ndet, (*isum)->getBinContent(ndet) + me->getEntries());
849  }
850  (*isum)->setAxisTitle("FED #", 1);
851  string title = " ";
852  title = sname.substr(7, (sname.find('_', 7) - 7)) + " per FED";
853  (*isum)->setAxisTitle(title, 2);
854  }
855  break;
856  }
857  }
858  }
859  iBooker.goUp();
860  iGetter.setCurrentFolder(iBooker.pwd());
861  }
862  } else {
863  vector<string> subdirs = iGetter.getSubdirs();
864  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
865  if ((*it).find("Endcap") != string::npos || (*it).find("Barrel") != string::npos)
866  continue;
867  iBooker.cd(*it);
868  iGetter.cd(*it);
869  fillFEDErrorSummary(iBooker, iGetter, dir_name, me_names);
870  iBooker.goUp();
871  iGetter.setCurrentFolder(iBooker.pwd());
872  }
873  }
874  // printing cout<<"...leaving
875  // SiPixelActionExecutor::fillFEDErrorSummary!"<<endl;
876 }
877 
878 //=============================================================================================================
879 void SiPixelActionExecutor::fillGrandBarrelSummaryHistos(DQMStore::IBooker &iBooker,
880  DQMStore::IGetter &iGetter,
881  vector<string> &me_names,
882  bool isUpgrade) {
883  // cout<<"Entering
884  // SiPixelActionExecutor::fillGrandBarrelSummaryHistos...:"<<me_names.size()<<endl;
885  vector<MonitorElement *> gsum_mes;
886  string currDir = iBooker.pwd();
887  string path_name = iBooker.pwd();
888  string dir_name = path_name.substr(path_name.find_last_of('/') + 1);
889  if ((dir_name.find("DQMData") == 0) || (dir_name.find("Pixel") == 0) ||
890  (dir_name.find("AdditionalPixelErrors") == 0) || (dir_name.find("Endcap") == 0) ||
891  (dir_name.find("HalfCylinder") == 0) || (dir_name.find("Disk") == 0) || (dir_name.find("Blade") == 0) ||
892  (dir_name.find("Panel") == 0))
893  return;
894  vector<string> subdirs = iGetter.getSubdirs();
895  int nDirs = subdirs.size();
896  int iDir = 0;
897  int nbin = 0;
898  int nbin_i = 0;
899  int nbin_subdir = 0;
900  int cnt = 0;
901  bool first_subdir = true;
902  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
903  cnt++;
904  iBooker.cd(*it);
905  iGetter.cd(*it);
906  vector<string> contents = iGetter.getMEs();
907 
908  iBooker.goUp();
909  iGetter.setCurrentFolder(iBooker.pwd());
910 
911  string prefix;
912  if (source_type_ == 0)
913  prefix = "SUMRAW";
914  else if (source_type_ == 1)
915  prefix = "SUMDIG";
916  else if (source_type_ == 2)
917  prefix = "SUMCLU";
918  else if (source_type_ == 3)
919  prefix = "SUMTRK";
920  else if (source_type_ == 4)
921  prefix = "SUMHIT";
922  else if (source_type_ >= 7 && source_type_ < 20)
923  prefix = "SUMCAL";
924  else if (source_type_ == 20)
925  prefix = "SUMOFF";
926 
927  for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++) {
928  for (vector<string>::const_iterator iv = me_names.begin(); iv != me_names.end(); iv++) {
929  string var = "_" + (*iv) + "_";
930  if ((*im).find(var) != string::npos) {
931  if ((var == "_charge_" || var == "_nclusters_" || var == "_size_" || var == "_sizeX_" || var == "_sizeY_") &&
932  (*im).find("Track_") != string::npos)
933  continue;
934  string full_path = (*it) + "/" + (*im);
935  MonitorElement *me = iGetter.get(full_path);
936  if (!me)
937  continue;
938  if (source_type_ == 5 || source_type_ == 6) {
939  if ((*iv) == "errorType" || (*iv) == "NErrors" || (*iv) == "fullType" || (*iv) == "chanNmbr" ||
940  (*iv) == "TBMType" || (*iv) == "EvtNbr" || (*iv) == "evtSize" || (*iv) == "linkId" ||
941  (*iv) == "ROCId" || (*iv) == "DCOLId" || (*iv) == "PXId" || (*iv) == "ROCNmbr" ||
942  (*iv) == "TBMMessage" || (*iv) == "Type36Hitmap")
943  prefix = "SUMRAW";
944  else if ((*iv) == "ndigis" || (*iv) == "adc" || (*iv) == "ndigisFREQ" || (*iv) == "adcCOMB")
945  prefix = "SUMDIG";
946  else if ((*iv) == "nclusters" || (*iv) == "x" || (*iv) == "y" || (*iv) == "charge" ||
947  (*iv) == "chargeCOMB" || (*iv) == "size" || (*iv) == "sizeX" || (*iv) == "sizeY" ||
948  (*iv) == "minrow" || (*iv) == "maxrow" || (*iv) == "mincol" || (*iv) == "maxcol")
949  prefix = "SUMCLU";
950  if (currDir.find("Track") != string::npos)
951  prefix = "SUMTRK";
952  else if ((*iv) == "residualX_mean" || (*iv) == "residualY_mean" || (*iv) == "residualX_RMS" ||
953  (*iv) == "residualY_RMS")
954  prefix = "SUMTRK";
955  else if ((*iv) == "ClustX" || (*iv) == "ClustY" || (*iv) == "nRecHits" || (*iv) == "ErrorX" ||
956  (*iv) == "ErrorY")
957  prefix = "SUMHIT";
958  else if ((*iv) == "Gain1d_mean" || (*iv) == "GainChi2NDF1d_mean" || (*iv) == "GainChi2Prob1d_mean" ||
959  (*iv) == "Pedestal1d_mean" || (*iv) == "ScurveChi2ProbSummary_mean" ||
960  (*iv) == "ScurveFitResultSummary_mean" || (*iv) == "ScurveSigmasSummary_mean" ||
961  (*iv) == "ScurveThresholdSummary_mean" || (*iv) == "Gain1d_RMS" || (*iv) == "GainChi2NDF1d_RMS" ||
962  (*iv) == "GainChi2Prob1d_RMS" || (*iv) == "Pedestal1d_RMS" || (*iv) == "GainNPoints1d_mean" ||
963  (*iv) == "GainNPoints1d_RMS" || (*iv) == "GainHighPoint1d_mean" ||
964  (*iv) == "GainHighPoint1d_RMS" || (*iv) == "GainLowPoint1d_mean" ||
965  (*iv) == "GainLowPoint1d_RMS" || (*iv) == "GainEndPoint1d_mean" || (*iv) == "GainEndPoint1d_RMS" ||
966  (*iv) == "GainFitResult2d_mean" || (*iv) == "GainFitResult2d_RMS" ||
967  (*iv) == "GainDynamicRange2d_mean" || (*iv) == "GainDynamicRange2d_RMS" ||
968  (*iv) == "GainSaturate2d_mean" || (*iv) == "GainSaturate2d_RMS" ||
969  (*iv) == "ScurveChi2ProbSummary_RMS" || (*iv) == "ScurveFitResultSummary_RMS" ||
970  (*iv) == "ScurveSigmasSummary_RMS" || (*iv) == "ScurveThresholdSummary_RMS" ||
971  (*iv) == "pixelAliveSummary_mean" || (*iv) == "pixelAliveSummary_FracOfPerfectPix" ||
972  (*iv) == "SiPixelErrorsCalibDigis_NCalibErrors")
973  prefix = "SUMCAL";
974  } // end source_type if
975 
976  if (first_subdir && !isUpgrade) {
977  nbin = me->getTH1F()->GetNbinsX();
978  string me_name = prefix + "_" + (*iv) + "_" + dir_name;
979  if ((*iv) == "adcCOMB" || (*iv) == "chargeCOMB")
980  me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
981  else if (prefix == "SUMOFF" && dir_name == "Barrel")
982  nbin = 192;
983  else if ((*iv) == "adcCOMB")
984  nbin = 256;
985  else if (dir_name == "Barrel")
986  nbin = 768;
987  else if (prefix == "SUMOFF" && dir_name.find("Shell") != string::npos)
988  nbin = 48;
989  else if (dir_name.find("Shell") != string::npos)
990  nbin = 192;
991  else
992  nbin = nbin * nDirs;
993 
994  getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
995  } else if (first_subdir && isUpgrade) {
996  nbin = me->getTH1F()->GetNbinsX();
997  string me_name = prefix + "_" + (*iv) + "_" + dir_name;
998  if ((*iv) == "adcCOMB" || (*iv) == "chargeCOMB")
999  me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
1000  else if (prefix == "SUMOFF" && dir_name == "Barrel")
1001  nbin = 296;
1002  else if ((*iv) == "adcCOMB")
1003  nbin = 256;
1004  else if (dir_name == "Barrel")
1005  nbin = 1184;
1006  else if (prefix == "SUMOFF" && dir_name.find("Shell") != string::npos)
1007  nbin = 74;
1008  else if (dir_name.find("Shell") != string::npos)
1009  nbin = 296;
1010  else
1011  nbin = nbin * nDirs;
1012 
1013  getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
1014  }
1015 
1016  for (vector<MonitorElement *>::const_iterator igm = gsum_mes.begin(); igm != gsum_mes.end(); igm++) {
1017  if ((*igm)->getName().find(var) != string::npos) {
1018  if (prefix == "SUMOFF")
1019  (*igm)->setAxisTitle("Ladders", 1);
1020  else if ((*igm)->getName().find("adcCOMB_") != string::npos)
1021  (*igm)->setAxisTitle("Digi charge [ADC]", 1);
1022  else if ((*igm)->getName().find("chargeCOMB_") != string::npos)
1023  (*igm)->setAxisTitle("Cluster charge [kilo electrons]", 1);
1024  else
1025  (*igm)->setAxisTitle("Modules", 1);
1026 
1027  // Setting title
1028 
1029  string title = "";
1030  if ((*igm)->getName().find("NErrors_") != string::npos && prefix == "SUMOFF")
1031  title = "Total number of errors per Ladder";
1032  else if ((*igm)->getName().find("NErrors_") != string::npos && prefix == "SUMRAW")
1033  title = "Total number of errors per Module";
1034  else if (prefix == "SUMOFF")
1035  title = "mean " + (*iv) + " per Ladder";
1036  else if ((*igm)->getName().find("FREQ_") != string::npos && prefix != "SUMOFF")
1037  title = "NEvents with digis per Module";
1038  else if ((*igm)->getName().find("FREQ_") != string::npos && prefix == "SUMOFF")
1039  title = "NEvents with digis per Ladder/Blade";
1040  else if ((*igm)->getName().find("adcCOMB_") != string::npos)
1041  title = "NDigis";
1042  else if ((*igm)->getName().find("chargeCOMB_") != string::npos)
1043  title = "NClusters";
1044  else
1045  title = "mean " + (*iv) + " per Module";
1046  (*igm)->setAxisTitle(title, 2);
1047 
1048  // Setting binning
1049  if (!isUpgrade) {
1050  if ((*igm)->getName().find("ALLMODS_adcCOMB_") != string::npos) {
1051  nbin_subdir = 128;
1052  } else if ((*igm)->getName().find("ALLMODS_chargeCOMB_") != string::npos) {
1053  nbin_subdir = 100;
1054  } else if ((*igm)->getName().find("Ladder") != string::npos) {
1055  nbin_i = 0;
1056  nbin_subdir = 4;
1057  } else if ((*igm)->getName().find("Layer") != string::npos) {
1058  nbin_i = (cnt - 1) * 4;
1059  nbin_subdir = 4;
1060  } else if ((*igm)->getName().find("Shell") != string::npos) {
1061  if (prefix != "SUMOFF") {
1062  if (iDir == 0) {
1063  nbin_i = 0;
1064  nbin_subdir = 40;
1065  } else if (iDir == 1) {
1066  nbin_i = 40;
1067  nbin_subdir = 64;
1068  } else if (iDir == 2) {
1069  nbin_i = 104;
1070  nbin_subdir = 88;
1071  }
1072  } else {
1073  if (iDir == 0) {
1074  nbin_i = 0;
1075  nbin_subdir = 10;
1076  } else if (iDir == 1) {
1077  nbin_i = 10;
1078  nbin_subdir = 16;
1079  } else if (iDir == 2) {
1080  nbin_i = 26;
1081  nbin_subdir = 22;
1082  }
1083  }
1084  } else if ((*igm)->getName().find("Barrel") != string::npos) {
1085  if (prefix != "SUMOFF") {
1086  if (iDir == 0) {
1087  nbin_i = 0;
1088  nbin_subdir = 192;
1089  } else if (iDir == 1) {
1090  nbin_i = 192;
1091  nbin_subdir = 192;
1092  } else if (iDir == 2) {
1093  nbin_i = 384;
1094  nbin_subdir = 192;
1095  } else if (iDir == 3) {
1096  nbin_i = 576;
1097  nbin_subdir = 192;
1098  }
1099  } else {
1100  if (iDir == 0) {
1101  nbin_i = 0;
1102  nbin_subdir = 48;
1103  } else if (iDir == 1) {
1104  nbin_i = 48;
1105  nbin_subdir = 48;
1106  } else if (iDir == 2) {
1107  nbin_i = 96;
1108  nbin_subdir = 48;
1109  } else if (iDir == 3) {
1110  nbin_i = 144;
1111  nbin_subdir = 48;
1112  }
1113  }
1114  }
1115  } else if (isUpgrade) {
1116  if ((*igm)->getName().find("ALLMODS_adcCOMB_") != string::npos) {
1117  nbin_subdir = 128;
1118  } else if ((*igm)->getName().find("ALLMODS_chargeCOMB_") != string::npos) {
1119  nbin_subdir = 100;
1120  } else if ((*igm)->getName().find("Ladder") != string::npos) {
1121  nbin_i = 0;
1122  nbin_subdir = 4;
1123  } else if ((*igm)->getName().find("Layer") != string::npos) {
1124  nbin_i = (cnt - 1) * 4;
1125  nbin_subdir = 4;
1126  } else if ((*igm)->getName().find("Shell") != string::npos) {
1127  if (prefix != "SUMOFF") {
1128  if (iDir == 0) {
1129  nbin_i = 0;
1130  nbin_subdir = 24;
1131  } // 40(2*20)-->24(2*12)
1132  else if (iDir == 1) {
1133  nbin_i = 24;
1134  nbin_subdir = 56;
1135  } // 64(32*2)-->56(2*28)
1136  else if (iDir == 2) {
1137  nbin_i = 80;
1138  nbin_subdir = 88;
1139  } // 88(44*2)-->same88(44*2)
1140  else if (iDir == 3) {
1141  nbin_i = 168;
1142  nbin_subdir = 128;
1143  }
1144  } else {
1145  if (iDir == 0) {
1146  nbin_i = 0;
1147  nbin_subdir = 6;
1148  } // 10-->6
1149  else if (iDir == 1) {
1150  nbin_i = 6;
1151  nbin_subdir = 14;
1152  } // 16-->14
1153  else if (iDir == 2) {
1154  nbin_i = 20;
1155  nbin_subdir = 22;
1156  } // 22-->same22
1157  else if (iDir == 3) {
1158  nbin_i = 42;
1159  nbin_subdir = 32;
1160  }
1161  }
1162  } else if ((*igm)->getName().find("Barrel") != string::npos) {
1163  if (prefix != "SUMOFF") {
1164  if (iDir == 0) {
1165  nbin_i = 0;
1166  nbin_subdir = 296;
1167  } // 192=76 8/4-->296=1184/4
1168  else if (iDir == 1) {
1169  nbin_i = 296;
1170  nbin_subdir = 296;
1171  } // 296*2,*3,*4=1184
1172  else if (iDir == 2) {
1173  nbin_i = 592;
1174  nbin_subdir = 296;
1175  } else if (iDir == 3) {
1176  nbin_i = 888;
1177  nbin_subdir = 296;
1178  } else if (iDir == 4) {
1179  nbin_i = 1184;
1180  nbin_subdir = 296;
1181  }
1182  } else {
1183  if (iDir == 0) {
1184  nbin_i = 0;
1185  nbin_subdir = 74;
1186  } // 48=192/4-->74=296/4
1187  else if (iDir == 1) {
1188  nbin_i = 74;
1189  nbin_subdir = 74;
1190  } // 74*2,...*4=296
1191  else if (iDir == 2) {
1192  nbin_i = 148;
1193  nbin_subdir = 74;
1194  } else if (iDir == 3) {
1195  nbin_i = 222;
1196  nbin_subdir = 74;
1197  } else if (iDir == 4) {
1198  nbin_i = 296;
1199  nbin_subdir = 74;
1200  }
1201  }
1202  }
1203  }
1204 
1205  if ((*igm)->getName().find("ndigisFREQ") == string::npos) {
1206  if (((*igm)->getName().find("adcCOMB") != string::npos &&
1207  me->getName().find("adcCOMB") != string::npos) ||
1208  ((*igm)->getName().find("chargeCOMB") != string::npos &&
1209  me->getName().find("chargeCOMB") != string::npos)) {
1210  (*igm)->getTH1F()->Add(me->getTH1F());
1211  } else if (((*igm)->getName().find("charge_") != string::npos &&
1212  (*igm)->getName().find("Track_") == string::npos &&
1213  me->getName().find("charge_") != string::npos &&
1214  me->getName().find("Track_") == string::npos) ||
1215  ((*igm)->getName().find("nclusters_") != string::npos &&
1216  (*igm)->getName().find("Track_") == string::npos &&
1217  me->getName().find("nclusters_") != string::npos &&
1218  me->getName().find("Track_") == string::npos) ||
1219  ((*igm)->getName().find("size_") != string::npos &&
1220  (*igm)->getName().find("Track_") == string::npos &&
1221  me->getName().find("size_") != string::npos &&
1222  me->getName().find("Track_") == string::npos) ||
1223  ((*igm)->getName().find("charge_OffTrack_") != string::npos &&
1224  me->getName().find("charge_OffTrack_") != string::npos) ||
1225  ((*igm)->getName().find("nclusters_OffTrack_") != string::npos &&
1226  me->getName().find("nclusters_OffTrack_") != string::npos) ||
1227  ((*igm)->getName().find("size_OffTrack_") != string::npos &&
1228  me->getName().find("size_OffTrack_") != string::npos) ||
1229  ((*igm)->getName().find("charge_OnTrack_") != string::npos &&
1230  me->getName().find("charge_OnTrack_") != string::npos) ||
1231  ((*igm)->getName().find("nclusters_OnTrack_") != string::npos &&
1232  me->getName().find("nclusters_OnTrack_") != string::npos) ||
1233  ((*igm)->getName().find("size_OnTrack_") != string::npos &&
1234  me->getName().find("size_OnTrack_") != string::npos) ||
1235  ((*igm)->getName().find("charge_") == string::npos &&
1236  (*igm)->getName().find("nclusters_") == string::npos &&
1237  (*igm)->getName().find("size_") == string::npos)) {
1238  for (int k = 1; k < nbin_subdir + 1; k++)
1239  if (me->getBinContent(k) > 0)
1240  (*igm)->setBinContent(k + nbin_i, me->getBinContent(k));
1241  }
1242  } else if (me->getName().find("ndigisFREQ") != string::npos) {
1243  for (int k = 1; k < nbin_subdir + 1; k++)
1244  if (me->getBinContent(k) > 0)
1245  (*igm)->setBinContent(k + nbin_i, me->getBinContent(k));
1246  }
1247  } // end var in igm (gsum_mes)
1248  } // end igm loop
1249  } // end var in im (contents)
1250  } // end of iv loop
1251  } // end of im loop
1252  iDir++;
1253  first_subdir = false; // We are done processing the first directory, we
1254  // don't add any new MEs in the future passes.
1255  } // end of it loop (subdirs)
1256  // cout<<"...leaving
1257  // SiPixelActionExecutor::fillGrandBarrelSummaryHistos!"<<endl;
1258 }
1259 
1260 //=============================================================================================================
1261 void SiPixelActionExecutor::fillGrandEndcapSummaryHistos(DQMStore::IBooker &iBooker,
1262  DQMStore::IGetter &iGetter,
1263  vector<string> &me_names,
1264  bool isUpgrade) {
1265  // printing cout<<"Entering
1266  // SiPixelActionExecutor::fillGrandEndcapSummaryHistos..."<<endl;
1267  vector<MonitorElement *> gsum_mes;
1268  string currDir = iBooker.pwd();
1269  string path_name = iBooker.pwd();
1270  string dir_name = path_name.substr(path_name.find_last_of('/') + 1);
1271  if ((dir_name.find("DQMData") == 0) || (dir_name.find("Pixel") == 0) ||
1272  (dir_name.find("AdditionalPixelErrors") == 0) || (dir_name.find("Barrel") == 0) ||
1273  (dir_name.find("Shell") == 0) || (dir_name.find("Layer") == 0) || (dir_name.find("Ladder") == 0))
1274  return;
1275  vector<string> subdirs = iGetter.getSubdirs();
1276  int iDir = 0;
1277  int nbin = 0;
1278  int nbin_i = 0;
1279  int nbin_subdir = 0;
1280  int cnt = 0;
1281  bool first_subdir = true;
1282  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
1283  cnt++;
1284  iBooker.cd(*it);
1285  iGetter.cd(*it);
1286  vector<string> contents = iGetter.getMEs();
1287  iBooker.goUp();
1288  iGetter.setCurrentFolder(iBooker.pwd());
1289 
1290  string prefix;
1291  if (source_type_ == 0)
1292  prefix = "SUMRAW";
1293  else if (source_type_ == 1)
1294  prefix = "SUMDIG";
1295  else if (source_type_ == 2)
1296  prefix = "SUMCLU";
1297  else if (source_type_ == 3)
1298  prefix = "SUMTRK";
1299  else if (source_type_ == 4)
1300  prefix = "SUMHIT";
1301  else if (source_type_ >= 7 && source_type_ < 20)
1302  prefix = "SUMCAL";
1303  else if (source_type_ == 20)
1304  prefix = "SUMOFF";
1305 
1306  for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++) {
1307  for (vector<string>::const_iterator iv = me_names.begin(); iv != me_names.end(); iv++) {
1308  string var = "_" + (*iv) + "_";
1309  if ((*im).find(var) != string::npos) {
1310  if ((var == "_charge_" || var == "_nclusters_" || var == "_size_" || var == "_sizeX_" || var == "_sizeY_") &&
1311  (*im).find("Track_") != string::npos)
1312  continue;
1313  string full_path = (*it) + "/" + (*im);
1314  MonitorElement *me = iGetter.get(full_path);
1315  if (!me)
1316  continue;
1317  if (source_type_ == 5 || source_type_ == 6) {
1318  if ((*iv) == "errorType" || (*iv) == "NErrors" || (*iv) == "fullType" || (*iv) == "chanNmbr" ||
1319  (*iv) == "TBMType" || (*iv) == "EvtNbr" || (*iv) == "evtSize" || (*iv) == "linkId" ||
1320  (*iv) == "ROCId" || (*iv) == "DCOLId" || (*iv) == "PXId" || (*iv) == "ROCNmbr" ||
1321  (*iv) == "TBMMessage" || (*iv) == "Type36Hitmap")
1322  prefix = "SUMRAW";
1323  else if ((*iv) == "ndigis" || (*iv) == "adc" || (*iv) == "ndigisFREQ" || (*iv) == "adcCOMB")
1324  prefix = "SUMDIG";
1325  else if ((*iv) == "nclusters" || (*iv) == "x" || (*iv) == "y" || (*iv) == "charge" ||
1326  (*iv) == "chargeCOMB" || (*iv) == "size" || (*iv) == "sizeX" || (*iv) == "sizeY" ||
1327  (*iv) == "minrow" || (*iv) == "maxrow" || (*iv) == "mincol" || (*iv) == "maxcol")
1328  prefix = "SUMCLU";
1329  if (currDir.find("Track") != string::npos)
1330  prefix = "SUMTRK";
1331  else if ((*iv) == "residualX_mean" || (*iv) == "residualY_mean" || (*iv) == "residualX_RMS" ||
1332  (*iv) == "residualY_RMS")
1333  prefix = "SUMTRK";
1334  else if ((*iv) == "ClustX" || (*iv) == "ClustY" || (*iv) == "nRecHits" || (*iv) == "ErrorX" ||
1335  (*iv) == "ErrorY")
1336  prefix = "SUMHIT";
1337  else if ((*iv) == "Gain1d_mean" || (*iv) == "GainChi2NDF1d_mean" || (*iv) == "GainChi2Prob1d_mean" ||
1338  (*iv) == "Pedestal1d_mean" || (*iv) == "ScurveChi2ProbSummary_mean" ||
1339  (*iv) == "ScurveFitResultSummary_mean" || (*iv) == "ScurveSigmasSummary_mean" ||
1340  (*iv) == "ScurveThresholdSummary_mean" || (*iv) == "Gain1d_RMS" || (*iv) == "GainChi2NDF1d_RMS" ||
1341  (*iv) == "GainChi2Prob1d_RMS" || (*iv) == "Pedestal1d_RMS" || (*iv) == "GainNPoints1d_mean" ||
1342  (*iv) == "GainNPoints1d_RMS" || (*iv) == "GainHighPoint1d_mean" ||
1343  (*iv) == "GainHighPoint1d_RMS" || (*iv) == "GainLowPoint1d_mean" ||
1344  (*iv) == "GainLowPoint1d_RMS" || (*iv) == "GainEndPoint1d_mean" || (*iv) == "GainEndPoint1d_RMS" ||
1345  (*iv) == "GainFitResult2d_mean" || (*iv) == "GainFitResult2d_RMS" ||
1346  (*iv) == "GainDynamicRange2d_mean" || (*iv) == "GainDynamicRange2d_RMS" ||
1347  (*iv) == "GainSaturate2d_mean" || (*iv) == "GainSaturate2d_RMS" ||
1348  (*iv) == "ScurveChi2ProbSummary_RMS" || (*iv) == "ScurveFitResultSummary_RMS" ||
1349  (*iv) == "ScurveSigmasSummary_RMS" || (*iv) == "ScurveThresholdSummary_RMS" ||
1350  (*iv) == "pixelAliveSummary_mean" || (*iv) == "pixelAliveSummary_FracOfPerfectPix" ||
1351  (*iv) == "SiPixelErrorsCalibDigis_NCalibErrors")
1352  prefix = "SUMCAL";
1353  }
1354 
1355  if (first_subdir && !isUpgrade) {
1356  nbin = me->getTH1F()->GetNbinsX();
1357  string me_name = prefix + "_" + (*iv) + "_" + dir_name;
1358  if ((*iv) == "adcCOMB" || (*iv) == "chargeCOMB")
1359  me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
1360  else if (prefix == "SUMOFF" && dir_name == "Endcap")
1361  nbin = 96;
1362  else if (dir_name == "Endcap")
1363  nbin = 672;
1364  else if (prefix == "SUMOFF" && dir_name.find("HalfCylinder") != string::npos)
1365  nbin = 24;
1366  else if (dir_name.find("HalfCylinder") != string::npos)
1367  nbin = 168;
1368  else if (prefix == "SUMOFF" && dir_name.find("Disk") != string::npos)
1369  nbin = 12;
1370  else if (dir_name.find("Disk") != string::npos)
1371  nbin = 84;
1372  else if (dir_name.find("Blade") != string::npos)
1373  nbin = 7;
1374  getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
1375  } else if (first_subdir && isUpgrade) {
1376  nbin = me->getTH1F()->GetNbinsX();
1377  string me_name = prefix + "_" + (*iv) + "_" + dir_name;
1378  if ((*iv) == "adcCOMB" || (*iv) == "chargeCOMB")
1379  me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
1380  else if (prefix == "SUMOFF" && dir_name == "Endcap")
1381  nbin = 336;
1382  else if (dir_name == "Endcap")
1383  nbin = 672;
1384  else if (prefix == "SUMOFF" && dir_name.find("HalfCylinder") != string::npos)
1385  nbin = 84;
1386  else if (dir_name.find("HalfCylinder") != string::npos)
1387  nbin = 168;
1388  else if (prefix == "SUMOFF" && dir_name.find("Disk") != string::npos)
1389  nbin = 28;
1390  else if (dir_name.find("Disk") != string::npos)
1391  nbin = 56;
1392  else if (dir_name.find("Blade") != string::npos)
1393  nbin = 2;
1394  getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
1395  }
1396 
1397  for (vector<MonitorElement *>::const_iterator igm = gsum_mes.begin(); igm != gsum_mes.end(); igm++) {
1398  if ((*igm)->getName().find(var) != string::npos) {
1399  if (prefix == "SUMOFF")
1400  (*igm)->setAxisTitle("Blades", 1);
1401  else if ((*igm)->getName().find("adcCOMB_") != string::npos)
1402  (*igm)->setAxisTitle("Digi charge [ADC]", 1);
1403  else if ((*igm)->getName().find("chargeCOMB_") != string::npos)
1404  (*igm)->setAxisTitle("Cluster charge [kilo electrons]", 1);
1405  else
1406  (*igm)->setAxisTitle("Modules", 1);
1407  string title = "";
1408  if ((*igm)->getName().find("NErrors_") != string::npos && prefix == "SUMOFF")
1409  title = "Total number of errors per Blade";
1410  else if ((*igm)->getName().find("NErrors_") != string::npos && prefix == "SUMRAW")
1411  title = "Total number of errors per Module";
1412  else if (prefix == "SUMOFF")
1413  title = "mean " + (*iv) + " per Blade";
1414  else if ((*igm)->getName().find("FREQ_") != string::npos)
1415  title = "NEvents with digis per Module";
1416  else if ((*igm)->getName().find("adcCOMB_") != string::npos)
1417  title = "NDigis";
1418  else if ((*igm)->getName().find("chargeCOMB_") != string::npos)
1419  title = "NClusters";
1420  else
1421  title = "mean " + (*iv) + " per Module";
1422  (*igm)->setAxisTitle(title, 2);
1423  nbin_i = 0;
1424  if (!isUpgrade) {
1425  if ((*igm)->getName().find("ALLMODS_adcCOMB_") != string::npos) {
1426  nbin_subdir = 128;
1427  } else if ((*igm)->getName().find("ALLMODS_chargeCOMB_") != string::npos) {
1428  nbin_subdir = 100;
1429  } else if ((*igm)->getName().find("Panel_") != string::npos) {
1430  nbin_subdir = 7;
1431  } else if ((*igm)->getName().find("Blade") != string::npos) {
1432  if ((*im).find("_1") != string::npos)
1433  nbin_subdir = 4;
1434  if ((*im).find("_2") != string::npos) {
1435  nbin_i = 4;
1436  nbin_subdir = 3;
1437  }
1438  } else if ((*igm)->getName().find("Disk") != string::npos) {
1439  nbin_i = ((cnt - 1) % 12) * 7;
1440  nbin_subdir = 7;
1441  } else if ((*igm)->getName().find("HalfCylinder") != string::npos) {
1442  if (prefix != "SUMOFF") {
1443  nbin_subdir = 84;
1444  if ((*im).find("_2") != string::npos)
1445  nbin_i = 84;
1446  } else {
1447  nbin_subdir = 12;
1448  if ((*im).find("_2") != string::npos)
1449  nbin_i = 12;
1450  }
1451  } else if ((*igm)->getName().find("Endcap") != string::npos) {
1452  if (prefix != "SUMOFF") {
1453  nbin_subdir = 168;
1454  if ((*im).find("_mO") != string::npos)
1455  nbin_i = 168;
1456  if ((*im).find("_pI") != string::npos)
1457  nbin_i = 336;
1458  if ((*im).find("_pO") != string::npos)
1459  nbin_i = 504;
1460  } else {
1461  nbin_subdir = 24;
1462  if ((*im).find("_mO") != string::npos)
1463  nbin_i = 24;
1464  if ((*im).find("_pI") != string::npos)
1465  nbin_i = 48;
1466  if ((*im).find("_pO") != string::npos)
1467  nbin_i = 72;
1468  }
1469  }
1470  } else if (isUpgrade) {
1471  if ((*igm)->getName().find("ALLMODS_adcCOMB_") != string::npos) {
1472  nbin_subdir = 128;
1473  } else if ((*igm)->getName().find("ALLMODS_chargeCOMB_") != string::npos) {
1474  nbin_subdir = 100;
1475  } else if ((*igm)->getName().find("Panel_") != string::npos) {
1476  nbin_subdir = 2;
1477  } else if ((*igm)->getName().find("Blade") != string::npos) {
1478  if ((*im).find("_1") != string::npos)
1479  nbin_subdir = 1;
1480  if ((*im).find("_2") != string::npos) {
1481  nbin_i = 1;
1482  nbin_subdir = 1;
1483  }
1484  } else if ((*igm)->getName().find("Disk") != string::npos) {
1485  nbin_i = ((cnt - 1) % 28) * 2;
1486  nbin_subdir = 2;
1487  } else if ((*igm)->getName().find("HalfCylinder") != string::npos) {
1488  if (prefix != "SUMOFF") {
1489  nbin_subdir = 56;
1490  if ((*im).find("_2") != string::npos)
1491  nbin_i = 56;
1492  if ((*im).find("_3") != string::npos)
1493  nbin_i = 112;
1494  } else {
1495  nbin_subdir = 28;
1496  if ((*im).find("_2") != string::npos)
1497  nbin_i = 28;
1498  if ((*im).find("_3") != string::npos)
1499  nbin_i = 56;
1500  }
1501  } else if ((*igm)->getName().find("Endcap") != string::npos) {
1502  if (prefix != "SUMOFF") {
1503  nbin_subdir = 168;
1504  if ((*im).find("_mO") != string::npos)
1505  nbin_i = 168;
1506  if ((*im).find("_pI") != string::npos)
1507  nbin_i = 336;
1508  if ((*im).find("_pO") != string::npos)
1509  nbin_i = 504;
1510  } else {
1511  nbin_subdir = 84;
1512  if ((*im).find("_mO") != string::npos)
1513  nbin_i = 84;
1514  if ((*im).find("_pI") != string::npos)
1515  nbin_i = 168;
1516  if ((*im).find("_pO") != string::npos)
1517  nbin_i = 252;
1518  }
1519  }
1520  }
1521 
1522  if ((*igm)->getName().find("ndigisFREQ") == string::npos) {
1523  if (((*igm)->getName().find("adcCOMB") != string::npos &&
1524  me->getName().find("adcCOMB") != string::npos) ||
1525  ((*igm)->getName().find("chargeCOMB") != string::npos &&
1526  me->getName().find("chargeCOMB") != string::npos)) {
1527  (*igm)->getTH1F()->Add(me->getTH1F());
1528  } else if (((*igm)->getName().find("charge_") != string::npos &&
1529  (*igm)->getName().find("Track_") == string::npos &&
1530  me->getName().find("charge_") != string::npos &&
1531  me->getName().find("Track_") == string::npos) ||
1532  ((*igm)->getName().find("nclusters_") != string::npos &&
1533  (*igm)->getName().find("Track_") == string::npos &&
1534  me->getName().find("nclusters_") != string::npos &&
1535  me->getName().find("Track_") == string::npos) ||
1536  ((*igm)->getName().find("size_") != string::npos &&
1537  (*igm)->getName().find("Track_") == string::npos &&
1538  me->getName().find("size_") != string::npos &&
1539  me->getName().find("Track_") == string::npos) ||
1540  ((*igm)->getName().find("charge_OffTrack_") != string::npos &&
1541  me->getName().find("charge_OffTrack_") != string::npos) ||
1542  ((*igm)->getName().find("nclusters_OffTrack_") != string::npos &&
1543  me->getName().find("nclusters_OffTrack_") != string::npos) ||
1544  ((*igm)->getName().find("size_OffTrack_") != string::npos &&
1545  me->getName().find("size_OffTrack_") != string::npos) ||
1546  ((*igm)->getName().find("charge_OnTrack_") != string::npos &&
1547  me->getName().find("charge_OnTrack_") != string::npos) ||
1548  ((*igm)->getName().find("nclusters_OnTrack_") != string::npos &&
1549  me->getName().find("nclusters_OnTrack_") != string::npos) ||
1550  ((*igm)->getName().find("size_OnTrack_") != string::npos &&
1551  me->getName().find("size_OnTrack_") != string::npos) ||
1552  ((*igm)->getName().find("charge_") == string::npos &&
1553  (*igm)->getName().find("nclusters_") == string::npos &&
1554  (*igm)->getName().find("size_") == string::npos)) {
1555  for (int k = 1; k < nbin_subdir + 1; k++)
1556  if (me->getBinContent(k) > 0)
1557  (*igm)->setBinContent(k + nbin_i, me->getBinContent(k));
1558  }
1559  } else if (me->getName().find("ndigisFREQ") != string::npos) {
1560  for (int k = 1; k < nbin_subdir + 1; k++)
1561  if (me->getBinContent(k) > 0)
1562  (*igm)->setBinContent(k + nbin_i, me->getBinContent(k));
1563  }
1564  // }// for
1565  }
1566  }
1567  }
1568  }
1569  }
1570 
1571  iDir++;
1572  first_subdir = false; // We are done processing the first directory, we
1573  // don't add any new MEs in the future passes.
1574  } // end for it (subdirs)
1575 }
1576 //=============================================================================================================
1577 //
1578 // -- Get Summary ME
1579 //
1580 void SiPixelActionExecutor::getGrandSummaryME(
1581  DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, int nbin, string &me_name, vector<MonitorElement *> &mes) {
1582  // printing cout<<"Entering SiPixelActionExecutor::getGrandSummaryME for:
1583  // "<<me_name<<endl;
1584  if ((iBooker.pwd()).find("Pixel") == string::npos)
1585  return; // If one doesn't find pixel
1586  vector<string> contents = iGetter.getMEs();
1587 
1588  for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) {
1589  // printing cout<<"in grand summary me: "<<me_name<<","<<(*it)<<endl;
1590  if ((*it).find(me_name) == 0) {
1591  string fullpathname = iBooker.pwd() + "/" + me_name;
1592  MonitorElement *me = iGetter.get(fullpathname);
1593 
1594  if (me) {
1595  me->Reset();
1596  mes.push_back(me);
1597  return;
1598  }
1599  }
1600  }
1601 
1602  MonitorElement *temp_me(nullptr);
1603  if (me_name.find("ALLMODS_adcCOMB_") != string::npos)
1604  temp_me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 128, 0, 256);
1605  else if (me_name.find("ALLMODS_chargeCOMB_") != string::npos)
1606  temp_me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 100, 0, 200);
1607  else
1608  temp_me = iBooker.book1D(me_name.c_str(), me_name.c_str(), nbin, 1., nbin + 1.);
1609  if (temp_me)
1610  mes.push_back(temp_me);
1611 
1612  // if(temp_me) cout<<"finally found grand ME: "<<me_name<<endl;
1613 }
1614 
1615 //=============================================================================================================
1616 //
1617 // -- Get Summary ME
1618 //
1619 SiPixelActionExecutor::MonitorElement *SiPixelActionExecutor::getSummaryME(DQMStore::IBooker &iBooker,
1620  DQMStore::IGetter &iGetter,
1621  string me_name,
1622  bool isUpgrade) {
1623  // printing cout<<"Entering SiPixelActionExecutor::getSummaryME for:
1624  // "<<me_name<<endl;
1625  MonitorElement *me = nullptr;
1626  if ((iBooker.pwd()).find("Pixel") == string::npos)
1627  return me;
1628  vector<string> contents = iGetter.getMEs();
1629 
1630  for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) {
1631  if ((*it).find(me_name) == 0) {
1632  string fullpathname = iBooker.pwd() + "/" + (*it);
1633  me = iGetter.get(fullpathname);
1634  if (me) {
1635  me->Reset();
1636  return me;
1637  }
1638  }
1639  }
1640  contents.clear();
1641  if (!isUpgrade) {
1642  if (me_name.find("SUMOFF") == string::npos) {
1643  if (me_name.find("Blade_") != string::npos)
1644  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 7, 1., 8.);
1645  else
1646  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 4, 1., 5.);
1647  } else if (me_name.find("Layer_1") != string::npos) {
1648  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 10, 1., 11.);
1649  } else if (me_name.find("Layer_2") != string::npos) {
1650  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 16, 1., 17.);
1651  } else if (me_name.find("Layer_3") != string::npos) {
1652  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 22, 1., 23.);
1653  } else if (me_name.find("Disk_") != string::npos) {
1654  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 12, 1., 13.);
1655  }
1656  } // endifNOTUpgrade
1657  else if (isUpgrade) {
1658  if (me_name.find("SUMOFF") == string::npos) {
1659  if (me_name.find("Blade_") != string::npos)
1660  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 2, 1., 3.);
1661  else
1662  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 1, 1., 2.);
1663  } else if (me_name.find("Layer_1") != string::npos) {
1664  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 6, 1., 7.);
1665  } else if (me_name.find("Layer_2") != string::npos) {
1666  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 14, 1., 15.);
1667  } else if (me_name.find("Layer_3") != string::npos) {
1668  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 22, 1., 23.);
1669  } else if (me_name.find("Layer_4") != string::npos) {
1670  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 32, 1., 33.);
1671  } else if (me_name.find("Disk_") != string::npos) {
1672  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 28, 1., 29.);
1673  }
1674  } // endifUpgrade
1675 
1676  return me;
1677 }
1678 
1679 //=============================================================================================================
1680 SiPixelActionExecutor::MonitorElement *SiPixelActionExecutor::getFEDSummaryME(DQMStore::IBooker &iBooker,
1681  DQMStore::IGetter &iGetter,
1682  string me_name) {
1683  // printing cout<<"Entering SiPixelActionExecutor::getFEDSummaryME..."<<endl;
1684  MonitorElement *me = nullptr;
1685  if ((iBooker.pwd()).find("Pixel") == string::npos)
1686  return me;
1687  vector<string> contents = iGetter.getMEs();
1688 
1689  for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) {
1690  if ((*it).find(me_name) == 0) {
1691  string fullpathname = iBooker.pwd() + "/" + (*it);
1692 
1693  me = iGetter.get(fullpathname);
1694 
1695  if (me) {
1696  me->Reset();
1697  return me;
1698  }
1699  }
1700  }
1701  contents.clear();
1702  me = iBooker.book1D(me_name.c_str(), me_name.c_str(), 40, -0.5, 39.5);
1703 
1704  return me;
1705 }
1706 
1707 //=============================================================================================================
1708 void SiPixelActionExecutor::bookOccupancyPlots(DQMStore::IBooker &iBooker,
1709  DQMStore::IGetter &iGetter,
1710  bool hiRes,
1711  bool isbarrel) // Polymorphism
1712 {
1713  if (Tier0Flag_)
1714  return;
1715  vector<string> subdirs = iGetter.getSubdirs();
1716  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
1717  if (isbarrel && (*it).find("Barrel") == string::npos)
1718  continue;
1719  if (!isbarrel && (*it).find("Endcap") == string::npos)
1720  continue;
1721 
1722  if ((*it).find("Module_") != string::npos)
1723  continue;
1724  if ((*it).find("Panel_") != string::npos)
1725  continue;
1726  if ((*it).find("Ladder_") != string::npos)
1727  continue;
1728  if ((*it).find("Blade_") != string::npos)
1729  continue;
1730  if ((*it).find("Layer_") != string::npos)
1731  continue;
1732  if ((*it).find("Disk_") != string::npos)
1733  continue;
1734  iBooker.cd(*it);
1735  iGetter.cd(*it);
1736  bookOccupancyPlots(iBooker, iGetter, hiRes, isbarrel);
1737  if (!hiRes) {
1738  // occupancyprinting cout<<"booking low res barrel occ plot now!"<<endl;
1739  OccupancyMap = iBooker.book2D((isbarrel ? "barrelOccupancyMap" : "endcapOccupancyMap"),
1740  "Barrel Digi Occupancy Map (4 pix per bin)",
1741  isbarrel ? 208 : 130,
1742  0.,
1743  isbarrel ? 416. : 260.,
1744  80,
1745  0.,
1746  160.);
1747  } else {
1748  // occupancyprinting cout<<"booking high res barrel occ plot now!"<<endl;
1749  OccupancyMap = iBooker.book2D((isbarrel ? "barrelOccupancyMap" : "endcapOccupancyMap"),
1750  "Barrel Digi Occupancy Map (1 pix per bin)",
1751  isbarrel ? 416 : 260,
1752  0.,
1753  isbarrel ? 416. : 260.,
1754  160,
1755  0.,
1756  160.);
1757  }
1758  OccupancyMap->setAxisTitle("Columns", 1);
1759  OccupancyMap->setAxisTitle("Rows", 2);
1760 
1761  iBooker.goUp();
1762  iGetter.setCurrentFolder(iBooker.pwd());
1763  }
1764 }
1765 //=============================================================================================================
1766 void SiPixelActionExecutor::bookOccupancyPlots(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool hiRes) {
1767  if (Tier0Flag_)
1768  return;
1769  // Barrel
1770  iGetter.cd();
1771  iBooker.cd();
1772  iGetter.setCurrentFolder("Pixel");
1773  iBooker.setCurrentFolder("Pixel");
1774  this->bookOccupancyPlots(iBooker, iGetter, hiRes, true);
1775 
1776  // Endcap
1777  iGetter.cd();
1778  iBooker.cd();
1779  iGetter.setCurrentFolder("Pixel");
1780  iBooker.setCurrentFolder("Pixel");
1781  this->bookOccupancyPlots(iBooker, iGetter, hiRes, false);
1782 }
1783 
1784 void SiPixelActionExecutor::createOccupancy(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) {
1785  // std::cout<<"entering SiPixelActionExecutor::createOccupancy..."<<std::endl;
1786  if (Tier0Flag_)
1787  return;
1788  iBooker.cd();
1789  iGetter.cd();
1790  fillOccupancy(iBooker, iGetter, true);
1791  iBooker.cd();
1792  iGetter.cd();
1793  fillOccupancy(iBooker, iGetter, false);
1794  iBooker.cd();
1795  iGetter.cd();
1796 
1797  // std::cout<<"leaving SiPixelActionExecutor::createOccupancy..."<<std::endl;
1798 }
1799 
1800 //=============================================================================================================
1801 
1802 void SiPixelActionExecutor::fillOccupancy(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isbarrel) {
1803  // occupancyprinting cout<<"entering
1804  // SiPixelActionExecutor::fillOccupancy..."<<std::endl;
1805  if (Tier0Flag_)
1806  return;
1807  string currDir = iBooker.pwd();
1808  string dname = currDir.substr(currDir.find_last_of('/') + 1);
1809 
1810  if (dname.find("Layer_") != string::npos || dname.find("Disk_") != string::npos) {
1811  vector<string> meVec = iGetter.getMEs();
1812  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
1813  string full_path = currDir + "/" + (*it);
1814  if (full_path.find("hitmap_siPixelDigis") != string::npos) { // If we have the hitmap ME
1815  MonitorElement *me = iGetter.get(full_path);
1816  if (!me)
1817  continue;
1818  string path = full_path;
1819  while (path.find_last_of('/') != 5) // Stop before Pixel/
1820  {
1821  path = path.substr(0, path.find_last_of('/'));
1822  // cout << "\t" <<
1823  // path
1824  //<< endl;
1825  OccupancyMap = iGetter.get(path + "/" + (isbarrel ? "barrel" : "endcap") + "OccupancyMap");
1826 
1827  if (OccupancyMap) {
1828  for (int i = 1; i != me->getNbinsX() + 1; i++)
1829  for (int j = 1; j != me->getNbinsY() + 1; j++) {
1830  float previous = OccupancyMap->getBinContent(i, j);
1831  OccupancyMap->setBinContent(i, j, previous + me->getBinContent(i, j));
1832  }
1833  OccupancyMap->getTH2F()->SetEntries(OccupancyMap->getTH2F()->Integral());
1834  }
1835  }
1836  }
1837  }
1838  } else {
1839  vector<string> subdirs = iGetter.getSubdirs();
1840  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
1841  iGetter.cd(*it);
1842  iBooker.cd(*it);
1843  if (*it != "Pixel" &&
1844  ((isbarrel && (*it).find("Barrel") == string::npos) || (!isbarrel && (*it).find("Endcap") == string::npos)))
1845  continue;
1846  fillOccupancy(iBooker, iGetter, isbarrel);
1847  iBooker.goUp();
1848  iGetter.setCurrentFolder(iBooker.pwd());
1849  }
1850  }
1851 
1852  // occupancyprinting cout<<"leaving
1853  // SiPixelActionExecutor::fillOccupancy..."<<std::endl;
1854 }
1855 
1856 //=============================================================================================================
1857 
1858 void SiPixelActionExecutor::normaliseAvDigiOcc(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) {
1859  // occupancyprinting cout<<"entering
1860  // SiPixelActionExecutor::normaliseAvDigiOcc..."<<std::endl;
1861 
1862  iGetter.cd();
1863 
1864  MonitorElement *roccupancyPlot = iGetter.get("Pixel/averageDigiOccupancy");
1865 
1866  float totalDigisBPIX = 0.;
1867  float totalDigisFPIX = 0.;
1868  for (int i = 1; i != 41; i++) {
1869  if (i < 33)
1870  totalDigisBPIX += roccupancyPlot->getBinContent(i);
1871  else
1872  totalDigisFPIX += roccupancyPlot->getBinContent(i);
1873  }
1874  float averageBPIXOcc = totalDigisBPIX / 32.;
1875  float averageFPIXOcc = totalDigisFPIX / 8.;
1876  for (int i = 1; i != 41; i++) {
1877  if (i < 33)
1878  roccupancyPlot->setBinContent(i, averageBPIXOcc != 0. ? roccupancyPlot->getBinContent(i) / averageBPIXOcc : 0.);
1879  else
1880  roccupancyPlot->setBinContent(i, averageFPIXOcc != 0. ? roccupancyPlot->getBinContent(i) / averageFPIXOcc : 0.);
1881  }
1882 
1883  iGetter.setCurrentFolder(iBooker.pwd());
1884 }
1885 
1886 //=============================================================================================================
1887 
1888 void SiPixelActionExecutor::normaliseAvDigiOccVsLumi(DQMStore::IBooker &iBooker,
1889  DQMStore::IGetter &iGetter,
1890  int lumisec) {
1891  iGetter.cd();
1892 
1893  MonitorElement *avgfedDigiOccvsLumi = iGetter.get("Pixel/avgfedDigiOccvsLumi");
1894 
1895  float totalDigisBPIX = 0.;
1896  float totalDigisFPIX = 0.;
1897  for (int i = 1; i != 41; i++) {
1898  if (i < 33)
1899  totalDigisBPIX += avgfedDigiOccvsLumi->getBinContent(lumisec, i);
1900  else
1901  totalDigisFPIX += avgfedDigiOccvsLumi->getBinContent(lumisec, i);
1902  }
1903  float averageBPIXOcc = totalDigisBPIX / 32.;
1904  float averageFPIXOcc = totalDigisFPIX / 8.;
1905  for (int i = 1; i != 41; i++) {
1906  if (i < 33) {
1907  avgfedDigiOccvsLumi->setBinContent(
1908  lumisec, i, averageBPIXOcc != 0. ? avgfedDigiOccvsLumi->getBinContent(lumisec, i) / averageBPIXOcc : 0.);
1909  } else {
1910  avgfedDigiOccvsLumi->setBinContent(
1911  lumisec, i, averageFPIXOcc != 0. ? avgfedDigiOccvsLumi->getBinContent(lumisec, i) / averageFPIXOcc : 0.);
1912  }
1913  }
1914 
1915  iGetter.setCurrentFolder(iBooker.pwd());
1916 }
1917 
1918 //=============================================================================================================
1919 
1920 void SiPixelActionExecutor::bookEfficiency(DQMStore::IBooker &iBooker, bool isUpgrade) {
1921  // Barrel
1922  iBooker.cd();
1923  iBooker.setCurrentFolder("Pixel/Barrel");
1924  if (!isUpgrade) {
1925  if (Tier0Flag_) {
1926  HitEfficiency_L1 = iBooker.book2D(
1927  "HitEfficiency_L1", "Hit Efficiency in Barrel_Layer1;Module;Ladder", 9, -4.5, 4.5, 21, -10.5, 10.5);
1928  HitEfficiency_L2 = iBooker.book2D(
1929  "HitEfficiency_L2", "Hit Efficiency in Barrel_Layer2;Module;Ladder", 9, -4.5, 4.5, 33, -16.5, 16.5);
1930  HitEfficiency_L3 = iBooker.book2D(
1931  "HitEfficiency_L3", "Hit Efficiency in Barrel_Layer3;Module;Ladder", 9, -4.5, 4.5, 45, -22.5, 22.5);
1932  } else {
1933  HitEfficiency_L1 = iBooker.book2D(
1934  "HitEfficiency_L1", "Hit Efficiency in Barrel_Layer1;Module;Ladder", 9, -4.5, 4.5, 21, -10.5, 10.5);
1935  HitEfficiency_L2 = iBooker.book2D(
1936  "HitEfficiency_L2", "Hit Efficiency in Barrel_Layer2;Module;Ladder", 9, -4.5, 4.5, 33, -16.5, 16.5);
1937  HitEfficiency_L3 = iBooker.book2D(
1938  "HitEfficiency_L3", "Hit Efficiency in Barrel_Layer3;Module;Ladder", 9, -4.5, 4.5, 45, -22.5, 22.5);
1939  }
1940  } // endifNOTUpgrade
1941  else if (isUpgrade) {
1942  if (Tier0Flag_) {
1943  HitEfficiency_L1 =
1944  iBooker.book2D("HitEfficiency_L1", "Hit Efficiency in Barrel_Layer1;z-side;Ladder", 2, -1., 1., 12, -6., 6.);
1945  HitEfficiency_L2 = iBooker.book2D(
1946  "HitEfficiency_L2", "Hit Efficiency in Barrel_Layer2;z-side;Ladder", 2, -1., 1., 28, -14., 14.);
1947  HitEfficiency_L3 = iBooker.book2D(
1948  "HitEfficiency_L3", "Hit Efficiency in Barrel_Layer3;z-side;Ladder", 2, -1., 1., 44, -22., 22.);
1949  HitEfficiency_L4 = iBooker.book2D(
1950  "HitEfficiency_L4", "Hit Efficiency in Barrel_Layer4;z-side;Ladder", 2, -1., 1., 64, -32., 32.);
1951  } else {
1952  HitEfficiency_L1 =
1953  iBooker.book2D("HitEfficiency_L1", "Hit Efficiency in Barrel_Layer1;Module;Ladder", 8, -4., 4., 12, -6., 6.);
1954  HitEfficiency_L2 = iBooker.book2D(
1955  "HitEfficiency_L2", "Hit Efficiency in Barrel_Layer2;Module;Ladder", 8, -4., 4., 28, -14., 14.);
1956  HitEfficiency_L3 = iBooker.book2D(
1957  "HitEfficiency_L3", "Hit Efficiency in Barrel_Layer3;Module;Ladder", 8, -4., 4., 44, -22., 22.);
1958  HitEfficiency_L4 = iBooker.book2D(
1959  "HitEfficiency_L4", "Hit Efficiency in Barrel_Layer4;Module;Ladder", 8, -4., 4., 64, -32., 32.);
1960  }
1961  } // endifUpgrade
1962  // Endcap
1963  iBooker.cd();
1964  iBooker.setCurrentFolder("Pixel/Endcap");
1965  if (!isUpgrade) {
1966  if (Tier0Flag_) {
1967  HitEfficiency_Dp1 = iBooker.book2D(
1968  "HitEfficiency_Dp1", "Hit Efficiency in Endcap_Disk_p1;Blade;Panel", 26, -13., 13., 2, 0.5, 2.5);
1969  HitEfficiency_Dp2 = iBooker.book2D(
1970  "HitEfficiency_Dp2", "Hit Efficiency in Endcap_Disk_p2;Blade;Panel", 26, -13., 13., 2, 0.5, 2.5);
1971  HitEfficiency_Dm1 = iBooker.book2D(
1972  "HitEfficiency_Dm1", "Hit Efficiency in Endcap_Disk_m1;Blade;Panel", 26, -13., 13., 2, 0.5, 2.5);
1973  HitEfficiency_Dm2 = iBooker.book2D(
1974  "HitEfficiency_Dm2", "Hit Efficiency in Endcap_Disk_m2;;Blade;Panel", 26, -13., 13., 2, 0.5, 2.5);
1975  } else {
1976  HitEfficiency_Dp1 = iBooker.book2D(
1977  "HitEfficiency_Dp1", "Hit Efficiency in Endcap_Disk_p1;Blades;Modules", 24, -12., 12., 7, 1., 8.);
1978  HitEfficiency_Dp2 = iBooker.book2D(
1979  "HitEfficiency_Dp2", "Hit Efficiency in Endcap_Disk_p2;Blades;Modules", 24, -12., 12., 7, 1., 8.);
1980  HitEfficiency_Dm1 = iBooker.book2D(
1981  "HitEfficiency_Dm1", "Hit Efficiency in Endcap_Disk_m1;Blades;Modules", 24, -12., 12., 7, 1., 8.);
1982  HitEfficiency_Dm2 = iBooker.book2D(
1983  "HitEfficiency_Dm2", "Hit Efficiency in Endcap_Disk_m2;Blades;Modules", 24, -12., 12., 7, 1., 8.);
1984  }
1985  } else if (isUpgrade) {
1986  if (Tier0Flag_) {
1987  HitEfficiency_Dp1 =
1988  iBooker.book2D("HitEfficiency_Dp1", "Hit Efficiency in Endcap_Disk_p1;Blades;", 28, -17., 11., 1, 0., 1.);
1989  HitEfficiency_Dp2 =
1990  iBooker.book2D("HitEfficiency_Dp2", "Hit Efficiency in Endcap_Disk_p2;Blades;", 28, -17., 11., 1, 0., 1.);
1991  HitEfficiency_Dp3 =
1992  iBooker.book2D("HitEfficiency_Dp3", "Hit Efficiency in Endcap_Disk_p3;Blades;", 28, -17., 11., 1, 0., 1.);
1993  HitEfficiency_Dm1 =
1994  iBooker.book2D("HitEfficiency_Dm1", "Hit Efficiency in Endcap_Disk_m1;Blades;", 28, -17., 11., 1, 0., 1.);
1995  HitEfficiency_Dm2 =
1996  iBooker.book2D("HitEfficiency_Dm2", "Hit Efficiency in Endcap_Disk_m2;Blades;", 28, -17., 11., 1, 0., 1.);
1997  HitEfficiency_Dm3 =
1998  iBooker.book2D("HitEfficiency_Dm3", "Hit Efficiency in Endcap_Disk_m3;Blades;", 28, -17., 11., 1, 0., 1.);
1999  } else {
2000  HitEfficiency_Dp1 = iBooker.book2D(
2001  "HitEfficiency_Dp1", "Hit Efficiency in Endcap_Disk_p1;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2002  HitEfficiency_Dp2 = iBooker.book2D(
2003  "HitEfficiency_Dp2", "Hit Efficiency in Endcap_Disk_p2;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2004  HitEfficiency_Dp3 = iBooker.book2D(
2005  "HitEfficiency_Dp3", "Hit Efficiency in Endcap_Disk_p3;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2006  HitEfficiency_Dm1 = iBooker.book2D(
2007  "HitEfficiency_Dm1", "Hit Efficiency in Endcap_Disk_m1;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2008  HitEfficiency_Dm2 = iBooker.book2D(
2009  "HitEfficiency_Dm2", "Hit Efficiency in Endcap_Disk_m2;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2010  HitEfficiency_Dm3 = iBooker.book2D(
2011  "HitEfficiency_Dm3", "Hit Efficiency in Endcap_Disk_m3;Blades;Modules", 28, -17., 11., 2, 1., 3.);
2012  }
2013  } // endif(isUpgrade)
2014  iBooker.cd();
2015  iBooker.cd("Pixel/");
2016  string bins[] = {"Layer1", "Layer2", "Layer3", "Disk1+", "Disk2+", "Disk1-", "Disk2-"};
2017  HitEfficiencySummary = iBooker.book1D("HitEfficiencySummary", "Hit efficiency per sub detector", 7, 0, 7);
2018  HitEfficiencySummary->setAxisTitle("Sub Detector", 1);
2019  HitEfficiencySummary->setAxisTitle("Efficiency (%)", 2);
2020  for (int i = 1; i < 8; i++) {
2021  HitEfficiencySummary->setBinLabel(i, bins[i - 1]);
2022  }
2023 }
2024 
2025 //=============================================================================================================
2026 
2027 void SiPixelActionExecutor::createEfficiency(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isUpgrade) {
2028  // std::cout<<"entering
2029  // SiPixelActionExecutor::createEfficiency..."<<std::endl;
2030  iGetter.cd();
2031  iBooker.cd();
2032  fillEfficiency(iBooker, iGetter, true, isUpgrade); // Barrel
2033  iGetter.cd();
2034  iBooker.cd();
2035  fillEfficiency(iBooker, iGetter, false, isUpgrade); // Endcap
2036  iGetter.cd();
2037  iBooker.cd();
2038  // std::cout<<"leaving SiPixelActionExecutor::createEfficiency..."<<std::endl;
2039 }
2040 
2041 //=============================================================================================================
2042 
2043 int SiPixelActionExecutor::getLadder(const std::string &dname_) {
2044  int biny_ = 0;
2045  string lad = dname_.substr(dname_.find("Ladder_") + 7, 2);
2046  if (dname_.find(lad) != string::npos) {
2047  biny_ = atoi(lad.c_str());
2048  }
2049  return biny_;
2050 }
2051 
2052 //=============================================================================================================
2053 
2054 int SiPixelActionExecutor::getBlade(const std::string &dname_) {
2055  int binx_ = 0;
2056  string blad = dname_.substr(dname_.find("Blade_") + 6, 2);
2057  if (dname_.find(blad) != string::npos) {
2058  binx_ = atoi(blad.c_str());
2059  }
2060  return binx_;
2061 }
2062 
2063 //=============================================================================================================
2064 
2065 void SiPixelActionExecutor::fillEfficiency(DQMStore::IBooker &iBooker,
2066  DQMStore::IGetter &iGetter,
2067  bool isbarrel,
2068  bool isUpgrade) {
2069  // cout<<"entering SiPixelActionExecutor::fillEfficiency..."<<std::endl;
2070  string currDir = iBooker.pwd();
2071  string dname = currDir.substr(currDir.find_last_of('/') + 1);
2072  // cout<<"currDir= "<<currDir<< " , dname= "<<dname<<std::endl;
2073 
2074  if (Tier0Flag_) { // Offline
2075  if (isbarrel && dname.find("Ladder_") != string::npos) {
2076  if (!isUpgrade) {
2077  vector<string> meVec = iGetter.getMEs();
2078  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
2079  string full_path = currDir + "/" + (*it);
2080 
2081  if (full_path.find("missingMod_") != string::npos) { // If we have missing hits ME
2082 
2083  // Get the MEs that contain missing and valid hits
2084  MonitorElement *missing = iGetter.get(full_path);
2085  if (!missing)
2086  continue;
2087  string new_path = full_path.replace(full_path.find("missing"), 7, "valid");
2088  MonitorElement *valid = iGetter.get(new_path);
2089  if (!valid)
2090  continue;
2091  int binx = 0;
2092  int biny = 0;
2093  // get the ladder number
2094  biny = getLadder(dname); // Current
2095  if (currDir.find("Shell_mO") != string::npos || currDir.find("Shell_pO") != string::npos) {
2096  biny = -biny;
2097  }
2098  const int nMod = 4;
2099  for (int i = 1; i < nMod + 1; i++) {
2100  float hitEfficiency = -1.0;
2101  float missingHits = 0;
2102  float validHits = 0;
2103  binx = i; // Module
2104  if (currDir.find("Shell_m") != string::npos)
2105  binx = -binx;
2106 
2107  missingHits = missing->getBinContent(i);
2108  validHits = valid->getBinContent(i);
2109  if (validHits + missingHits > 0.)
2110  hitEfficiency = validHits / (validHits + missingHits);
2111 
2112  if (currDir.find("Layer_1") != string::npos) {
2113  HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
2114  if (HitEfficiency_L1)
2115  HitEfficiency_L1->setBinContent(HitEfficiency_L1->getTH2F()->FindBin(binx, biny),
2116  (float)hitEfficiency);
2117  } else if (currDir.find("Layer_2") != string::npos) {
2118  HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
2119  if (HitEfficiency_L2)
2120  HitEfficiency_L2->setBinContent(HitEfficiency_L2->getTH2F()->FindBin(binx, biny),
2121  (float)hitEfficiency);
2122  } else if (currDir.find("Layer_3") != string::npos) {
2123  HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");
2124  if (HitEfficiency_L3)
2125  HitEfficiency_L3->setBinContent(HitEfficiency_L3->getTH2F()->FindBin(binx, biny),
2126  (float)hitEfficiency);
2127  }
2128  }
2129  }
2130  }
2131  } // endifNOTUpgradeInBPix
2132  else if (isUpgrade) {
2133  vector<string> meVec = iGetter.getMEs();
2134  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
2135  string full_path = currDir + "/" + (*it);
2136  if (full_path.find("missing_") != string::npos) { // If we have missing hits ME
2137  MonitorElement *me = iGetter.get(full_path);
2138  if (!me)
2139  continue;
2140  float missingHits = me->getEntries();
2141  string new_path = full_path.replace(full_path.find("missing"), 7, "valid");
2142  me = iGetter.get(new_path);
2143  if (!me)
2144  continue;
2145  float validHits = me->getEntries();
2146  float hitEfficiency = -1.;
2147  if (validHits + missingHits > 0.)
2148  hitEfficiency = validHits / (validHits + missingHits);
2149  int binx = 0;
2150  int biny = 0;
2151  biny = getLadder(dname);
2152  if (currDir.find("Shell_mO") != string::npos || currDir.find("Shell_pO") != string::npos) {
2153  biny = -biny;
2154  }
2155  if (currDir.find("Shell_m") != string::npos) {
2156  binx = 1;
2157  } else {
2158  binx = 2;
2159  } // x-axis: z-side
2160  if (currDir.find("Layer_1") != string::npos) {
2161  HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
2162  if (HitEfficiency_L1)
2163  HitEfficiency_L1->setBinContent(binx, biny, (float)hitEfficiency);
2164  } else if (currDir.find("Layer_2") != string::npos) {
2165  HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
2166  if (HitEfficiency_L2)
2167  HitEfficiency_L2->setBinContent(binx, biny, (float)hitEfficiency);
2168  } else if (currDir.find("Layer_3") != string::npos) {
2169  HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");
2170  if (HitEfficiency_L3)
2171  HitEfficiency_L3->setBinContent(binx, biny, (float)hitEfficiency);
2172  } else if (currDir.find("Layer_4") != string::npos) {
2173  HitEfficiency_L4 = iGetter.get("Pixel/Barrel/HitEfficiency_L4");
2174  if (HitEfficiency_L4)
2175  HitEfficiency_L4->setBinContent(binx, biny, (float)hitEfficiency);
2176  }
2177  }
2178  }
2179  } // endifUpgradeInBPix
2180  } else if (!isbarrel && dname.find("Blade_") != string::npos && !isUpgrade) {
2181  vector<string> meVec = iGetter.getMEs();
2182  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
2183  string full_path = currDir + "/" + (*it);
2184  if (full_path.find("missing_") != string::npos) { // If we have missing hits ME
2185  MonitorElement *missing = iGetter.get(full_path);
2186  if (!missing)
2187  continue;
2188  // float missingHits = missing->getEntries();
2189  string new_path = full_path.replace(full_path.find("missing"), 7, "valid");
2190  MonitorElement *valid = iGetter.get(new_path);
2191  if (!valid)
2192  continue;
2193  // float validHits = valid->getEntries();
2194  int binx = 0;
2195  int biny = 0;
2196  binx = getBlade(dname);
2197  if (currDir.find("HalfCylinder_mI") != string::npos || currDir.find("HalfCylinder_pI") != string::npos) {
2198  binx = binx + 14;
2199  } else {
2200  binx = 13 - binx;
2201  }
2202  const int nPanel = 2;
2203  for (int i = 1; i < nPanel + 1; i++) {
2204  float hitEfficiency = -1.;
2205  float missingHits = 0;
2206  float validHits = 0;
2207  biny = i;
2208  missingHits = missing->getBinContent(i);
2209  validHits = valid->getBinContent(i);
2210  if (validHits + missingHits > 0.)
2211  hitEfficiency = validHits / (validHits + missingHits);
2212  if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2213  HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
2214  if (HitEfficiency_Dm1)
2215  HitEfficiency_Dm1->setBinContent(binx, biny, (float)hitEfficiency);
2216  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2217  HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
2218  if (HitEfficiency_Dm2)
2219  HitEfficiency_Dm2->setBinContent(binx, biny, (float)hitEfficiency);
2220  } else if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2221  HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
2222  if (HitEfficiency_Dp1)
2223  HitEfficiency_Dp1->setBinContent(binx, biny, (float)hitEfficiency);
2224  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2225  HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
2226  if (HitEfficiency_Dp2)
2227  HitEfficiency_Dp2->setBinContent(binx, biny, (float)hitEfficiency);
2228  }
2229  } // EndOfFor
2230  }
2231  }
2232  } else if (!isbarrel && dname.find("Blade_") != string::npos && isUpgrade) {
2233  vector<string> meVec = iGetter.getMEs();
2234  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
2235  string full_path = currDir + "/" + (*it);
2236  if (full_path.find("missing_") != string::npos) { // If we have missing hits ME
2237  MonitorElement *me = iGetter.get(full_path);
2238  if (!me)
2239  continue;
2240  float missingHits = me->getEntries();
2241  string new_path = full_path.replace(full_path.find("missing"), 7, "valid");
2242  me = iGetter.get(new_path);
2243  if (!me)
2244  continue;
2245  float validHits = me->getEntries();
2246  float hitEfficiency = -1.;
2247  if (validHits + missingHits > 0.)
2248  hitEfficiency = validHits / (validHits + missingHits);
2249  int binx = 0;
2250  int biny = 1;
2251  binx = getBlade(dname);
2252  if (currDir.find("HalfCylinder_mI") != string::npos || currDir.find("HalfCylinder_pI") != string::npos) {
2253  binx = binx + 12;
2254  } else {
2255  if (binx == 1)
2256  binx = 17;
2257  else if (binx == 2)
2258  binx = 16;
2259  else if (binx == 3)
2260  binx = 15;
2261  else if (binx == 4)
2262  binx = 14;
2263  else if (binx == 5)
2264  binx = 13;
2265  else if (binx == 6)
2266  binx = 12;
2267  else if (binx == 7)
2268  binx = 11;
2269  else if (binx == 8)
2270  binx = 10;
2271  else if (binx == 9)
2272  binx = 9;
2273  else if (binx == 10)
2274  binx = 8;
2275  else if (binx == 11)
2276  binx = 7;
2277  else if (binx == 12)
2278  binx = 6;
2279  else if (binx == 13)
2280  binx = 5;
2281  else if (binx == 14)
2282  binx = 4;
2283  else if (binx == 15)
2284  binx = 3;
2285  else if (binx == 16)
2286  binx = 2;
2287  else if (binx == 17)
2288  binx = 1;
2289  }
2290  if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2291  HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
2292  if (HitEfficiency_Dm1)
2293  HitEfficiency_Dm1->setBinContent(binx, biny, (float)hitEfficiency);
2294  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2295  HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
2296  if (HitEfficiency_Dm2)
2297  HitEfficiency_Dm2->setBinContent(binx, biny, (float)hitEfficiency);
2298  } else if (currDir.find("Disk_3") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2299  HitEfficiency_Dm3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm3");
2300  if (HitEfficiency_Dm3)
2301  HitEfficiency_Dm3->setBinContent(binx, biny, (float)hitEfficiency);
2302  } else if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2303  HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
2304  if (HitEfficiency_Dp1)
2305  HitEfficiency_Dp1->setBinContent(binx, biny, (float)hitEfficiency);
2306  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2307  HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
2308  if (HitEfficiency_Dp2)
2309  HitEfficiency_Dp2->setBinContent(binx, biny, (float)hitEfficiency);
2310  } else if (currDir.find("Disk_3") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2311  HitEfficiency_Dp3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp3");
2312  if (HitEfficiency_Dp3)
2313  HitEfficiency_Dp3->setBinContent(binx, biny, (float)hitEfficiency);
2314  }
2315  // std::cout<<"EFFI: "<<currDir<<" , x: "<<binx<<" , y:
2316  // "<<biny<<std::endl;
2317  }
2318  }
2319  } else {
2320  vector<string> subdirs = iGetter.getSubdirs();
2321  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
2322  iBooker.cd(*it);
2323  iGetter.cd(*it);
2324  if (*it != "Pixel" &&
2325  ((isbarrel && (*it).find("Barrel") == string::npos) || (!isbarrel && (*it).find("Endcap") == string::npos)))
2326  continue;
2327  fillEfficiency(iBooker, iGetter, isbarrel, isUpgrade);
2328  iBooker.goUp();
2329  iGetter.setCurrentFolder(iBooker.pwd());
2330  }
2331  }
2332  } else { // Online
2333  if (dname.find("Module_") != string::npos) {
2334  vector<string> meVec = iGetter.getMEs();
2335  for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
2336  string full_path = currDir + "/" + (*it);
2337  if (full_path.find("missing_") != string::npos) { // If we have missing hits ME
2338  MonitorElement *me = iGetter.get(full_path);
2339  if (!me)
2340  continue;
2341  float missingHits = me->getEntries();
2342  string new_path = full_path.replace(full_path.find("missing"), 7, "valid");
2343  me = iGetter.get(new_path);
2344  if (!me)
2345  continue;
2346  float validHits = me->getEntries();
2347  float hitEfficiency = -1.;
2348  if (validHits + missingHits > 0.)
2349  hitEfficiency = validHits / (validHits + missingHits);
2350  int binx = 0;
2351  int biny = 0;
2352  if (isbarrel) {
2353  if (currDir.find("Shell_m") != string::npos) {
2354  if (currDir.find("Module_4") != string::npos) {
2355  binx = 1;
2356  } else if (currDir.find("Module_3") != string::npos) {
2357  binx = 2;
2358  }
2359  if (currDir.find("Module_2") != string::npos) {
2360  binx = 3;
2361  } else if (currDir.find("Module_1") != string::npos) {
2362  binx = 4;
2363  }
2364  } else if (currDir.find("Shell_p") != string::npos) {
2365  if (currDir.find("Module_1") != string::npos) {
2366  binx = 5;
2367  } else if (currDir.find("Module_2") != string::npos) {
2368  binx = 6;
2369  }
2370  if (currDir.find("Module_3") != string::npos) {
2371  binx = 7;
2372  } else if (currDir.find("Module_4") != string::npos) {
2373  binx = 8;
2374  }
2375  }
2376  if (!isUpgrade) {
2377  if (currDir.find("01") != string::npos) {
2378  biny = 1;
2379  } else if (currDir.find("02") != string::npos) {
2380  biny = 2;
2381  } else if (currDir.find("03") != string::npos) {
2382  biny = 3;
2383  } else if (currDir.find("04") != string::npos) {
2384  biny = 4;
2385  } else if (currDir.find("05") != string::npos) {
2386  biny = 5;
2387  } else if (currDir.find("06") != string::npos) {
2388  biny = 6;
2389  } else if (currDir.find("07") != string::npos) {
2390  biny = 7;
2391  } else if (currDir.find("08") != string::npos) {
2392  biny = 8;
2393  } else if (currDir.find("09") != string::npos) {
2394  biny = 9;
2395  } else if (currDir.find("10") != string::npos) {
2396  biny = 10;
2397  } else if (currDir.find("11") != string::npos) {
2398  biny = 11;
2399  } else if (currDir.find("12") != string::npos) {
2400  biny = 12;
2401  } else if (currDir.find("13") != string::npos) {
2402  biny = 13;
2403  } else if (currDir.find("14") != string::npos) {
2404  biny = 14;
2405  } else if (currDir.find("15") != string::npos) {
2406  biny = 15;
2407  } else if (currDir.find("16") != string::npos) {
2408  biny = 16;
2409  } else if (currDir.find("17") != string::npos) {
2410  biny = 17;
2411  } else if (currDir.find("18") != string::npos) {
2412  biny = 18;
2413  } else if (currDir.find("19") != string::npos) {
2414  biny = 19;
2415  } else if (currDir.find("20") != string::npos) {
2416  biny = 20;
2417  } else if (currDir.find("21") != string::npos) {
2418  biny = 21;
2419  } else if (currDir.find("22") != string::npos) {
2420  biny = 22;
2421  }
2422  if (currDir.find("Shell_mO") != string::npos || currDir.find("Shell_pO") != string::npos) {
2423  if (currDir.find("Layer_1") != string::npos) {
2424  biny = biny + 10;
2425  } else if (currDir.find("Layer_2") != string::npos) {
2426  biny = biny + 16;
2427  } else if (currDir.find("Layer_3") != string::npos) {
2428  biny = biny + 22;
2429  }
2430  }
2431  } else if (isUpgrade) {
2432  if (currDir.find("01") != string::npos) {
2433  biny = 1;
2434  } else if (currDir.find("02") != string::npos) {
2435  biny = 2;
2436  } else if (currDir.find("03") != string::npos) {
2437  biny = 3;
2438  } else if (currDir.find("04") != string::npos) {
2439  biny = 4;
2440  } else if (currDir.find("05") != string::npos) {
2441  biny = 5;
2442  } else if (currDir.find("06") != string::npos) {
2443  biny = 6;
2444  } else if (currDir.find("07") != string::npos) {
2445  biny = 7;
2446  } else if (currDir.find("08") != string::npos) {
2447  biny = 8;
2448  } else if (currDir.find("09") != string::npos) {
2449  biny = 9;
2450  } else if (currDir.find("10") != string::npos) {
2451  biny = 10;
2452  } else if (currDir.find("11") != string::npos) {
2453  biny = 11;
2454  } else if (currDir.find("12") != string::npos) {
2455  biny = 12;
2456  } else if (currDir.find("13") != string::npos) {
2457  biny = 13;
2458  } else if (currDir.find("14") != string::npos) {
2459  biny = 14;
2460  } else if (currDir.find("15") != string::npos) {
2461  biny = 15;
2462  } else if (currDir.find("16") != string::npos) {
2463  biny = 16;
2464  } else if (currDir.find("17") != string::npos) {
2465  biny = 17;
2466  } else if (currDir.find("18") != string::npos) {
2467  biny = 18;
2468  } else if (currDir.find("19") != string::npos) {
2469  biny = 19;
2470  } else if (currDir.find("20") != string::npos) {
2471  biny = 20;
2472  } else if (currDir.find("21") != string::npos) {
2473  biny = 21;
2474  } else if (currDir.find("22") != string::npos) {
2475  biny = 22;
2476  } else if (currDir.find("23") != string::npos) {
2477  biny = 23;
2478  } else if (currDir.find("24") != string::npos) {
2479  biny = 24;
2480  } else if (currDir.find("25") != string::npos) {
2481  biny = 25;
2482  } else if (currDir.find("25") != string::npos) {
2483  biny = 25;
2484  } else if (currDir.find("26") != string::npos) {
2485  biny = 26;
2486  } else if (currDir.find("27") != string::npos) {
2487  biny = 27;
2488  } else if (currDir.find("28") != string::npos) {
2489  biny = 28;
2490  } else if (currDir.find("29") != string::npos) {
2491  biny = 29;
2492  } else if (currDir.find("30") != string::npos) {
2493  biny = 30;
2494  } else if (currDir.find("31") != string::npos) {
2495  biny = 31;
2496  } else if (currDir.find("32") != string::npos) {
2497  biny = 32;
2498  }
2499  if (currDir.find("Shell_mO") != string::npos || currDir.find("Shell_pO") != string::npos) {
2500  if (currDir.find("Layer_1") != string::npos) {
2501  biny = biny + 6;
2502  } else if (currDir.find("Layer_2") != string::npos) {
2503  biny = biny + 14;
2504  } else if (currDir.find("Layer_3") != string::npos) {
2505  biny = biny + 22;
2506  } else if (currDir.find("Layer_4") != string::npos) {
2507  biny = biny + 32;
2508  }
2509  }
2510  }
2511  } else { // endcap
2512  if (!isUpgrade) {
2513  if (currDir.find("01") != string::npos) {
2514  binx = 1;
2515  } else if (currDir.find("02") != string::npos) {
2516  binx = 2;
2517  } else if (currDir.find("03") != string::npos) {
2518  binx = 3;
2519  } else if (currDir.find("04") != string::npos) {
2520  binx = 4;
2521  } else if (currDir.find("05") != string::npos) {
2522  binx = 5;
2523  } else if (currDir.find("06") != string::npos) {
2524  binx = 6;
2525  } else if (currDir.find("07") != string::npos) {
2526  binx = 7;
2527  } else if (currDir.find("08") != string::npos) {
2528  binx = 8;
2529  } else if (currDir.find("09") != string::npos) {
2530  binx = 9;
2531  } else if (currDir.find("10") != string::npos) {
2532  binx = 10;
2533  } else if (currDir.find("11") != string::npos) {
2534  binx = 11;
2535  } else if (currDir.find("12") != string::npos) {
2536  binx = 12;
2537  }
2538  if (currDir.find("HalfCylinder_mO") != string::npos || currDir.find("HalfCylinder_pO") != string::npos) {
2539  binx = binx + 12;
2540  }
2541  if (currDir.find("Panel_1/Module_1") != string::npos) {
2542  biny = 1;
2543  } else if (currDir.find("Panel_2/Module_1") != string::npos) {
2544  biny = 2;
2545  } else if (currDir.find("Panel_1/Module_2") != string::npos) {
2546  biny = 3;
2547  } else if (currDir.find("Panel_2/Module_2") != string::npos) {
2548  biny = 4;
2549  } else if (currDir.find("Panel_1/Module_3") != string::npos) {
2550  biny = 5;
2551  } else if (currDir.find("Panel_2/Module_3") != string::npos) {
2552  biny = 6;
2553  } else if (currDir.find("Panel_1/Module_4") != string::npos) {
2554  biny = 7;
2555  }
2556  } else if (isUpgrade) {
2557  if (currDir.find("01") != string::npos) {
2558  binx = 1;
2559  } else if (currDir.find("02") != string::npos) {
2560  binx = 2;
2561  } else if (currDir.find("03") != string::npos) {
2562  binx = 3;
2563  } else if (currDir.find("04") != string::npos) {
2564  binx = 4;
2565  } else if (currDir.find("05") != string::npos) {
2566  binx = 5;
2567  } else if (currDir.find("06") != string::npos) {
2568  binx = 6;
2569  } else if (currDir.find("07") != string::npos) {
2570  binx = 7;
2571  } else if (currDir.find("08") != string::npos) {
2572  binx = 8;
2573  } else if (currDir.find("09") != string::npos) {
2574  binx = 9;
2575  } else if (currDir.find("10") != string::npos) {
2576  binx = 10;
2577  } else if (currDir.find("11") != string::npos) {
2578  binx = 11;
2579  } else if (currDir.find("12") != string::npos) {
2580  binx = 12;
2581  } else if (currDir.find("13") != string::npos) {
2582  binx = 13;
2583  } else if (currDir.find("14") != string::npos) {
2584  binx = 14;
2585  } else if (currDir.find("15") != string::npos) {
2586  binx = 15;
2587  } else if (currDir.find("16") != string::npos) {
2588  binx = 16;
2589  } else if (currDir.find("17") != string::npos) {
2590  binx = 17;
2591  }
2592  if (currDir.find("HalfCylinder_mO") != string::npos || currDir.find("HalfCylinder_pO") != string::npos) {
2593  binx = binx + 17;
2594  }
2595  if (currDir.find("Panel_1/Module_1") != string::npos) {
2596  biny = 1;
2597  } else if (currDir.find("Panel_2/Module_1") != string::npos) {
2598  biny = 2;
2599  }
2600  } // endif(isUpgrade)
2601  }
2602 
2603  if (currDir.find("Layer_1") != string::npos) {
2604  HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
2605  if (HitEfficiency_L1)
2606  HitEfficiency_L1->setBinContent(binx, biny, (float)hitEfficiency);
2607  } else if (currDir.find("Layer_2") != string::npos) {
2608  HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
2609  if (HitEfficiency_L2)
2610  HitEfficiency_L2->setBinContent(binx, biny, (float)hitEfficiency);
2611  } else if (currDir.find("Layer_3") != string::npos) {
2612  HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");
2613  if (HitEfficiency_L3)
2614  HitEfficiency_L3->setBinContent(binx, biny, (float)hitEfficiency);
2615  } else if (isUpgrade && (currDir.find("Layer_4") != string::npos)) {
2616  HitEfficiency_L4 = iGetter.get("Pixel/Barrel/HitEfficiency_L4");
2617  if (HitEfficiency_L4)
2618  HitEfficiency_L4->setBinContent(binx, biny, (float)hitEfficiency);
2619  } else if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2620  HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
2621  if (HitEfficiency_Dm1)
2622  HitEfficiency_Dm1->setBinContent(binx, biny, (float)hitEfficiency);
2623  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2624  HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
2625  if (HitEfficiency_Dm2)
2626  HitEfficiency_Dm2->setBinContent(binx, biny, (float)hitEfficiency);
2627  } else if (currDir.find("Disk_3") != string::npos && currDir.find("HalfCylinder_m") != string::npos) {
2628  HitEfficiency_Dm3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm3");
2629  if (HitEfficiency_Dm3)
2630  HitEfficiency_Dm3->setBinContent(binx, biny, (float)hitEfficiency);
2631  } else if (currDir.find("Disk_1") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2632  HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
2633  if (HitEfficiency_Dp1)
2634  HitEfficiency_Dp1->setBinContent(binx, biny, (float)hitEfficiency);
2635  } else if (currDir.find("Disk_2") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2636  HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
2637  if (HitEfficiency_Dp2)
2638  HitEfficiency_Dp2->setBinContent(binx, biny, (float)hitEfficiency);
2639  } else if (currDir.find("Disk_3") != string::npos && currDir.find("HalfCylinder_p") != string::npos) {
2640  HitEfficiency_Dp3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp3");
2641  if (HitEfficiency_Dp3)
2642  HitEfficiency_Dp3->setBinContent(binx, biny, (float)hitEfficiency);
2643  }
2644  }
2645  }
2646  } else {
2647  // cout<<"finding subdirs now"<<std::endl;
2648  vector<string> subdirs = iGetter.getSubdirs();
2649  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
2650  iBooker.cd(*it);
2651  iGetter.cd(*it);
2652  if (*it != "Pixel" &&
2653  ((isbarrel && (*it).find("Barrel") == string::npos) || (!isbarrel && (*it).find("Endcap") == string::npos)))
2654  continue;
2655  fillEfficiency(iBooker, iGetter, isbarrel, isUpgrade);
2656  iBooker.goUp();
2657  iGetter.setCurrentFolder(iBooker.pwd());
2658  }
2659  }
2660  } // end online/offline
2661  // cout<<"leaving SiPixelActionExecutor::fillEfficiency..."<<std::endl;
2662 }
2663 
2664 //=============================================================================================================
2665 
2666 void SiPixelActionExecutor::fillEfficiencySummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) {
2667  // cout<<"entering
2668  // SiPixelActionExecutor::fillEfficiencySummary..."<<std::endl; First we get
2669  // the summary plot"
2670  if (!Tier0Flag_)
2671  return;
2672  HitEfficiencySummary = iGetter.get("Pixel/HitEfficiencySummary");
2673  // Now we will loop over the hit efficiency plots and fill it"
2674  string hitEfficiencyPostfix[] = {"L1", "L2", "L3", "Dp1", "Dp2", "Dm1", "Dm2"};
2675  std::vector<std::vector<float>> ignoreXBins = {
2676  {-4, 2}, {4, 4, -1, -3, 3, -4, -3, -2, -1, -4, -3, -2, -1, 1, -4}, {1, -4, 1}, {}, {}, {}, {}};
2677  std::vector<std::vector<float>> ignoreYBins = {
2678  {-9, -3}, {1, 16, 1, -13, -13, -5, -5, -5, -5, -6, -6, -6, -6, -8, -8}, {3, 14, 6}, {}, {}, {}, {}};
2679 
2680  for (int i = 0; i < 7; i++) {
2681  string subdetName = "Endcap/";
2682  if (i < 3)
2683  subdetName = "Barrel/";
2684  char meName[50];
2685  sprintf(meName, "Pixel/%sHitEfficiency_%s", subdetName.c_str(), hitEfficiencyPostfix[i].c_str());
2686  MonitorElement *tempHitEffMap = iGetter.get(meName);
2687  float totalEff = 0.;
2688  int totalBins = 0;
2689  TH1 *hitEffMap = tempHitEffMap->getTH1();
2690  for (int xBin = 1; xBin < tempHitEffMap->getNbinsX() + 1; xBin++) {
2691  if (fabs(hitEffMap->GetXaxis()->GetBinCenter(xBin)) < 1.)
2692  continue;
2693  for (int yBin = 1; yBin < tempHitEffMap->getNbinsY() + 1; yBin++) {
2694  if (fabs(hitEffMap->GetYaxis()->GetBinCenter(yBin)) < 0.5)
2695  continue;
2696  bool ignoreBin = false;
2697  for (unsigned int j = 0; j < ignoreXBins[i].size(); j++) {
2698  if (hitEffMap->GetXaxis()->GetBinCenter(xBin) == ignoreXBins[i][j] &&
2699  hitEffMap->GetYaxis()->GetBinCenter(yBin) == ignoreYBins[i][j]) {
2700  ignoreBin = true;
2701  break;
2702  }
2703  }
2704  if (ignoreBin)
2705  continue;
2706  if (!(tempHitEffMap->getBinContent(xBin, yBin) < 0.))
2707  totalEff += tempHitEffMap->getBinContent(xBin, yBin);
2708  totalBins++;
2709  }
2710  }
2711  float overalEff = 0.;
2712  // std::cout << i << " " << totalEff << " " << totalBins << std::endl;
2713  if (totalBins > 0)
2714  overalEff = totalEff / float(totalBins);
2715  HitEfficiencySummary->setBinContent(i + 1, overalEff);
2716  }
2717 }
compare_using_db.dir_name
dir_name
Definition: compare_using_db.py:32
SiPixelActionExecutor::normaliseAvDigiOccVsLumi
void normaliseAvDigiOccVsLumi(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, int lumisec)
Definition: SiPixelActionExecutor.cc:1888
relmon_rootfiles_spy.contents
contents
Definition: relmon_rootfiles_spy.py:129
SiPixelOfflineDQM_source_cff.hiRes
hiRes
Definition: SiPixelOfflineDQM_source_cff.py:43
mps_fire.i
i
Definition: mps_fire.py:429
hcallasereventfilter2012_cfi.prefix
prefix
Definition: hcallasereventfilter2012_cfi.py:10
SiPixelActionExecutor::getGrandSummaryME
void getGrandSummaryME(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, int nbin, std::string &me_name, std::vector< MonitorElement *> &mes)
Definition: SiPixelActionExecutor.cc:1580
SiPixelOfflineDQM_client_cff.Tier0Flag
Tier0Flag
Definition: SiPixelOfflineDQM_client_cff.py:17
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksByDensity.h:46
SiPixelActionExecutor::HitEfficiency_Dp3
MonitorElement * HitEfficiency_Dp3
Definition: SiPixelActionExecutor.h:130
trigObjTnPSource_cfi.var
var
Definition: trigObjTnPSource_cfi.py:21
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
castor_dqm_sourceclient_file_cfg.path
path
Definition: castor_dqm_sourceclient_file_cfg.py:37
hlt_dqm_clientPB-live_cfg.me1
me1
Definition: hlt_dqm_clientPB-live_cfg.py:71
SiPixelActionExecutor::createEfficiency
void createEfficiency(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:2027
dqm::implementation::IGetter::getMEs
virtual std::vector< std::string > getMEs() const
Definition: DQMStore.cc:759
dqm::implementation::NavigatorBase::pwd
virtual std::string pwd()
Definition: DQMStore.cc:20
SiPixelActionExecutor::HitEfficiency_L4
MonitorElement * HitEfficiency_L4
Definition: SiPixelActionExecutor.h:127
SiPixelActionExecutor::bookOccupancyPlots
void bookOccupancyPlots(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool hiRes, bool isbarrel)
Definition: SiPixelActionExecutor.cc:1708
SiPixelActionExecutor::DEV_size_Barrel
MonitorElement * DEV_size_Barrel
Definition: SiPixelActionExecutor.h:139
SiPixelActionExecutor::HitEfficiency_L3
MonitorElement * HitEfficiency_L3
Definition: SiPixelActionExecutor.h:126
SiPixelActionExecutor::configParser_
SiPixelConfigParser * configParser_
Definition: SiPixelActionExecutor.h:109
SiPixelActionExecutor::DEV_nclusters_Barrel
MonitorElement * DEV_nclusters_Barrel
Definition: SiPixelActionExecutor.h:138
SiPixelActionExecutor::fillGrandEndcapSummaryHistos
void fillGrandEndcapSummaryHistos(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::vector< std::string > &me_names, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:1261
isum
double isum
Definition: MuScleFitUtils.cc:77
SiPixelConfigParser::getMENamesForEndcapSummary
bool getMENamesForEndcapSummary(std::string &structure_name, std::vector< std::string > &me_names)
Definition: SiPixelConfigParser.cc:78
SiPixelActionExecutor::HitEfficiency_Dm2
MonitorElement * HitEfficiency_Dm2
Definition: SiPixelActionExecutor.h:132
SiPixelActionExecutor::HitEfficiency_L2
MonitorElement * HitEfficiency_L2
Definition: SiPixelActionExecutor.h:125
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
SiPixelActionExecutor::getBlade
int getBlade(const std::string &dname)
Definition: SiPixelActionExecutor.cc:2054
photonAnalyzer_cfi.xBin
xBin
Definition: photonAnalyzer_cfi.py:81
combine.missing
missing
Definition: combine.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
SiPixelActionExecutor::normaliseAvDigiOcc
void normaliseAvDigiOcc(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter)
Definition: SiPixelActionExecutor.cc:1858
SiPixelConfigParser::getCalibType
bool getCalibType(int &u_freq)
Definition: SiPixelConfigParser.cc:154
SiPixelActionExecutor::fillEfficiency
void fillEfficiency(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isbarrel, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:2065
trigObjTnPSource_cfi.bins
bins
Definition: trigObjTnPSource_cfi.py:20
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
dqm::legacy::MonitorElement::getTH2F
virtual TH2F * getTH2F() const
Definition: MonitorElement.h:499
photonAnalyzer_cfi.yBin
yBin
Definition: photonAnalyzer_cfi.py:85
SiPixelConfigParser::getMENamesForGrandEndcapSummary
bool getMENamesForGrandEndcapSummary(std::string &structure_name, std::vector< std::string > &me_names)
Definition: SiPixelConfigParser.cc:118
SiPixelActionExecutor::configWriter_
SiPixelConfigWriter * configWriter_
Definition: SiPixelActionExecutor.h:110
SiPixelConfigParser::getFrequencyForGrandBarrelSummary
bool getFrequencyForGrandBarrelSummary(int &u_freq)
Definition: SiPixelConfigParser.cc:126
LogMessageMonitor_cff.modules
modules
Definition: LogMessageMonitor_cff.py:7
SiPixelConfigParser::getFrequencyForEndcapSummary
bool getFrequencyForEndcapSummary(int &u_freq)
Definition: SiPixelConfigParser.cc:103
SiPixelConfigParser::getMENamesForFEDErrorSummary
bool getMENamesForFEDErrorSummary(std::string &structure_name, std::vector< std::string > &me_names)
Definition: SiPixelConfigParser.cc:86
validateGeometry_cfg.valid
valid
Definition: validateGeometry_cfg.py:21
makeGlobalPositionRcd_cfg.tag
tag
Definition: makeGlobalPositionRcd_cfg.py:6
SiPixelActionExecutor::bookDeviations
void bookDeviations(DQMStore::IBooker &iBooker, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:230
SiPixelOfflineDQM_client_cff.isUpgrade
isUpgrade
Definition: SiPixelOfflineDQM_client_cff.py:19
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
SiPixelActionExecutor::createSummary
void createSummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:161
SiPixelActionExecutor::DEV_adc_Endcap
MonitorElement * DEV_adc_Endcap
Definition: SiPixelActionExecutor.h:140
hlt_dqm_clientPB-live_cfg.me
me
Definition: hlt_dqm_clientPB-live_cfg.py:62
SiPixelConfigParser::getMENamesForGrandBarrelSummary
bool getMENamesForGrandBarrelSummary(std::string &structure_name, std::vector< std::string > &me_names)
Definition: SiPixelConfigParser.cc:110
SiPixelConfigParser::getFrequencyForBarrelSummary
bool getFrequencyForBarrelSummary(int &u_freq)
Definition: SiPixelConfigParser.cc:96
SiPixelActionExecutor::fillFEDErrorSummary
void fillFEDErrorSummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::string dir_name, std::vector< std::string > &me_names)
Definition: SiPixelActionExecutor.cc:687
SiPixelConfigParser::getFrequencyForGrandEndcapSummary
bool getFrequencyForGrandEndcapSummary(int &u_freq)
Definition: SiPixelConfigParser.cc:133
SiPixelActionExecutor::DEV_charge_Endcap
MonitorElement * DEV_charge_Endcap
Definition: SiPixelActionExecutor.h:142
SiPixelConfigParser::getMENamesForBarrelSummary
bool getMENamesForBarrelSummary(std::string &structure_name, std::vector< std::string > &me_names)
Definition: SiPixelConfigParser.cc:71
dqm::impl::MonitorElement::setBinLabel
virtual 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:772
SiPixelActionExecutor::getLadder
int getLadder(const std::string &dname)
Definition: SiPixelActionExecutor.cc:2043
SiPixelActionExecutor::getSummaryME
MonitorElement * getSummaryME(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::string me_name, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:1619
SiPixelActionExecutor::fillOccupancy
void fillOccupancy(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, bool isbarrel)
Definition: SiPixelActionExecutor.cc:1802
SiPixelActionExecutor::getFEDSummaryME
MonitorElement * getFEDSummaryME(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::string me_name)
Definition: SiPixelActionExecutor.cc:1680
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
SiPixelActionExecutor::SiPixelActionExecutor
SiPixelActionExecutor(bool offlineXMLfile, bool Tier0Flag)
Definition: SiPixelActionExecutor.cc:27
SiPixelActionExecutor::DEV_charge_Barrel
MonitorElement * DEV_charge_Barrel
Definition: SiPixelActionExecutor.h:137
PixelSLinkDataInputSource_cfi.fedid
fedid
Definition: PixelSLinkDataInputSource_cfi.py:6
SiPixelActionExecutor::HitEfficiency_Dm3
MonitorElement * HitEfficiency_Dm3
Definition: SiPixelActionExecutor.h:133
dqm::impl::MonitorElement::setBinContent
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:692
SiPixelConfigParser::getMessageLimitForQTests
bool getMessageLimitForQTests(int &u_freq)
Definition: SiPixelConfigParser.cc:140
SiPixelActionExecutor::HitEfficiency_L1
MonitorElement * HitEfficiency_L1
Definition: SiPixelActionExecutor.h:124
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
SiPixelActionExecutor::DEV_nclusters_Endcap
MonitorElement * DEV_nclusters_Endcap
Definition: SiPixelActionExecutor.h:143
SiPixelActionExecutor::HitEfficiency_Dp1
MonitorElement * HitEfficiency_Dp1
Definition: SiPixelActionExecutor.h:128
L1TdeStage2BMTF_cfi.ignoreBin
ignoreBin
Definition: L1TdeStage2BMTF_cfi.py:18
dqm::implementation::IGetter::get
virtual MonitorElement * get(std::string const &fullpath) const
Definition: DQMStore.cc:712
SiPixelConfigParser::getSourceType
bool getSourceType(int &u_freq)
Definition: SiPixelConfigParser.cc:147
dqm::legacy::MonitorElement::getTH1
virtual TH1 * getTH1() const
Definition: MonitorElement.h:479
SiPixelActionExecutor::DEV_size_Endcap
MonitorElement * DEV_size_Endcap
Definition: SiPixelActionExecutor.h:144
nlast
int nlast
Definition: AMPTWrapper.h:57
SiPixelActionExecutor::HitEfficiencySummary
MonitorElement * HitEfficiencySummary
Definition: SiPixelActionExecutor.h:134
tname
std::string tname(const std::string &tableName, const std::string &schemaVersion)
Definition: CredentialStore.cc:150
SiPixelActionExecutor::fillEfficiencySummary
void fillEfficiencySummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter)
Definition: SiPixelActionExecutor.cc:2666
SiPixelActionExecutor::createOccupancy
void createOccupancy(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter)
Definition: SiPixelActionExecutor.cc:1784
SiPixelActionExecutor::HitEfficiency_Dp2
MonitorElement * HitEfficiency_Dp2
Definition: SiPixelActionExecutor.h:129
SiPixelConfigParser::getDocument
void getDocument(std::string filename)
Definition: SiPixelConfigParser.cc:16
SiPixelConfigParser::getFrequencyForTrackerMap
bool getFrequencyForTrackerMap(int &u_freq)
Definition: SiPixelConfigParser.cc:52
gather_cfg.cout
cout
Definition: gather_cfg.py:144
SiPixelActionExecutor::fillGrandBarrelSummaryHistos
void fillGrandBarrelSummaryHistos(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::vector< std::string > &me_names, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:879
SiPixelActionExecutor::fillSummary
void fillSummary(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::string dir_name, std::vector< std::string > &me_names, bool isbarrel, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:287
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
SiPixelActionExecutor::DEV_ndigis_Endcap
MonitorElement * DEV_ndigis_Endcap
Definition: SiPixelActionExecutor.h:141
SiPixelActionExecutor::GetBladeSubdirs
void GetBladeSubdirs(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, std::vector< std::string > &blade_subdirs)
Definition: SiPixelActionExecutor.cc:269
mysort.rows
rows
Definition: mysort.py:12
SiPixelActionExecutor::HitEfficiency_Dm1
MonitorElement * HitEfficiency_Dm1
Definition: SiPixelActionExecutor.h:131
SiPixelActionExecutor::DEV_ndigis_Barrel
MonitorElement * DEV_ndigis_Barrel
Definition: SiPixelActionExecutor.h:136
hlt_dqm_clientPB-live_cfg.me2
me2
Definition: hlt_dqm_clientPB-live_cfg.py:76
SiPixelActionExecutor::bookEfficiency
void bookEfficiency(DQMStore::IBooker &iBooker, bool isUpgrade)
Definition: SiPixelActionExecutor.cc:1920
SiPixelActionExecutor::DEV_adc_Barrel
MonitorElement * DEV_adc_Barrel
Definition: SiPixelActionExecutor.h:135
dqm::impl::MonitorElement::getBinContent
virtual double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:593
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:801
dqm::implementation::IGetter::getSubdirs
virtual DQM_DEPRECATED std::vector< std::string > getSubdirs() const
Definition: DQMStore.cc:739
Generated for CMSSW Reference Manual by doxygen 1.8.14