Classes
class	card

class	Cluster

class	clusterInfo

class	crystal

class	crystalMax

class	ecalRegion_t

class	ecaltp_t

class	etaStrip_t

class	etaStripPeak_t

class	GCTcard_t

class	GCTcluster_t

class	GCTCorrfiber_t

class	GCTinternal_t

class	GCTintTowers_t

class	GCTtoCorr_t

class	GCTtower_t

class	linkECAL

class	RCTcard_t

class	RCTcluster_t

class	RCTtoGCTfiber_t

class	RCTtower_t

class	region3x4

class	SimpleCaloHit

class	tower_t

class	towerHCAL

class	towers3x4

Functions
void	algo_top (const GCTcard_t &GCTcard, GCTtoCorr_t &GCTtoCorr, unsigned int nGCTCard, std::unique_ptr< l1tp2::CaloCrystalClusterCollection > const &gctClusters, std::unique_ptr< l1tp2::CaloTowerCollection > const &gctTowers, std::unique_ptr< l1tp2::CaloTowerCollection > const &gctFullTowers, std::unique_ptr< l1t::EGammaBxCollection > const &gctEGammas, std::unique_ptr< l1tp2::DigitizedClusterCorrelatorCollection > const &gctDigitizedClustersCorrelator, std::unique_ptr< l1tp2::DigitizedTowerCorrelatorCollection > const &gctDigitizedTowersCorrelator, std::unique_ptr< l1tp2::DigitizedClusterGTCollection > const &gctDigitizedClustersGT, l1tp2::ParametricCalibration calib_)

ecaltp_t	bestOf2 (const ecaltp_t ecaltp0, const ecaltp_t ecaltp1)

bool	compareClusterET (const Cluster &lhs, const Cluster &rhs)

void	doProximityAndBremsStitching (const RCTcard_t(&inputCards)[N_RCTCARDS_PHI], RCTcard_t(&outputCards)[N_RCTCARDS_PHI], int iStartingCard, bool isPositiveEta)

int	getAbsID_iEta_fromFirmwareCardTowerLink (int nCard, int nTower, int nLink)

int	getAbsID_iPhi_fromFirmwareCardTowerLink (int nCard, int nTower, int nLink)

clusterInfo	getBremsValuesNeg (crystal tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI], ap_uint< 5 > seed_eta, ap_uint< 5 > seed_phi)

clusterInfo	getBremsValuesPos (crystal tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI], ap_uint< 5 > seed_eta, ap_uint< 5 > seed_phi)

int	getCard_iEtaMax (int cc)

int	getCard_iEtaMin (int cc)

int	getCard_iPhiMax (int cc)

int	getCard_iPhiMin (int cc)

int	getCard_refCrystal_iEta (int cc)

int	getCard_refCrystal_iPhi (int cc)

Cluster	getClusterFromRegion3x4 (crystal temp[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI])

clusterInfo	getClusterPosition (const ecalRegion_t ecalRegion)

GCTcard_t	getClustersCombined (const GCTcard_t &GCTcard, unsigned int nGCTCard)

GCTinternal_t	getClustersTowers (const GCTcard_t &GCTcard, unsigned int nGCTCard)

clusterInfo	getClusterValues (crystal tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI], ap_uint< 5 > seed_eta, ap_uint< 5 > seed_phi)

void	getECALTowersEt (crystal tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI], ap_uint< 12 > towerEt[12])

GCTintTowers_t	getFullTowers (const GCTinternal_t &GCTinternal)

crystalMax	getPeakBin15N (const etaStripPeak_t etaStrip)

ecaltp_t	getPeakBin20N (const etaStrip_t etaStrip)

int	getRegionNumber (const int local_iEta)

int	getTower_absEtaID (float eta)

int	getTower_absPhiID (float phi)

float	getTowerEta_fromAbsID (int id)

float	getTowerPhi_fromAbsID (int id)

ecalRegion_t	initStructure (crystal temporary[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI])

bool	isValidCard (int cc)

Cluster	packCluster (ap_uint< 15 > &clusterEt, ap_uint< 5 > &etaMax_t, ap_uint< 5 > &phiMax_t)

bool	passes_iso (float pt, float iso)

bool	passes_looseTkiso (float pt, float iso)

bool	passes_looseTkss (float pt, float eta, float ss)

bool	passes_ss (float pt, float eta, float ss)

void	printl1tp2TowerInfo (l1tp2::CaloTower thisTower, std::string description="")

void	removeClusterFromCrystal (crystal temp[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI], ap_uint< 5 > seed_eta, ap_uint< 5 > seed_phi, ap_uint< 2 > brems)

void	stitchClusterOverRegionBoundary (std::vector< Cluster > &cluster_list, int towerEtaUpper, int towerEtaLower, int cc)

void	writeToCorrelatorAndGTOutputs (const GCTinternal_t &GCTinternal, GCTtoCorr_t &GCTtoCorrOutput, std::unique_ptr< l1tp2::CaloCrystalClusterCollection > const &gctClustersOutput, std::unique_ptr< l1tp2::CaloTowerCollection > const &gctTowersOutput, std::unique_ptr< l1t::EGammaBxCollection > const &gctEGammas, std::unique_ptr< l1tp2::DigitizedClusterCorrelatorCollection > const &gctDigitizedClustersCorrelator, std::unique_ptr< l1tp2::DigitizedTowerCorrelatorCollection > const &gctDigitizedTowersCorrelator, std::unique_ptr< l1tp2::DigitizedClusterGTCollection > const &gctDigitizedClustersGT, int nGCTCard, int fiberStart, int fiberEnd, int corrFiberIndexOffset, int corrTowPhiOffset)

Variables
static constexpr float	a0 = 0.21

static constexpr float	a0_80 = 0.85

static constexpr float	a1_80 = 0.0080

static constexpr float	b0 = 0.38

static constexpr float	b1 = 1.9

static constexpr float	b2 = 0.05

static constexpr float	c0_ss = 0.94

static constexpr float	c1_ss = 0.052

static constexpr float	c2_ss = 0.044

static constexpr int	CRYSTAL_IN_ETA = 15

static constexpr int	CRYSTAL_IN_PHI = 20

static constexpr int	CRYSTALS_IN_TOWER_ETA = 5

static constexpr int	CRYSTALS_IN_TOWER_PHI = 5

static constexpr float	cut_500_MeV = 0.5

static constexpr float	d0 = 0.96

static constexpr float	d1 = 0.0003

static constexpr float	e0_looseTkss = 0.944

static constexpr float	e1_looseTkss = 0.65

static constexpr float	e2_looseTkss = 0.4

static constexpr float	ECAL_eta_range = 1.4841

static constexpr float	ECAL_LSB = 0.5

static constexpr float	eta0_loose_ss = 0.4

static constexpr float	eta1_loose_ss = 0.5

static constexpr float	eta2_loose_ss = 0.75

static constexpr float	eta3_loose_ss = 0.85

static constexpr int	GCTCARD_0_TOWER_IPHI_OFFSET = 20

static constexpr int	GCTCARD_1_TOWER_IPHI_OFFSET = 44

static constexpr int	GCTCARD_2_TOWER_IPHI_OFFSET = 68

static const unsigned int	GCTcardtoRCTcardnumber [N_GCTCARDS][N_RCTCARDS_PHI *2]

static constexpr float	half_crystal_size = 0.00873

static constexpr float	HCAL_LSB = 0.5

static constexpr float	high_pt_threshold = 130.

static constexpr float	loose_ss_offset = 0.89

static constexpr int	n_clusters_4link = 8

static constexpr int	n_clusters_link = 2

static constexpr int	N_CLUSTERS_PER_REGION = 4

static constexpr int	n_crystals_cardEta = (n_towers_Eta * n_towers_cardEta)

static constexpr int	n_crystals_cardPhi = (n_towers_Phi * n_towers_cardPhi)

static constexpr int	N_GCTCARDS = 3

static constexpr int	N_GCTCLUSTERS_FIBER = 2

static constexpr int	N_GCTCORR_FIBERS = 48

static constexpr int	N_GCTETA = 34

static constexpr int	N_GCTINTERNAL_FIBERS = 64

static constexpr int	N_GCTPHI = 32

static constexpr int	N_GCTPOSITIVE_FIBERS = 32

static constexpr int	N_GCTTOWERS_CLUSTER_ISO_ONESIDE = 5

static constexpr int	N_GCTTOWERS_FIBER = 17

static constexpr int	n_links_card = 4

static constexpr int	N_RCTCARDS_PHI = 8

static constexpr int	N_RCTCLUSTERS_FIBER = 2

static constexpr int	N_RCTGCT_FIBERS = 4

static constexpr int	N_RCTTOWERS_FIBER = 17

static constexpr int	N_REGIONS_PER_CARD = 6

static constexpr int	n_towers_cardEta = 17

static constexpr int	n_towers_cardPhi = 4

static constexpr int	n_towers_Eta = 34

static constexpr int	n_towers_halfPhi = 36

static constexpr int	n_towers_per_link = 17

static constexpr int	n_towers_Phi = 72

static constexpr float	slideIsoPtThreshold = 80

static constexpr int	TOWER_IN_ETA = 3

static constexpr int	TOWER_IN_PHI = 4

Function Documentation

◆ algo_top()

void p2eg::algo_top	(	const GCTcard_t &	GCTcard,
		p2eg::GCTtoCorr_t &	GCTtoCorr,
		unsigned int	nGCTCard,
		std::unique_ptr< l1tp2::CaloCrystalClusterCollection > const &	gctClusters,
		std::unique_ptr< l1tp2::CaloTowerCollection > const &	gctTowers,
		std::unique_ptr< l1tp2::CaloTowerCollection > const &	gctFullTowers,
		std::unique_ptr< l1t::EGammaBxCollection > const &	gctEGammas,
		std::unique_ptr< l1tp2::DigitizedClusterCorrelatorCollection > const &	gctDigitizedClustersCorrelator,
		std::unique_ptr< l1tp2::DigitizedTowerCorrelatorCollection > const &	gctDigitizedTowersCorrelator,
		std::unique_ptr< l1tp2::DigitizedClusterGTCollection > const &	gctDigitizedClustersGT,
		l1tp2::ParametricCalibration	calib_
	)

inline

Definition at line 363 of file Phase2L1GCT.h.

