|
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
|
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "CondFormats/SiStripObjects/interface/SiStripDetSummary.h"
void SiStripDetSummary::add(DetId detid, float value) {
int layer = 0;
int stereo = 0;
int detNum = 0;
// Using the operator[] if the element does not exist it is created with the default value. That is 0 for integral types.
switch (detid.subdetId()) {
case StripSubdetector::TIB:
layer = trackerTopo_->tibLayer(detid);
stereo = trackerTopo_->tibStereo(detid);
detNum = 1000;
break;
case StripSubdetector::TOB:
layer = trackerTopo_->tobLayer(detid);
stereo = trackerTopo_->tobStereo(detid);
detNum = 2000;
break;
case StripSubdetector::TEC:
// is this module in TEC+ or TEC-?
layer = trackerTopo_->tecWheel(detid);
stereo = trackerTopo_->tecStereo(detid);
detNum = 3000;
break;
case StripSubdetector::TID:
// is this module in TID+ or TID-?
layer = trackerTopo_->tidWheel(detid);
stereo = trackerTopo_->tidStereo(detid);
detNum = 4000;
break;
}
detNum += layer * 10 + stereo;
valueMap_[detNum].mean += value;
valueMap_[detNum].rms += value * value;
valueMap_[detNum].count += 1;
}
void SiStripDetSummary::print(std::stringstream& ss, const bool mean) const {
// Compute the mean for each detector and for each layer.
std::map<unsigned int, Values>::const_iterator valueMapIt = valueMap_.begin();
ss << "subDet" << std::setw(15) << "layer" << std::setw(16) << "mono/stereo" << std::setw(20);
if (mean)
ss << "mean +- rms" << std::endl;
else
ss << "count" << std::endl;
std::string detector;
std::string oldDetector;
for (; valueMapIt != valueMap_.end(); ++valueMapIt) {
int count = valueMapIt->second.count;
double mean = 0.;
double rms = 0.;
if (computeMean_ && count != 0) {
mean = (valueMapIt->second.mean) / count;
rms = (valueMapIt->second.rms) / count - mean * mean;
if (rms <= 0)
rms = 0;
else
rms = sqrt(rms);
}
// Detector type
switch ((valueMapIt->first) / 1000) {
case 1:
detector = "TIB ";
break;
case 2:
detector = "TOB ";
break;
case 3:
detector = "TEC ";
break;
case 4:
detector = "TID ";
break;
}
if (detector != oldDetector) {
ss << std::endl << detector;
oldDetector = detector;
} else
ss << " ";
// Layer number
int layer = (valueMapIt->first) / 10 - (valueMapIt->first) / 1000 * 100;
int stereo = valueMapIt->first - layer * 10 - (valueMapIt->first) / 1000 * 1000;
ss << std::setw(15) << layer << std::setw(13) << stereo << std::setw(18);
if (computeMean_)
ss << mean << " +- " << rms << std::endl;
else
ss << count << std::endl;
}
}
|