|
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
|
#include "DQM/SiStripCommissioningSources/interface/VpspScanTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripFecKey.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <algorithm>
using namespace sistrip;
// -----------------------------------------------------------------------------
//
VpspScanTask::VpspScanTask(DQMStore* dqm, const FedChannelConnection& conn)
: CommissioningTask(dqm, conn, "VpspScanTask"), vpsp_() {}
// -----------------------------------------------------------------------------
//
VpspScanTask::~VpspScanTask() {}
// -----------------------------------------------------------------------------
//
void VpspScanTask::book() {
uint16_t nbins = 60;
std::string title;
vpsp_.resize(2);
for (uint16_t iapv = 0; iapv < 2; iapv++) {
if (connection().i2cAddr(iapv)) {
std::stringstream extra_info;
extra_info << sistrip::apv_ << iapv;
title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::VPSP_SCAN,
sistrip::FED_KEY,
fedKey(),
sistrip::LLD_CHAN,
connection().lldChannel(),
extra_info.str())
.title();
vpsp_[iapv].histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));
vpsp_[iapv].vNumOfEntries_.resize(nbins, 0);
vpsp_[iapv].vSumOfContents_.resize(nbins, 0);
vpsp_[iapv].vSumOfSquares_.resize(nbins, 0);
}
}
}
// -----------------------------------------------------------------------------
//
void VpspScanTask::fill(const SiStripEventSummary& summary, const edm::DetSet<SiStripRawDigi>& digis) {
// Retrieve VPSP setting and CCU channel
uint32_t vpsp = summary.vpsp();
uint32_t ccu_chan = summary.vpspCcuChan();
// Check CCU channel from EventSummary is consistent with this module
if (SiStripFecKey(fecKey()).ccuChan() != ccu_chan) {
return;
}
if (digis.data.size() != 256) {
edm::LogWarning(mlDqmSource_) << "[VpspScanTask::" << __func__ << "]"
<< " Unexpected number of digis! " << digis.data.size();
return;
}
// Fill histo with baseline(calc'ed from median value of data)
for (uint16_t iapv = 0; iapv < 2; iapv++) {
if (vpsp >= vpsp_[iapv].vNumOfEntries_.size()) {
edm::LogWarning(mlDqmSource_) << "[VpspScanTask::" << __func__ << "]"
<< " Unexpected VPSP value! " << vpsp;
return;
}
std::vector<uint16_t> baseline;
baseline.reserve(128);
for (uint16_t idigi = 128 * iapv; idigi < 128 * (iapv + 1); idigi++) {
baseline.push_back(digis.data[idigi].adc());
}
sort(baseline.begin(), baseline.end());
uint16_t index = baseline.size() % 2 ? baseline.size() / 2 : baseline.size() / 2 - 1;
if (!baseline.empty()) {
updateHistoSet(vpsp_[iapv], vpsp, baseline[index]);
}
}
}
// -----------------------------------------------------------------------------
//
void VpspScanTask::update() {
for (uint32_t iapv = 0; iapv < vpsp_.size(); iapv++) {
updateHistoSet(vpsp_[iapv]);
}
}
|