|
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
|
#ifndef DataFormats_L1TCalorimeterPhase2_GCTHadDigiCluster_h
#define DataFormats_L1TCalorimeterPhase2_GCTHadDigiCluster_h
#include <ap_int.h>
#include <vector>
#include "DataFormats/L1TCalorimeterPhase2/interface/CaloPFCluster.h"
namespace l1tp2 {
class GCTHadDigiCluster {
private:
// Data
unsigned long long int clusterData;
// Constants
static constexpr float LSB_PT = 0.5; // 0.5 GeV
// start of the unused bits
static constexpr int n_bits_unused_start = 31;
// reference to corresponding float cluster
edm::Ref<l1tp2::CaloPFClusterCollection> clusterRef_;
public:
GCTHadDigiCluster() { clusterData = 0; }
GCTHadDigiCluster(ap_uint<64> data) { clusterData = data; }
// Note types of the constructor
GCTHadDigiCluster(ap_uint<12> pt, int etaCr, int phiCr, ap_uint<4> hoe) {
// To use .range() we need an ap class member
ap_uint<64> temp_data;
ap_uint<7> etaCrDigitized = abs(etaCr);
ap_int<7> phiCrDigitized = phiCr;
temp_data.range(11, 0) = pt.range();
temp_data.range(18, 12) = etaCrDigitized.range();
temp_data.range(25, 19) = phiCrDigitized.range();
clusterData = temp_data;
}
// Setters
void setRef(const edm::Ref<l1tp2::CaloPFClusterCollection> &clusterRef) { clusterRef_ = clusterRef; }
// Getters
ap_uint<64> data() const { return clusterData; }
// Other getters
float ptLSB() const { return LSB_PT; }
ap_uint<12> pt() const { return data().range(11, 0); }
float ptFloat() const { return pt() * ptLSB(); }
// crystal eta (unsigned 7 bits)
int eta() const { return (ap_uint<7>)data().range(18, 12); }
// crystal phi (signed 7 bits)
int phi() const { return (ap_int<7>)data().range(25, 19); }
// HoE value
ap_uint<4> hoe() const { return data().range(30, 26); }
// Check that unused bits are zero
const int unusedBitsStart() const { return n_bits_unused_start; }
bool passNullBitsCheck(void) const { return ((data() >> unusedBitsStart()) == 0); }
// Get the underlying ref
edm::Ref<l1tp2::CaloPFClusterCollection> clusterRef() const { return clusterRef_; }
};
// Collection typedefs
// This represents the 36 GCTHadDigiClusters in one link (one link spans 4 RCT cards, each RCT card sends 9 clusters (zero-padded and sorted by decreasing pT)
// The ordering of the 4 RCT cards in this std::vector is, e.g. for GCT1.SLR3, real phi -50 to -20 degrees, then real phi -20 to 10 degrees, then real phi 10 to 40 degrees, and lastly real phi 40 to 70 degrees
typedef std::vector<l1tp2::GCTHadDigiCluster> GCTHadDigiClusterLink;
// This represents the 12 links sending GCTHadDigiClusters in the full barrel: there are 12 links = (3 GCT cards) * (two SLRs per GCT) * (one positive eta link and one negative eta link)
// The ordering of the links in this std::vector is (GCT1.SLR1 negEta, GCT.SLR1 posEta, GCT1.SLR3 negEta, GCT1.SLR3 posEta, then analogously for GCT2 and GCT3)
typedef std::vector<l1tp2::GCTHadDigiClusterLink> GCTHadDigiClusterCollection;
} // namespace l1tp2
#endif
|