|
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
|
#include "DQM/SiStripCommissioningSources/interface/Averages.h"
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
// ----------------------------------------------------------------------------
//
Averages::Averages() : n_(0), s_(0.), x_(0.), xx_(0.), median_(), mode_() { ; }
// ----------------------------------------------------------------------------
//
void Averages::add(const uint32_t& x, const uint32_t& e) {
mode_[x]++;
add(static_cast<float>(x), static_cast<float>(e));
}
// ----------------------------------------------------------------------------
//
void Averages::add(const uint32_t& x) {
mode_[x]++;
add(static_cast<float>(x), -1.);
}
// ----------------------------------------------------------------------------
//
void Averages::add(const float& x, const float& e) {
n_++;
if (e > 0.) {
float wt = 1. / sqrt(e);
s_ += wt;
x_ += x * wt;
xx_ += x * x * wt;
} else {
s_++;
x_ += x;
xx_ += x * x;
}
median_.push_back(x);
}
// ----------------------------------------------------------------------------
//
void Averages::add(const float& x) { add(x, -1.); }
// ----------------------------------------------------------------------------
//
void Averages::calc(Params& params) {
params.num_ = n_;
if (s_ > 0.) {
float m = x_ / s_;
float t = xx_ / s_ - m * m;
if (t > 0.) {
t = sqrt(t);
} else {
t = 0.;
}
params.mean_ = m;
params.rms_ = t;
params.weight_ = s_;
}
if (!median_.empty()) {
sort(median_.begin(), median_.end());
uint16_t index = median_.size() % 2 ? median_.size() / 2 : median_.size() / 2 - 1;
params.median_ = median_[index];
params.max_ = median_.back();
params.min_ = median_.front();
}
if (!mode_.empty()) {
uint32_t max = 0;
std::map<uint32_t, uint32_t>::const_iterator imap = mode_.begin();
for (; imap != mode_.end(); imap++) {
if (imap->second > max) {
max = imap->second;
params.mode_ = imap->first;
}
}
}
}
|