ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
LVL1::gFEXaltMetAlgo Class Reference

#include <gFEXaltMetAlgo.h>

Inheritance diagram for LVL1::gFEXaltMetAlgo:
Collaboration diagram for LVL1::gFEXaltMetAlgo:

Public Member Functions

 gFEXaltMetAlgo (const std::string &type, const std::string &name, const IInterface *parent)
 Constructor.
virtual StatusCode initialize () override
 standard Athena-Algorithm method
virtual void setAlgoConstant (std::vector< int > &&A_thr, std::vector< int > &&B_thr, const int rhoPlusThr) override
virtual void altMetAlgo (const gTowersCentral &Atwr, const gTowersCentral &Btwr, std::array< uint32_t, 4 > &outTOB) const override

Private Member Functions

void metFPGA (const gTowersCentral &twrs, int &MET_x, int &MET_y, const unsigned short FPGA_NO) const
void metTotal (const int A_MET_x, const int A_MET_y, const int B_MET_x, const int B_MET_y, int &MET_x, int &MET_y, int &MET) const
int get_rho (const gTowersCentral &twrs) const
int get_sigma (const gTowersCentral &twrs) const
void rho_MET (const gTowersCentral &twrs, int &MET_x, int &MET_y, const int rho, const int sigma) const
int sumEtFPGAnc (const gTowersCentral &twrs, const unsigned short FPGA_NO) const
int sumEtFPGArms (const gTowersCentral &twrs, const int sigma) const
int sumEt (const int A_sumEt, const int B_sumEt) const
float sinLUT (const unsigned int phiIDX, const unsigned int aw) const
float cosLUT (const unsigned int phiIDX, const unsigned int aw) const

Private Attributes

std::array< std::vector< int >, 2 > m_etaThr
int m_rhoPlusThr {}

Detailed Description

Definition at line 24 of file gFEXaltMetAlgo.h.

Constructor & Destructor Documentation

◆ gFEXaltMetAlgo()

LVL1::gFEXaltMetAlgo::gFEXaltMetAlgo ( const std::string & type,
const std::string & name,
const IInterface * parent )

Constructor.

Definition at line 20 of file gFEXaltMetAlgo.cxx.

20 :
21base_class(type, name, parent)
22{}

Member Function Documentation

◆ altMetAlgo()

void LVL1::gFEXaltMetAlgo::altMetAlgo ( const gTowersCentral & Atwr,
const gTowersCentral & Btwr,
std::array< uint32_t, 4 > & outTOB ) const
overridevirtual

Definition at line 41 of file gFEXaltMetAlgo.cxx.