References p2eg::GCTinternal_t::computeClusterIsolationInPlace(), l1tPhase2CaloJetEmulator_cfi::gctFullTowers, getClustersTowers(), getFullTowers(), N_GCTINTERNAL_FIBERS, N_GCTPOSITIVE_FIBERS, N_RCTGCT_FIBERS, p2eg::GCTinternal_t::setIsolationInfo(), writeToCorrelatorAndGTOutputs(), and p2eg::GCTintTowers_t::writeToPFOutput().

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                        {
   //-------------------------//
   // Initialize the GCT area
   //-------------------------//
   p2eg::GCTinternal_t GCTinternal = p2eg::getClustersTowers(GCTcard, nGCTCard);
 
   //------------------------------------------------//
   // Combine towers and clusters to get full towers
   //------------------------------------------------//
   p2eg::GCTintTowers_t GCTintTowers = p2eg::getFullTowers(GCTinternal);
 
   //---------------------------//
   // Compute cluster isolation
   //--------------------------//
   GCTinternal.computeClusterIsolationInPlace(nGCTCard);
   GCTinternal.setIsolationInfo();
 
   //-----------------------------------------------------------------------------------------------------------------------//
   // Output to correlator and CMSSW collections.
   // For positive eta, skip overlap region, i.e. fibers i = 0, 1, 2, 3, and i = 28, 29, 30, 31.
   // For negative eta, skip overlap region, i.e. fibers 32, 33, 34, 35, and 61, 62, 63, 64.
   //-----------------------------------------------------------------------------------------------------------------------//
   int posEtaFiberStart = p2eg::N_RCTGCT_FIBERS;  // 4, since there are 4 fibers in one RCT card
   int posEtaFiberEnd = (p2eg::N_GCTPOSITIVE_FIBERS - p2eg::N_RCTGCT_FIBERS);
   int negEtaFiberStart = (p2eg::N_GCTPOSITIVE_FIBERS + p2eg::N_RCTGCT_FIBERS);
   int negEtaFiberEnd =
       (p2eg::N_GCTINTERNAL_FIBERS - p2eg::N_RCTGCT_FIBERS);  // first term is number of gct internal fibers
 
   // When indexing into the correlator output, note that the output to correlator does NOT include the overlap region,
   // so fiber number "i" in GCT is not fiber "i" in the correlator output, it's reduced by 4 in positive eta, and 12 in negative eta.
   // (4 because we are skipping one RCT card in positive eta, and 12 because we are skipping three RCT cards in negative eta)
   int posEtaCorrelatorFiberIndexOffset = 4;
   int negEtaCorrelatorFiberIndexOffset = 12;
 
   // The offset in the actual towPhi value is going to be the same in pos/neg eta; shifted down by 4 due to no overlap region
   int correlatorTowPhiOffset = 4;
 
   // Positive eta
   p2eg::writeToCorrelatorAndGTOutputs(GCTinternal,
                                       GCTtoCorr,
                                       gctClusters,
                                       gctTowers,
                                       gctEGammas,
                                       gctDigitizedClustersCorrelator,
                                       gctDigitizedTowersCorrelator,
                                       gctDigitizedClustersGT,
                                       nGCTCard,
                                       posEtaFiberStart,
                                       posEtaFiberEnd,
                                       posEtaCorrelatorFiberIndexOffset,
                                       correlatorTowPhiOffset);
   // Negative eta
   p2eg::writeToCorrelatorAndGTOutputs(GCTinternal,
                                       GCTtoCorr,
                                       gctClusters,
                                       gctTowers,
                                       gctEGammas,
                                       gctDigitizedClustersCorrelator,
                                       gctDigitizedTowersCorrelator,
                                       gctDigitizedClustersGT,
                                       nGCTCard,
                                       negEtaFiberStart,
                                       negEtaFiberEnd,
                                       negEtaCorrelatorFiberIndexOffset,
                                       correlatorTowPhiOffset);
 
   //-----------------------------------------------------------------------------------------------------------------------//
   // CMSSW outputs for GCT Full Towers (clusters + towers) output for PFClusters.
   //-----------------------------------------------------------------------------------------------------------------------//
   GCTintTowers.writeToPFOutput(nGCTCard, gctFullTowers);
 }

◆ bestOf2()

p2eg::ecaltp_t p2eg::bestOf2	(	const ecaltp_t	ecaltp0,
		const ecaltp_t	ecaltp1
	)

inline

Definition at line 940 of file Phase2L1RCT.h.

References p2eg::ecaltp_t::energy, and x.

Referenced by getPeakBin15N(), and getPeakBin20N().

                                                                                         {
   p2eg::ecaltp_t x;
   x = (ecaltp0.energy > ecaltp1.energy) ? ecaltp0 : ecaltp1;
 
   return x;
 }

◆ compareClusterET()

bool p2eg::compareClusterET	(	const Cluster &	lhs,
		const Cluster &	rhs
	)

inline

Definition at line 873 of file Phase2L1CaloEGammaUtils.h.

References p2eg::Cluster::clusterEnergy().

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                                        {
     return (lhs.clusterEnergy() > rhs.clusterEnergy());
   }

◆ doProximityAndBremsStitching()

void p2eg::doProximityAndBremsStitching	(	const RCTcard_t(&)	inputCards[N_RCTCARDS_PHI],
		RCTcard_t(&)	outputCards[N_RCTCARDS_PHI],
		int	iStartingCard,
		bool	isPositiveEta
	)

Referenced by getClustersCombined().

◆ getAbsID_iEta_fromFirmwareCardTowerLink()

int p2eg::getAbsID_iEta_fromFirmwareCardTowerLink	(	int	nCard,
		int	nTower,
		int	nLink
	)

inline

Definition at line 191 of file Phase2L1CaloEGammaUtils.h.

References n_towers_per_link, and funct::void.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                                                        {
     // iEta only depends on the tower position in the link
     (void)nCard;
     (void)nLink;
     if ((nCard % 2) == 1) {  // if cc is odd (positive eta), e.g. nTower = 0 will correspond to absolute iEta ID 17.
       return n_towers_per_link + nTower;
     } else {  // if cc is even (negative eta): e.g. nTower = 0 will correspond to absolute iEta ID 16.
       return (16 - nTower);
     }
   }

◆ getAbsID_iPhi_fromFirmwareCardTowerLink()

int p2eg::getAbsID_iPhi_fromFirmwareCardTowerLink	(	int	nCard,
		int	nTower,
		int	nLink
	)

inline

Definition at line 205 of file Phase2L1CaloEGammaUtils.h.

References TOWER_IN_PHI, and funct::void.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                                                        {
     // iPhi only depends on the card and link number
     (void)nTower;
     if ((nCard % 2) == 1) {  // if cc is odd (positive eta),
       // e.g. cc=3, link #2, global iPhi = int(3/2) * 4 + 2 = 6
       return (int(nCard / 2) * TOWER_IN_PHI) + nLink;
     } else {  // if cc is even (negative eta)
       // e.g. cc=4, link #2, global iPhi = int(4/2) * 4 + (4 - 2 - 1)
       //                                 = 2*4 + 1
       //                                 = 9
       // minus one is because TOWER_IN_PHI is 4
       return (int(nCard / 2) * TOWER_IN_PHI) + (TOWER_IN_PHI - nLink - 1);
     }
   }

◆ getBremsValuesNeg()

clusterInfo p2eg::getBremsValuesNeg	(	crystal	tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI],
		ap_uint< 5 >	seed_eta,
		ap_uint< 5 >	seed_phi
	)

◆ getBremsValuesPos()

clusterInfo p2eg::getBremsValuesPos	(	crystal	tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI],
		ap_uint< 5 >	seed_eta,
		ap_uint< 5 >	seed_phi
	)

◆ getCard_iEtaMax()

int p2eg::getCard_iEtaMax ( int cc )

inline

Definition at line 108 of file Phase2L1CaloEGammaUtils.h.

References cms::cuda::assert(), gpuPixelDoublets::cc, CRYSTALS_IN_TOWER_ETA, TrackingDataMCValidation_Standalone_cff::etamax, isValidCard(), n_towers_cardEta, and n_towers_Eta.

Referenced by p2eg::SimpleCaloHit::isInCard().

                                      {
     assert(isValidCard(cc));
 
     int etamax = 0;
     if (cc % 2 == 0)                                            // Even card: negative eta
       etamax = (n_towers_cardEta * CRYSTALS_IN_TOWER_ETA - 1);  // First eta half. 5 crystals in eta in 1 tower.
     else                                                        // Odd card: positive eta
       etamax = (n_towers_Eta * CRYSTALS_IN_TOWER_ETA - 1);
     return etamax;
   }

◆ getCard_iEtaMin()

int p2eg::getCard_iEtaMin ( int cc )

inline

Definition at line 120 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc, CRYSTALS_IN_TOWER_ETA, TrackingDataMCValidation_Standalone_cff::etamin, and n_towers_cardEta.

Referenced by p2eg::SimpleCaloHit::isInCard().

                                      {
     int etamin = 0;
     if (cc % 2 == 0)  // Even card: negative eta
       etamin = (0);
     else  // Odd card: positive eta
       etamin = (n_towers_cardEta * CRYSTALS_IN_TOWER_ETA);
     return etamin;
   }

◆ getCard_iPhiMax()

int p2eg::getCard_iPhiMax ( int cc )

inline

Definition at line 130 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc, CRYSTALS_IN_TOWER_PHI, and phimax.

Referenced by p2eg::SimpleCaloHit::isInCard().

                                      {
     int phimax = ((cc / 2) + 1) * 4 * CRYSTALS_IN_TOWER_PHI - 1;
     return phimax;
   }

◆ getCard_iPhiMin()

int p2eg::getCard_iPhiMin ( int cc )

inline

Definition at line 136 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc, CRYSTALS_IN_TOWER_PHI, and phimin.

Referenced by p2eg::SimpleCaloHit::isInCard().

                                      {
     int phimin = (cc / 2) * 4 * CRYSTALS_IN_TOWER_PHI;
     return phimin;
   }

◆ getCard_refCrystal_iEta()

int p2eg::getCard_refCrystal_iEta ( int cc )

inline

Definition at line 142 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc, and CRYSTALS_IN_TOWER_ETA.

Referenced by p2eg::Cluster::crystaliEtaFromCardRegionInfo(), and p2eg::SimpleCaloHit::crystalLocaliEta().

                                              {
     if ((cc % 2) == 1) {  // if cc is odd (positive eta)
       return (17 * CRYSTALS_IN_TOWER_ETA);
     } else {  // if cc is even (negative eta) the bottom left corner is further in eta, hence +4
       return ((16 * CRYSTALS_IN_TOWER_ETA) + 4);
     }
   }

◆ getCard_refCrystal_iPhi()

