|
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
|
#include "DQM/RPCMonitorDigi/interface/RPCLinkSynchroHistoMaker.h"
#include "DataFormats/RPCDigi/interface/RPCRawSynchro.h"
#include <algorithm>
#include <utility>
struct OrderLbSpread {
bool operator()(const std::pair<double, unsigned int>& lb1, const std::pair<double, unsigned int>& lb2) {
return lb1.first > lb2.first;
}
};
struct OrderLbOccup {
bool operator()(const std::pair<unsigned int, unsigned int>& lb1, const std::pair<unsigned int, unsigned int>& lb2) {
return lb1.first > lb2.first;
}
};
void RPCLinkSynchroHistoMaker::fill(TH1F* hDelay, TH2F* hDelaySpread, TH2F* hTopOccup, TH2F* hTopSpread) const {
hDelay->Reset();
hDelaySpread->Reset();
hTopOccup->Reset();
hTopSpread->Reset();
typedef std::vector<std::pair<unsigned int, unsigned int> > TopOccup;
typedef std::vector<std::pair<double, unsigned int> > TopSpread;
TopOccup topOccup(10, std::make_pair(0, 0));
TopSpread topSpread(10, std::make_pair(0., 0));
for (unsigned int idx = 0; idx < theLinkStat.theLinkStatMap.size(); ++idx) {
const RPCLinkSynchroStat::BoardAndCounts& bc = theLinkStat.theLinkStatMap[idx];
int sum = bc.second.sum();
double rms = bc.second.rms();
hDelaySpread->Fill(bc.second.mean() - 3., bc.second.rms());
if (sum == 0)
continue;
for (int i = 0; i <= 7; ++i)
hDelay->Fill(i - 3, bc.second.counts()[i]);
std::pair<unsigned int, unsigned int> canOccup = std::make_pair(sum, idx);
std::pair<double, unsigned int> canSpread = std::make_pair(rms, idx);
TopOccup::iterator io = upper_bound(topOccup.begin(), topOccup.end(), canOccup, OrderLbOccup());
TopSpread::iterator is = upper_bound(topSpread.begin(), topSpread.end(), canSpread, OrderLbSpread());
if (io != topOccup.end()) {
topOccup.insert(io, canOccup);
topOccup.erase(topOccup.end() - 1);
}
if (is != topSpread.end()) {
topSpread.insert(is, canSpread);
topSpread.erase(topSpread.end() - 1);
}
}
for (int itop = 0; itop < 10; itop++) {
const RPCLinkSynchroStat::BoardAndCounts& occup = theLinkStat.theLinkStatMap[topOccup[itop].second];
const RPCLinkSynchroStat::BoardAndCounts& spread = theLinkStat.theLinkStatMap[topSpread[itop].second];
hTopOccup->GetYaxis()->SetBinLabel(itop + 1, occup.first.name().c_str());
hTopSpread->GetYaxis()->SetBinLabel(itop + 1, spread.first.name().c_str());
for (unsigned int icount = 0; icount < occup.second.counts().size(); icount++) {
hTopOccup->SetBinContent(icount + 1, itop + 1, float(occup.second.counts()[icount]));
hTopSpread->SetBinContent(icount + 1, itop + 1, float(spread.second.counts()[icount]));
}
}
}
|