42 {
43
44 //FPGA A observables
45 int A_MET_x_nc = 0x0;
46 int A_MET_y_nc = 0x0;
47 int A_MET_x_rms = 0x0;
48 int A_MET_y_rms = 0x0;
49
50 int A_sumEt_nc = 0x0;
51 int A_sumEt_rms = 0x0;
52
53 //FPGA B observables
54 int B_MET_x_nc = 0x0;
55 int B_MET_y_nc = 0x0;
56 int B_MET_x_rms = 0x0;
57 int B_MET_y_rms = 0x0;
58
59 int B_sumEt_nc = 0x0;
60 int B_sumEt_rms = 0x0;
61
62 //Global observables
63 int MET_x_nc = 0x0;
64 int MET_y_nc = 0x0;
65 int MET_nc = 0x0;
66 int MET_x_rms = 0x0;
67 int MET_y_rms = 0x0;
68 int MET_rms = 0x0;
69
70
71 int total_sumEt_nc = 0x0;
72 int total_sumEt_rms = 0x0;
73
74 metFPGA(Atwr, A_MET_x_nc, A_MET_y_nc, 0);
75 metFPGA(Btwr, B_MET_x_nc, B_MET_y_nc, 1);
76
77 metTotal(A_MET_x_nc, A_MET_y_nc, B_MET_x_nc, B_MET_y_nc, MET_x_nc, MET_y_nc, MET_nc);
78
79 int A_rho{get_rho(Atwr)};
80 int B_rho{get_rho(Btwr)};
81 int A_sigma{3*get_sigma(Atwr)};
82 int B_sigma{3*get_sigma(Btwr)};
83
84 rho_MET(Atwr, A_MET_x_rms, A_MET_y_rms, A_rho, A_sigma);
85 rho_MET(Btwr, B_MET_x_rms, B_MET_y_rms, B_rho, B_sigma);
86
87 metTotal(A_MET_x_rms, A_MET_y_rms, B_MET_x_rms, B_MET_y_rms, MET_x_rms, MET_y_rms, MET_rms);
88
89 A_sumEt_nc = sumEtFPGAnc(Atwr, 0);
90 B_sumEt_nc = sumEtFPGAnc(Btwr, 1);
91 total_sumEt_nc = sumEt(A_sumEt_nc, B_sumEt_nc);
92 total_sumEt_nc = total_sumEt_nc/4;
93
94 A_sumEt_rms = sumEtFPGArms(Atwr, A_sigma);
95 B_sumEt_rms = sumEtFPGArms(Btwr, B_sigma);
96 total_sumEt_rms = sumEt(A_sumEt_rms, B_sumEt_rms);
97 total_sumEt_rms = total_sumEt_rms/4;
98 //Define a vector to be filled with all the TOBs of one event
99
100 //TOB order
101 // 1) MET_x | MET_y <- ncMET
102 // 2) MET_x | MET_y <- rms
103 // 3) MET | sumET <- ncMET
104 // 4) MET | sumET <- rms
105
106 // fill in TOBs
107 // The order of the TOBs is given according to the TOB ID (TODO: check how it's done in fw)
108
109 // First TOB is (MET, SumEt)
110 outTOB[0] = (MET_y_nc& 0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
111 outTOB[0] = outTOB[0] | (MET_x_nc & 0x00000FFF) << 12;//Quantity 1 (in bit number 12)
112 if (MET_y_nc != 0) outTOB[0] = outTOB[0] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
113 if (MET_x_nc != 0) outTOB[0] = outTOB[0] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
114 outTOB[0] = outTOB[0] | (2 & 0x0000001F) << 26;//TOB ID temporary set to 2 according to JwoJ convention (need updates in EDM)
115
116// Second TOB is (MET_x, MET_y)
117 outTOB[1] = (MET_y_rms& 0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
118 outTOB[1] = outTOB[1] | (MET_x_rms & 0x00000FFF) << 12;//Quantity 1 (in bit number 12)
119 if (MET_y_rms != 0) outTOB[1] = outTOB[1] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
120 if (MET_x_rms != 0) outTOB[1] = outTOB[1] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
121 outTOB[1] = outTOB[1] | (2 & 0x0000001F) << 26;//TOB ID temporary set to 2 according to JwoJ convention (need updates in EDM)
122
123// Third TOB is hard components (MHT_x, MHT_y)
124 outTOB[2] = (total_sumEt_nc& 0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
125 outTOB[2] = outTOB[2] | (MET_nc & 0x00000FFF) << 12;//Quantity 1 (in bit number 12)
126 if (total_sumEt_nc != 0) outTOB[2] = outTOB[2] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
127 if (MET_nc != 0) outTOB[2] = outTOB[2] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
128 outTOB[2] = outTOB[2] | (1 & 0x0000001F) << 26;//TOB ID temporary set to 1 according to JwoJ convention (need updates in EDM)
129
130 // Fourth TOB is hard components (MST_x, MST_y)
131 outTOB[3] = (total_sumEt_rms& 0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
132 outTOB[3] = outTOB[3] | (MET_rms & 0x00000FFF) << 12;//Quantity 1 (in bit number 12)
133 if (total_sumEt_rms != 0) outTOB[3] = outTOB[3] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
134 if (MET_rms != 0) outTOB[3] = outTOB[3] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
135 outTOB[3] = outTOB[3] | (1 & 0x0000001F) << 26;//TOB ID temporary set to 1 according to JwoJ convention (need updates in EDM)
136
137
138}
LVL1::gFEXaltMetAlgo::metFPGA
void metFPGA(const gTowersCentral &twrs, int &MET_x, int &MET_y, const unsigned short FPGA_NO) const
Definition gFEXaltMetAlgo.cxx:142
LVL1::gFEXaltMetAlgo::metTotal
void metTotal(const int A_MET_x, const int A_MET_y, const int B_MET_x, const int B_MET_y, int &MET_x, int &MET_y, int &MET) const
Definition gFEXaltMetAlgo.cxx:157
LVL1::gFEXaltMetAlgo::sumEtFPGArms
int sumEtFPGArms(const gTowersCentral &twrs, const int sigma) const
Definition gFEXaltMetAlgo.cxx:239
LVL1::gFEXaltMetAlgo::get_sigma
int get_sigma(const gTowersCentral &twrs) const
Definition gFEXaltMetAlgo.cxx:193
LVL1::gFEXaltMetAlgo::sumEtFPGAnc
int sumEtFPGAnc(const gTowersCentral &twrs, const unsigned short FPGA_NO) const
Definition gFEXaltMetAlgo.cxx:226
LVL1::gFEXaltMetAlgo::rho_MET
void rho_MET(const gTowersCentral &twrs, int &MET_x, int &MET_y, const int rho, const int sigma) const
Definition gFEXaltMetAlgo.cxx:211
LVL1::gFEXaltMetAlgo::sumEt
int sumEt(const int A_sumEt, const int B_sumEt) const
Definition gFEXaltMetAlgo.cxx:253
LVL1::gFEXaltMetAlgo::get_rho
int get_rho(const gTowersCentral &twrs) const
Definition gFEXaltMetAlgo.cxx:179

◆ cosLUT()

float LVL1::gFEXaltMetAlgo::cosLUT ( const unsigned int phiIDX,
const unsigned int aw ) const
private

Definition at line 274 of file gFEXaltMetAlgo.cxx.

275{
276 float c = ((float)phiIDX)/std::pow(2,aw);
277 float rad = (2*M_PI) *c;
278 float rcos = std::cos(rad);
279 return rcos;
280}
M_PI
#define M_PI
Definition ActiveFraction.h:11

◆ get_rho()

int LVL1::gFEXaltMetAlgo::get_rho ( const gTowersCentral & twrs) const
private

Definition at line 179 of file gFEXaltMetAlgo.cxx.

179 {
180 const int rows = twrs.size();
181 const int cols = twrs[0].size();
182 const int n{rows*cols};
183 float rho = 0;
184 for(int i = 0; i < rows; i++) {
185 for(int j = 0; j < cols; j++) {
186 rho += twrs[i][j] < m_rhoPlusThr ? twrs[i][j] : 0;
187 }
188 }
189 return rho/n;
190}

◆ get_sigma()

int LVL1::gFEXaltMetAlgo::get_sigma ( const gTowersCentral & twrs) const
private

Definition at line 193 of file gFEXaltMetAlgo.cxx.

193 {
194
195 int rows = twrs.size();
196 int cols = twrs[0].size();
197 const int n{rows*cols};
198 int sigma = 0;
199 for(int i = 0; i < rows; ++i) {
200 for(int j = 0; j < cols; ++j) {
201 const int towers{twrs[i][j]};
202 sigma += twrs[i][j] < m_rhoPlusThr ? towers*towers: 0;
203 }
204 }
205
206 return static_cast<int>(std::sqrt(sigma * 1. / n));
207
208}
pdg_comparison.sigma
sigma
Definition pdg_comparison.py:324

◆ initialize()

StatusCode LVL1::gFEXaltMetAlgo::initialize ( )
overridevirtual

standard Athena-Algorithm method

Definition at line 24 of file gFEXaltMetAlgo.cxx.

24 {
25
26 return StatusCode::SUCCESS;
27
28}

◆ metFPGA()

void LVL1::gFEXaltMetAlgo::metFPGA ( const gTowersCentral & twrs,
int & MET_x,
int & MET_y,
const unsigned short FPGA_NO ) const
private

Definition at line 142 of file gFEXaltMetAlgo.cxx.

142 {
143
144 int rows = twrs.size();
145 int cols = twrs[0].size();
146 for( int irow = 0; irow < rows; irow++ ){
147 for(int jcolumn = 0; jcolumn<cols; jcolumn++){
148 bool filter{twrs[irow][jcolumn] > m_etaThr[FPGA_NO][jcolumn]};
149 MET_x += filter ? (twrs[irow][jcolumn])*cosLUT(irow, 5) : 0;
150 MET_y += filter ? (twrs[irow][jcolumn])*sinLUT(irow, 5) : 0;
151
152 }
153 }
154}
LVL1::gFEXaltMetAlgo::cosLUT
float cosLUT(const unsigned int phiIDX, const unsigned int aw) const
Definition gFEXaltMetAlgo.cxx:274
LVL1::gFEXaltMetAlgo::sinLUT
float sinLUT(const unsigned int phiIDX, const unsigned int aw) const
Definition gFEXaltMetAlgo.cxx:262
LVL1::gFEXaltMetAlgo::m_etaThr
std::array< std::vector< int >, 2 > m_etaThr
Definition gFEXaltMetAlgo.h:46

◆ metTotal()

void LVL1::gFEXaltMetAlgo::metTotal ( const int A_MET_x,
const int A_MET_y,
const int B_MET_x,
const int B_MET_y,
int & MET_x,
int & MET_y,
int & MET ) const
inlineprivate

Definition at line 157 of file gFEXaltMetAlgo.cxx.

159 {
160
161 MET_x = A_MET_x + B_MET_x;
162 MET_y = A_MET_y + B_MET_y;
163
164 if (MET_x < -0x0007FF) MET_x = -0x0007FF;
165 if (MET_y < -0x0007FF) MET_y = -0x0007FF;
166
167 if (MET_x > 0x0007FF) MET_x = 0x0007FF;
168 if (MET_y > 0x0007FF) MET_y = 0x0007FF;
169
170 int MET2 = MET_x * MET_x + MET_y * MET_y;
171
172 if (MET2 > 0x000FFF) MET = 0x000FFF;
173 else if (MET2 < 0) MET = 0x000FFF;
174 else MET = std::sqrt(MET2);
175
176}

◆ rho_MET()

void LVL1::gFEXaltMetAlgo::rho_MET ( const gTowersCentral & twrs,
int & MET_x,
int & MET_y,
const int rho,
const int sigma ) const
private

Definition at line 211 of file gFEXaltMetAlgo.cxx.

211 {
212
213 int rows = twrs.size();
214 int cols = twrs[0].size();
215 for( int irow = 0; irow < rows; irow++ ){
216 for(int jcolumn = 0; jcolumn < cols; jcolumn++){
217 const int ET_gTower_sub{(twrs[irow][jcolumn] - rho) & 0xFFF};
218 const bool filter{ET_gTower_sub > sigma && !(ET_gTower_sub & 0x800)};
219 MET_x += filter ? (ET_gTower_sub)*cosLUT(irow, 5) : 0;
220 MET_y += filter ? (ET_gTower_sub)*sinLUT(irow, 5) : 0;
221
222 }
223 }
224}

◆ setAlgoConstant()

void LVL1::gFEXaltMetAlgo::setAlgoConstant ( std::vector< int > && A_thr,
std::vector< int > && B_thr,
const int rhoPlusThr )
overridevirtual

Definition at line 31 of file gFEXaltMetAlgo.cxx.

33 {
34 m_etaThr[0] = std::move(A_thr);
35 m_etaThr[1] = std::move(B_thr);
36 m_rhoPlusThr = rhoPlusThr;
37}

◆ sinLUT()

float LVL1::gFEXaltMetAlgo::sinLUT ( const unsigned int phiIDX,
const unsigned int aw ) const
private

Definition at line 262 of file gFEXaltMetAlgo.cxx.

263{
264 float c = ((float)phiIDX)/std::pow(2,aw);
265 float rad = (2*M_PI) *c;
266 float rsin = std::sin(rad);
267 return rsin;
268}

◆ sumEt()

int LVL1::gFEXaltMetAlgo::sumEt ( const int A_sumEt,
const int B_sumEt ) const
private

Definition at line 253 of file gFEXaltMetAlgo.cxx.

253 {
254
255 return A_sumEt + B_sumEt;
256}

◆ sumEtFPGAnc()

int LVL1::gFEXaltMetAlgo::sumEtFPGAnc ( const gTowersCentral & twrs,
const unsigned short FPGA_NO ) const
private

Definition at line 226 of file gFEXaltMetAlgo.cxx.

226 {
227
228 int partial_sumEt = 0;
229 const int rows = twrs.size();
230 const int cols = twrs[0].size();
231 for( int irow = 0; irow < rows; irow++ ){
232 for(int jcolumn = 0; jcolumn<cols; jcolumn++){
233 partial_sumEt += twrs[irow][jcolumn] > m_etaThr[FPGA_NO][jcolumn] ? twrs[irow][jcolumn] : 0;
234 }
235 }
236 return partial_sumEt;
237}

◆ sumEtFPGArms()

int LVL1::gFEXaltMetAlgo::sumEtFPGArms ( const gTowersCentral & twrs,
const int sigma ) const
private

Definition at line 239 of file gFEXaltMetAlgo.cxx.

239 {
240
241 int partial_sumEt = 0;
242 const int rows = twrs.size();
243 const int cols = twrs[0].size();
244 for(int i{0}; i < rows; ++i){
245 for(int j{0}; j < cols; ++j) {
246 partial_sumEt += twrs[i][j] > sigma ? twrs[i][j] : 0;
247 }
248 }
249 return partial_sumEt;
250}

Member Data Documentation

◆ m_etaThr

std::array<std::vector<int>, 2> LVL1::gFEXaltMetAlgo::m_etaThr
private

Definition at line 46 of file gFEXaltMetAlgo.h.

◆ m_rhoPlusThr

int LVL1::gFEXaltMetAlgo::m_rhoPlusThr {}
private

Definition at line 47 of file gFEXaltMetAlgo.h.

47{};
The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0