|
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
|
#include "DQM/SiStripCommissioningSources/interface/FedTimingTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace sistrip;
// -----------------------------------------------------------------------------
//
FedTimingTask::FedTimingTask(DQMStore* dqm, const FedChannelConnection& conn)
: CommissioningTask(dqm, conn, "FedTimingTask"),
timing_(),
nBins_(
40) //@@ this should be from number of scope mode samples (mean booking in event loop and putting scope mode length in trigger fed)
{
LogDebug("Commissioning") << "[FedTimingTask::FedTimingTask] Constructing object...";
}
// -----------------------------------------------------------------------------
//
FedTimingTask::~FedTimingTask() { LogDebug("Commissioning") << "[FedTimingTask::FedTimingTask] Destructing object..."; }
// -----------------------------------------------------------------------------
//
void FedTimingTask::book() {
LogDebug("Commissioning") << "[FedTimingTask::book]";
uint16_t nbins = 24 * nBins_; // 24 "fine" pll skews possible
std::string title;
title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::FED_TIMING,
sistrip::FED_KEY,
fedKey(),
sistrip::LLD_CHAN,
connection().lldChannel())
.title();
timing_.histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));
timing_.vNumOfEntries_.resize(nbins, 0);
timing_.vSumOfContents_.resize(nbins, 0);
timing_.vSumOfSquares_.resize(nbins, 0);
}
// -----------------------------------------------------------------------------
//
/*
Some notes:
- use all samples
- extract number of samples from trigger fed
- need to book histos in event loop?
- why only use fine skew setting when filling histos? should use coarse setting as well?
- why do different settings every 100 events - change more freq?
*/
void FedTimingTask::fill(const SiStripEventSummary& summary, const edm::DetSet<SiStripRawDigi>& digis) {
LogDebug("Commissioning") << "[FedTimingTask::fill]";
//@@ if scope mode length is in trigger fed, then
//@@ can add check here on number of digis
if (digis.data.size() < nBins_) {
edm::LogWarning("Commissioning") << "[FedTimingTask::fill]"
<< " Unexpected number of digis! " << digis.data.size();
} else {
uint32_t pll_fine = summary.pllFine();
for (uint16_t coarse = 0; coarse < nBins_ /*digis.data.size()*/; coarse++) {
uint16_t fine = (coarse + 1) * 24 - (pll_fine + 1);
updateHistoSet(timing_, fine, digis.data[coarse].adc());
}
}
}
// -----------------------------------------------------------------------------
//
void FedTimingTask::update() {
LogDebug("Commissioning") << "[FedTimingTask::update]";
updateHistoSet(timing_);
}
|