int p2eg::getCard_refCrystal_iPhi ( int cc )

inline

Definition at line 151 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc, CRYSTALS_IN_TOWER_PHI, createfilelist::int, and TOWER_IN_PHI.

Referenced by p2eg::Cluster::crystaliPhiFromCardRegionInfo(), and p2eg::SimpleCaloHit::crystalLocaliPhi().

                                              {
     if ((cc % 2) == 1) {
       // if cc is odd: positive eta
       return int(cc / 2) * TOWER_IN_PHI * CRYSTALS_IN_TOWER_PHI;
     } else {
       // if cc is even, the bottom left corner is further in phi, hence the +4 and -1
       return (((int(cc / 2) * TOWER_IN_PHI) + 4) * CRYSTALS_IN_TOWER_PHI) - 1;
     }
   }

◆ getClusterFromRegion3x4()

Cluster p2eg::getClusterFromRegion3x4 ( crystal temp[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI] )

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ getClusterPosition()

p2eg::clusterInfo p2eg::getClusterPosition ( const ecalRegion_t ecalRegion )

inline

Definition at line 1057 of file Phase2L1RCT.h.

References p2eg::clusterInfo::brems, p2eg::crystalMax::energy, p2eg::clusterInfo::energy, p2eg::clusterInfo::et2x5, p2eg::clusterInfo::et5x5, p2eg::crystalMax::etaMax, p2eg::clusterInfo::etaMax, p2eg::ecalRegion_t::etaStrip0, p2eg::ecalRegion_t::etaStrip1, p2eg::ecalRegion_t::etaStrip10, p2eg::ecalRegion_t::etaStrip11, p2eg::ecalRegion_t::etaStrip12, p2eg::ecalRegion_t::etaStrip13, p2eg::ecalRegion_t::etaStrip14, p2eg::ecalRegion_t::etaStrip2, p2eg::ecalRegion_t::etaStrip3, p2eg::ecalRegion_t::etaStrip4, p2eg::ecalRegion_t::etaStrip5, p2eg::ecalRegion_t::etaStrip6, p2eg::ecalRegion_t::etaStrip7, p2eg::ecalRegion_t::etaStrip8, p2eg::ecalRegion_t::etaStrip9, getPeakBin15N(), getPeakBin20N(), p2eg::crystalMax::phiMax, p2eg::clusterInfo::phiMax, p2eg::etaStripPeak_t::pk0, p2eg::etaStripPeak_t::pk1, p2eg::etaStripPeak_t::pk10, p2eg::etaStripPeak_t::pk11, p2eg::etaStripPeak_t::pk12, p2eg::etaStripPeak_t::pk13, p2eg::etaStripPeak_t::pk14, p2eg::etaStripPeak_t::pk2, p2eg::etaStripPeak_t::pk3, p2eg::etaStripPeak_t::pk4, p2eg::etaStripPeak_t::pk5, p2eg::etaStripPeak_t::pk6, p2eg::etaStripPeak_t::pk7, p2eg::etaStripPeak_t::pk8, p2eg::etaStripPeak_t::pk9, and p2eg::clusterInfo::seedEnergy.

                                                                                  {
   p2eg::etaStripPeak_t etaStripPeak;
   p2eg::clusterInfo cluster;
 
   etaStripPeak.pk0 = p2eg::getPeakBin20N(ecalRegion.etaStrip0);
   etaStripPeak.pk1 = p2eg::getPeakBin20N(ecalRegion.etaStrip1);
   etaStripPeak.pk2 = p2eg::getPeakBin20N(ecalRegion.etaStrip2);
   etaStripPeak.pk3 = p2eg::getPeakBin20N(ecalRegion.etaStrip3);
   etaStripPeak.pk4 = p2eg::getPeakBin20N(ecalRegion.etaStrip4);
   etaStripPeak.pk5 = p2eg::getPeakBin20N(ecalRegion.etaStrip5);
   etaStripPeak.pk6 = p2eg::getPeakBin20N(ecalRegion.etaStrip6);
   etaStripPeak.pk7 = p2eg::getPeakBin20N(ecalRegion.etaStrip7);
   etaStripPeak.pk8 = p2eg::getPeakBin20N(ecalRegion.etaStrip8);
   etaStripPeak.pk9 = p2eg::getPeakBin20N(ecalRegion.etaStrip9);
   etaStripPeak.pk10 = p2eg::getPeakBin20N(ecalRegion.etaStrip10);
   etaStripPeak.pk11 = p2eg::getPeakBin20N(ecalRegion.etaStrip11);
   etaStripPeak.pk12 = p2eg::getPeakBin20N(ecalRegion.etaStrip12);
   etaStripPeak.pk13 = p2eg::getPeakBin20N(ecalRegion.etaStrip13);
   etaStripPeak.pk14 = p2eg::getPeakBin20N(ecalRegion.etaStrip14);
 
   p2eg::crystalMax peakIn15;
   peakIn15 = p2eg::getPeakBin15N(etaStripPeak);
 
   cluster.seedEnergy = peakIn15.energy;
   cluster.energy = 0;
   cluster.etaMax = peakIn15.etaMax;
   cluster.phiMax = peakIn15.phiMax;
   cluster.brems = 0;
   cluster.et5x5 = 0;
   cluster.et2x5 = 0;
 
   return cluster;
 }

◆ getClustersCombined()

p2eg::GCTcard_t p2eg::getClustersCombined	(	const GCTcard_t &	GCTcard,
		unsigned int	nGCTCard
	)

inline

Definition at line 91 of file Phase2L1GCT.h.

References p2eg::RCTcluster_t::crEta, p2eg::RCTcluster_t::crPhi, doProximityAndBremsStitching(), HLT_2024v10_cff::dPhi, p2eg::RCTcluster_t::et, HLT_2024v10_cff::eta1, HLT_2024v10_cff::eta2, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), dqmiolumiharvest::j, dqmdumpme::k, N_RCTCARDS_PHI, N_RCTCLUSTERS_FIBER, N_RCTGCT_FIBERS, p2eg::GCTcard_t::RCTcardEtaNeg, p2eg::GCTcard_t::RCTcardEtaPos, p2eg::RCTtoGCTfiber_t::RCTclusters, p2eg::RCTcard_t::RCTtoGCTfiber, p2eg::RCTcluster_t::towEta, p2eg::RCTcluster_t::towPhi, and SiPixelPI::two.

Referenced by getClustersTowers().

                                                                                                   {
   p2eg::GCTcard_t GCTcombinedClusters;
 
   // Initialize the output
   for (int i = 0; i < p2eg::N_RCTCARDS_PHI; i++) {
     for (int j = 0; j < p2eg::N_RCTGCT_FIBERS; j++) {
       for (int k = 0; k < p2eg::N_RCTCLUSTERS_FIBER; k++) {
         GCTcombinedClusters.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k] =
             GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k];
         GCTcombinedClusters.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTclusters[k] =
             GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTclusters[k];
       }
     }
   }
   bool isPositiveEta;
   int iStartingCard;
 
   // we will store new et in the GCTcombinedClusters, 0'ing lower clusters after stitching, dont need to care about other variables they stay the
   // same as input for now at least
   // we combine even phi boundaries positive eta. Start at card 0 (third argument), and tell the function this is positive eta (fourth argument)
   isPositiveEta = true;
   iStartingCard = 0;
   p2eg::doProximityAndBremsStitching(
       GCTcard.RCTcardEtaPos, GCTcombinedClusters.RCTcardEtaPos, iStartingCard, isPositiveEta);
 
   // now we combine odd phi boundaries positive eta
   isPositiveEta = true;
   iStartingCard = 1;
   p2eg::doProximityAndBremsStitching(
       GCTcard.RCTcardEtaPos, GCTcombinedClusters.RCTcardEtaPos, iStartingCard, isPositiveEta);
 
   // repeat above steps for NEGATIVE eta, even phi boundaries
   isPositiveEta = false;
   iStartingCard = 0;
   p2eg::doProximityAndBremsStitching(
       GCTcard.RCTcardEtaNeg, GCTcombinedClusters.RCTcardEtaNeg, iStartingCard, isPositiveEta);
 
   // lastly, NEGATIVE eta, odd phi boundaries
   isPositiveEta = false;
   iStartingCard = 1;
   p2eg::doProximityAndBremsStitching(
       GCTcard.RCTcardEtaNeg, GCTcombinedClusters.RCTcardEtaNeg, iStartingCard, isPositiveEta);
 
   // we need to store what we did before we start phi stitching
   p2eg::GCTcard_t GCTout;
   for (int i = 0; i < p2eg::N_RCTCARDS_PHI; i++) {
     for (int j = 0; j < p2eg::N_RCTGCT_FIBERS; j++) {
       for (int k = 0; k < p2eg::N_RCTCLUSTERS_FIBER; k++) {
         GCTout.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k] =
             GCTcombinedClusters.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k];
         GCTout.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTclusters[k] =
             GCTcombinedClusters.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTclusters[k];
       }
     }
   }
 
   // now we combine eta boundaries, just positive and negative eta
   // Uses RCTcardEtaPos and RCTcardEtaNeg
   for (int i = 0; i < p2eg::N_RCTCARDS_PHI; i++) {
     for (int j = 0; j < p2eg::N_RCTGCT_FIBERS; j++) {
       for (int k = 0; k < p2eg::N_RCTCLUSTERS_FIBER; k++) {
         ap_uint<15> phi1 = (i * 4 + GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].towPhi) * 5 +
                            GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].crPhi;
         ap_uint<15> eta1 = GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].crEta;
         if (GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].towEta == 0 && eta1 == 0) {
           for (int j1 = 0; j1 < p2eg::N_RCTGCT_FIBERS; j1++) {
             for (int k1 = 0; k1 < p2eg::N_RCTCLUSTERS_FIBER; k1++) {
               ap_uint<15> phi2 = (i * 4 + (3 - GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].towPhi)) * 5 +
                                  (4 - GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].crPhi);
               ap_uint<15> eta2 = GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].crEta;
               if (GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].towEta == 0 && eta2 == 0) {
                 ap_uint<15> dPhi;
                 dPhi = (phi1 > phi2) ? (phi1 - phi2) : (phi2 - phi1);
                 if (dPhi < 2) {
                   ap_uint<12> one = GCTcombinedClusters.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].et;
                   ap_uint<12> two = GCTcombinedClusters.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].et;
                   if (one > two) {
                     GCTout.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].et = one + two;
                     GCTout.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].et = 0;
                   } else {
                     GCTout.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k].et = 0;
                     GCTout.RCTcardEtaNeg[i].RCTtoGCTfiber[j1].RCTclusters[k1].et = one + two;
                   }
                 }
               }
             }
           }
         }
       }
     }
   }
   return GCTout;
 }

