#include <Phase2L1CaloEGammaUtils.h>

Public Member Functions
void	computeClusterIsolationInPlace (int nGCTCard)

void	setIsolationInfo (void)

Public Attributes
GCTCorrfiber_t	GCTCorrfiber [N_GCTINTERNAL_FIBERS]

Detailed Description

Definition at line 1546 of file Phase2L1CaloEGammaUtils.h.

Member Function Documentation

◆ computeClusterIsolationInPlace()

void p2eg::GCTinternal_t::computeClusterIsolationInPlace ( int nGCTCard )

inline

Definition at line 1550 of file Phase2L1CaloEGammaUtils.h.

References p2eg::GCTtower_t::ecalEt, EgHLTOffHistBins_cfi::et, p2eg::GCTcluster_t::et, nano_mu_digi_cff::float, p2eg::GCTCorrfiber_t::GCTclusters, GCTCorrfiber, p2eg::GCTCorrfiber_t::GCTtowers, l1tTowerCalibrationProducer_cfi::iEta, p2eg::GCTcluster_t::inCardToweriEta(), p2eg::GCTcluster_t::inCardToweriPhi(), p2eg::GCTcluster_t::iso, SiStripPI::max, SiStripPI::min, p2eg::N_GCTCLUSTERS_FIBER, p2eg::N_GCTETA, p2eg::N_GCTINTERNAL_FIBERS, p2eg::N_GCTPHI, p2eg::N_GCTPOSITIVE_FIBERS, p2eg::N_GCTTOWERS_CLUSTER_ISO_ONESIDE, p2eg::N_GCTTOWERS_FIBER, and scaleFactor.

