|
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
|
#include <iostream>
#include <map>
#include <cmath>
#include <vector>
#include <TMath.h>
#include <TH1.h>
#include "CalibTracker/SiPixelQuality/interface/SiPixelModuleStatus.h"
// ----------------------------------------------------------------------
SiPixelModuleStatus::SiPixelModuleStatus(int detId, int nrocs) : fDetid_(detId), fNrocs_(nrocs) {
for (int i = 0; i < fNrocs_; ++i) {
SiPixelRocStatus a;
fRocs_.push_back(a);
}
};
// ----------------------------------------------------------------------
SiPixelModuleStatus::~SiPixelModuleStatus() {}
// ----------------------------------------------------------------------
void SiPixelModuleStatus::fillDIGI(int iroc) {
if (iroc < fNrocs_)
fRocs_[iroc].fillDIGI();
}
// ----------------------------------------------------------------------
void SiPixelModuleStatus::fillFEDerror25(PixelFEDChannel ch) {
int roc_first = int(ch.roc_first);
int roc_last = int(ch.roc_last);
for (int iroc = 0; iroc < fNrocs_; iroc++) {
if (iroc >= roc_first && iroc <= roc_last) {
fRocs_[iroc].fillFEDerror25();
}
}
}
// ----------------------------------------------------------------------
int SiPixelModuleStatus::detid() { return fDetid_; }
// ----------------------------------------------------------------------
int SiPixelModuleStatus::nrocs() { return fNrocs_; }
// ----------------------------------------------------------------------
void SiPixelModuleStatus::setDetId(int detid) { fDetid_ = detid; }
// ----------------------------------------------------------------------
void SiPixelModuleStatus::setNrocs(int nRoc) { fNrocs_ = nRoc; }
// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateDIGI(int iroc, unsigned int nhit) {
if (iroc < fNrocs_)
fRocs_[iroc].updateDIGI(nhit);
}
// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateFEDerror25(int iroc, bool fedError25) {
if (iroc < fNrocs_)
fRocs_[iroc].updateFEDerror25(fedError25);
}
// ----------------------------------------------------------------------
unsigned int SiPixelModuleStatus::digiOccROC(int iroc) { return (iroc < fNrocs_ ? fRocs_[iroc].digiOccROC() : 0); }
// ----------------------------------------------------------------------
bool SiPixelModuleStatus::fedError25(int iroc) { return (iroc < fNrocs_ ? fRocs_[iroc].isFEDerror25() : false); }
// ----------------------------------------------------------------------
unsigned int SiPixelModuleStatus::digiOccMOD() {
unsigned int count(0);
for (int iroc = 0; iroc < fNrocs_; ++iroc) {
count += digiOccROC(iroc);
}
return count;
}
// ----------------------------------------------------------------------
double SiPixelModuleStatus::perRocDigiOcc() {
unsigned int ave(0);
for (int iroc = 0; iroc < fNrocs_; ++iroc) {
unsigned int inc = digiOccROC(iroc);
ave += inc;
}
return (1.0 * ave) / fNrocs_;
}
double SiPixelModuleStatus::perRocDigiOccVar() {
double fModAverage = SiPixelModuleStatus::perRocDigiOcc();
double sig = 1.0;
for (int iroc = 0; iroc < fNrocs_; ++iroc) {
unsigned int inc = digiOccROC(iroc);
sig += (fModAverage - inc) * (fModAverage - inc);
}
double fModSigma = sig / (fNrocs_ - 1);
return TMath::Sqrt(fModSigma);
}
// ----------------------------------------------------------------------
// Return the address not the value of ROC status
SiPixelRocStatus* SiPixelModuleStatus::getRoc(int iroc) { return (iroc < fNrocs_ ? &fRocs_[iroc] : nullptr); }
// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateModuleDIGI(int iroc, unsigned int nhits) {
if (iroc < fNrocs_)
fRocs_[iroc].updateDIGI(nhits);
}
// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateModuleStatus(SiPixelModuleStatus newData) {
bool isSameModule = true;
if (fDetid_ != newData.detid() || fNrocs_ != newData.nrocs()) {
isSameModule = false;
}
if (isSameModule) {
for (int iroc = 0; iroc < fNrocs_; ++iroc) { // loop over rocs
//update occupancy
fRocs_[iroc].updateDIGI(newData.digiOccROC(iroc));
//update FEDerror25
fRocs_[iroc].updateFEDerror25(newData.fedError25(iroc));
} // loop over rocs
} // if same module
}
|