|
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
|
#include "DQM/SiStripCommissioningSources/interface/FineDelayTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/SiStripCommon/interface/SiStripDetKey.h"
#define NBINS (500)
#define LOWBIN (-125)
#define HIGHBIN (125)
// -----------------------------------------------------------------------------
//
CommissioningTask::HistoSet FineDelayTask::timing_;
CommissioningTask::MonitorElement* FineDelayTask::mode_ = nullptr;
// -----------------------------------------------------------------------------
//
FineDelayTask::FineDelayTask(DQMStore* dqm, const FedChannelConnection& conn)
: CommissioningTask(dqm, conn, "FineDelayTask") {
LogDebug("Commissioning") << "[FineDelayTask::FineDelayTask] Constructing object...";
}
// -----------------------------------------------------------------------------
//
FineDelayTask::~FineDelayTask() { LogDebug("Commissioning") << "[FineDelayTask::FineDelayTask] Destructing object..."; }
// -----------------------------------------------------------------------------
//
void FineDelayTask::book() {
LogDebug("Commissioning") << "[FineDelayTask::book]";
std::string title;
int nBins = NBINS;
SiStripDetKey detkeytracker((uint32_t)0);
// see if the global timing histogram is already booked
if (timing_.histo()) {
// if already booked, use it
LogDebug("Commissioning") << "[FineDelayTask::book] using existing histogram.";
} else {
// make a new histo on the tracker level if not existing yet
LogDebug("Commissioning") << "[LatencyTask::book] booking a new histogram.";
// construct the histo title
title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::FINE_DELAY,
sistrip::DET_KEY,
detkeytracker.key(),
sistrip::TRACKER,
0,
sistrip::extrainfo::clusterCharge_)
.title();
dqm()->setCurrentFolder(detkeytracker.path());
timing_.histo(dqm()->bookProfile(title,
title, // name and title
nBins,
LOWBIN,
HIGHBIN, // binning + range
100,
0.,
-1.,
"s")); // Y range : automatic
timing_.vNumOfEntries_.resize(nBins, 0);
timing_.vSumOfContents_.resize(nBins, 0);
timing_.vSumOfSquares_.resize(nBins, 0);
}
LogDebug("Commissioning") << "Binning is " << timing_.vNumOfEntries_.size();
LogDebug("Commissioning") << "[FineDelayTask::book] done";
if (!mode_) {
std::string pwd = dqm()->pwd();
std::string rootDir = pwd.substr(0, pwd.find(std::string(sistrip::root_) + "/") + (sizeof(sistrip::root_)));
dqm()->setCurrentFolder(rootDir);
mode_ = dqm()->bookInt("latencyCode");
}
}
// -----------------------------------------------------------------------------
//
void FineDelayTask::fill(const SiStripEventSummary& summary, const edm::DetSet<SiStripRawDigi>& digis) {
LogDebug("Commissioning") << "[FineDelayTask::fill]";
// retrieve the delay from the EventSummary
float delay = summary.ttcrx();
uint32_t latencyCode = (summary.layerScanned() >> 24) & 0xff;
LogDebug("Commissioning") << "[FineDelayTask::fill]: layerScanned() is " << summary.layerScanned();
int latencyShift =
latencyCode & 0x3f; // number of bunch crossings between current value and start of scan... must be positive
if (latencyShift > 32)
latencyShift -= 64; // allow negative values: we cover [-32,32].. should not be needed.
if ((latencyCode >> 6) == 2)
latencyShift -= 3; // layer in deconv, rest in peak
else if ((latencyCode >> 6) == 1)
latencyShift += 3; // layer in peak, rest in deconv
float correctedDelayBase =
delay - (latencyShift * 25.); // shifts the delay so that 0 corresponds to the current settings.
LogDebug("Commissioning") << "[FineDelayTask::fill]; the delay is " << delay;
// loop on the strips to find the (maybe) non-zero digi
for (unsigned int strip = 0; strip < digis.data.size(); strip++) {
if (digis.data[strip].adc() != 0) {
if ((digis.data[strip].adc() >> 8) == 255)
continue; // skip hit if TOF is in overflow
// apply the TOF correction
float tof = (digis.data[strip].adc() >> 8) / 10.;
float correctedDelay = correctedDelayBase - tof;
// compute the bin
float nbins = NBINS;
float lowbin = LOWBIN;
float highbin = HIGHBIN;
int bin = int((correctedDelay - lowbin) / ((highbin - lowbin) / nbins));
LogDebug("Commissioning") << "[FineDelayTask::fill]; using a hit with value " << (digis.data[strip].adc() & 0xff)
<< " at corrected delay of " << correctedDelay << " in bin " << bin << " (tof is "
<< tof << "( since adc = " << digis.data[strip].adc() << "))";
updateHistoSet(timing_, bin, digis.data[strip].adc() & 0xff);
if (mode_)
mode_->Fill(latencyCode);
}
}
}
// -----------------------------------------------------------------------------
//
void FineDelayTask::update() {
LogDebug("Commissioning") << "[FineDelayTask::update]";
updateHistoSet(timing_);
}
|