|
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
|
#include "DQM/L1TMonitorClient/interface/L1TRPCTFClient.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "TRandom.h"
#include <TF1.h>
#include <cstdio>
#include <sstream>
#include <cmath>
#include <TProfile.h>
#include <TProfile2D.h>
using namespace edm;
using namespace std;
L1TRPCTFClient::L1TRPCTFClient(const edm::ParameterSet &ps) {
parameters_ = ps;
initialize();
}
L1TRPCTFClient::~L1TRPCTFClient() { LogInfo("TriggerDQM") << "[TriggerDQM]: ending... "; }
//--------------------------------------------------------
void L1TRPCTFClient::initialize() {
counterLS_ = 0;
counterEvt_ = 0;
// base folder for the contents of this job
monitorName_ = parameters_.getUntrackedParameter<string>("monitorName", "");
// cout << "Monitor name = " << monitorName_ << endl;
prescaleLS_ = parameters_.getUntrackedParameter<int>("prescaleLS", -1);
// cout << "DQM lumi section prescale = " << prescaleLS_ << " lumi section(s)"<< endl;
prescaleEvt_ = parameters_.getUntrackedParameter<int>("prescaleEvt", -1);
// cout << "DQM event prescale = " << prescaleEvt_ << " events(s)"<< endl;
output_dir_ = parameters_.getUntrackedParameter<string>("output_dir", "");
// cout << "DQM output dir = " << output_dir_ << endl;
input_dir_ = parameters_.getUntrackedParameter<string>("input_dir", "");
// cout << "DQM input dir = " << input_dir_ << endl;
verbose_ = parameters_.getUntrackedParameter<bool>("verbose", false);
m_runInEventLoop = parameters_.getUntrackedParameter<bool>("runInEventLoop", false);
m_runInEndLumi = parameters_.getUntrackedParameter<bool>("runInEndLumi", false);
m_runInEndRun = parameters_.getUntrackedParameter<bool>("runInEndRun", false);
m_runInEndJob = parameters_.getUntrackedParameter<bool>("runInEndJob", false);
LogInfo("TriggerDQM");
}
//--------------------------------------------------------
void L1TRPCTFClient::book(DQMStore::IBooker &ibooker) {
LogInfo("TriggerDQM") << "[TriggerDQM]: Begin Job";
ibooker.setCurrentFolder(output_dir_);
m_deadChannels = ibooker.book2D("RPCTF_deadchannels", "RPCTF deadchannels", 33, -16.5, 16.5, 144, -0.5, 143.5);
m_noisyChannels = ibooker.book2D("RPCTF_noisychannels", "RPCTF noisy channels", 33, -16.5, 16.5, 144, -0.5, 143.5);
}
//--------------------------------------------------------
void L1TRPCTFClient::dqmEndLuminosityBlock(DQMStore::IBooker &ibooker,
DQMStore::IGetter &igetter,
const edm::LuminosityBlock &lumiSeg,
const edm::EventSetup &c) {
if (verbose_)
std::cout << "L1TRPCTFClient::endLuminosityBlock" << std::endl;
if (m_runInEndLumi) {
book(ibooker);
processHistograms(igetter);
}
}
void L1TRPCTFClient::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
book(ibooker);
processHistograms(igetter);
}
//--------------------------------------------------------
void L1TRPCTFClient::processHistograms(DQMStore::IGetter &igetter) {
igetter.setCurrentFolder(input_dir_);
{
MonitorElement *me = igetter.get(input_dir_ + "/RPCTF_muons_eta_phi_bx0");
if (me) {
const QReport *qreport;
qreport = me->getQReport("DeadChannels_RPCTF_2D");
if (qreport) {
vector<dqm::me_util::Channel> badChannels = qreport->getBadChannels();
for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
++channel) {
m_deadChannels->setBinContent((*channel).getBinX(), (*channel).getBinY(), 100);
} // for(badchannels)
} //if (qreport)
qreport = me->getQReport("HotChannels_RPCTF_2D");
if (qreport) {
vector<dqm::me_util::Channel> badChannels = qreport->getBadChannels();
for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
++channel) {
// (*channel).getBinY() == 0 for NoisyChannels QTEST
m_noisyChannels->setBinContent((*channel).getBinX(), 100);
} // for(badchannels)
} //if (qreport)
// else std::cout << "dupa" << std::endl;
} // if (me)
}
if (verbose_) {
std::vector<string> meVec = igetter.getMEs();
for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
std::string full_path = input_dir_ + "/" + (*it);
MonitorElement *me = igetter.get(full_path);
// for this MEs, get list of associated QTs
std::vector<QReport *> Qtest_map = me->getQReports();
if (!Qtest_map.empty()) {
std::cout << "Test: " << full_path << std::endl;
for (std::vector<QReport *>::const_iterator it = Qtest_map.begin(); it != Qtest_map.end(); ++it) {
std::cout << " Name " << (*it)->getQRName() << " Status " << (*it)->getStatus() << std::endl;
std::vector<dqm::me_util::Channel> badChannels = (*it)->getBadChannels();
vector<dqm::me_util::Channel>::iterator badchsit = badChannels.begin();
while (badchsit != badChannels.end()) {
int ix = (*badchsit).getBinX();
int iy = (*badchsit).getBinY();
std::cout << "(" << ix << "," << iy << ") ";
++badchsit;
}
std::cout << std::endl;
}
}
} //
}
}
|