◆ getClustersTowers()

p2eg::GCTinternal_t p2eg::getClustersTowers	(	const GCTcard_t &	GCTcard,
		unsigned int	nGCTCard
	)

inline

Definition at line 189 of file Phase2L1GCT.h.

References p2eg::GCTCorrfiber_t::GCTclusters, p2eg::GCTinternal_t::GCTCorrfiber, p2eg::GCTCorrfiber_t::GCTtowers, getClustersCombined(), mps_fire::i, p2eg::GCTcluster_t::initFromRCTCluster(), p2eg::GCTtower_t::initFromRCTTower(), dqmiolumiharvest::j, dqmdumpme::k, N_GCTPOSITIVE_FIBERS, N_RCTCARDS_PHI, N_RCTCLUSTERS_FIBER, N_RCTGCT_FIBERS, N_RCTTOWERS_FIBER, p2eg::GCTcard_t::RCTcardEtaNeg, p2eg::GCTcard_t::RCTcardEtaPos, p2eg::RCTtoGCTfiber_t::RCTclusters, p2eg::RCTcard_t::RCTtoGCTfiber, and p2eg::RCTtoGCTfiber_t::RCTtowers.

Referenced by algo_top().

                                                                                                     {
   p2eg::GCTcard_t GCTcombinedClusters;
   p2eg::GCTinternal_t GCTout;
 
   // here we will stitch the clusters in phi and eta
   GCTcombinedClusters = p2eg::getClustersCombined(GCTcard, nGCTCard);
 
   // create internal structure of GCT card
   // we start from RCT card 0 - it is overlap with other GCT card and fill structure that we will use to send data to Correlator
   // we only need to care about clusters et in combinrdClusters, since the rest remains unchanged wrt input, the cluster that we set to 0
   // remain in the data at the same place , it will just get 0 et now
   // we need to code Positive and Negative Eta differently !  For negative Eta link 0 for each RCT
   // region becomes 3 in GCT output, the RCT card is rotated around 0:0 point of the card
   // First 16 fibers - positive Eta , second 16 - negative. Eta coded 0...16 and towEtaNeg = 0 or 1 for clusters ;
   // Phi is coded 0...15 , in case if whole card 0...33 and subdevision 1/5 in crPhi and crEta 0...4 for
   // position in tower
   //
   // towers are put in link starting from eta=0, the link number defines Eta negative or positive and Phi position of tower.
   for (int i = 0; i < p2eg::N_RCTCARDS_PHI; i++) {
     for (int j = 0; j < p2eg::N_RCTGCT_FIBERS; j++) {
       for (int k = 0; k < p2eg::N_RCTCLUSTERS_FIBER; k++) {
         bool isPositiveEta;
         // positive eta: initialize from RCT clusters in pos object
         isPositiveEta = true;
         GCTout.GCTCorrfiber[i * 4 + j].GCTclusters[k].initFromRCTCluster(
             i, isPositiveEta, GCTcombinedClusters.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTclusters[k]);
         // negative eta:  initialize from RCT clusters in neg object
         isPositiveEta = false;
         GCTout.GCTCorrfiber[i * 4 + (3 - j) + p2eg::N_GCTPOSITIVE_FIBERS].GCTclusters[k].initFromRCTCluster(
             i, isPositiveEta, GCTcombinedClusters.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTclusters[k]);
       }
       for (int k = 0; k < N_RCTTOWERS_FIBER; k++) {
         GCTout.GCTCorrfiber[i * 4 + j].GCTtowers[k].initFromRCTTower(
             GCTcard.RCTcardEtaPos[i].RCTtoGCTfiber[j].RCTtowers[k]);  // pos eta
         GCTout.GCTCorrfiber[i * 4 + (3 - j) + p2eg::N_GCTPOSITIVE_FIBERS].GCTtowers[k].initFromRCTTower(
             GCTcard.RCTcardEtaNeg[i].RCTtoGCTfiber[j].RCTtowers[k]);  // neg eta
       }
     }
   }
   return GCTout;
 }

◆ getClusterValues()

clusterInfo p2eg::getClusterValues	(	crystal	tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI],
		ap_uint< 5 >	seed_eta,
		ap_uint< 5 >	seed_phi
	)

◆ getECALTowersEt()

void p2eg::getECALTowersEt	(	crystal	tempX[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI],
		ap_uint< 12 >	towerEt[12]
	)

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ getFullTowers()

p2eg::GCTintTowers_t p2eg::getFullTowers ( const GCTinternal_t & GCTinternal )

inline

Definition at line 234 of file Phase2L1GCT.h.

References p2eg::GCTtower_t::ecalEt, p2eg::GCTcluster_t::et, p2eg::GCTtower_t::et, PVValHelper::eta, HLT_2024v10_cff::eta1, p2eg::GCTCorrfiber_t::GCTclusters, p2eg::GCTinternal_t::GCTCorrfiber, p2eg::GCTintTowers_t::GCTtower, p2eg::GCTCorrfiber_t::GCTtowers, p2eg::GCTtower_t::hcalEt, p2eg::GCTtower_t::hoe, mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, dqmdumpme::k, N_GCTCLUSTERS_FIBER, N_GCTETA, N_GCTINTERNAL_FIBERS, N_GCTPOSITIVE_FIBERS, N_GCTTOWERS_FIBER, PVValHelper::phi, p2eg::GCTcluster_t::towEta, and p2eg::GCTcluster_t::towPhi.

Referenced by algo_top().

                                                                                   {
   p2eg::GCTintTowers_t GCTintTowers;
   // Positive eta
   for (int i = 0; i < p2eg::N_GCTPOSITIVE_FIBERS; i = i + 4) {
     for (int i1 = 0; i1 < 4; i1++) {
       for (int k = 0; k < p2eg::N_GCTTOWERS_FIBER; k++) {
         ap_uint<15> phi = i + i1;
         ap_uint<15> eta = p2eg::N_GCTETA / 2 + k;
         GCTintTowers.GCTtower[eta][phi].et = GCTinternal.GCTCorrfiber[phi].GCTtowers[k].et;
         GCTintTowers.GCTtower[eta][phi].hoe = GCTinternal.GCTCorrfiber[phi].GCTtowers[k].hoe;
         GCTintTowers.GCTtower[eta][phi].ecalEt = GCTinternal.GCTCorrfiber[phi].GCTtowers[k].ecalEt;
         GCTintTowers.GCTtower[eta][phi].hcalEt = GCTinternal.GCTCorrfiber[phi].GCTtowers[k].hcalEt;
         for (int ic1 = 0; ic1 < 4; ic1++) {
           for (int jc = 0; jc < p2eg::N_GCTCLUSTERS_FIBER; jc++) {
             ap_uint<15> eta1 = p2eg::N_GCTETA / 2 + GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].towEta;
             ap_uint<15> phi1 = GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].towPhi;
             if (eta == eta1 && phi == phi1) {
               GCTintTowers.GCTtower[eta][phi].et =
                   (GCTintTowers.GCTtower[eta][phi].et + GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].et);
             }
           }
         }
       }
     }
   }
 
   // Negative eta
   for (int i = p2eg::N_GCTPOSITIVE_FIBERS; i < p2eg::N_GCTINTERNAL_FIBERS; i = i + 4) {
     for (int i1 = 0; i1 < 4; i1++) {
       for (int k = 0; k < p2eg::N_GCTTOWERS_FIBER; k++) {
         ap_uint<15> eta = p2eg::N_GCTETA / 2 - k - 1;
         ap_uint<15> phi = i + i1 - p2eg::N_GCTPOSITIVE_FIBERS;
         GCTintTowers.GCTtower[eta][phi].et = GCTinternal.GCTCorrfiber[i + i1].GCTtowers[k].et;
         GCTintTowers.GCTtower[eta][phi].hoe = GCTinternal.GCTCorrfiber[i + i1].GCTtowers[k].hoe;
         GCTintTowers.GCTtower[eta][phi].ecalEt = GCTinternal.GCTCorrfiber[i + i1].GCTtowers[k].ecalEt;
         GCTintTowers.GCTtower[eta][phi].hcalEt = GCTinternal.GCTCorrfiber[i + i1].GCTtowers[k].hcalEt;
         for (int ic1 = 0; ic1 < 4; ic1++) {
           for (int jc = 0; jc < p2eg::N_GCTCLUSTERS_FIBER; jc++) {
             ap_uint<15> eta1 = p2eg::N_GCTETA / 2 - 1 - GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].towEta;
             ap_uint<15> phi1 = GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].towPhi;
             if (eta == eta1 && phi == phi1) {
               GCTintTowers.GCTtower[eta][phi].et =
                   (GCTintTowers.GCTtower[eta][phi].et + GCTinternal.GCTCorrfiber[i + ic1].GCTclusters[jc].et);
             }
           }
         }
       }
     }
   }
 
   return GCTintTowers;
 }

◆ getPeakBin15N()

p2eg::crystalMax p2eg::getPeakBin15N ( const etaStripPeak_t etaStrip )

inline

Definition at line 982 of file Phase2L1RCT.h.

Referenced by getClusterPosition().

                                                                            {
   p2eg::crystalMax x;
 
   p2eg::ecaltp_t best01 = p2eg::bestOf2(etaStrip.pk0, etaStrip.pk1);
   p2eg::ecaltp_t best23 = p2eg::bestOf2(etaStrip.pk2, etaStrip.pk3);
   p2eg::ecaltp_t best45 = p2eg::bestOf2(etaStrip.pk4, etaStrip.pk5);
   p2eg::ecaltp_t best67 = p2eg::bestOf2(etaStrip.pk6, etaStrip.pk7);
   p2eg::ecaltp_t best89 = p2eg::bestOf2(etaStrip.pk8, etaStrip.pk9);
   p2eg::ecaltp_t best1011 = p2eg::bestOf2(etaStrip.pk10, etaStrip.pk11);
   p2eg::ecaltp_t best1213 = p2eg::bestOf2(etaStrip.pk12, etaStrip.pk13);
 
   p2eg::ecaltp_t best0123 = p2eg::bestOf2(best01, best23);
   p2eg::ecaltp_t best4567 = p2eg::bestOf2(best45, best67);
   p2eg::ecaltp_t best891011 = p2eg::bestOf2(best89, best1011);
   p2eg::ecaltp_t best121314 = p2eg::bestOf2(best1213, etaStrip.pk14);
 
   p2eg::ecaltp_t best01234567 = p2eg::bestOf2(best0123, best4567);
   p2eg::ecaltp_t best891011121314 = p2eg::bestOf2(best891011, best121314);
 
   p2eg::ecaltp_t bestOf15 = p2eg::bestOf2(best01234567, best891011121314);
 
   x.energy = bestOf15.energy;
   x.etaMax = bestOf15.eta;
   x.phiMax = bestOf15.phi;
 
   return x;
 }

