|
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
|
#ifndef SiStripDetSummary_h
#define SiStripDetSummary_h
#include "DataFormats/DetId/interface/DetId.h"
class TrackerTopology;
#include <sstream>
#include <map>
#include <cmath>
#include <iomanip>
#include <iostream>
/**
* @class SiStripDetSummary
* @author M. De Mattia
* @date 26/3/2009
* Class to compute and print summary information.
*
* If values are passed together with DetIds (method add( detId, value)), it computes the mean value and
* rms of a given quantity and is able to print a summary divided by layer/disk for each subdetector. <br>
* If instead only DetIds are passed (method add( detId )), it prints the count divided by layer/disk for
* each subdetector. <br>
* <br>
* Note: consider the possibility to move this class inside SiStripBaseObject as a protected member class.
*
*/
class SiStripDetSummary {
public:
explicit SiStripDetSummary(const TrackerTopology* tTopo) : computeMean_(true), trackerTopo_(tTopo) {
// Initialize valueMap_ with zeros
// WARNING: this initialization is strongly connected with how the map is filled in the add method
// TIB: layers = 4, stereo = the first 2
// TOB: layers = 6, stereo = the first 2
// TEC: wheels = 9, stereo = 9
// TID: wheels = 3, stereo = 3
unsigned int layers[] = {4, 6, 9, 3};
unsigned int stereo[] = {2, 2, 9, 3};
Values initValues;
for (unsigned int subDet = 0; subDet < 4; ++subDet) {
// Layers start from 1
for (unsigned int layer = 1; layer <= layers[subDet]; ++layer) {
valueMap_[1000 * (subDet + 1) + layer * 10] = initValues;
if (layer <= stereo[subDet])
valueMap_[1000 * (subDet + 1) + layer * 10 + 1] = initValues;
}
}
}
/// Used to compute the mean value of the value variable divided by subdetector, layer and mono/stereo
void add(DetId detid, float value);
/// Used to compute the number of entries divided by subdetector, layer and mono/stereo
inline void add(DetId detid) {
computeMean_ = false;
add(detid, 0);
}
/**
* Method used to write the output. By default mean == true and it writes the mean value. If mean == false
* it will write the count.
*/
void print(std::stringstream& ss, const bool mean = true) const;
inline void clear() { valueMap_.clear(); }
struct Values {
Values() : mean(0.), rms(0.), count(0) {}
double mean;
double rms;
unsigned int count;
};
std::map<unsigned int, Values> getCounts() { return valueMap_; }
protected:
// Maps to store the value and the counts
std::map<unsigned int, Values> valueMap_;
bool computeMean_;
private:
const TrackerTopology* trackerTopo_;
};
#endif
|