|
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
|
#include "CalibTracker/SiStripAPVAnalysis/interface/ApvAnalysis.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/TkApvMask.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/TkNoiseCalculator.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/TkPedestalCalculator.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/TkCommonModeCalculator.h"
#include <algorithm>
using namespace std;
ApvAnalysis::ApvAnalysis(int nEvForUpdate) {
theTkCommonModeCalculator = nullptr;
theTkPedestalCalculator = nullptr;
theTkNoiseCalculator = nullptr;
theTkApvMask = nullptr;
nEventsForNoiseCalibration_ = 0;
eventsRequiredToUpdate_ = nEvForUpdate;
}
void ApvAnalysis::newEvent() const {
theTkPedestalCalculator->newEvent();
theTkNoiseCalculator->newEvent();
theTkCommonModeCalculator->newEvent();
}
void ApvAnalysis::updateCalibration(edm::DetSet<SiStripRawDigi>& in) {
theTkPedestalCalculator->updatePedestal(in);
PedestalType noise;
if (theTkPedestalCalculator->status()->isUpdating()) {
nEventsForNoiseCalibration_++;
if (theTkNoiseCalculator->noise().empty()) {
noise = theTkPedestalCalculator->rawNoise();
theTkNoiseCalculator->setStripNoise(noise);
theTkApvMask->calculateMask(noise);
}
PedestalType pedestal = theTkPedestalCalculator->pedestal();
PedestalType tmp;
tmp.clear();
edm::DetSet<SiStripRawDigi>::const_iterator it = in.data.begin();
int i = 0;
for (; it != in.data.end(); it++) {
tmp.push_back((*it).adc() - pedestal[i]);
i++;
}
PedestalType tmp2 = theTkCommonModeCalculator->doIt(tmp);
if (!tmp2.empty()) {
theTkNoiseCalculator->updateNoise(tmp2);
}
if (nEventsForNoiseCalibration_ % eventsRequiredToUpdate_ == 1 && nEventsForNoiseCalibration_ > 1) {
noise = theTkNoiseCalculator->noise();
theTkApvMask->calculateMask(noise);
}
}
}
|