◆ getPeakBin20N()

p2eg::ecaltp_t p2eg::getPeakBin20N ( const etaStrip_t etaStrip )

inline

Definition at line 950 of file Phase2L1RCT.h.

References bestOf2(), p2eg::etaStrip_t::cr0, p2eg::etaStrip_t::cr1, p2eg::etaStrip_t::cr10, p2eg::etaStrip_t::cr11, p2eg::etaStrip_t::cr12, p2eg::etaStrip_t::cr13, p2eg::etaStrip_t::cr14, p2eg::etaStrip_t::cr15, p2eg::etaStrip_t::cr16, p2eg::etaStrip_t::cr17, p2eg::etaStrip_t::cr18, p2eg::etaStrip_t::cr19, p2eg::etaStrip_t::cr2, p2eg::etaStrip_t::cr3, p2eg::etaStrip_t::cr4, p2eg::etaStrip_t::cr5, p2eg::etaStrip_t::cr6, p2eg::etaStrip_t::cr7, p2eg::etaStrip_t::cr8, and p2eg::etaStrip_t::cr9.

Referenced by getClusterPosition().

                                                                      {
   p2eg::ecaltp_t best01 = p2eg::bestOf2(etaStrip.cr0, etaStrip.cr1);
   p2eg::ecaltp_t best23 = p2eg::bestOf2(etaStrip.cr2, etaStrip.cr3);
   p2eg::ecaltp_t best45 = p2eg::bestOf2(etaStrip.cr4, etaStrip.cr5);
   p2eg::ecaltp_t best67 = p2eg::bestOf2(etaStrip.cr6, etaStrip.cr7);
   p2eg::ecaltp_t best89 = p2eg::bestOf2(etaStrip.cr8, etaStrip.cr9);
   p2eg::ecaltp_t best1011 = p2eg::bestOf2(etaStrip.cr10, etaStrip.cr11);
   p2eg::ecaltp_t best1213 = p2eg::bestOf2(etaStrip.cr12, etaStrip.cr13);
   p2eg::ecaltp_t best1415 = p2eg::bestOf2(etaStrip.cr14, etaStrip.cr15);
   p2eg::ecaltp_t best1617 = p2eg::bestOf2(etaStrip.cr16, etaStrip.cr17);
   p2eg::ecaltp_t best1819 = p2eg::bestOf2(etaStrip.cr18, etaStrip.cr19);
 
   p2eg::ecaltp_t best0123 = p2eg::bestOf2(best01, best23);
   p2eg::ecaltp_t best4567 = p2eg::bestOf2(best45, best67);
   p2eg::ecaltp_t best891011 = p2eg::bestOf2(best89, best1011);
   p2eg::ecaltp_t best12131415 = p2eg::bestOf2(best1213, best1415);
   p2eg::ecaltp_t best16171819 = p2eg::bestOf2(best1617, best1819);
 
   p2eg::ecaltp_t best01234567 = p2eg::bestOf2(best0123, best4567);
   p2eg::ecaltp_t best89101112131415 = p2eg::bestOf2(best891011, best12131415);
 
   p2eg::ecaltp_t best0to15 = p2eg::bestOf2(best01234567, best89101112131415);
   p2eg::ecaltp_t bestOf20 = p2eg::bestOf2(best0to15, best16171819);
 
   return bestOf20;
 }

◆ getRegionNumber()

int p2eg::getRegionNumber ( const int local_iEta )

inline

Definition at line 246 of file Phase2L1CaloEGammaUtils.h.

References cms::cuda::assert(), CRYSTALS_IN_TOWER_ETA, createfilelist::int, and TOWER_IN_ETA.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                    {
     int no = int(local_iEta / (TOWER_IN_ETA * CRYSTALS_IN_TOWER_ETA));
     assert(no < 6);
     return no;
   }

◆ getTower_absEtaID()

int p2eg::getTower_absEtaID ( float eta )

inline

Definition at line 167 of file Phase2L1CaloEGammaUtils.h.

References ECAL_eta_range, PVValHelper::eta, createfilelist::int, and n_towers_Eta.

                                           {
     float size_cell = 2 * ECAL_eta_range / n_towers_Eta;
     int etaID = int((eta + ECAL_eta_range) / size_cell);
     return etaID;
   }

◆ getTower_absPhiID()

int p2eg::getTower_absPhiID ( float phi )

inline

Definition at line 177 of file Phase2L1CaloEGammaUtils.h.

References createfilelist::int, M_PI, and n_towers_Phi.

                                           {
     float size_cell = 2 * M_PI / n_towers_Phi;
     int phiID = int((phi + M_PI) / size_cell);
     return phiID;
   }

◆ getTowerEta_fromAbsID()

float p2eg::getTowerEta_fromAbsID ( int id )

inline

Definition at line 226 of file Phase2L1CaloEGammaUtils.h.

References ECAL_eta_range, PVValHelper::eta, and n_towers_Eta.

Referenced by p2eg::GCTtower_t::createCaloTowerFromFiberIdx(), p2eg::GCTtower_t::createFullTowerFromCardIdx(), p2eg::GCTtower_t::printGCTTowerInfoFromGlobalIdx(), and Phase2L1CaloEGammaEmulator::produce().

                                              {
     float size_cell = 2 * ECAL_eta_range / n_towers_Eta;
     float eta = (id * size_cell) - ECAL_eta_range + 0.5 * size_cell;
     return eta;
   }

◆ getTowerPhi_fromAbsID()

float p2eg::getTowerPhi_fromAbsID ( int id )

inline

Definition at line 236 of file Phase2L1CaloEGammaUtils.h.

References M_PI, and n_towers_Phi.

Referenced by p2eg::GCTtower_t::createCaloTowerFromFiberIdx(), p2eg::GCTtower_t::createFullTowerFromCardIdx(), p2eg::GCTtower_t::printGCTTowerInfoFromGlobalIdx(), and Phase2L1CaloEGammaEmulator::produce().

                                              {
     float size_cell = 2 * M_PI / n_towers_Phi;
     float phi = (id * size_cell) - M_PI + 0.5 * size_cell;
     return phi;
   }

◆ initStructure()

ecalRegion_t p2eg::initStructure ( crystal temporary[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI] )

◆ isValidCard()

bool p2eg::isValidCard ( int cc )

inline

Definition at line 102 of file Phase2L1CaloEGammaUtils.h.

References gpuPixelDoublets::cc.

Referenced by getCard_iEtaMax().

102 { return ((cc > -1) && (cc < 36)); }

gpuPixelDoublets::cc

uint32_t cc[maxCellsPerHit]

Definition: gpuFishbone.h:49

◆ packCluster()

p2eg::Cluster p2eg::packCluster	(	ap_uint< 15 > &	clusterEt,
		ap_uint< 5 > &	etaMax_t,
		ap_uint< 5 > &	phiMax_t
	)

inline

Definition at line 1096 of file Phase2L1RCT.h.

References timeUnitHelper::pack().

                                                                                                        {
   ap_uint<12> peggedEt;
   p2eg::Cluster pack;
 
   ap_uint<5> towerEta = (etaMax_t) / 5;
   ap_uint<2> towerPhi = (phiMax_t) / 5;
   ap_uint<3> clusterEta = etaMax_t - 5 * towerEta;
   ap_uint<3> clusterPhi = phiMax_t - 5 * towerPhi;
 
   peggedEt = (clusterEt > 0xFFF) ? (ap_uint<12>)0xFFF : (ap_uint<12>)clusterEt;
 
   pack = p2eg::Cluster(peggedEt, towerEta, towerPhi, clusterEta, clusterPhi, 0);
 
   return pack;
 }

◆ passes_iso()

bool p2eg::passes_iso	(	float	pt,
		float	iso
	)

inline

Definition at line 904 of file Phase2L1CaloEGammaUtils.h.

References a0, a0_80, a1_80, high_pt_threshold, DiDispStaMuonMonitor_cfi::pt, and slideIsoPtThreshold.

Referenced by p2eg::GCTcluster_t::setRelIsoAndFlags().

                                               {
     bool is_iso = true;
     if (pt > high_pt_threshold)
       is_iso = true;
     else if (pt < slideIsoPtThreshold) {
       if (!((a0_80 - a1_80 * pt) > iso))
         is_iso = false;
     } else {
       if (iso > a0)
         is_iso = false;
     }
     return is_iso;
   }

◆ passes_looseTkiso()

bool p2eg::passes_looseTkiso	(	float	pt,
		float	iso
	)

inline

Definition at line 918 of file Phase2L1CaloEGammaUtils.h.

References b0, b1, b2, JetChargeProducer_cfi::exp, high_pt_threshold, and DiDispStaMuonMonitor_cfi::pt.

Referenced by p2eg::GCTcluster_t::setRelIsoAndFlags().

                                                      {
     bool is_iso;
     if (pt > high_pt_threshold)
       is_iso = true;
     else
       is_iso = (b0 + b1 * std::exp(-b2 * pt) > iso);
     return is_iso;
   }

◆ passes_looseTkss()

bool p2eg::passes_looseTkss	(	float	pt,
		float	eta,
		float	ss
	)

inline

Definition at line 939 of file Phase2L1CaloEGammaUtils.h.

References funct::abs(), e0_looseTkss, e1_looseTkss, e2_looseTkss, PVValHelper::eta, eta0_loose_ss, eta1_loose_ss, eta2_loose_ss, eta3_loose_ss, JetChargeProducer_cfi::exp, high_pt_threshold, loose_ss_offset, DiDispStaMuonMonitor_cfi::pt, and contentValuesCheck::ss.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                               {
     bool is_ss;
     if (pt > high_pt_threshold)
       is_ss = true;
     else if ((abs(eta) > eta0_loose_ss && abs(eta) < eta1_loose_ss) ||
              (abs(eta) > eta2_loose_ss && abs(eta) < eta3_loose_ss))  // temporary adjustment
       is_ss = ((loose_ss_offset - e1_looseTkss * std::exp(-e2_looseTkss * pt)) <= ss);
     else
       is_ss = ((e0_looseTkss - e1_looseTkss * std::exp(-e2_looseTkss * pt)) <= ss);
     return is_ss;
   }

