|
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 "CalibTracker/SiStripAPVAnalysis/interface/ApvAnalysisFactory.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/ApvFactoryService.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQM/SiStripCommon/interface/ExtractTObject.h"
#include "DQM/SiStripCommon/interface/UpdateTProfile.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQM/SiStripCommissioningSources/interface/PedsOnlyTask.h"
using namespace sistrip;
// -----------------------------------------------------------------------------
//
PedsOnlyTask::PedsOnlyTask(DQMStore *dqm, const FedChannelConnection &rCONN)
: CommissioningTask(dqm, rCONN, "PedsOnlyTask") {
//@@ NOT GUARANTEED TO BE THREAD SAFE!
pApvFactory_ = edm::Service<ApvFactoryService>().operator->()->getApvFactory();
LogTrace(mlDqmSource_) << "[PedsOnlyTask::" << __func__ << "]"
<< " Constructing object...";
}
// -----------------------------------------------------------------------------
//
PedsOnlyTask::~PedsOnlyTask() {
LogTrace(mlDqmSource_) << "[PedsOnlyTask::" << __func__ << "]"
<< " Destructing object...";
// Have to delete pApvFactory_ manually even though we didn't create it
// ourself. :(
if (pApvFactory_) {
delete pApvFactory_;
}
}
// -----------------------------------------------------------------------------
//
void PedsOnlyTask::book() {
LogTrace(mlDqmSource_) << "[PedsOnlyTask::" << __func__ << "]";
const uint16_t nBINS = 256;
{
// Pedestals
std::string title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::PEDS_ONLY,
sistrip::FED_KEY,
fedKey(),
sistrip::LLD_CHAN,
connection().lldChannel(),
sistrip::extrainfo::pedestals_)
.title();
HistoSet oHSet;
oHSet.isProfile_ = true;
oHSet.vNumOfEntries_.resize(nBINS, 0);
oHSet.vSumOfContents_.resize(nBINS, 0);
oHSet.vSumOfSquares_.resize(nBINS, 0);
oHSet.histo(dqm()->bookProfile(title, title, nBINS, -0.5, nBINS * 1. - 0.5, 1025, 0., 1025.));
peds_.push_back(oHSet);
}
{
// Raw Noise
std::string title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
sistrip::PEDS_ONLY,
sistrip::FED_KEY,
fedKey(),
sistrip::LLD_CHAN,
connection().lldChannel(),
sistrip::extrainfo::rawNoise_)
.title();
HistoSet oHSet;
oHSet.isProfile_ = true;
oHSet.vNumOfEntries_.resize(nBINS, 0);
oHSet.vSumOfContents_.resize(nBINS, 0);
oHSet.vSumOfSquares_.resize(nBINS, 0);
oHSet.histo(dqm()->bookProfile(title, title, nBINS, -0.5, nBINS * 1. - 0.5, 1025, 0., 1025.));
peds_.push_back(oHSet);
}
pApvFactory_->instantiateApvs(connection().detId(), connection().nApvs());
}
// -----------------------------------------------------------------------------
//
void PedsOnlyTask::fill(const SiStripEventSummary &rSummary, const edm::DetSet<SiStripRawDigi> &rDigis) {
pApvFactory_->updatePair(connection().detId(), connection().apvPairNumber(), rDigis);
}
// -----------------------------------------------------------------------------
//
void PedsOnlyTask::update() {
static UpdateTProfile updateTProfile;
TProfile *pedsProf = ExtractTObject<TProfile>().extract(peds_[0].histo());
TProfile *noiseProf = ExtractTObject<TProfile>().extract(peds_[1].histo());
for (uint16_t nLclApv = 2 * connection().apvPairNumber(), nMaxLclApv = nLclApv + 2; nMaxLclApv > nLclApv; ++nLclApv) {
ApvAnalysis::PedestalType lclPedestals;
ApvAnalysis::PedestalType lclNoises;
pApvFactory_->getPedestal(connection().detId(), nLclApv, lclPedestals);
pApvFactory_->getRawNoise(connection().detId(), nLclApv, lclNoises);
const uint16_t nSTART_BIN = 128 * (nLclApv % 2);
for (uint16_t nBin = 0, nAbsBin = nSTART_BIN + nBin + 1; 128 > nBin; ++nBin, ++nAbsBin) {
updateTProfile.setBinContent(pedsProf, nAbsBin, 5, lclPedestals[nBin], lclNoises[nBin]);
updateTProfile.setBinContent(noiseProf, nAbsBin, 5, lclNoises[nBin], 0);
} // End loop over BINs
} // End loop over Local Apvs
}
// -----------------------------------------------------------------------------
|