157 if(
verbose_)
LogInfo (
"HLTMuonVal") <<
">>> EndRun (HLTOverallSummary) <<<" << std::endl;
160 LogInfo (
"HLTMuonVal") <<
"No dqmstore... skipping processing step" << endl;
164 std::vector<string> histoNameVector;
169 float defaultValueIfNotFound = 1.0;
178 muonQualityBit =
dbe_->
get(
"HLT_Muon");
180 if (!muonQualityBit) {
181 LogInfo (
"HLTMuonVal") <<
"Can't find muonQuality bit... making a bit, setting it to zero" << endl;
184 muonQualityBit->
Fill(defaultValueIfNotFound);
189 eleQualityBit =
dbe_->
get(
"HLT_Electron");
191 if (!eleQualityBit) {
192 LogInfo (
"HLTMuonVal") <<
"Can't find eleQuality bit... making a bit, setting it to zero" << endl;
195 eleQualityBit->
Fill(defaultValueIfNotFound);
199 photonQualityBit =
dbe_->
get(
"HLT_Photon");
201 if (!photonQualityBit) {
202 LogInfo (
"HLTMuonVal") <<
"Can't find photonQuality bit... making a bit, setting it to zero" << endl;
205 photonQualityBit->
Fill(defaultValueIfNotFound);
218 MonitorElement* CertificationSummaryMap =
dbe_->
book2D(
"certificationSummaryMap",
"HLT: CertificationSummaryMap",3,-0.5,2.5,1,-0.5,0.5);
220 TH2 * reportSummaryMapTH2 = reportSummaryMap->
getTH2F();
222 reportSummaryMapTH2->GetXaxis()->SetBinLabel(1,
"Muon");
223 reportSummaryMapTH2->GetXaxis()->SetBinLabel(2,
"Electron");
224 reportSummaryMapTH2->GetXaxis()->SetBinLabel(3,
"Photon");
226 reportSummaryMapTH2->GetYaxis()->SetBinLabel(1,
"Quality");
229 TH2 * CertificationSummaryMapTH2 = CertificationSummaryMap->getTH2F();
231 CertificationSummaryMapTH2->GetXaxis()->SetBinLabel(1,
"Muon");
232 CertificationSummaryMapTH2->GetXaxis()->SetBinLabel(2,
"Electron");
233 CertificationSummaryMapTH2->GetXaxis()->SetBinLabel(3,
"Photon");
234 CertificationSummaryMapTH2->GetYaxis()->SetBinLabel(1,
"Quality");
240 float photonValue = photonQualityBit->getFloatValue();
242 float electronValue = eleQualityBit->getFloatValue();
244 float muonValue = muonQualityBit->getFloatValue();
246 float hltOverallValue = 1.0;
248 if ( (photonValue > 0.99)
249 && (electronValue > 0.99)
250 && (muonValue > 0.99) ) {
252 hltOverallValue = 1.0;
256 hltOverallValue = 0.0;
260 hltQualityBit->Fill(hltOverallValue);
262 unsigned int hltSummaryValue = 0x0;
264 unsigned int ELECTRON_MASK = 0x1;
265 unsigned int PHOTON_MASK = 0x2;
266 unsigned int MUON_MASK = 0x4;
268 if (electronValue > 0.99) hltSummaryValue = hltSummaryValue | ELECTRON_MASK;
269 if (photonValue > 0.99) hltSummaryValue = hltSummaryValue | PHOTON_MASK;
270 if (muonValue > 0.99) hltSummaryValue = hltSummaryValue | MUON_MASK;
272 hltQualitySummaryWord->Fill(hltSummaryValue);
274 reportSummaryMapTH2->SetBinContent(reportSummaryMapTH2->GetBin(1,1), muonValue);
275 reportSummaryMapTH2->SetBinContent(reportSummaryMapTH2->GetBin(2,1), electronValue);
276 reportSummaryMapTH2->SetBinContent(reportSummaryMapTH2->GetBin(3,1), photonValue);
278 CertificationSummaryMapTH2->SetBinContent(CertificationSummaryMapTH2->GetBin(1,1), muonValue);
279 CertificationSummaryMapTH2->SetBinContent(CertificationSummaryMapTH2->GetBin(2,1), electronValue);
280 CertificationSummaryMapTH2->SetBinContent(CertificationSummaryMapTH2->GetBin(3,1), photonValue);
MonitorElement * bookFloat(const char *name)
Book float.
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. "my/long/dir/my_histo")
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 setCurrentFolder(const std::string &fullpath)