◆ passes_ss()

bool p2eg::passes_ss	(	float	pt,
		float	eta,
		float	ss
	)

inline

Definition at line 927 of file Phase2L1CaloEGammaUtils.h.

References funct::abs(), c0_ss, c1_ss, c2_ss, PVValHelper::eta, eta0_loose_ss, eta1_loose_ss, eta2_loose_ss, eta3_loose_ss, JetChargeProducer_cfi::exp, high_pt_threshold, loose_ss_offset, DiDispStaMuonMonitor_cfi::pt, and contentValuesCheck::ss.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                        {
     bool is_ss;
     if (pt > high_pt_threshold)
       is_ss = true;
     else if ((abs(eta) > eta0_loose_ss && abs(eta) < eta1_loose_ss) ||
              (abs(eta) > eta2_loose_ss && abs(eta) < eta3_loose_ss))  // temporary adjustment
       is_ss = ((loose_ss_offset + c1_ss * std::exp(-c2_ss * pt)) <= ss);
     else
       is_ss = ((c0_ss + c1_ss * std::exp(-c2_ss * pt)) <= ss);
     return is_ss;
   }

◆ printl1tp2TowerInfo()

void p2eg::printl1tp2TowerInfo	(	l1tp2::CaloTower	thisTower,
		std::string	description = `""`
	)

inline

Definition at line 1616 of file Phase2L1CaloEGammaUtils.h.

References gather_cfg::cout, makeListRunsInFiles::description, l1tp2::CaloTower::ecalTowerEt(), l1tp2::CaloTower::hcalTowerEt(), l1tp2::CaloTower::towerEta(), l1tp2::CaloTower::towerIEta(), l1tp2::CaloTower::towerIPhi(), and l1tp2::CaloTower::towerPhi().

                                                                                         {
     std::cout << "[Print l1tp2::CaloTower info:] [" << description << "]: "
               << ".ecalTowerEta() (float): " << thisTower.ecalTowerEt() << ", "
               << ".hcalTowerEta() (float): " << thisTower.hcalTowerEt() << ", "
               << ".towerIEta(): " << thisTower.towerIEta() << ", "
               << ".towerIPhi(): " << thisTower.towerIPhi() << ", "
               << ".towerEta() " << thisTower.towerEta() << ", "
               << ".towerPhi() " << thisTower.towerPhi() << std::endl;
   }

◆ removeClusterFromCrystal()

void p2eg::removeClusterFromCrystal	(	crystal	temp[CRYSTAL_IN_ETA][CRYSTAL_IN_PHI],
		ap_uint< 5 >	seed_eta,
		ap_uint< 5 >	seed_phi,
		ap_uint< 2 >	brems
	)

◆ stitchClusterOverRegionBoundary()

void p2eg::stitchClusterOverRegionBoundary	(	std::vector< Cluster > &	cluster_list,
		int	towerEtaUpper,
		int	towerEtaLower,
		int	cc
	)

inline

Definition at line 1475 of file Phase2L1RCT.h.

References alignmentValidation::c1, gpuPixelDoublets::cc, p2eg::Cluster::clusterEnergy(), p2eg::Cluster::clusterEta(), p2eg::Cluster::clusterPhi(), HLT_2024v10_cff::dPhi, p2eg::Cluster::getBrems(), p2eg::Cluster::getCalib(), p2eg::Cluster::getIsLooseTkss(), p2eg::Cluster::getIsSS(), mps_fire::i, dqmiolumiharvest::j, p2eg::Cluster::region(), p2eg::Cluster::satur(), p2eg::Cluster::towerEta(), p2eg::Cluster::towerEtaInCard(), p2eg::Cluster::towerPhi(), p2eg::Cluster::uint_et2x5(), p2eg::Cluster::uint_et5x5(), and funct::void.

Referenced by Phase2L1CaloEGammaEmulator::produce().

                                                           {
   (void)cc;  // for printout statements
 
   int crystalEtaUpper = 0;
   int crystalEtaLower = 4;
 
   for (size_t i = 0; i < cluster_list.size(); i++) {
     for (size_t j = 0; j < cluster_list.size(); j++) {
       // Do not double-count
       if (i == j)
         continue;
 
       p2eg::Cluster c1 = cluster_list[i];
       p2eg::Cluster c2 = cluster_list[j];
 
       p2eg::Cluster newc1;
       p2eg::Cluster newc2;
 
       // Use the .towerEtaInCard() method to get the tower eta in the entire RCT card
       if ((c1.clusterEnergy() > 0) && (c1.towerEtaInCard() == towerEtaUpper) && (c1.clusterEta() == crystalEtaUpper)) {
         if ((c2.clusterEnergy() > 0) && (c2.towerEtaInCard() == towerEtaLower) &&
             (c2.clusterEta() == crystalEtaLower)) {
           ap_uint<5> phi1 = c1.towerPhi() * 5 + c1.clusterPhi();
           ap_uint<5> phi2 = c2.towerPhi() * 5 + c2.clusterPhi();
           ap_uint<5> dPhi;
           dPhi = (phi1 > phi2) ? (phi1 - phi2) : (phi2 - phi1);
 
           if (dPhi < 2) {
             ap_uint<12> totalEnergy = c1.clusterEnergy() + c2.clusterEnergy();
             ap_uint<15> totalEt2x5 = c1.uint_et2x5() + c2.uint_et2x5();
             ap_uint<15> totalEt5x5 = c1.uint_et5x5() + c2.uint_et5x5();
 
             bool rct_is_iso = false;         // RCT has no isolation information
             bool rct_is_looseTkiso = false;  // RCT has no isolation information
 
             // Initialize a cluster with the larger cluster's position and total energy
             if (c1.clusterEnergy() > c2.clusterEnergy()) {
               newc1 = p2eg::Cluster(totalEnergy,
                                     c1.towerEta(),
                                     c1.towerPhi(),
                                     c1.clusterEta(),
                                     c1.clusterPhi(),
                                     c1.satur(),
                                     totalEt5x5,
                                     totalEt2x5,
                                     c1.getBrems(),
                                     c1.getCalib(),
                                     c1.getIsSS(),
                                     c1.getIsLooseTkss(),
                                     rct_is_iso,
                                     rct_is_looseTkiso,
                                     c1.region());
               newc2 = p2eg::Cluster(0,
                                     c2.towerEta(),
                                     c2.towerPhi(),
                                     c2.clusterEta(),
                                     c2.clusterPhi(),
                                     c2.satur(),
                                     0,
                                     0,
                                     0,
                                     1.0,
                                     false,
                                     false,
                                     rct_is_iso,
                                     rct_is_looseTkiso,
                                     c2.region());
               cluster_list[i] = newc1;
               cluster_list[j] = newc2;
             } else {
               // Analogous to above portion
               newc1 = p2eg::Cluster(0,
                                     c1.towerEta(),
                                     c1.towerPhi(),
                                     c1.clusterEta(),
                                     c1.clusterPhi(),
                                     c1.satur(),
                                     0,
                                     0,
                                     0,
                                     1.0,
                                     false,
                                     false,
                                     rct_is_iso,
                                     rct_is_looseTkiso,
                                     c1.region());
               newc2 = p2eg::Cluster(totalEnergy,
                                     c2.towerEta(),
                                     c2.towerPhi(),
                                     c2.clusterEta(),
                                     c2.clusterPhi(),
                                     c2.satur(),
                                     totalEt5x5,
                                     totalEt2x5,
                                     c2.getBrems(),
                                     c2.getCalib(),
                                     c2.getIsSS(),
                                     c2.getIsLooseTkss(),
                                     rct_is_iso,
                                     rct_is_looseTkiso,
                                     c2.region());
               cluster_list[i] = newc1;
               cluster_list[j] = newc2;
             }
           }
         }
       }
     }
   }
 }

◆ writeToCorrelatorAndGTOutputs()

void p2eg::writeToCorrelatorAndGTOutputs	(	const GCTinternal_t &	GCTinternal,
		p2eg::GCTtoCorr_t &	GCTtoCorrOutput,
		std::unique_ptr< l1tp2::CaloCrystalClusterCollection > const &	gctClustersOutput,
		std::unique_ptr< l1tp2::CaloTowerCollection > const &	gctTowersOutput,
		std::unique_ptr< l1t::EGammaBxCollection > const &	gctEGammas,
		std::unique_ptr< l1tp2::DigitizedClusterCorrelatorCollection > const &	gctDigitizedClustersCorrelator,
		std::unique_ptr< l1tp2::DigitizedTowerCorrelatorCollection > const &	gctDigitizedTowersCorrelator,
		std::unique_ptr< l1tp2::DigitizedClusterGTCollection > const &	gctDigitizedClustersGT,
		int	nGCTCard,
		int	fiberStart,
		int	fiberEnd,
		int	corrFiberIndexOffset,
		int	corrTowPhiOffset = `4`
	)

inline

Definition at line 290 of file Phase2L1GCT.h.

References nano_mu_digi_cff::bx, p2eg::GCTcluster_t::createCaloCrystalCluster(), p2eg::GCTtower_t::createCaloTowerFromFiberIdx(), p2eg::GCTcluster_t::createDigitizedClusterCorrelator(), p2eg::GCTcluster_t::createDigitizedClusterGT(), p2eg::GCTtower_t::createDigitizedTowerCorrelator(), p2eg::GCTcluster_t::createL1TEGamma(), p2eg::GCTcluster_t::etFloat(), p2eg::GCTCorrfiber_t::GCTclusters, p2eg::GCTtoCorr_t::GCTCorrfiber, p2eg::GCTinternal_t::GCTCorrfiber, p2eg::GCTCorrfiber_t::GCTtowers, mps_fire::i, sistrip::SpyUtilities::isValid(), dqmdumpme::k, N_GCTCLUSTERS_FIBER, N_GCTTOWERS_FIBER, and p2eg::GCTcluster_t::towPhi.

