30 " Creating SiStripQualityChecker " <<
"\n" ;
35 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TIB",
"TIB"));
36 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TOB",
"TOB"));
37 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TECF",
"TEC/side_2"));
38 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TECB",
"TEC/side_1"));
39 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TIDF",
"TID/side_2"));
40 SubDetFolderMap.insert(std::pair<std::string, std::string>(
"TIDB",
"TID/side_1"));
45 "\n------------------------------------------"
46 "\nUnAvailable Service TkHistoMap: please insert in the configuration file an instance like"
47 "\n\tprocess.TkDetMap = cms.Service(\"TkDetMap\")"
48 "\n------------------------------------------";
60 TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>(
"Rate", tracking_mes));
67 TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>(
"Chi2", tracking_mes));
74 TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>(
"RecHits", tracking_mes));
82 " Deleting SiStripQualityChecker " <<
"\n" ;
91 std::string strip_dir =
"";
93 if (strip_dir.size() == 0) strip_dir =
"SiStrip";
99 std::string hname, htitle;
100 hname =
"detFractionReportMap";
101 htitle =
"SiStrip Report for Good Detector Fraction";
105 hname =
"sToNReportMap";
106 htitle =
"SiStrip Report for Signal-to-Noise";
114 hname =
"reportSummaryMap";
115 htitle =
"SiStrip Report Summary Map";
124 for (std::map<std::string, std::string>::const_iterator it =
SubDetFolderMap.begin();
127 std::string det = it->first;
135 else if (det ==
"TECB") local_mes.
detectorTag =
"TEC-";
136 else if (det ==
"TIDF") local_mes.
detectorTag =
"TID+";
137 else if (det ==
"TIDB") local_mes.
detectorTag =
"TID-";
141 me_name =
"SiStrip_" + det;
144 me_name =
"SiStrip_DetFraction_" + det;
147 me_name =
"SiStrip_SToNFlag_" + det;
149 SubDetMEsMap.insert(std::pair<std::string, SubDetMEs>(det, local_mes));
155 std::string tracking_dir =
"";
157 if (tracking_dir.size() == 0) tracking_dir =
"Tracking";
162 std::string hname, htitle;
163 hname =
"reportSummaryMap";
164 htitle =
"Tracking Report Summary Map";
173 dqm_store->
setCurrentFolder(tracking_dir+
"/EventInfo/reportSummaryContents");
175 for (std::map<std::string, TrackingMEs>::iterator it =
TrackingMEsMap.begin();
178 std::string
name = it->first;
179 it->second.TrackingFlag = dqm_store->
bookFloat(
"Track"+name);
193 for (std::map<std::string, SubDetMEs>::const_iterator it =
SubDetMEsMap.begin();
217 for (std::map<std::string, TrackingMEs>::iterator it =
TrackingMEsMap.begin();
219 it->second.TrackingFlag->Fill(-1.0);
228 for (std::map<std::string, SubDetMEs>::const_iterator it =
SubDetMEsMap.begin();
245 for (std::map<std::string, TrackingMEs>::iterator it =
TrackingMEsMap.begin();
247 it->second.TrackingFlag->Reset();
268 unsigned int xbin = 0;
269 float global_flag = 0;
271 std::string mdir =
"MechanicalView";
273 std::string mechanicalview_dir = dqm_store->
pwd();
276 for (std::map<std::string, SubDetMEs>::const_iterator it =
SubDetMEsMap.begin();
278 std::string det = it->first;
279 std::map<std::string, std::string>::const_iterator cPos =
SubDetFolderMap.find(det);
281 std::string dname = mechanicalview_dir +
"/" + cPos->second;
282 if (!dqm_store->
dirExists(dname))
continue;
283 dqm_store->
cd(dname);
290 global_flag = global_flag/xbin*1.0;
300 std::string
dir =
"Tracking";
302 dir =
"TrackParameters";
305 std::vector<MonitorElement*> meVec1 = dqm_store->
getContents(dqm_store->
pwd()+
"/GeneralProperties");
306 std::vector<MonitorElement*> meVec2 = dqm_store->
getContents(dqm_store->
pwd()+
"/HitProperties");
307 std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size());
308 std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
311 for (std::vector<MonitorElement*>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
314 std::vector<QReport *> qt_reports = me->
getQReports();
315 if (qt_reports.size() == 0)
continue;
321 for (std::map<std::string, TrackingMEs>::const_iterator it =
TrackingMEsMap.begin();
324 std::string hname = it->second.HistoName;
325 if (name.find(hname) != std::string::npos) {
326 status = qt_reports[0]->getQTresult();
327 it->second.TrackingFlag->Fill(status);
332 gstatus = gstatus *
status;
342 SubDetMEs& mes,
unsigned int xbin,
float& gflag) {
345 if (status_flag < 1)
return;
347 std::vector<std::string> subDirVec = dqm_store->
getSubdirs();
349 unsigned int ybin = 0;
352 float tot_ston_stat = 0;
354 for (std::vector<std::string>::const_iterator ic = subDirVec.begin();
355 ic != subDirVec.end(); ic++) {
356 std::string dname = (*ic);
357 if (dname.find(
"BadModuleList") != std::string::npos)
continue;
358 std::vector<MonitorElement*> meVec;
360 dqm_store->
cd((*ic));
366 int lnum = atoi(dname.substr(dname.find_last_of(
"_")+1).c_str());
367 ndet = cabling->connectedNumber(mes.
detectorTag, lnum);
371 for (std::vector<MonitorElement*>::const_iterator it = meVec.begin();
372 it != meVec.end(); it++) {
375 std::vector<QReport *> reports = me->
getQReports();
377 if (reports.size() == 0)
continue;
380 if( name.find(
"Summary_ClusterStoNCorr__OnTrack") != std::string::npos){
382 if (reports[0]->getQTresult() == -1) {
391 float eff_fac = 1 - (errdet*1.0/ndet);
398 tot_errdet += errdet;
399 tot_ston_stat += ston_stat;
401 dqm_store->
cd((*ic));
404 float tot_eff_fac = 1 - (tot_errdet*1.0/tot_ndet);
406 float tot_ston_fac = tot_ston_stat/ybin;
408 gflag =
std::min(tot_eff_fac,tot_ston_fac);
416 std::ostringstream det_summary_str;
417 for (std::map<std::string, SubDetMEs>::const_iterator it =
SubDetMEsMap.begin();
419 std::string det = it->first;
420 det_summary_str << std::setprecision(4);
421 det_summary_str << std::setiosflags(std::ios::fixed);
423 det_summary_str <<
" Printing Status for " << det <<
" : " << std::endl;
427 float fval1, fval2, fval3;
428 fval1 = fval2 = fval3 = -1.0;
432 if (sval.size() > 0) fval1 = atof(sval.c_str());
434 if (sval.size() > 0) fval2 = atof(sval.c_str());
436 if (sval.size() > 0) fval3 = atof(sval.c_str());
438 det_summary_str << std::setw(7) <<
" % of good detectors " << fval1
439 <<
" SToN Flag " << fval2
440 <<
" Summary Flag " << fval3 << std::endl;
442 std::cout << det_summary_str.str() << std::endl;
450 std::map<uint32_t,uint16_t> bad_modules;
451 for (std::vector<MonitorElement*>::const_iterator it = layer_mes.begin();
452 it != layer_mes.end(); it++) {
455 std::vector<QReport *> qreports = me->
getQReports();
456 if (qreports.size() == 0)
continue;
458 std::vector<DQMChannel> bad_channels_me;
460 bad_channels_me = qreports[0]->getBadChannels();
463 bad_channels_me = qreports[0]->getBadChannels();
464 lname = name.substr(name.find(
"TkHMap_")+7);
465 lname = lname.substr(lname.find(
"_T")+1);
468 for (std::vector<DQMChannel>::iterator it = bad_channels_me.begin(); it != bad_channels_me.end(); it++){
469 int xval = (*it).getBinX();
470 int yval = (*it).getBinY();
472 std::map<uint32_t,uint16_t>::iterator iPos = bad_modules.find(detId);
474 if (iPos != bad_modules.end()){
481 bad_modules.insert(std::pair<uint32_t,uint16_t>(detId,flag));
485 for(std::map<uint32_t,uint16_t>::const_iterator it = bad_modules.begin();
486 it != bad_modules.end(); it++) {
487 uint32_t detId = it->first;
488 uint16_t
flag = it->second;
489 std::map<uint32_t,uint16_t>::iterator iPos =
badModuleList.find(detId);
493 badModuleList.insert(std::pair<uint32_t,uint16_t>(detId,flag));
496 errdet = bad_modules.size();
504 th2d->SetBinContent(xbin, ybin, val);
513 std::string mdir =
"MechanicalView";
515 std::string mechanical_dir = dqm_store->
pwd();
519 uint32_t detId = it->first;
520 std::string subdet_folder ;
522 if (!dqm_store->
dirExists(subdet_folder)) {
523 subdet_folder = mechanical_dir + subdet_folder.substr(subdet_folder.find(
"MechanicalView")+14);
524 if (!dqm_store->
dirExists(subdet_folder))
continue;
526 std::string bad_module_folder = subdet_folder +
"/" +
"BadModuleList";
529 std::ostringstream detid_str;
531 std::string full_path = bad_module_folder +
"/" + detid_str.str();
534 else me = dqm_store->
bookInt(detid_str.str());
535 me->
Fill(it->second);
561 std::string rdir =
"ReadoutView";
563 std::string fullpath = dqm_store->
pwd()
564 +
"/FedMonitoringSummary/"
565 +
"lumiErrorFraction";
569 float global_fraction = 0.0;
571 for (
int ibin = 1; ibin <= th1->GetNbinsX(); ibin++) {
572 std::string
label = th1->GetXaxis()->GetBinLabel(ibin);
573 std::map<std::string, SubDetMEs>::iterator iPos =
SubDetMEsMap.find(label);
575 float fraction = 1.0 - th1->GetBinContent(ibin);
576 global_fraction += fraction;
578 iPos->second.DetFraction->Fill(fraction);
579 iPos->second.SToNFlag->Fill(-1.0);
580 iPos->second.SummaryFlag->Fill(fraction);
583 global_fraction = global_fraction/dets;
593 std::string
dir =
"Tracking";
595 dir =
"TrackParameters";
598 std::vector<MonitorElement*> meVec1 = dqm_store->
getContents(dqm_store->
pwd()+
"/GeneralProperties");
599 std::vector<MonitorElement*> meVec2 = dqm_store->
getContents(dqm_store->
pwd()+
"/HitProperties");
600 std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size());
601 std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
603 for (std::vector<MonitorElement*>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
610 for (std::map<std::string, TrackingMEs>::const_iterator it =
TrackingMEsMap.begin();
613 std::string hname = it->second.HistoName;
614 float lower_cut = it->second.LowerCut;
615 float upper_cut = it->second.UpperCut;
616 if (name.find(hname) != std::string::npos) {
617 if (me->
getMean() <= lower_cut || me->
getMean() > upper_cut) status = 0.0;
619 it->second.TrackingFlag->Fill(status);
624 if (status == -1.0) gstatus = -1.0;
625 else gstatus = gstatus *
status;
static void setBadModuleFlag(std::string &hname, uint16_t &flg)
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
const std::string & getName(void) const
get name of ME
MonitorElement * SToNFlag
uint32_t getDetFromBin(int layer, int ix, int iy)
bool bookedTrackingStatus_
void getSubDetFolder(const uint32_t &detid, std::string &folder_name)
void bookStatus(DQMStore *dqm_store)
void fillDetectorStatusAtLumi(DQMStore *dqm_store)
std::vector< std::string > getSubdirs(void) const
void fillDetectorStatus(DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling)
static void getMEValue(MonitorElement *me, std::string &val)
std::map< std::string, std::string > SubDetFolderMap
void cd(void)
go to top directory (ie. root)
MonitorElement * SummaryFlag
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)
void fillTrackingStatusAtLumi(DQMStore *dqm_store)
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
void initialiseBadModuleList()
MonitorElement * bookFloat(const char *name)
Book float.
SiStripQualityChecker(edm::ParameterSet const &ps)
int getNbinsY(void) const
get # of bins in Y-axis
std::map< std::string, SubDetMEs > SubDetMEsMap
virtual ~SiStripQualityChecker()
void getModuleStatus(std::vector< MonitorElement * > &layer_mes, int &errdet)
void fillStatus(DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling)
Kind kind(void) const
Get the type of the monitor element.
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. "my/long/dir/my_histo")
std::vector< MonitorElement * > getContents(const std::string &path) const
MonitorElement * DetFractionReportMap
MonitorElement * TrackingFlag
bool dirExists(const std::string &path) const
true if directory exists
std::vector< QReport * > getQReports(void) const
get map of QReports
static int getMEStatus(MonitorElement *me)
MonitorElement * TrackSummaryReportGlobal
TH1F * getTH1F(void) const
void fillStatusAtLumi(DQMStore *dqm_store)
MonitorElement * SummaryReportMap
MonitorElement * SummaryReportGlobal
MonitorElement * SToNReportMap
void fillTrackingStatus(DQMStore *dqm_store)
int getNbinsX(void) const
get # of bins in X-axis
static const int STATUS_OK
MonitorElement * TrackSummaryReportMap
std::map< uint32_t, uint16_t > badModuleList
std::map< std::string, TrackingMEs > TrackingMEsMap
TH2F * getTH2F(void) const
MonitorElement * bookInt(const char *name)
Book int.
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void Reset(void)
reset ME (ie. contents, errors, etc)
void fillSubDetStatus(DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling, SubDetMEs &mes, unsigned int xbin, float &gflag)
static bool goToDir(DQMStore *dqm_store, std::string name)
MonitorElement * DetFraction
void setCurrentFolder(const std::string &fullpath)
void fillFaultyModuleStatus(DQMStore *dqm_store)
void fillStatusHistogram(MonitorElement *, int xbin, int ybin, float val)
const std::string & pwd(void) const
static void getTopFolderPath(DQMStore *dqm_store, std::string type, std::string &path)