|
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
|
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctJetCounts.h"
using std::endl;
using std::ostream;
using std::vector;
/// static maximum number of jet counts
/// This can be up to 12 but we use some of the
/// available bandwidth for other information.
const unsigned L1GctJetCounts::MAX_TOTAL_COUNTS = 12;
/// MAX_TRUE_COUNTS specifies the bandwidth remaining
/// for real jet count information.
const unsigned L1GctJetCounts::MAX_TRUE_COUNTS = 6;
// default constructor
L1GctJetCounts::L1GctJetCounts() : m_data0(0), m_data1(0), m_bx(0) {}
// constructor for unpacking
L1GctJetCounts::L1GctJetCounts(uint32_t data0, uint32_t data1)
: m_data0(data0 & 0x7fff7fff), // Mask off bits 15 and 31 for better compression and consistency
m_data1(data1 & 0x7fff7fff), // with emulator constructor - these bits are not jet count data!
m_bx(0) {}
// constructor for unpacking
L1GctJetCounts::L1GctJetCounts(uint32_t data0, uint32_t data1, int16_t bx)
: m_data0(data0 & 0x7fff7fff), // Mask off bits 15 and 31 for better compression and consistency
m_data1(data1 & 0x7fff7fff), // with emulator constructor - these bits are not jet count data!
m_bx(bx) {}
// constructor for emulator
L1GctJetCounts::L1GctJetCounts(const std::vector<unsigned>& counts) : m_data0(0), m_data1(0), m_bx(0) {
// Assumes all required output data has been packed
// into 12 5-bit fields
if (counts.size() != MAX_TOTAL_COUNTS) {
} else {
for (unsigned int i = 0; i < 3; ++i) {
m_data0 += (counts[i] << (5 * i));
m_data0 += (counts[i + 3] << (5 * i + 16));
m_data1 += (counts[i + 6] << (5 * i));
m_data1 += (counts[i + 9] << (5 * i + 16));
}
}
}
// constructor for emulator
L1GctJetCounts::L1GctJetCounts(const std::vector<unsigned>& counts, int16_t bx) : m_data0(0), m_data1(0), m_bx(bx) {
if (counts.size() != MAX_TOTAL_COUNTS) {
} else {
for (unsigned int i = 0; i < 3; ++i) {
m_data0 += (counts[i] << (5 * i));
m_data0 += (counts[i + 3] << (5 * i + 16));
m_data1 += (counts[i + 6] << (5 * i));
m_data1 += (counts[i + 9] << (5 * i + 16));
}
}
}
// destructor
L1GctJetCounts::~L1GctJetCounts() {}
// return counts by index
unsigned L1GctJetCounts::count(unsigned i) const {
if (i < 6) {
return ((m_data0 >> (i < 3 ? (5 * i) : ((5 * i) + 1))) & 0x1f);
} else if (i < MAX_TOTAL_COUNTS) {
return ((m_data1 >> (i < 9 ? ((5 * i) - 30) : ((5 * i) - 29))) & 0x1f);
} else {
return 0;
}
}
// pretty print
ostream& operator<<(ostream& s, const L1GctJetCounts& c) {
s << "L1GctJetCounts : ";
for (unsigned int i = 0; i < 12; ++i) {
s << "\n count " << i << "=" << c.count(i);
}
return s;
}
|