|
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
151
152
|
#include "DQM/SiStripMonitorSummary/interface/SiStripCablingDQM.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "TCanvas.h"
SiStripCablingDQM::SiStripCablingDQM(edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> token,
edm::RunNumber_t iRun,
edm::ParameterSet const &hPSet,
edm::ParameterSet const &fPSet,
const TrackerTopology *tTopo,
const TkDetMap *tkDetMap)
: SiStripBaseCondObjDQMGet<SiStripDetCabling, SiStripDetCablingRcd>{token, iRun, hPSet, fPSet, tTopo} {
if (HistoMaps_On_) {
Tk_HM_ = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/Histo_Map", "Cabling_TkMap", 0.);
}
}
SiStripCablingDQM::~SiStripCablingDQM() {}
void SiStripCablingDQM::getActiveDetIds(const edm::EventSetup &eSetup) {
// Get active and total detIds
getConditionObject(eSetup);
if (!condObj_) {
edm::LogError("InvalidCablingHandle") << "Invalid Cabling Handle";
return;
}
condObj_->addActiveDetectorsRawIds(activeDetIds);
condObj_->addAllDetectorsRawIds(activeDetIds);
// Initialize arrays for counting:
int counterTIB[4];
for (int i = 0; i < 4; i++)
counterTIB[i] = 0;
int counterTID[2][3];
for (int i = 0; i < 2; i++) {
for (int j = 0; j < 3; j++)
counterTID[i][j] = 0;
}
int counterTOB[6];
for (int i = 0; i < 6; i++)
counterTOB[i] = 0;
int counterTEC[2][9];
for (int i = 0; i < 2; i++) {
for (int j = 0; j < 9; j++)
counterTEC[i][j] = 0;
}
// fill arrays for counting and fill Histo_Map with value for connected :
for (const auto detId : activeDetIds) {
StripSubdetector subdet(detId);
if (HistoMaps_On_) {
Tk_HM_->fill(detId, condObj_->nApvPairs(detId) * 2);
}
if (fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")) {
int32_t n_conn = 0;
for (uint32_t connDet_i = 0; connDet_i < condObj_->getConnections(detId).size(); connDet_i++) {
if (condObj_->getConnections(detId)[connDet_i] != nullptr &&
condObj_->getConnections(detId)[connDet_i]->isConnected() != 0)
n_conn++;
}
fillTkMap(detId, n_conn * 2.);
}
switch (subdet.subdetId()) {
case StripSubdetector::TIB: {
int i = tTopo_->tibLayer(detId) - 1;
counterTIB[i]++;
break;
}
case StripSubdetector::TID: {
int j = tTopo_->tidWheel(detId) - 1;
int side = tTopo_->tidSide(detId);
if (side == 2) {
counterTID[0][j]++;
} else if (side == 1) {
counterTID[1][j]++;
}
break;
}
case StripSubdetector::TOB: {
int i = tTopo_->tobLayer(detId) - 1;
counterTOB[i]++;
break;
}
case StripSubdetector::TEC: {
int j = tTopo_->tecWheel(detId) - 1;
int side = tTopo_->tecSide(detId);
if (side == 2) {
counterTEC[0][j]++;
} else if (side == 1) {
counterTEC[1][j]++;
}
break;
}
}
} // idet
// obtained from tracker.dat and hard-coded
int TIBDetIds[4] = {672, 864, 540, 648};
int TIDDetIds[2][3] = {{136, 136, 136}, {136, 136, 136}};
int TOBDetIds[6] = {1008, 1152, 648, 720, 792, 888};
int TECDetIds[2][9] = {{408, 408, 408, 360, 360, 360, 312, 312, 272}, {408, 408, 408, 360, 360, 360, 312, 312, 272}};
DQMStore *dqmStore_ = edm::Service<DQMStore>().operator->();
std::string FolderName = fPSet_.getParameter<std::string>("FolderName_For_QualityAndCabling_SummaryHistos");
dqmStore_->setCurrentFolder(FolderName);
// dqmStore_->cd("SiStrip/MechanicalView/");
MonitorElement *ME;
ME = dqmStore_->book2D("SummaryOfCabling", "SummaryOfCabling", 6, 0.5, 6.5, 9, 0.5, 9.5);
ME->setAxisTitle("Sub Det", 1);
ME->setAxisTitle("Layer", 2);
ME->setBinLabel(1, "TIB");
ME->setBinLabel(2, "TID F");
ME->setBinLabel(3, "TID B");
ME->setBinLabel(4, "TOB");
ME->setBinLabel(5, "TEC F");
ME->setBinLabel(6, "TEC B");
for (int i = 0; i < 4; i++) {
ME->Fill(1, i + 1, float(counterTIB[i]) / TIBDetIds[i]);
}
for (int i = 0; i < 2; i++) {
for (int j = 0; j < 3; j++) {
ME->Fill(i + 2, j + 1, float(counterTID[i][j]) / TIDDetIds[i][j]);
}
}
for (int i = 0; i < 6; i++) {
ME->Fill(4, i + 1, float(counterTOB[i]) / TOBDetIds[i]);
}
for (int i = 0; i < 2; i++) {
for (int j = 0; j < 9; j++) {
ME->Fill(i + 5, j + 1, float(counterTEC[i][j]) / TECDetIds[i][j]);
}
}
if (fPSet_.getParameter<bool>("OutputSummaryAtLayerLevelAsImage")) {
TCanvas c1("c1");
ME->getTH1()->Draw("TEXT");
ME->getTH1()->SetStats(kFALSE);
std::string name(ME->getTitle());
name += ".png";
c1.Print(name.c_str());
}
}
|