|
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
|
// test
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctHFBitCounts.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctHFRingEtSums.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctInternHFData.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctInternEtSum.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctInternJetData.h"
#include <iostream>
#include <cstdlib>
using namespace std;
int main() {
// test HF bit counts set methods
L1GctHFBitCounts b;
unsigned c[4];
for (unsigned i = 0; i < 4; ++i)
c[i] = 0;
for (unsigned i = 0; i < 4; ++i) {
for (unsigned j = 0; j < 0x7; ++j) {
b.setBitCount(i, j);
c[i] = j;
// check all counts are as expected
for (unsigned k = 0; k < 4; ++k) {
bool check = (b.bitCount(k) == c[k]);
if (!check) {
std::cout << "L1GctHFBitCounts failed : ";
std::cout << " bitCount(" << std::dec << k << ") = " << std::hex << b.bitCount(k);
std::cout << ", expected " << std::hex << c[k] << std::endl;
exit(1);
}
}
}
}
// test HF ring sums set methods
L1GctHFRingEtSums s;
for (unsigned i = 0; i < 4; ++i)
c[i] = 0;
for (unsigned i = 0; i < 4; ++i) {
for (unsigned j = 0; j < 0x7; ++j) {
s.setEtSum(i, j);
c[i] = j;
// check all counts are as expected
for (unsigned k = 0; k < 4; ++k) {
bool check = (s.etSum(k) == c[k]);
if (!check) {
std::cout << "L1GctHFBitCounts failed : ";
std::cout << " bitCount(" << std::dec << k << ") = " << std::hex << s.etSum(k);
std::cout << ", expected " << std::hex << c[k] << std::endl;
exit(1);
}
}
}
}
// test intern HF data
L1GctInternHFData d;
for (unsigned i = 0; i < 4; ++i)
c[i] = 0;
for (unsigned i = 0; i < 4; ++i) {
for (unsigned j = 0; j < 0xff; ++j) {
d.setValue(i, j);
c[i] = j;
// check all counts are as expected
for (unsigned k = 0; k < 4; ++k) {
bool check = (d.value(k) == c[k]);
if (!check) {
std::cout << "L1GctHFBitCounts failed : ";
std::cout << " bitCount(" << std::dec << k << ") = " << std::hex << d.value(k);
std::cout << ", expected " << std::hex << c[k] << std::endl;
exit(1);
}
}
}
}
// test intern Et Sum
L1GctInternJetData jd;
for (unsigned rank = 0; rank < 0x3f; ++rank) {
for (unsigned tauVeto = 0; tauVeto < 2; ++tauVeto) {
for (unsigned phi = 0; phi < 18; ++phi) {
for (unsigned eta = 0; eta < 11; ++eta) {
for (unsigned et = 0; et < 0xfff; ++et) {
for (unsigned oflow = 0; oflow < 2; ++oflow) {
for (unsigned sgnEta = 0; sgnEta < 2; ++sgnEta) {
jd.setData(sgnEta, oflow, et, eta, phi, tauVeto, rank);
if (jd.sgnEta() != sgnEta || jd.oflow() != oflow || jd.et() != et || jd.eta() != eta ||
jd.phi() != phi || jd.tauVeto() != tauVeto || jd.rank() != rank) {
std::cout << "L1GctInternEtSum failed : " << std::endl;
std::cout << "Expected sgnEta=" << sgnEta;
std::cout << " oflow=" << oflow;
std::cout << " et=" << et;
std::cout << " eta=" << eta;
std::cout << " phi=" << phi;
std::cout << " tauVeto=" << tauVeto;
std::cout << " rank=" << rank << std::endl;
std::cout << "Got " << jd << std::endl;
exit(1);
}
}
}
}
}
}
}
}
exit(0);
}
|