Referenced by algo_top().

                               {
   for (int i = fiberStart; i < fiberEnd; i++) {
     // In each fiber, first do clusters
     for (int k = 0; k < p2eg::N_GCTCLUSTERS_FIBER; k++) {
       // First do CMSSW cluster outputs
       p2eg::GCTcluster_t thisCluster = GCTinternal.GCTCorrfiber[i].GCTclusters[k];
       if (thisCluster.etFloat() > 0.0) {
         // Make l1tp2::CaloCrystalCluster
         gctClustersOutput->push_back(thisCluster.createCaloCrystalCluster());
 
         // Make l1t::EGamma
         int bx = 0;
         l1t::EGamma thisEGamma = thisCluster.createL1TEGamma();
         gctEGammas->push_back(bx, thisEGamma);
       }
 
       // Then the clusters to the correlator: all fields are the same with the exception of towPhi, which
       // needs to be subtracted by 4 because the output to correlator does NOT include the overlap region.
       GCTtoCorrOutput.GCTCorrfiber[i - corrFiberIndexOffset].GCTclusters[k] = thisCluster;
       GCTtoCorrOutput.GCTCorrfiber[i - corrFiberIndexOffset].GCTclusters[k].towPhi =
           (thisCluster.towPhi - corrTowPhiOffset);
 
       // Make l1tp2::DigitizedClusterCorrelator. The function needs corrTowPhiOffset to know the towPhi in the card excluding the overlap region.
       // The correlator clusters don't need to know the fiber offset.
       if (thisCluster.etFloat() > 0.0) {
         gctDigitizedClustersCorrelator->push_back(thisCluster.createDigitizedClusterCorrelator(corrTowPhiOffset));
       }
 
       // Make l1tp2::DigitizedClusterGT.
       if (thisCluster.etFloat() > 0.0) {
         bool isValid = true;
         gctDigitizedClustersGT->push_back(thisCluster.createDigitizedClusterGT(isValid));
       }
     }
 
     // Next do tower outputs
     for (int k = 0; k < p2eg::N_GCTTOWERS_FIBER; k++) {
       // First do CMSSW tower outputs
       p2eg::GCTtower_t thisTower = GCTinternal.GCTCorrfiber[i].GCTtowers[k];
       l1tp2::CaloTower thisL1CaloTower = thisTower.createCaloTowerFromFiberIdx(nGCTCard, i, k);
       gctTowersOutput->push_back(thisL1CaloTower);
 
       // Then the towers to the correlator. Note the same corrFiberIndexOffset as was done for the clusters
       GCTtoCorrOutput.GCTCorrfiber[i - corrFiberIndexOffset].GCTtowers[k] = thisTower;
 
       // For the collection, the three arguments are (1) the GCT card, (2) the fiber index in the GCT card (excluding the overlap region), and (3) the tower index in the fiber
       l1tp2::DigitizedTowerCorrelator thisDigitizedTowerCorrelator =
           thisTower.createDigitizedTowerCorrelator(nGCTCard, i - corrFiberIndexOffset, k);
       gctDigitizedTowersCorrelator->push_back(thisDigitizedTowerCorrelator);
     }
   }
 }

Variable Documentation

◆ a0

constexpr float p2eg::a0 = 0.21

static

Definition at line 57 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_iso().

◆ a0_80

constexpr float p2eg::a0_80 = 0.85

static

Definition at line 57 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_iso().

◆ a1_80

constexpr float p2eg::a1_80 = 0.0080

static

Definition at line 57 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_iso().

◆ b0

constexpr float p2eg::b0 = 0.38

static

Definition at line 58 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkiso().

◆ b1

constexpr float p2eg::b1 = 1.9

static

Definition at line 58 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkiso().

◆ b2

constexpr float p2eg::b2 = 0.05

static

Definition at line 58 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkiso().

◆ c0_ss

constexpr float p2eg::c0_ss = 0.94

static

Definition at line 59 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_ss().

◆ c1_ss

constexpr float p2eg::c1_ss = 0.052

static

Definition at line 59 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_ss().

◆ c2_ss

constexpr float p2eg::c2_ss = 0.044

static

Definition at line 59 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_ss().

◆ CRYSTAL_IN_ETA

constexpr int p2eg::CRYSTAL_IN_ETA = 15

static

Definition at line 50 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ CRYSTAL_IN_PHI

constexpr int p2eg::CRYSTAL_IN_PHI = 20

static

Definition at line 51 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ CRYSTALS_IN_TOWER_ETA

constexpr int p2eg::CRYSTALS_IN_TOWER_ETA = 5

static

Definition at line 44 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::linkECAL::addCrystalE(), p2eg::SimpleCaloHit::crystaliEta(), p2eg::Cluster::crystaliEtaFromCardRegionInfo(), getCard_iEtaMax(), getCard_iEtaMin(), getCard_refCrystal_iEta(), getRegionNumber(), p2eg::GCTcluster_t::globalClusteriEta(), Phase2L1CaloEGammaEmulator::produce(), p2eg::Cluster::realEta(), p2eg::GCTcluster_t::realEta(), p2eg::linkECAL::setCrystalE(), and p2eg::linkECAL::zeroOut().

◆ CRYSTALS_IN_TOWER_PHI

constexpr int p2eg::CRYSTALS_IN_TOWER_PHI = 5

static

Definition at line 45 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::linkECAL::addCrystalE(), p2eg::GCTcluster_t::createDigitizedClusterCorrelator(), p2eg::SimpleCaloHit::crystaliPhi(), p2eg::Cluster::crystaliPhiFromCardRegionInfo(), getCard_iPhiMax(), getCard_iPhiMin(), getCard_refCrystal_iPhi(), p2eg::GCTcluster_t::globalClusteriPhi(), Phase2L1CaloEGammaEmulator::produce(), p2eg::Cluster::realPhi(), p2eg::GCTcluster_t::realPhi(), p2eg::linkECAL::setCrystalE(), and p2eg::linkECAL::zeroOut().

◆ cut_500_MeV

constexpr float p2eg::cut_500_MeV = 0.5

static

Definition at line 65 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ d0

constexpr float p2eg::d0 = 0.96

static

Definition at line 60 of file Phase2L1CaloEGammaUtils.h.

◆ d1

constexpr float p2eg::d1 = 0.0003

static

Definition at line 60 of file Phase2L1CaloEGammaUtils.h.

◆ e0_looseTkss

constexpr float p2eg::e0_looseTkss = 0.944

static

Definition at line 61 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss().

◆ e1_looseTkss

constexpr float p2eg::e1_looseTkss = 0.65

static

Definition at line 61 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss().

◆ e2_looseTkss

constexpr float p2eg::e2_looseTkss = 0.4

static

Definition at line 61 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss().

◆ ECAL_eta_range

constexpr float p2eg::ECAL_eta_range = 1.4841

static

Definition at line 53 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::SimpleCaloHit::crystaliEta(), getTower_absEtaID(), getTowerEta_fromAbsID(), p2eg::Cluster::realEta(), and p2eg::GCTcluster_t::realEta().

◆ ECAL_LSB

constexpr float p2eg::ECAL_LSB = 0.5

static

Definition at line 67 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::tower_t::applyCalibration(), p2eg::Cluster::applyCalibration(), p2eg::GCTtower_t::ecalEtFloat(), p2eg::GCTcluster_t::et2x5Float(), p2eg::GCTcluster_t::et5x5Float(), p2eg::GCTcluster_t::etFloat(), p2eg::tower_t::getEt(), p2eg::Cluster::getEt2x5(), p2eg::Cluster::getEt5x5(), p2eg::Cluster::getPt(), p2eg::GCTcluster_t::isoFloat(), Phase2L1CaloEGammaEmulator::produce(), and p2eg::GCTtower_t::totalEtFloat().

◆ eta0_loose_ss

constexpr float p2eg::eta0_loose_ss = 0.4

static

Definition at line 63 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss(), and passes_ss().

◆ eta1_loose_ss

constexpr float p2eg::eta1_loose_ss = 0.5

static

Definition at line 63 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss(), and passes_ss().

◆ eta2_loose_ss

constexpr float p2eg::eta2_loose_ss = 0.75

static

Definition at line 63 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss(), and passes_ss().

◆ eta3_loose_ss

constexpr float p2eg::eta3_loose_ss = 0.85

static

Definition at line 63 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss(), and passes_ss().

◆ GCTCARD_0_TOWER_IPHI_OFFSET

constexpr int p2eg::GCTCARD_0_TOWER_IPHI_OFFSET = 20

static

Definition at line 90 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTcluster_t::globalClusteriPhi(), p2eg::GCTtower_t::globalToweriPhi(), and p2eg::GCTtower_t::globalToweriPhiFromGCTcardiPhi().

◆ GCTCARD_1_TOWER_IPHI_OFFSET

constexpr int p2eg::GCTCARD_1_TOWER_IPHI_OFFSET = 44

static

Definition at line 92 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTcluster_t::globalClusteriPhi(), p2eg::GCTtower_t::globalToweriPhi(), and p2eg::GCTtower_t::globalToweriPhiFromGCTcardiPhi().

◆ GCTCARD_2_TOWER_IPHI_OFFSET

constexpr int p2eg::GCTCARD_2_TOWER_IPHI_OFFSET = 68

static

Definition at line 93 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTcluster_t::globalClusteriPhi(), p2eg::GCTtower_t::globalToweriPhi(), and p2eg::GCTtower_t::globalToweriPhiFromGCTcardiPhi().

◆ GCTcardtoRCTcardnumber

const unsigned int p2eg::GCTcardtoRCTcardnumber[N_GCTCARDS][N_RCTCARDS_PHI *2]

static

Initial value:

= {
      
      {11, 13, 15, 17, 19, 21, 23, 25, 10, 12, 14, 16, 18, 20, 22, 24},
      
      {23, 25, 27, 29, 31, 33, 35, 1, 22, 24, 26, 28, 30, 32, 34, 0},
      
      {35, 1, 3, 5, 7, 9, 11, 13, 34, 0, 2, 4, 6, 8, 10, 12}}

Definition at line 1603 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ half_crystal_size

constexpr float p2eg::half_crystal_size = 0.00873

static

Definition at line 54 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::Cluster::realEta(), p2eg::GCTcluster_t::realEta(), p2eg::Cluster::realPhi(), and p2eg::GCTcluster_t::realPhi().

◆ HCAL_LSB

constexpr float p2eg::HCAL_LSB = 0.5

static

Definition at line 68 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTtower_t::hcalEtFloat(), and Phase2L1CaloEGammaEmulator::produce().

◆ high_pt_threshold

constexpr float p2eg::high_pt_threshold = 130.

static

Definition at line 64 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_iso(), passes_looseTkiso(), passes_looseTkss(), and passes_ss().

◆ loose_ss_offset

constexpr float p2eg::loose_ss_offset = 0.89

static

Definition at line 62 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_looseTkss(), and passes_ss().

◆ n_clusters_4link

