p2eg Namespace Reference

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_2024v13_cff::dPhi, p2eg::RCTcluster_t::et, HLT_2024v13_cff::eta1, HLT_2024v13_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_2024v13_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.

References bestOf2(), p2eg::ecaltp_t::energy, p2eg::ecaltp_t::eta, p2eg::ecaltp_t::phi, 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 x.

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().

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

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, genparticles_cff::iso, 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, genparticles_cff::iso, 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 1613 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_2024v13_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::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 1600 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().