|
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
|
#include "DQM/SiStripCommissioningSources/interface/FedCablingTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <algorithm>
#include <sstream>
#include <iomanip>
using namespace sistrip;
// -----------------------------------------------------------------------------
//
FedCablingTask::FedCablingTask(DQMStore* dqm, const FedChannelConnection& conn)
: CommissioningTask(dqm, conn, "FedCablingTask"), histos_() {}
// -----------------------------------------------------------------------------
//
FedCablingTask::~FedCablingTask() {}
// -----------------------------------------------------------------------------
//
void FedCablingTask::book() {
histos_.resize(2);
std::string title;
uint16_t nbins = 0;
std::string extra_info = "";
for (uint16_t iter = 0; iter < 2; iter++) {
// Define number of histo bins and title
if (iter == 0) {
nbins = 1024;
extra_info = sistrip::feDriver_;
} else if (iter == 1) {
nbins = 96;
extra_info = sistrip::fedChannel_;
} else {
edm::LogWarning(mlDqmSource_) << "[FedCablingTask::" << __func__ << "]"
<< " Unexpected number of HistoSets: " << iter;
}
title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::FED_CABLING,
sistrip::FED_KEY,
fedKey(),
sistrip::LLD_CHAN,
connection().lldChannel(),
extra_info)
.title();
histos_[iter].histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));
histos_[iter].vNumOfEntries_.resize(nbins, 0);
histos_[iter].vSumOfContents_.resize(nbins, 0);
histos_[iter].vSumOfSquares_.resize(nbins, 0);
}
}
// -----------------------------------------------------------------------------
//
void FedCablingTask::fill(const SiStripEventSummary& summary,
const uint16_t& fed_id,
const std::map<uint16_t, float>& fed_ch) {
if (fed_ch.empty()) {
edm::LogWarning(mlDqmSource_) << "[FedCablingTask::" << __func__ << "]"
<< " No FED channels with high signal!";
return;
} else {
LogTrace(mlDqmSource_) << "[FedCablingTask::" << __func__ << "]"
<< " Found " << fed_ch.size() << " FED channels with high signal!";
}
std::map<uint16_t, float>::const_iterator ichan = fed_ch.begin();
for (; ichan != fed_ch.end(); ichan++) {
updateHistoSet(histos_[0], fed_id, ichan->second);
updateHistoSet(histos_[1], ichan->first, ichan->second);
}
}
// -----------------------------------------------------------------------------
//
void FedCablingTask::update() {
for (uint32_t iter = 0; iter < histos_.size(); iter++) {
updateHistoSet(histos_[iter]);
}
}
|