|
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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
#ifndef L1TMUONPHASE2_PHASE2GMT_CONSTANTS
#define L1TMUONPHASE2_PHASE2GMT_CONSTANTS
#include "ap_int.h"
namespace Phase2L1GMT {
// INPUT Dataformat
// Track
const int BITSTTCURV = 15;
const int BITSTTCURV2 = 27;
const int BITSTTPHI = 12;
const int BITSTTTANL = 16;
// Not used in the emulator now, but will in future
const int BITSTTZ0 = 12;
const int BITSTTD0 = 13;
const int BITSTTCHI2 = 4;
const int BITSTTBENDCHI2 = 3;
const int BITSTTMASK = 7;
const int BITSTTTRACKMVA = 3;
const int BITSTTOTHERMVA = 6;
// Bitwidth for common internal variables in GMT
const int BITSPT = 13;
const int BITSPHI = 13;
const int BITSETA = 13;
const int BITSZ0 = 10;
const int BITSD0 = 12;
//Muon ROI
const int BITSSTUBCOORD = 10;
const int BITSSTUBETA = 8;
const int BITSSTUBID = 12;
const int BITSSTUBPHIQUALITY = 4;
const int BITSSTUBETAQUALITY = 4;
const int BITSSTUBTIME = 8;
const int BITSSTAMUONQUALITY = 6;
const int BITSMUONBX = 4;
//PreTrackMatherdMuon
const int BITSMATCHQUALITY = 9;
const int BITSMUONBETA = 4;
//Track Muon Match
const int BITSSIGMAETA = 4;
const int BITSSIGMACOORD = 4;
const int BITSPROPCOORD = 6; //lowered this to make room for another prop coefficient
const int BITSPROPSIGMACOORD_A = 6;
const int BITSPROPSIGMACOORD_B = 6;
const int BITSPROPSIGMAETA_A = 5;
const int BITSPROPSIGMAETA_B = 5;
// OUTPUT Dataformat
// Bitwidth for standalone muons to CL1 and GT
const int BITSSAZ0 = 5;
const int BITSSAD0 = 7;
const int BITSSAQUAL = 4;
// Bitwidth for dataformat to GT
const int BITSGTPT = 16;
const int BITSGTPHI = 13;
const int BITSGTETA = 14;
const int BITSGTZ0 = 10;
const int BITSGTD0 = 10;
const int BITSGTQUAL = 6;
const int BITSGTQUALSC = 7;
const int BITSGTISO = 6;
const int BITSGTBETA = 4;
// Bitwidth for Tau->3mu object
const int BITSTMPT = 8;
const int BITSTMPHI = 8;
const int BITSTMETA = 8;
const int BITSTMMASS2 = 8;
const int BITSTMTYPE = 6;
const int BITSTMIDX = 4;
const int BITSTMQUAL = 4;
const float maxCurv_ = 0.00855; // 2 GeV pT Rinv is in cm
const float maxPhi_ = 1.026; // relative to the center of the sector
const float maxTanl_ = 8.0;
const float maxZ0_ = 25.6;
const float maxD0_ = 15.36;
// Updated barrelLimit according to Karol, https://indico.cern.ch/event/1113802/#1-phase2-gmt-performance-and-i
const int barrelLimit0_ = 1.4 / 0.00076699039 / 8;
const int barrelLimit1_ = 1.1 / 0.00076699039 / 8;
const int barrelLimit2_ = 0.95 / 0.00076699039 / 8;
const int barrelLimit3_ = 0.95 / 0.00076699039 / 8;
const int barrelLimit4_ = 0;
// LSB
const float LSBpt = 0.03125;
const float LSBphi = 2. * M_PI / pow(2, BITSPHI);
const float LSBeta = 2. * M_PI / pow(2, BITSETA);
const float LSBGTiso = 0.25; // 0.25GeV, LSBphi * 8
const float LSBGTz0 = 0.05; // 0.5mm, in sync with GTT and Correlator
const float LSBGTd0 = 0.03; // from GT interface doc
const float LSBSAz0 = 1.6; // 0.05 * 32 cm, with range +- 25.6
const float LSBSAd0 = 3.84; // 0.03 * 128 cm, with range +- 245.76
typedef ap_uint<64> wordtype;
typedef ap_uint<1> valid_gt_t; //valid
typedef ap_uint<1> q_gt_t; //charge
typedef ap_uint<BITSGTPT> pt_gt_t; //pt of tracker muon
typedef ap_int<BITSGTPHI> phi_gt_t; //phi of tracker muon
typedef ap_int<BITSGTETA> eta_gt_t; //eta of tracker muon
typedef ap_int<BITSGTZ0> z0_gt_t; //z0 of tracker muon
typedef ap_int<BITSGTD0> d0_gt_t; //d0 of tracker muon
typedef ap_uint<BITSGTISO> iso_gt_t; //isolation of tracker muon
typedef ap_uint<BITSGTBETA> beta_gt_t; //beta of tracker muon
typedef ap_uint<BITSGTQUAL> qual_gt_t; //quality of tracker muon
typedef ap_int<BITSGTQUALSC> qualsc_gt_t; //quality score of tracker muon (qual/q)
//Standalone muon datatype
typedef ap_uint<1> valid_sa_t; //valid
typedef ap_uint<BITSGTPT> pt_sa_t; //pt of standalone muon
typedef ap_int<BITSGTPHI> phi_sa_t; //phi of standalone muon
typedef ap_int<BITSGTETA> eta_sa_t; //eta of standalone muon
typedef ap_int<BITSSAZ0> z0_sa_t; //z0 of standalone muon
typedef ap_int<BITSSAD0> d0_sa_t; //d0 of standalone muon
typedef ap_uint<1> q_sa_t; //charge of standalone muon
typedef ap_uint<BITSSAQUAL> qual_sa_t; //quality of standalone muon
inline uint64_t twos_complement(long long int v, uint bits) {
uint64_t mask = (1 << bits) - 1;
if (v >= 0)
return v & mask;
else
return (~(-v) + 1) & mask;
}
template <class T>
inline int wordconcat(T& word, int bstart, long int input, int bitsize) {
int bend = bstart + bitsize - 1;
word.range(bend, bstart) = twos_complement(input, bitsize);
return bend + 1;
}
} // namespace Phase2L1GMT
#endif
|