|
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
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
|
#include "DQM/SiStripCommissioningSources/interface/CalibrationTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <CondFormats/SiStripObjects/interface/SiStripPedestals.h>
#include <arpa/inet.h>
#include <cstdio>
#include <fstream>
#include <netdb.h>
#include <sys/socket.h>
#include <sys/unistd.h>
// -----------------------------------------------------------------------------
//
CalibrationTask::CalibrationTask(DQMStore* dqm,
const FedChannelConnection& conn,
const sistrip::RunType& rtype,
const char* filename,
uint32_t run,
const SiStripPedestals& pedestals)
: CommissioningTask(dqm, conn, "CalibrationTask"),
runType_(rtype),
nBins_(0),
lastCalChan_(1000),
lastCalSel_(1000),
lastLatency_(1000),
extrainfo_(),
run_(run) {
// each latency point is 25ns, each calSel is 25/8 --> 64 points = 8 calsel + 8 latency --> as set in the DAQ configuration
if (runType_ == sistrip::CALIBRATION)
nBins_ = 64;
// each latency point is 25ns, each calSel is 25/8 --> 80 points = 8 calsel + 10 latency --> as set in the DAQ configuration
else if (runType_ == sistrip::CALIBRATION_DECO)
nBins_ = 80;
LogDebug("Commissioning") << "[CalibrationTask::CalibrationTask] Constructing object...";
// load the pedestals
SiStripPedestals::Range detPedRange = pedestals.getRange(conn.detId());
int start = conn.apvPairNumber() * 256;
int stop = start + 256;
int value = 0;
ped.reserve(256);
for (int strip = start; strip < stop; ++strip) {
value = int(pedestals.getPed(strip, detPedRange));
if (value > 895)
value -= 1024;
ped.push_back(value);
}
}
// -----------------------------------------------------------------------------
//
CalibrationTask::~CalibrationTask() {
LogDebug("Commissioning") << "[CalibrationTask::CalibrationTask] Destructing object...";
}
// -----------------------------------------------------------------------------
//
void CalibrationTask::book() {
LogDebug("Commissioning") << "[CalibrationTask::book]";
// book 16 histograms, one for each strip in a calibration group --> APV1 --> 16 = dimension of one calChan
if (calib1_.find(extrainfo_) == calib1_.end()) {
dqm()->setCurrentFolder(directory_);
calib1_[extrainfo_].resize(16);
for (int i = 0; i < 16; ++i) {
std::string postfix = extrainfo_ + "_istrip_" + std::to_string(i);
std::string title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
runType_,
sistrip::FED_KEY,
fedKey(),
sistrip::APV,
connection().i2cAddr(0),
postfix)
.title();
if (runType_ == sistrip::CALIBRATION)
calib1_[extrainfo_][i].histo(dqm()->book1D(title, title, nBins_, 0, 200));
else if (runType_ == sistrip::CALIBRATION_DECO)
calib1_[extrainfo_][i].histo(dqm()->book1D(title, title, nBins_, 0, 250));
calib1_[extrainfo_][i].isProfile_ = false;
calib1_[extrainfo_][i].vNumOfEntries_.resize(nBins_, 0);
}
}
// book 16 histograms, one for each strip in a calibration group --> APV2 --> 16 = dimension of one calChan
if (calib2_.find(extrainfo_) == calib2_.end()) {
dqm()->setCurrentFolder(directory_);
calib2_[extrainfo_].resize(16);
for (int i = 0; i < 16; ++i) {
std::string postfix = extrainfo_ + "_istrip_" + std::to_string(i);
std::string title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
runType_,
sistrip::FED_KEY,
fedKey(),
sistrip::APV,
connection().i2cAddr(1),
postfix)
.title();
if (runType_ == sistrip::CALIBRATION)
calib2_[extrainfo_][i].histo(dqm()->book1D(title, title, nBins_, 0, 200));
else if (runType_ == sistrip::CALIBRATION_DECO)
calib2_[extrainfo_][i].histo(dqm()->book1D(title, title, nBins_, 0, 250));
calib2_[extrainfo_][i].isProfile_ = false;
calib2_[extrainfo_][i].vNumOfEntries_.resize(nBins_, 0);
}
}
}
// -----------------------------------------------------------------------------
//
void CalibrationTask::fill(const SiStripEventSummary& summary, const edm::DetSet<SiStripRawDigi>& digis) {
LogDebug("Commissioning") << "[CalibrationTask::fill]";
if (lastCalChan_ != summary.calChan()) { // change in the calChan value
lastCalChan_ = summary.calChan();
lastLatency_ = summary.latency();
lastCalSel_ = summary.calSel();
extrainfo_ = "calChan_" + std::to_string(lastCalChan_);
book(); // book histograms and load the right one
} else {
lastCalChan_ = summary.calChan();
lastLatency_ = summary.latency();
lastCalSel_ = summary.calSel();
extrainfo_ = "calChan_" + std::to_string(lastCalChan_);
}
// Check if CalChan changed. In that case, save, reset histo, change title, and continue
int isub = lastCalChan_ < 4 ? lastCalChan_ + 4 : lastCalChan_ - 4;
int bin = 0;
if (runType_ == sistrip::CALIBRATION)
bin = (100 - summary.latency()) * 8 + (7 - summary.calSel());
else if (runType_ == sistrip::CALIBRATION_DECO)
bin = (102 - summary.latency()) * 8 + (7 - summary.calSel());
// Fill the histograms data-ped -(data-ped)_isub, the second term corresponds to the common mode substraction, looking at a strip far away.
for (int k = 0; k < 16; ++k) {
updateHistoSet(calib1_[extrainfo_][k],
bin,
digis.data[lastCalChan_ + k * 8].adc() - ped[lastCalChan_ + k * 8] -
(digis.data[isub + k * 8].adc() - ped[isub + k * 8]));
updateHistoSet(calib2_[extrainfo_][k],
bin,
digis.data[128 + lastCalChan_ + k * 8].adc() - ped[128 + lastCalChan_ + k * 8] -
(digis.data[128 + isub + k * 8].adc() - ped[128 + isub + k * 8]));
}
update(); //TODO: temporary: find a better solution later
}
// -----------------------------------------------------------------------------
//
void CalibrationTask::update() {
LogDebug("Commissioning") << "[CalibrationTask::update]"; // huge output
for (const auto& element : calib1_) { // all pulse for different calChan
for (auto vecelement : element.second) // all strips in a calCan
updateHistoSet(vecelement);
}
for (const auto& element : calib2_) { // all pulse for different calChan
for (auto vecelement : element.second) // all strips in a calCan
updateHistoSet(vecelement);
}
}
void CalibrationTask::setCurrentFolder(const std::string& dir) { directory_ = dir; }
|