|
|
|
#include <SiStripQualityChecker.h>
Definition at line 20 of file SiStripQualityChecker.h.
| SiStripQualityChecker::SiStripQualityChecker | ( | edm::ParameterSet const & | ps | ) |
Definition at line 28 of file SiStripQualityChecker.cc.
References badModuleList, bookedStripStatus_, bookedTrackingStatus_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SiStripQualityChecker::TrackingMEs::HistoName, SiStripQualityChecker::TrackingMEs::LowerCut, cmsCodeRules::cppFunctionSkipper::operator, pSet_, SubDetFolderMap, tkDetMap_, SiStripQualityChecker::TrackingMEs::TrackingFlag, TrackingMEsMap, SiStripQualityChecker::TrackingMEs::UpperCut, and useGoodTracks_.
:pSet_(ps) { edm::LogInfo("SiStripQualityChecker") << " Creating SiStripQualityChecker " << "\n" ; bookedStripStatus_ = false; bookedTrackingStatus_ = false; SubDetFolderMap.insert(std::pair<std::string, std::string>("TIB", "TIB")); SubDetFolderMap.insert(std::pair<std::string, std::string>("TOB", "TOB")); SubDetFolderMap.insert(std::pair<std::string, std::string>("TECF", "TEC/side_2")); SubDetFolderMap.insert(std::pair<std::string, std::string>("TECB", "TEC/side_1")); SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDF", "TID/side_2")); SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDB", "TID/side_1")); badModuleList.clear(); if(!edm::Service<TkDetMap>().isAvailable()){ edm::LogError("TkHistoMap") << "\n------------------------------------------" "\nUnAvailable Service TkHistoMap: please insert in the configuration file an instance like" "\n\tprocess.TkDetMap = cms.Service(\"TkDetMap\")" "\n------------------------------------------"; } tkDetMap_=edm::Service<TkDetMap>().operator->(); TrackingMEs tracking_mes; edm::ParameterSet TkPSet; TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackRatePSet"); tracking_mes.TrackingFlag = 0; tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("Rate", tracking_mes)); TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackChi2PSet"); tracking_mes.TrackingFlag = 0; tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("Chi2", tracking_mes)); TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackHitPSet"); tracking_mes.TrackingFlag = 0; tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("RecHits", tracking_mes)); // // LS analysis // TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackRateLSPSet"); // tracking_mes.TrackingFlag = 0; // tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); // tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); // tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); // TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("Rate", tracking_mes)); // // TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackChi2LSPSet"); // tracking_mes.TrackingFlag = 0; // tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); // tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); // tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); // TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("Chi2", tracking_mes)); // // TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackHitLSPSet"); // tracking_mes.TrackingFlag = 0; // tracking_mes.HistoName = TkPSet.getParameter<std::string>("Name"); // tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut"); // tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut"); // TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("RecHits", tracking_mes)); useGoodTracks_ = pSet_.getUntrackedParameter<bool>("UseGoodTracks", false); if (useGoodTracks_) edm::LogInfo("SiStripQualityChecker") << " use GoodTrack histograms for certification " << "\n" ; }
| SiStripQualityChecker::~SiStripQualityChecker | ( | ) | [virtual] |
Definition at line 104 of file SiStripQualityChecker.cc.
{
edm::LogInfo("SiStripQualityChecker") <<
" Deleting SiStripQualityChecker " << "\n" ;
}
| void SiStripQualityChecker::bookStatus | ( | DQMStore * | dqm_store | ) |
Definition at line 111 of file SiStripQualityChecker.cc.
References DQMStore::book2D(), bookedStripStatus_, bookedTrackingStatus_, DQMStore::bookFloat(), DQMStore::cd(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, SiStripUtility::getTopFolderPath(), mergeVDriftHistosByStation::name, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetFolderMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.
Referenced by SiStripActionExecutor::createStatus(), fillStatus(), and fillStatusAtLumi().
{
if (!bookedStripStatus_) {
dqm_store->cd();
std::string strip_dir = "";
SiStripUtility::getTopFolderPath(dqm_store, "SiStrip", strip_dir);
if (strip_dir.size() == 0) strip_dir = "SiStrip";
// Non Standard Plots and should be put outside EventInfo folder
dqm_store->setCurrentFolder(strip_dir+"/MechanicalView");
std::string hname, htitle;
hname = "detFractionReportMap";
htitle = "SiStrip Report for Good Detector Fraction";
DetFractionReportMap = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
DetFractionReportMap->setAxisTitle("Sub Detector Type", 1);
DetFractionReportMap->setAxisTitle("Layer/Disc Number", 2);
hname = "sToNReportMap";
htitle = "SiStrip Report for Signal-to-Noise";
SToNReportMap = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
SToNReportMap->setAxisTitle("Sub Detector Type", 1);
SToNReportMap->setAxisTitle("Layer/Disc Number", 2);
// this is the main reportSummary 2D plot and should be in EventInfo
dqm_store->setCurrentFolder(strip_dir+"/EventInfo");
hname = "reportSummaryMap";
htitle = "SiStrip Report Summary Map";
SummaryReportMap = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
SummaryReportMap->setAxisTitle("Sub Detector Type", 1);
SummaryReportMap->setAxisTitle("Layer/Disc Number", 2);
SummaryReportGlobal = dqm_store->bookFloat("reportSummary");
int ibin = 0;
dqm_store->setCurrentFolder(strip_dir+"/EventInfo/reportSummaryContents");
for (std::map<std::string, std::string>::const_iterator it = SubDetFolderMap.begin();
it != SubDetFolderMap.end(); it++) {
ibin++;
std::string det = it->first;
DetFractionReportMap->setBinLabel(ibin,det);
SToNReportMap->setBinLabel(ibin,det);
SummaryReportMap->setBinLabel(ibin,det);
SubDetMEs local_mes;
if (det == "TECF") local_mes.detectorTag = "TEC+";
else if (det == "TECB") local_mes.detectorTag = "TEC-";
else if (det == "TIDF") local_mes.detectorTag = "TID+";
else if (det == "TIDB") local_mes.detectorTag = "TID-";
else local_mes.detectorTag = det;
std::string me_name;
me_name = "SiStrip_" + det;
local_mes.SummaryFlag = dqm_store->bookFloat(me_name);
me_name = "SiStrip_DetFraction_" + det;
local_mes.DetFraction = dqm_store->bookFloat(me_name);
me_name = "SiStrip_SToNFlag_" + det;
local_mes.SToNFlag = dqm_store->bookFloat(me_name);
SubDetMEsMap.insert(std::pair<std::string, SubDetMEs>(det, local_mes));
}
bookedStripStatus_ = true;
}
if (!bookedTrackingStatus_) {
dqm_store->cd();
std::string tracking_dir = "";
SiStripUtility::getTopFolderPath(dqm_store, "Tracking", tracking_dir);
if (tracking_dir.size() == 0) tracking_dir = "Tracking";
dqm_store->setCurrentFolder(tracking_dir+"/EventInfo");
TrackSummaryReportGlobal = dqm_store->bookFloat("reportSummary");
std::string hname, htitle;
hname = "reportSummaryMap";
htitle = "Tracking Report Summary Map";
TrackSummaryReportMap = dqm_store->book2D(hname, htitle, 3,0.5,3.5,1,0.5,1.5);
TrackSummaryReportMap->setAxisTitle("Track Quality Type", 1);
TrackSummaryReportMap->setAxisTitle("QTest Flag", 2);
TrackSummaryReportMap->setBinLabel(1,"Rate");
TrackSummaryReportMap->setBinLabel(2,"Chi2");
TrackSummaryReportMap->setBinLabel(3,"RecHits");
dqm_store->setCurrentFolder(tracking_dir+"/EventInfo/reportSummaryContents");
int ibin = 0;
for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
it != TrackingMEsMap.end(); it++) {
ibin++;
std::string name = it->first;
it->second.TrackingFlag = dqm_store->bookFloat("Track"+name);
TrackSummaryReportMap->setBinLabel(ibin,name);
}
bookedTrackingStatus_ = true;
dqm_store->cd();
}
}
| void SiStripQualityChecker::fillDetectorStatus | ( | DQMStore * | dqm_store, |
| const edm::ESHandle< SiStripDetCabling > & | cabling | ||
| ) | [private] |
Definition at line 303 of file SiStripQualityChecker.cc.
References DQMStore::cd(), DQMStore::dirExists(), MonitorElement::Fill(), fillSubDetStatus(), SiStripUtility::goToDir(), initialiseBadModuleList(), DQMStore::pwd(), SubDetFolderMap, SubDetMEsMap, and SummaryReportGlobal.
Referenced by fillStatus().
{
unsigned int xbin = 0;
float global_flag = 0;
dqm_store->cd();
std::string mdir = "MechanicalView";
if (!SiStripUtility::goToDir(dqm_store, mdir)) return;
std::string mechanicalview_dir = dqm_store->pwd();
initialiseBadModuleList();
for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin();
it != SubDetMEsMap.end(); it++) {
std::string det = it->first;
std::map<std::string, std::string>::const_iterator cPos = SubDetFolderMap.find(det);
if (cPos == SubDetFolderMap.end()) continue;
std::string dname = mechanicalview_dir + "/" + cPos->second;
if (!dqm_store->dirExists(dname)) continue;
dqm_store->cd(dname);
SubDetMEs local_mes = it->second;
xbin++;
float flag;
fillSubDetStatus(dqm_store, cabling, local_mes, xbin,flag);
global_flag += flag;
}
global_flag = global_flag/xbin*1.0;
if (SummaryReportGlobal) SummaryReportGlobal->Fill(global_flag);
dqm_store->cd();
}
| void SiStripQualityChecker::fillDetectorStatusAtLumi | ( | DQMStore * | dqm_store | ) | [private] |
Definition at line 605 of file SiStripQualityChecker.cc.
References DQMStore::cd(), MonitorElement::DQM_KIND_TH1F, MonitorElement::Fill(), DQMStore::get(), MonitorElement::getTH1F(), SiStripUtility::goToDir(), MonitorElement::kind(), label, DQMStore::pwd(), SubDetMEsMap, and SummaryReportGlobal.
Referenced by fillStatusAtLumi().
{
dqm_store->cd();
std::string rdir = "ReadoutView";
if (!SiStripUtility::goToDir(dqm_store, rdir)) return;
std::string fullpath = dqm_store->pwd()
+ "/FedSummary/PerLumiSection/"
+ "lumiErrorFraction";
MonitorElement* me = dqm_store->get(fullpath);
if (me && me->kind() == MonitorElement::DQM_KIND_TH1F) {
TH1F* th1 = me->getTH1F();
float global_fraction = 0.0;
float dets = 0.0;
for (int ibin = 1; ibin <= th1->GetNbinsX(); ibin++) {
std::string label = th1->GetXaxis()->GetBinLabel(ibin);
std::map<std::string, SubDetMEs>::iterator iPos = SubDetMEsMap.find(label);
if (iPos != SubDetMEsMap.end()) {
float fraction = 1.0 - th1->GetBinContent(ibin);
global_fraction += fraction;
dets++;
iPos->second.DetFraction->Fill(fraction);
iPos->second.SToNFlag->Fill(-1.0);
iPos->second.SummaryFlag->Fill(fraction);
}
}
global_fraction = global_fraction/dets;
if (SummaryReportGlobal) SummaryReportGlobal->Fill(global_fraction);
}
dqm_store->cd();
}
| void SiStripQualityChecker::fillDummyStatus | ( | ) |
Definition at line 213 of file SiStripQualityChecker.cc.
References bookedStripStatus_, bookedTrackingStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, MonitorElement::Fill(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), resetStatus(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.
Referenced by SiStripActionExecutor::fillDummyStatus(), fillStatus(), and fillStatusAtLumi().
{
resetStatus();
if (bookedStripStatus_) {
for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin();
it != SubDetMEsMap.end(); it++) {
SubDetMEs local_mes = it->second;
local_mes.DetFraction->Fill(-1.0);
local_mes.SToNFlag->Fill(-1.0);
local_mes.SummaryFlag->Fill(-1.0);
}
for (int xbin = 1; xbin < SummaryReportMap->getNbinsX()+1; xbin++) {
for (int ybin = 1; ybin < SummaryReportMap->getNbinsY()+1; ybin++) {
DetFractionReportMap->Fill(xbin, ybin, -1.0);
SToNReportMap->Fill(xbin, ybin, -1.0);
SummaryReportMap->Fill(xbin, ybin, -1.0);
}
}
SummaryReportGlobal->Fill(-1.0);
}
if (bookedTrackingStatus_) {
TrackSummaryReportGlobal->Fill(-1.0);
for (int xbin = 1; xbin < TrackSummaryReportMap->getNbinsX()+1; xbin++) {
for (int ybin = 1; ybin < TrackSummaryReportMap->getNbinsY()+1; ybin++) {
TrackSummaryReportMap->Fill(xbin, ybin, -1.0);
}
}
for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
it != TrackingMEsMap.end(); it++) {
it->second.TrackingFlag->Fill(-1.0);
}
/*
for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsLSMap.begin();
it != TrackingMEsLSMap.end(); it++) {
it->second.TrackingFlag->Fill(-1.0);
}
*/
}
}
| void SiStripQualityChecker::fillFaultyModuleStatus | ( | DQMStore * | dqm_store | ) |
Definition at line 556 of file SiStripQualityChecker.cc.
References badModuleList, DQMStore::bookInt(), DQMStore::cd(), DQMStore::dirExists(), MonitorElement::Fill(), DQMStore::get(), SiStripFolderOrganizer::getSubDetFolder(), SiStripUtility::goToDir(), DQMStore::pwd(), MonitorElement::Reset(), and DQMStore::setCurrentFolder().
Referenced by fillStatus().
{
if (badModuleList.size() == 0) return;
dqm_store->cd();
std::string mdir = "MechanicalView";
if (!SiStripUtility::goToDir(dqm_store, mdir)) return;
std::string mechanical_dir = dqm_store->pwd();
SiStripFolderOrganizer folder_organizer;
for (std::map<uint32_t,uint16_t>::const_iterator it = badModuleList.begin() ; it != badModuleList.end(); it++) {
uint32_t detId = it->first;
std::string subdet_folder ;
folder_organizer.getSubDetFolder(detId,subdet_folder);
if (!dqm_store->dirExists(subdet_folder)) {
subdet_folder = mechanical_dir + subdet_folder.substr(subdet_folder.find("MechanicalView")+14);
if (!dqm_store->dirExists(subdet_folder)) continue;
}
std::string bad_module_folder = subdet_folder + "/" + "BadModuleList";
dqm_store->setCurrentFolder(bad_module_folder);
std::ostringstream detid_str;
detid_str << detId;
std::string full_path = bad_module_folder + "/" + detid_str.str();
MonitorElement* me = dqm_store->get(full_path);
if (me) me->Reset();
else me = dqm_store->bookInt(detid_str.str());
me->Fill(it->second);
}
dqm_store->cd();
}
| void SiStripQualityChecker::fillStatus | ( | DQMStore * | dqm_store, |
| const edm::ESHandle< SiStripDetCabling > & | cabling | ||
| ) |
Definition at line 290 of file SiStripQualityChecker.cc.
References bookedStripStatus_, bookedTrackingStatus_, bookStatus(), fillDetectorStatus(), fillDummyStatus(), fillFaultyModuleStatus(), fillTrackingStatus(), edm::ParameterSet::getUntrackedParameter(), and pSet_.
Referenced by SiStripActionExecutor::fillStatus().
{
if (!bookedStripStatus_ || !bookedTrackingStatus_) bookStatus(dqm_store);
fillDummyStatus();
fillDetectorStatus(dqm_store, cabling);
fillTrackingStatus(dqm_store);
int faulty_moduleflag = pSet_.getUntrackedParameter<bool>("PrintFaultyModuleList", false);
if (faulty_moduleflag) fillFaultyModuleStatus(dqm_store);
}
| void SiStripQualityChecker::fillStatusAtLumi | ( | DQMStore * | dqm_store | ) |
Definition at line 596 of file SiStripQualityChecker.cc.
References bookedStripStatus_, bookedTrackingStatus_, bookStatus(), fillDetectorStatusAtLumi(), fillDummyStatus(), and fillTrackingStatusAtLumi().
Referenced by SiStripActionExecutor::fillStatusAtLumi().
{
if (!bookedStripStatus_ || !bookedTrackingStatus_) bookStatus(dqm_store);
fillDummyStatus();
fillDetectorStatusAtLumi(dqm_store);
fillTrackingStatusAtLumi(dqm_store);
}
| void SiStripQualityChecker::fillStatusHistogram | ( | MonitorElement * | me, |
| int | xbin, | ||
| int | ybin, | ||
| float | val | ||
| ) | [private] |
Definition at line 547 of file SiStripQualityChecker.cc.
References MonitorElement::DQM_KIND_TH2F, MonitorElement::getTH2F(), and MonitorElement::kind().
Referenced by fillSubDetStatus(), fillTrackingStatus(), and fillTrackingStatusAtLumi().
{
if (me && me->kind() == MonitorElement::DQM_KIND_TH2F) {
TH2F* th2d = me->getTH2F();
th2d->SetBinContent(xbin, ybin, val);
}
}
| void SiStripQualityChecker::fillSubDetStatus | ( | DQMStore * | dqm_store, |
| const edm::ESHandle< SiStripDetCabling > & | cabling, | ||
| SubDetMEs & | mes, | ||
| unsigned int | xbin, | ||
| float & | gflag | ||
| ) | [private] |
Definition at line 383 of file SiStripQualityChecker.cc.
References DQMStore::cd(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, dqm::qstatus::ERROR, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), SiStripUtility::getMEStatus(), getModuleStatus(), MonitorElement::getName(), MonitorElement::getQReports(), DQMStore::getSubdirs(), edm::ParameterSet::getUntrackedParameter(), min, mergeVDriftHistosByStation::name, pSet_, dqm::qstatus::STATUS_OK, SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, and SummaryReportMap.
Referenced by fillDetectorStatus().
{
int status_flag = pSet_.getUntrackedParameter<int>("GlobalStatusFilling", 1);
if (status_flag < 1) return;
std::vector<std::string> subDirVec = dqm_store->getSubdirs();
unsigned int ybin = 0;
int tot_ndet = 0;
int tot_errdet = 0;
float tot_ston_stat = 0;
for (std::vector<std::string>::const_iterator ic = subDirVec.begin();
ic != subDirVec.end(); ic++) {
std::string dname = (*ic);
if (dname.find("BadModuleList") != std::string::npos) continue;
std::vector<MonitorElement*> meVec;
ybin++;
dqm_store->cd((*ic));
meVec = dqm_store->getContents((*ic));
uint16_t ndet = 100;
int errdet = 0;
int ston_stat = 1;
int lnum = atoi(dname.substr(dname.find_last_of("_")+1).c_str());
ndet = cabling->connectedNumber(mes.detectorTag, lnum);
getModuleStatus(meVec, errdet);
for (std::vector<MonitorElement*>::const_iterator it = meVec.begin();
it != meVec.end(); it++) {
MonitorElement * me = (*it);
if (!me) continue;
std::vector<QReport *> reports = me->getQReports();
if (reports.size() == 0) continue;
std::string name = me->getName();
if( name.find("Summary_ClusterStoNCorr__OnTrack") != std::string::npos){
int istat = SiStripUtility::getMEStatus((*it));
if (reports[0]->getQTresult() == -1) {
ston_stat =-1;
} else {
if (istat == dqm::qstatus::ERROR) ston_stat = 0;
else if (istat == dqm::qstatus::STATUS_OK) ston_stat = 1;
}
}
}
if (ndet > 0) {
float eff_fac = 1 - (errdet*1.0/ndet);
fillStatusHistogram(SToNReportMap, xbin, ybin, ston_stat);
fillStatusHistogram(DetFractionReportMap, xbin, ybin, eff_fac);
if (ston_stat < 0) fillStatusHistogram(SummaryReportMap, xbin, ybin, eff_fac);
else fillStatusHistogram(SummaryReportMap, xbin, ybin, ston_stat*eff_fac);
tot_ndet += ndet;
tot_errdet += errdet;
tot_ston_stat += ston_stat;
}
dqm_store->cd((*ic));
}
if (tot_ndet > 0) {
float tot_eff_fac = 1 - (tot_errdet*1.0/tot_ndet);
if (mes.DetFraction) mes.DetFraction->Fill(tot_eff_fac);
float tot_ston_fac = tot_ston_stat/ybin;
if (mes.SToNFlag) mes.SToNFlag->Fill(tot_ston_fac);
if (tot_ston_fac < 0){
gflag = tot_eff_fac;
}else{
gflag = std::min(tot_eff_fac,tot_ston_fac);
}
if (mes.SummaryFlag) mes.SummaryFlag->Fill(gflag);
}
}
| void SiStripQualityChecker::fillTrackingStatus | ( | DQMStore * | dqm_store | ) | [private] |
Definition at line 333 of file SiStripQualityChecker.cc.
References DQMStore::cd(), dir, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), MonitorElement::getName(), MonitorElement::getQReports(), SiStripUtility::goToDir(), relval_steps::merge(), mergeVDriftHistosByStation::name, DQMStore::pwd(), ntuplemaker::status, TrackingMEsMap, TrackSummaryReportGlobal, TrackSummaryReportMap, and useGoodTracks_.
Referenced by fillStatus().
{
dqm_store->cd();
std::string dir = "Tracking";
if (!SiStripUtility::goToDir(dqm_store, dir)) return;
dir = "TrackParameters";
if (!SiStripUtility::goToDir(dqm_store, dir)) return;
std::vector<MonitorElement*> meVec1;
std::vector<MonitorElement*> meVec2;
if (useGoodTracks_){
meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties/GoodTracks");
meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties/GoodTracks");
}else{
meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties");
meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties");
}
std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size());
std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
float gstatus = 1.0;
for (std::vector<MonitorElement*>::const_iterator itME = meVec.begin(); itME != meVec.end(); itME++) {
MonitorElement * me = (*itME);
if (!me) continue;
std::vector<QReport *> qt_reports = me->getQReports();
if (qt_reports.size() == 0) continue;
std::string name = me->getName();
float status = 1.0;
int ibin = 0;
for (std::map<std::string, TrackingMEs>::const_iterator it = TrackingMEsMap.begin();
it != TrackingMEsMap.end(); it++) {
ibin++;
std::string hname = it->second.HistoName;
if (name.find(hname) != std::string::npos) {
status = qt_reports[0]->getQTresult();
it->second.TrackingFlag->Fill(status);
fillStatusHistogram(TrackSummaryReportMap, ibin, 1, status);
break;
}
}
gstatus = gstatus * status;
}
TrackSummaryReportGlobal->Fill(gstatus);
dqm_store->cd();
}
| void SiStripQualityChecker::fillTrackingStatusAtLumi | ( | DQMStore * | dqm_store | ) | [private] |
Definition at line 637 of file SiStripQualityChecker.cc.
References DQMStore::cd(), dir, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), MonitorElement::getMean(), MonitorElement::getName(), SiStripUtility::goToDir(), relval_steps::merge(), mergeVDriftHistosByStation::name, DQMStore::pwd(), ntuplemaker::status, TrackingMEsMap, TrackSummaryReportGlobal, TrackSummaryReportMap, and useGoodTracks_.
Referenced by fillStatusAtLumi().
{
dqm_store->cd();
std::string dir = "Tracking";
if (!SiStripUtility::goToDir(dqm_store, dir)) return;
dir = "TrackParameters";
if (!SiStripUtility::goToDir(dqm_store, dir)) return;
std::vector<MonitorElement*> meVec1;
std::vector<MonitorElement*> meVec2;
if (useGoodTracks_){
meVec1 = dqm_store->getContents(dqm_store->pwd()+"/LSanalysis");
// meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties/LSanalysis");
}else{
meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties");
meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties");
}
std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size());
std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
float gstatus = 1.0;
for (std::vector<MonitorElement*>::const_iterator itME = meVec.begin(); itME != meVec.end(); itME++) {
MonitorElement * me = (*itME);
if (!me) continue;
std::string name = me->getName();
float status = -1.0;
int ibin = 0;
for (std::map<std::string, TrackingMEs>::const_iterator it = TrackingMEsMap.begin();
it != TrackingMEsMap.end(); it++) {
ibin++;
std::string hname = it->second.HistoName+"lumiFlag_";
float lower_cut = it->second.LowerCut;
float upper_cut = it->second.UpperCut;
if (name.find(hname) != std::string::npos) {
if (me->getMean() <= lower_cut || me->getMean() > upper_cut) status = 0.0;
else status = 1.0;
it->second.TrackingFlag->Fill(status);
fillStatusHistogram(TrackSummaryReportMap, ibin, 1, status);
break;
} else {
}
}
if (status == -1.0) gstatus = -1.0;
else gstatus = gstatus * status;
}
TrackSummaryReportGlobal->Fill(gstatus);
dqm_store->cd();
}
| void SiStripQualityChecker::getModuleStatus | ( | std::vector< MonitorElement * > & | layer_mes, |
| int & | errdet | ||
| ) | [private] |
Definition at line 493 of file SiStripQualityChecker.cc.
References badModuleList, MonitorElement::DQM_KIND_TPROFILE, MonitorElement::DQM_KIND_TPROFILE2D, TkDetMap::getDetFromBin(), MonitorElement::getName(), MonitorElement::getQReports(), MonitorElement::kind(), mergeVDriftHistosByStation::name, SiStripUtility::setBadModuleFlag(), and tkDetMap_.
Referenced by fillSubDetStatus().
{
std::string lname;
std::map<uint32_t,uint16_t> bad_modules;
for (std::vector<MonitorElement*>::const_iterator it = layer_mes.begin();
it != layer_mes.end(); it++) {
MonitorElement * me = (*it);
if (!me) continue;
std::vector<QReport *> qreports = me->getQReports();
if (qreports.size() == 0) continue;
std::string name = me->getName();
std::vector<DQMChannel> bad_channels_me;
if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) {
bad_channels_me = qreports[0]->getBadChannels();
lname = "";
} else if (me->kind() == MonitorElement::DQM_KIND_TPROFILE2D && name.find("TkHMap") != std::string::npos) {
bad_channels_me = qreports[0]->getBadChannels();
lname = name.substr(name.find("TkHMap_")+7);
lname = lname.substr(lname.find("_T")+1);
}
for (std::vector<DQMChannel>::iterator it = bad_channels_me.begin(); it != bad_channels_me.end(); it++){
int xval = (*it).getBinX();
int yval = (*it).getBinY();
uint32_t detId = tkDetMap_->getDetFromBin(lname, xval, yval);
std::map<uint32_t,uint16_t>::iterator iPos = bad_modules.find(detId);
uint16_t flag;
if (iPos != bad_modules.end()){
flag = iPos->second;
SiStripUtility::setBadModuleFlag(name,flag);
iPos->second = flag;
} else {
flag = 0;
SiStripUtility::setBadModuleFlag(name,flag);
bad_modules.insert(std::pair<uint32_t,uint16_t>(detId,flag));
}
}
}
for(std::map<uint32_t,uint16_t>::const_iterator it = bad_modules.begin();
it != bad_modules.end(); it++) {
uint32_t detId = it->first;
uint16_t flag = it->second;
std::map<uint32_t,uint16_t>::iterator iPos = badModuleList.find(detId);
if (iPos != badModuleList.end()){
iPos->second = flag;
} else {
badModuleList.insert(std::pair<uint32_t,uint16_t>(detId,flag));
}
}
errdet = bad_modules.size();
}
| void SiStripQualityChecker::initialiseBadModuleList | ( | ) | [private] |
Definition at line 588 of file SiStripQualityChecker.cc.
References badModuleList.
Referenced by fillDetectorStatus().
{
for (std::map<uint32_t,uint16_t>::iterator it=badModuleList.begin(); it!=badModuleList.end(); it++) {
it->second = 0;
}
}
| void SiStripQualityChecker::printStatusReport | ( | ) |
Definition at line 462 of file SiStripQualityChecker.cc.
References SiStripQualityChecker::SubDetMEs::DetFraction, SiStripUtility::getMEValue(), SiStripQualityChecker::SubDetMEs::SToNFlag, SubDetMEsMap, and SiStripQualityChecker::SubDetMEs::SummaryFlag.
{
std::ostringstream det_summary_str;
for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin();
it != SubDetMEsMap.end(); it++) {
std::string det = it->first;
det_summary_str << std::setprecision(4);
det_summary_str << std::setiosflags(std::ios::fixed);
det_summary_str << " Printing Status for " << det << " : " << std::endl;
SubDetMEs local_mes = it->second;
std::string sval;
float fval1, fval2, fval3;
fval1 = fval2 = fval3 = -1.0;
SiStripUtility::getMEValue(local_mes.DetFraction, sval);
if (sval.size() > 0) fval1 = atof(sval.c_str());
SiStripUtility::getMEValue(local_mes.SToNFlag, sval);
if (sval.size() > 0) fval2 = atof(sval.c_str());
SiStripUtility::getMEValue(local_mes.SummaryFlag, sval);
if (sval.size() > 0) fval3 = atof(sval.c_str());
det_summary_str << std::setw(7) << " % of good detectors " << fval1
<< " SToN Flag " << fval2
<< " Summary Flag " << fval3 << std::endl;
}
}
| void SiStripQualityChecker::resetStatus | ( | ) |
Definition at line 256 of file SiStripQualityChecker.cc.
References bookedStripStatus_, bookedTrackingStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, MonitorElement::Reset(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.
Referenced by fillDummyStatus().
{
if (bookedStripStatus_) {
for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin();
it != SubDetMEsMap.end(); it++) {
SubDetMEs local_mes = it->second;
local_mes.DetFraction->Reset();
local_mes.SToNFlag->Reset();
local_mes.SummaryFlag->Reset();
}
DetFractionReportMap->Reset();
SToNReportMap->Reset();
SummaryReportMap->Reset();
SummaryReportGlobal->Reset();
}
if (bookedTrackingStatus_) {
TrackSummaryReportGlobal->Reset();
TrackSummaryReportMap->Reset();
for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
it != TrackingMEsMap.end(); it++) {
it->second.TrackingFlag->Reset();
}
/*
for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsLSMap.begin();
it != TrackingMEsLSMap.end(); it++) {
it->second.TrackingFlag->Reset();
}
*/
}
}
std::map<uint32_t,uint16_t> SiStripQualityChecker::badModuleList [private] |
Definition at line 78 of file SiStripQualityChecker.h.
Referenced by fillFaultyModuleStatus(), getModuleStatus(), initialiseBadModuleList(), and SiStripQualityChecker().
bool SiStripQualityChecker::bookedStripStatus_ [private] |
Definition at line 82 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillStatus(), fillStatusAtLumi(), resetStatus(), and SiStripQualityChecker().
bool SiStripQualityChecker::bookedTrackingStatus_ [private] |
Definition at line 83 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillStatus(), fillStatusAtLumi(), resetStatus(), and SiStripQualityChecker().
float SiStripQualityChecker::cutoffChi2overDoF_ [private] |
Definition at line 90 of file SiStripQualityChecker.h.
float SiStripQualityChecker::cutoffRecHits_ [private] |
Definition at line 91 of file SiStripQualityChecker.h.
float SiStripQualityChecker::cutoffTrackRate_ [private] |
Definition at line 89 of file SiStripQualityChecker.h.
Definition at line 68 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().
int SiStripQualityChecker::globalStatusFilling_ [private] |
Definition at line 84 of file SiStripQualityChecker.h.
Definition at line 80 of file SiStripQualityChecker.h.
Referenced by fillStatus(), fillSubDetStatus(), and SiStripQualityChecker().
Definition at line 69 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().
std::map<std::string, std::string> SiStripQualityChecker::SubDetFolderMap [private] |
Definition at line 65 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDetectorStatus(), and SiStripQualityChecker().
std::map<std::string, SubDetMEs> SiStripQualityChecker::SubDetMEsMap [private] |
Definition at line 64 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDetectorStatus(), fillDetectorStatusAtLumi(), fillDummyStatus(), printStatusReport(), and resetStatus().
Definition at line 72 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDetectorStatus(), fillDetectorStatusAtLumi(), fillDummyStatus(), and resetStatus().
Definition at line 70 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().
TkDetMap* SiStripQualityChecker::tkDetMap_ [private] |
Definition at line 87 of file SiStripQualityChecker.h.
Referenced by getModuleStatus(), and SiStripQualityChecker().
std::map<std::string, TrackingMEs> SiStripQualityChecker::TrackingMEsMap [private] |
Definition at line 66 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillTrackingStatus(), fillTrackingStatusAtLumi(), resetStatus(), and SiStripQualityChecker().
Definition at line 76 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillTrackingStatus(), fillTrackingStatusAtLumi(), and resetStatus().
Definition at line 74 of file SiStripQualityChecker.h.
Referenced by bookStatus(), fillDummyStatus(), fillTrackingStatus(), fillTrackingStatusAtLumi(), and resetStatus().
bool SiStripQualityChecker::useGoodTracks_ [private] |
Definition at line 85 of file SiStripQualityChecker.h.
Referenced by fillTrackingStatus(), fillTrackingStatusAtLumi(), and SiStripQualityChecker().
1.7.3