32 m_tower_id(id_modifier + iphi + (nphi * ieta)),
89 if (SCpass){
m_et = outET; }
186 return getFWID(iPhiFW, iEtaFW);
213 bool is_central =
true;
214 if (iEta <= 7 || iEta >= 32) is_central =
false;
222 gFEXtowerID = (
iEta - 8) + (
iPhi * 12);
226 gFEXtowerID = 10000 + (
iEta - 20) + (
iPhi * 12);
237 gFEXtowerID = gFEXtowerID + (
iPhi*24);
241 else if (
iEta == 1 ){
242 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 12);
246 else if (
iEta == 2 ){
247 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 1);
251 else if (
iEta == 3 ){
252 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 13);
256 else if (
iEta >= 4 and
iEta <= 7 ){
257 gFEXtowerID = gFEXtowerID + ((
iPhi*12) + (
iEta -2));
265 gFEXtowerID = gFEXtowerID + (
iPhi*12) + (
iEta -32 +6);
267 else if (
iEta == 36 ){
268 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 22);
272 else if (
iEta == 37 ){
273 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 10);
277 else if (
iEta == 38 ){
278 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 23);
282 else if (
iEta == 39 ){
283 gFEXtowerID = gFEXtowerID + ((
iPhi*24) + 11);
297 float s_centralPhiWidth = (2*
M_PI)/32;
298 float s_forwardPhiWidth = (2*
M_PI)/16;
300 constexpr std::array<float, 40> s_EtaCenter = { -4.7, -4.2, -3.7, -3.4, -3.2, -3,
301 -2.8, -2.6, -2.35, -2.1, -1.9, -1.7, -1.5, -1.3, -1.1, -0.9,
302 -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1,
303 1.3, 1.5, 1.7, 1.9, 2.1, 2.35, 2.6, 2.8, 3.0,
304 3.2, 3.4, 3.7, 4.2, 4.7};
308 float Phi_gFex = -99;
309 if ((
iEta <= 3 ) || ( (
iEta >= 36) )){
310 Phi_gFex = ( (
iPhi * s_forwardPhiWidth) + s_forwardPhiWidth/2);
313 Phi_gFex = ( (
iPhi * s_centralPhiWidth) + s_centralPhiWidth/2);
316 if (Phi_gFex <
M_PI) {