21 declareInterface<IjFEXSmallRJetAlgo>(
this);
30 ATH_CHECK(m_jTowerContainerKey.initialize());
32 return StatusCode::SUCCESS;
40 if(! m_jTowerContainer.isValid()) {
41 ATH_MSG_ERROR(
"Could not retrieve jTowerContainer " << m_jTowerContainerKey.key());
42 return StatusCode::FAILURE;
45 return StatusCode::SUCCESS;
50 for(
int phi=0; phi<7; phi++) {
51 for (
int eta=0; eta<7; eta++) {
52 m_jFEXalgoTowerID[phi][eta] = inputTable[6-phi][eta];
56 for(
int phi=0; phi<7; phi++) {
57 for (
int eta=0; eta<7; eta++) {
58 m_jFEXalgoTowerID_displaced[phi][eta] = inputTableDisplaced[6-phi][eta];
69 if(m_map_Etvalues.find(TTID) != m_map_Etvalues.end()) {
70 return m_map_Etvalues.at(TTID).at(0);
81 for(
int mphi = 1; mphi < 6; mphi++) {
82 for(
int meta = 1; meta< 6; meta++) {
85 int seedTotalET_displaced = 0;
86 for(
int iphi = -1; iphi < 2; iphi++) {
87 for(
int ieta = -1; ieta < 2; ieta++) {
90 seedTotalET += getTTowerET(m_jFEXalgoTowerID [mphi + iphi][meta + ieta]);
91 seedTotalET_displaced += getTTowerET(m_jFEXalgoTowerID_displaced[mphi + iphi][meta + ieta]);
94 m_jFEXalgoSearchWindowSeedET[mphi -1][meta -1] = seedTotalET;
95 m_jFEXalgoSearchWindowSeedET_displaced[mphi -1][meta -1] = seedTotalET_displaced;
104 int central_seed = mymatrix[2][2];
105 for (
int iphi = 0; iphi < 5; iphi++) {
106 for (
int ieta = 0; ieta < 5; ieta++) {
108 if ((ieta == 2) && (iphi == 2)) {
112 if( (iphi > ieta) || (iphi == 0 && ieta == 0) || (iphi == 1 && ieta == 1) ) {
113 if(central_seed < mymatrix[iphi][ieta]) {
118 if((iphi < ieta) || (iphi == 3 && ieta == 3) || (iphi == 4 && ieta == 4)) {
119 if(central_seed <= mymatrix[iphi][ieta]) {
132 bool isCentralLM = CalculateLM(m_jFEXalgoSearchWindowSeedET) &&
133 ( getTTowerET(m_jFEXalgoTowerID[3][3]) >= getTTowerET(m_jFEXalgoTowerID[4][2]) || m_jFEXalgoSearchWindowSeedET[2][2] > m_jFEXalgoSearchWindowSeedET[3][1]);
135 bool isDisplacedLM = CalculateLM(m_jFEXalgoSearchWindowSeedET_displaced) &&
136 ( getTTowerET(m_jFEXalgoTowerID[3][3]) > getTTowerET(m_jFEXalgoTowerID[2][4]) && m_jFEXalgoSearchWindowSeedET[2][2] == m_jFEXalgoSearchWindowSeedET[1][3]);
138 if(isCentralLM || isDisplacedLM ){
140 return m_jFEXalgoSearchWindowSeedET[2][2] > seedThreshold;
147 return m_jFEXalgoSearchWindowSeedET[2][2];
154 int SRJetClusterET = 0;
155 for(
int nphi = -3; nphi< 4; nphi++) {
156 for(
int neta = -3; neta< 4; neta++) {
158 if(DeltaRSquared < 16) {
159 SRJetClusterET += getTTowerET(m_jFEXalgoTowerID[3+nphi][3+neta]);
163 return SRJetClusterET;
168 m_JetSaturation =
false;
169 for(
int nphi = -3; nphi< 4; nphi++) {
170 for(
int neta = -3; neta< 4; neta++) {
172 if(DeltaRSquared < 16) {
173 m_JetSaturation = m_JetSaturation || getTTowerSat(m_jFEXalgoTowerID[3+nphi][3+neta]);
180 return m_JetSaturation;
184 int SmallETRing = getSmallClusterET() - m_jFEXalgoSearchWindowSeedET[3][3];
189 return m_jFEXalgoTowerID[3][3];
194 m_map_Etvalues=et_map;
203 const LVL1::jTower * tmpTower = m_jTowerContainer->findTower(TTID);