|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
/*
* See header file for a description of this class.
*
* \author G. Mila - INFN Torino
*
* threadsafe version (//-) oct/nov 2014 - WATWanAbdullah -ncpp-um-my
*
*/
#include "DQM/DTMonitorClient/src/DTSummaryClients.h"
// Framework
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <string>
using namespace edm;
using namespace std;
DTSummaryClients::DTSummaryClients(const ParameterSet& ps) : nevents(0) {
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients") << "[DTSummaryClients]: Constructor";
bookingdone = false;
}
DTSummaryClients::~DTSummaryClients() {
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients") << "DTSummaryClients: analyzed " << nevents << " events";
}
void DTSummaryClients::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients") << "[DTSummaryClients]: endJob";
}
void DTSummaryClients::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker,
DQMStore::IGetter& igetter,
edm::LuminosityBlock const& lumiSeg,
edm::EventSetup const& context) {
if (!bookingdone) {
// book the summary histos
ibooker.setCurrentFolder("DT/EventInfo");
summaryReport = ibooker.bookFloat("reportSummary");
// Initialize to 1 so that no alarms are thrown at the beginning of the run
summaryReport->Fill(1.);
summaryReportMap = ibooker.book2D("reportSummaryMap", "DT Report Summary Map", 12, 1, 13, 5, -2, 3);
summaryReportMap->setAxisTitle("sector", 1);
summaryReportMap->setAxisTitle("wheel", 2);
ibooker.setCurrentFolder("DT/EventInfo/reportSummaryContents");
for (int wheel = -2; wheel != 3; ++wheel) {
stringstream streams;
streams << "DT_Wheel" << wheel;
string meName = streams.str();
theSummaryContents.push_back(ibooker.bookFloat(meName));
// Initialize to 1 so that no alarms are thrown at the beginning of the run
theSummaryContents[wheel + 2]->Fill(1.);
}
}
bookingdone = true;
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients")
<< "[DTSummaryClients]: End of LS transition, performing the DQM client operation" << endl;
// reset the monitor elements
summaryReportMap->Reset();
summaryReport->Reset();
for (int ii = 0; ii != 5; ++ii) {
theSummaryContents[ii]->Reset();
}
bool noDTData = false;
// Check if DT data in each ROS have been read out and set the SummaryContents and the ErrorSummary
// accordignly
MonitorElement* dataIntegritySummary = igetter.get("DT/00-DataIntegrity/DataIntegritySummary");
if (dataIntegritySummary != nullptr) {
int nDisabledFED = 0;
for (int wheel = 1; wheel != 6; ++wheel) { // loop over the wheels
int nDisablesROS = 0;
for (int sect = 1; sect != 13; ++sect) { // loop over sectors
if (dataIntegritySummary->getBinContent(sect, wheel) == 0) {
nDisablesROS++;
}
}
if (nDisablesROS == 12) {
nDisabledFED++;
theSummaryContents[wheel - 1]->Fill(0);
}
}
if (nDisabledFED == 5) {
noDTData = true;
summaryReport->Fill(-1);
}
} else {
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients")
<< "Data Integrity Summary not found with name: DT/00-DataIntegrity/DataIntegritySummary" << endl;
}
double totalStatus = 0;
// protection
bool occupancyFound = true;
// Fill the map using, at the moment, only the information from DT occupancy
// problems at a granularity smaller than the chamber are ignored
for (int wheel = -2; wheel <= 2; wheel++) { // loop over wheels
// retrieve the occupancy summary
stringstream str;
str << "DT/01-Digi/OccupancySummary_W" << wheel;
MonitorElement* wheelOccupancySummary = igetter.get(str.str());
if (wheelOccupancySummary != nullptr) {
int nFailingChambers = 0;
for (int sector = 1; sector <= 12; sector++) { // loop over sectors
for (int station = 1; station != 5; ++station) { // loop over stations
double chamberStatus = wheelOccupancySummary->getBinContent(sector, station);
LogTrace("DTDQM|DTMonitorClient|DTSummaryClients")
<< "Wheel: " << wheel << " Stat: " << station << " Sect: " << sector << " status: " << chamberStatus
<< endl;
if (chamberStatus != 4) {
summaryReportMap->Fill(sector, wheel, 0.25);
} else {
nFailingChambers++;
}
LogTrace("DTDQM|DTMonitorClient|DTSummaryClients")
<< " sector (" << sector
<< ") status on the map is: " << summaryReportMap->getBinContent(sector, wheel + 3) << endl;
}
}
theSummaryContents[wheel + 2]->Fill((48. - nFailingChambers) / 48.);
totalStatus += (48. - nFailingChambers) / 48.;
} else {
occupancyFound = false;
LogVerbatim("DTDQM|DTMonitorClient|DTSummaryClients")
<< " Wheel Occupancy Summary not found with name: " << str.str() << endl;
}
}
if (occupancyFound && !noDTData)
summaryReport->Fill(totalStatus / 5.);
}
|