Referenced by p2eg::algo_top().

                                                       {
       for (unsigned int iFiber = 0; iFiber < N_GCTINTERNAL_FIBERS; iFiber++) {
         for (unsigned int iCluster = 0; iCluster < N_GCTCLUSTERS_FIBER; iCluster++) {
           // We will only save clusters with > 0 GeV, so only need to do this for clusters with >0 energy
           if (GCTCorrfiber[iFiber].GCTclusters[iCluster].et == 0) {
             GCTCorrfiber[iFiber].GCTclusters[iCluster].iso = 0;
             continue;
           }
 
           ap_uint<12> uint_isolation = 0;
 
           // do not add the GCT card off-set, so we remain in the gct local card iEta/iPhi
           int toweriEta_in_GCT_card = GCTCorrfiber[iFiber].GCTclusters[iCluster].inCardToweriEta();
           int toweriPhi_in_GCT_card = GCTCorrfiber[iFiber].GCTclusters[iCluster].inCardToweriPhi();
 
           // If cluster is in the overlap region, do not compute isolation
           bool inOverlapWithAnotherGCTCard = (((toweriPhi_in_GCT_card >= 0) && (toweriPhi_in_GCT_card < 4)) ||
                                               ((toweriPhi_in_GCT_card >= 28) && (toweriPhi_in_GCT_card < 32)));
           if (inOverlapWithAnotherGCTCard) {
             GCTCorrfiber[iFiber].GCTclusters[iCluster].iso = 0;
             continue;
           }
 
           // Size 5x5 in towers: include the overlap-region-between-GCT-cards-if-applicable. In eta direction, the min and max towers (inclusive!) are:
           int isoWindow_toweriEta_in_GCT_card_min = std::max(0, toweriEta_in_GCT_card - 2);
           int isoWindow_toweriEta_in_GCT_card_max = std::min(toweriEta_in_GCT_card + 2, N_GCTETA - 1);  // N_GCTETA = 34
           // e.g. if our window is centered at tower_iEta = 5, we want to sum towers_iEta 3, 4, (5), 6, 7, inclusive
           // e.g. if our window is near the boundary, tower_iEta = 32, we want to sum towers_iEta 30, 31, (32), 33
           // inclusive (but there are only N_GCTETA = 34 towers, so we stop at tower_iEta = 33)
 
           // in phi direction, the min and max towers (inclusive!) are:
           int isoWindow_toweriPhi_in_GCT_card_min = std::max(0, toweriPhi_in_GCT_card - 2);
           int isoWindow_toweriPhi_in_GCT_card_max = std::min(toweriPhi_in_GCT_card + 2, N_GCTPHI - 1);
 
           // Keep track of the number of towers we summed over
           int nTowersSummed = 0;
 
           // First add any nearby clusters to the isolation
           for (unsigned int candFiber = 0; candFiber < N_GCTINTERNAL_FIBERS; candFiber++) {
             for (unsigned int candCluster = 0; candCluster < N_GCTCLUSTERS_FIBER; candCluster++) {
               // Do not double-count the cluster we are calculating the isolation for
               if (!((candFiber == iFiber) && (candCluster == iCluster))) {
                 // Only consider clusters with et > 0 for isolation sum
                 if (GCTCorrfiber[candFiber].GCTclusters[candCluster].et > 0) {
                   // Get the candidate cluster's tower iEta and iPhi in GCT card
                   int candidate_toweriEta = GCTCorrfiber[candFiber].GCTclusters[candCluster].inCardToweriEta();
                   int candidate_toweriPhi = GCTCorrfiber[candFiber].GCTclusters[candCluster].inCardToweriPhi();
 
                   // If the tower that the candidate cluster is in, is within a 5x5 window, add the candidate cluster energy's to the isolation as a proxy for the ECAL energy
                   if (((candidate_toweriEta >= isoWindow_toweriEta_in_GCT_card_min) &&
                        (candidate_toweriEta <= isoWindow_toweriEta_in_GCT_card_max)) &&
                       ((candidate_toweriPhi >= isoWindow_toweriPhi_in_GCT_card_min) &&
                        (candidate_toweriPhi <= isoWindow_toweriPhi_in_GCT_card_max))) {
                     uint_isolation += GCTCorrfiber[candFiber].GCTclusters[candCluster].et;
                   }
                 }
               }
             }
           }
 
           //  From "tower index in GCT card", get which fiber it is in (out of 64 fibers), and which tower it is inside the fiber (out of 17 towers)
           for (int iEta = isoWindow_toweriEta_in_GCT_card_min; iEta <= isoWindow_toweriEta_in_GCT_card_max; iEta++) {
             for (int iPhi = isoWindow_toweriPhi_in_GCT_card_min; iPhi <= isoWindow_toweriPhi_in_GCT_card_max; iPhi++) {
               nTowersSummed += 1;
 
               int indexInto64Fibers;
               int indexInto17TowersInFiber;
 
               bool isTowerInPositiveEta = (iEta >= N_GCTTOWERS_FIBER);
               if (isTowerInPositiveEta) {
                 // phi index is simple (e.g. if real phi = +80 degrees, iPhi in GCT = 31)
                 indexInto64Fibers = iPhi;
                 // if real eta = 1.47, iEta in GCT card = 33. If real eta = 0.0, iEta in GCT = 17, so iEta in fiber = 17%17 = 0.
                 indexInto17TowersInFiber = (iEta % 17);
               } else {
                 // add offset (e.g. if real phi = +80 degrees, iPhi in GCT = 31, and my index into GCT fibers 31 + 32 = 63)
                 indexInto64Fibers = (iPhi + N_GCTPOSITIVE_FIBERS);
                 // e.g.  if real eta = 0, iEta innew GCT card = 16, i.e. our index into the GCT fiber is 16-16 = 0
                 indexInto17TowersInFiber = (16 - iEta);
               }
 
               ap_uint<12> ecalEt = GCTCorrfiber[indexInto64Fibers].GCTtowers[indexInto17TowersInFiber].ecalEt;
               uint_isolation += ecalEt;
             }
           }
 
           // Scale the isolation sum up if we summed over fewer than (5x5) = 25 towers
           float scaleFactor =
               ((float)(N_GCTTOWERS_CLUSTER_ISO_ONESIDE * N_GCTTOWERS_CLUSTER_ISO_ONESIDE) / (float)nTowersSummed);
 
           uint_isolation = (ap_uint<12>)(((float)uint_isolation) * scaleFactor);
 
           // Set the iso in the cluster
           GCTCorrfiber[iFiber].GCTclusters[iCluster].iso = uint_isolation;
         }
       }
     }

◆ setIsolationInfo()

void p2eg::GCTinternal_t::setIsolationInfo ( void )

inline

Definition at line 1648 of file Phase2L1CaloEGammaUtils.h.

References p2eg::GCTCorrfiber_t::GCTclusters, GCTCorrfiber, p2eg::N_GCTCLUSTERS_FIBER, p2eg::N_GCTINTERNAL_FIBERS, and p2eg::GCTcluster_t::setRelIsoAndFlags().

Referenced by p2eg::algo_top().

                                 {
       for (unsigned int iFiber = 0; iFiber < N_GCTINTERNAL_FIBERS; iFiber++) {
         for (unsigned int iCluster = 0; iCluster < N_GCTCLUSTERS_FIBER; iCluster++) {
           // update the cluster's isolation information
           GCTCorrfiber[iFiber].GCTclusters[iCluster].setRelIsoAndFlags();
         }
       }
     }

Member Data Documentation

◆ GCTCorrfiber

GCTCorrfiber_t p2eg::GCTinternal_t::GCTCorrfiber[N_GCTINTERNAL_FIBERS]

Definition at line 1548 of file Phase2L1CaloEGammaUtils.h.

Referenced by computeClusterIsolationInPlace(), p2eg::getClustersTowers(), p2eg::getFullTowers(), setIsolationInfo(), and p2eg::writeToCorrelatorAndGTOutputs().