constexpr int p2eg::n_clusters_4link = 8

static

Definition at line 41 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ n_clusters_link

constexpr int p2eg::n_clusters_link = 2

static

Definition at line 40 of file Phase2L1CaloEGammaUtils.h.

◆ N_CLUSTERS_PER_REGION

constexpr int p2eg::N_CLUSTERS_PER_REGION = 4

static

Definition at line 70 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ n_crystals_cardEta

constexpr int p2eg::n_crystals_cardEta = (n_towers_Eta * n_towers_cardEta)

static

Definition at line 35 of file Phase2L1CaloEGammaUtils.h.

◆ n_crystals_cardPhi

constexpr int p2eg::n_crystals_cardPhi = (n_towers_Phi * n_towers_cardPhi)

static

Definition at line 36 of file Phase2L1CaloEGammaUtils.h.

◆ N_GCTCARDS

constexpr int p2eg::N_GCTCARDS = 3

static

Definition at line 79 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ N_GCTCLUSTERS_FIBER

constexpr int p2eg::N_GCTCLUSTERS_FIBER = 2

static

Definition at line 82 of file Phase2L1CaloEGammaUtils.h.

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

◆ N_GCTCORR_FIBERS

constexpr int p2eg::N_GCTCORR_FIBERS = 48

static

Definition at line 80 of file Phase2L1CaloEGammaUtils.h.

◆ N_GCTETA

constexpr int p2eg::N_GCTETA = 34

static

Definition at line 86 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTinternal_t::computeClusterIsolationInPlace(), getFullTowers(), and p2eg::GCTintTowers_t::writeToPFOutput().

◆ N_GCTINTERNAL_FIBERS

constexpr int p2eg::N_GCTINTERNAL_FIBERS = 64

static

Definition at line 84 of file Phase2L1CaloEGammaUtils.h.

Referenced by algo_top(), p2eg::GCTinternal_t::computeClusterIsolationInPlace(), getFullTowers(), and p2eg::GCTinternal_t::setIsolationInfo().

◆ N_GCTPHI

constexpr int p2eg::N_GCTPHI = 32

static

Definition at line 87 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTinternal_t::computeClusterIsolationInPlace(), and p2eg::GCTintTowers_t::writeToPFOutput().

◆ N_GCTPOSITIVE_FIBERS

constexpr int p2eg::N_GCTPOSITIVE_FIBERS = 32

static

Definition at line 85 of file Phase2L1CaloEGammaUtils.h.

Referenced by algo_top(), p2eg::GCTinternal_t::computeClusterIsolationInPlace(), getClustersTowers(), getFullTowers(), p2eg::GCTtower_t::globalToweriEta(), and p2eg::GCTtower_t::globalToweriPhi().

◆ N_GCTTOWERS_CLUSTER_ISO_ONESIDE

constexpr int p2eg::N_GCTTOWERS_CLUSTER_ISO_ONESIDE = 5

static

Definition at line 95 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTinternal_t::computeClusterIsolationInPlace().

◆ N_GCTTOWERS_FIBER

constexpr int p2eg::N_GCTTOWERS_FIBER = 17

static

Definition at line 81 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::GCTinternal_t::computeClusterIsolationInPlace(), getFullTowers(), p2eg::GCTcluster_t::globalClusteriEta(), p2eg::GCTtower_t::globalToweriEta(), Phase2L1CaloEGammaEmulator::produce(), and writeToCorrelatorAndGTOutputs().

◆ n_links_card

constexpr int p2eg::n_links_card = 4

static

Definition at line 39 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ N_RCTCARDS_PHI

constexpr int p2eg::N_RCTCARDS_PHI = 8

static

Definition at line 74 of file Phase2L1CaloEGammaUtils.h.

Referenced by getClustersCombined(), getClustersTowers(), and Phase2L1CaloEGammaEmulator::produce().

◆ N_RCTCLUSTERS_FIBER

constexpr int p2eg::N_RCTCLUSTERS_FIBER = 2

static

Definition at line 77 of file Phase2L1CaloEGammaUtils.h.

Referenced by getClustersCombined(), getClustersTowers(), and Phase2L1CaloEGammaEmulator::produce().

◆ N_RCTGCT_FIBERS

constexpr int p2eg::N_RCTGCT_FIBERS = 4

static

Definition at line 75 of file Phase2L1CaloEGammaUtils.h.

Referenced by algo_top(), getClustersCombined(), getClustersTowers(), and Phase2L1CaloEGammaEmulator::produce().

◆ N_RCTTOWERS_FIBER

constexpr int p2eg::N_RCTTOWERS_FIBER = 17

static

Definition at line 76 of file Phase2L1CaloEGammaUtils.h.

Referenced by getClustersTowers().

◆ N_REGIONS_PER_CARD

constexpr int p2eg::N_REGIONS_PER_CARD = 6

static

Definition at line 71 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::card::card(), p2eg::card::getRegion3x4(), p2eg::card::getTowers3x4(), Phase2L1CaloEGammaEmulator::produce(), and p2eg::card::zeroOut().

◆ n_towers_cardEta

constexpr int p2eg::n_towers_cardEta = 17

static

Definition at line 33 of file Phase2L1CaloEGammaUtils.h.

Referenced by getCard_iEtaMax(), getCard_iEtaMin(), and Phase2L1CaloEGammaEmulator::produce().

◆ n_towers_cardPhi

constexpr int p2eg::n_towers_cardPhi = 4

static

Definition at line 34 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ n_towers_Eta

constexpr int p2eg::n_towers_Eta = 34

static

Definition at line 30 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::SimpleCaloHit::crystaliEta(), getCard_iEtaMax(), getTower_absEtaID(), getTowerEta_fromAbsID(), p2eg::Cluster::realEta(), and p2eg::GCTcluster_t::realEta().

◆ n_towers_halfPhi

constexpr int p2eg::n_towers_halfPhi = 36

static

Definition at line 32 of file Phase2L1CaloEGammaUtils.h.

Referenced by Phase2L1CaloEGammaEmulator::produce().

◆ n_towers_per_link

constexpr int p2eg::n_towers_per_link = 17

static

Definition at line 42 of file Phase2L1CaloEGammaUtils.h.

Referenced by getAbsID_iEta_fromFirmwareCardTowerLink().

◆ n_towers_Phi

constexpr int p2eg::n_towers_Phi = 72

static

Definition at line 31 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::SimpleCaloHit::crystaliPhi(), getTower_absPhiID(), getTowerPhi_fromAbsID(), p2eg::GCTcluster_t::globalClusteriPhi(), p2eg::GCTtower_t::globalToweriPhi(), p2eg::GCTtower_t::globalToweriPhiFromGCTcardiPhi(), p2eg::Cluster::realPhi(), and p2eg::GCTcluster_t::realPhi().

◆ slideIsoPtThreshold

constexpr float p2eg::slideIsoPtThreshold = 80

static

Definition at line 56 of file Phase2L1CaloEGammaUtils.h.

Referenced by passes_iso().

◆ TOWER_IN_ETA

constexpr int p2eg::TOWER_IN_ETA = 3

static

Definition at line 47 of file Phase2L1CaloEGammaUtils.h.

Referenced by p2eg::Cluster::crystaliEtaFromCardRegionInfo(), getRegionNumber(), Phase2L1CaloEGammaEmulator::produce(), p2eg::Cluster::towerEtaInCard(), p2eg::region3x4::zeroOut(), and p2eg::towers3x4::zeroOut().

◆ TOWER_IN_PHI

constexpr int p2eg::TOWER_IN_PHI = 4

static

Definition at line 48 of file Phase2L1CaloEGammaUtils.h.

Referenced by getAbsID_iPhi_fromFirmwareCardTowerLink(), getCard_refCrystal_iPhi(), Phase2L1CaloEGammaEmulator::produce(), p2eg::region3x4::zeroOut(), and p2eg::towers3x4::zeroOut().

Classes

Functions

Variables

Function Documentation

◆ algo_top()

◆ bestOf2()

◆ compareClusterET()

◆ doProximityAndBremsStitching()

◆ getAbsID_iEta_fromFirmwareCardTowerLink()

◆ getAbsID_iPhi_fromFirmwareCardTowerLink()

◆ getBremsValuesNeg()

◆ getBremsValuesPos()

◆ getCard_iEtaMax()

◆ getCard_iEtaMin()

◆ getCard_iPhiMax()

◆ getCard_iPhiMin()

◆ getCard_refCrystal_iEta()

◆ getCard_refCrystal_iPhi()

◆ getClusterFromRegion3x4()

◆ getClusterPosition()

◆ getClustersCombined()

◆ getClustersTowers()

◆ getClusterValues()

◆ getECALTowersEt()

◆ getFullTowers()

◆ getPeakBin15N()

◆ getPeakBin20N()

◆ getRegionNumber()

◆ getTower_absEtaID()

◆ getTower_absPhiID()

◆ getTowerEta_fromAbsID()

◆ getTowerPhi_fromAbsID()

◆ initStructure()

◆ isValidCard()

◆ packCluster()

◆ passes_iso()

◆ passes_looseTkiso()

◆ passes_looseTkss()

◆ passes_ss()

◆ printl1tp2TowerInfo()

◆ removeClusterFromCrystal()

◆ stitchClusterOverRegionBoundary()

◆ writeToCorrelatorAndGTOutputs()

Variable Documentation

◆ a0

◆ a0_80

◆ a1_80

◆ b0

◆ b1

◆ b2

◆ c0_ss

◆ c1_ss

◆ c2_ss

◆ CRYSTAL_IN_ETA

◆ CRYSTAL_IN_PHI

◆ CRYSTALS_IN_TOWER_ETA

◆ CRYSTALS_IN_TOWER_PHI

◆ cut_500_MeV

◆ d0

◆ d1

◆ e0_looseTkss

◆ e1_looseTkss

◆ e2_looseTkss

◆ ECAL_eta_range

◆ ECAL_LSB

◆ eta0_loose_ss

◆ eta1_loose_ss

◆ eta2_loose_ss

◆ eta3_loose_ss

◆ GCTCARD_0_TOWER_IPHI_OFFSET

◆ GCTCARD_1_TOWER_IPHI_OFFSET

◆ GCTCARD_2_TOWER_IPHI_OFFSET

◆ GCTcardtoRCTcardnumber

◆ half_crystal_size

◆ HCAL_LSB

◆ high_pt_threshold

◆ loose_ss_offset

◆ n_clusters_4link

◆ n_clusters_link

◆ N_CLUSTERS_PER_REGION