l1t::CaloTools Class Reference

#include <CaloTools.h>

Public Types
enum	SubDet { ECAL = 0x1, HCAL = 0x2, CALO = 0x3 }

Static Public Member Functions
static int	bin16Eta (int ieta)
static int	calHwEtSum (int iEta, int iPhi, const std::vector< l1t::CaloTower > &towers, int localEtaMin, int localEtaMax, int localPhiMin, int localPhiMax, SubDet etMode=CALO)
static int	calHwEtSum (int iEta, int iPhi, const std::vector< l1t::CaloTower > &towers, int localEtaMin, int localEtaMax, int localPhiMin, int localPhiMax, int iEtaAbsMax, SubDet etMode=CALO)
static size_t	calNrTowers (int iEtaMin, int iEtaMax, int iPhiMin, int iPhiMax, const std::vector< l1t::CaloTower > &towers, int minHwEt, int maxHwEt, SubDet etMode=CALO)
static int	caloEta (int ietaMP)
static size_t	caloTowerHash (int iEta, int iPhi)
static size_t	caloTowerHashMax ()
static l1t::EGamma	egP4Demux (l1t::EGamma &)
static l1t::EGamma	egP4MP (l1t::EGamma &)
static l1t::EtSum	etSumP4Demux (l1t::EtSum &)
static l1t::EtSum	etSumP4MP (l1t::EtSum &)
static const l1t::CaloCluster &	getCluster (const std::vector< l1t::CaloCluster > &clusters, int iEta, int iPhi)
static const l1t::CaloTower &	getTower (const std::vector< l1t::CaloTower > &towers, int iEta, int iPhi)
static unsigned int	gloriousDivision (uint32_t aNumerator, uint32_t aDenominator)
static int	gtEta (int ieta)
static int	gtPhi (int ieta, int iphi)
static bool	insertTower (std::vector< l1t::CaloTower > &towers, const l1t::CaloTower &tower)
static bool	isValidIEtaIPhi (int iEta, int iPhi)
static l1t::Jet	jetP4Demux (l1t::Jet &)
static l1t::Jet	jetP4MP (l1t::Jet &)
static int	mpEta (int ieta)
static math::PtEtaPhiMLorentzVector	p4Demux (l1t::L1Candidate *)
static math::PtEtaPhiMLorentzVector	p4MP (l1t::L1Candidate *)
static int	regionEta (int ieta)
static l1t::Tau	tauP4Demux (l1t::Tau &)
static l1t::Tau	tauP4MP (l1t::Tau &)
static float	towerEta (int ieta)
static float	towerEtaSize (int ieta)
static float	towerPhi (int ieta, int iphi)
static float	towerPhiSize (int ieta)

Static Public Attributes
static const int64_t	cos_coeff [72]
static const int	emul_to_data_sum_index_map [31]
static const int	kHBHEEnd = 28
static const int	kHBHENrPhi = 72
static const int	kHFBegin = 29
static const int	kHFEnd = 41
static const int	kHFNrPhi = 72 / kHFPhiSeg
static const int	kHFPhiSeg = 1
static const int	kNPhi = 72
static const int	kNrHBHETowers = kHBHEEnd * kHBHENrPhi * 2
static const int	kNrTowers = ((kHFEnd - kHFBegin + 1) * kHFNrPhi + kHBHEEnd * kHBHENrPhi) * 2
static const int	kSatEcal = 510
static const int	kSatHcal = 509
static const int	kSatJet = 65535
static const int	kSatTower = 511
static const int64_t	sin_coeff [72]

Private Member Functions
	CaloTools ()
	~CaloTools ()

Static Private Member Functions
static std::pair< float, float >	towerEtaBounds (int ieta)

Static Private Attributes
static const float	kGTEtaLSB = 0.0435
static const float	kGTEtLSB = 0.5
static const float	kGTPhiLSB = 0.0435
static const l1t::CaloCluster	nullCluster_
static const l1t::CaloTower	nullTower_

Detailed Description

Description: A collection of useful functions for the Calorimeter that are of generic interest

Implementation: currently implimented as a static class rather than a namespace, open to re-writing it as namespace

Author

: Sam Harper - RAL

Definition at line 29 of file CaloTools.h.

Member Enumeration Documentation

enum l1t::CaloTools::SubDet

Enumerator
ECAL
HCAL
CALO

Definition at line 60 of file CaloTools.h.

{ ECAL = 0x1, HCAL = 0x2, CALO = 0x3 };  //CALO is a short cut for ECAL|HCAL

Constructor & Destructor Documentation

l1t::CaloTools::CaloTools ( )

inlineprivate

Definition at line 32 of file CaloTools.h.

{}

l1t::CaloTools::~CaloTools ( )

inlineprivate

Definition at line 33 of file CaloTools.h.

{}

Member Function Documentation

int l1t::CaloTools::bin16Eta ( int  ieta )

static

Definition at line 265 of file CaloTools.cc.

References funct::abs().

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibrate().

                                   {
  int absIEta = abs(ieta);

  if (absIEta > 0 && absIEta <= 5)
    return 0;
  else if (absIEta <= 9)
    return 1;
  else if (absIEta <= 13)
    return 2;
  else if (absIEta <= 15)
    return 3;
  else if (absIEta <= 17)
    return 4;
  else if (absIEta <= 19)
    return 5;
  else if (absIEta <= 21)
    return 6;
  else if (absIEta == 22)
    return 7;
  else if (absIEta == 23)
    return 8;
  else if (absIEta == 24)
    return 9;
  else if (absIEta == 25)
    return 10;
  else if (absIEta == 26)
    return 11;
  else if (absIEta <= 28)
    return 12;
  else if (absIEta <= 32)
    return 13;
  else if (absIEta <= 36)
    return 14;
  else if (absIEta <= 41)
    return 15;
  else
    return -1;  // error
}

int l1t::CaloTools::calHwEtSum ( int  iEta, int  iPhi, const std::vector< l1t::CaloTower > &  towers, int  localEtaMin, int  localEtaMax, int  localPhiMin, int  localPhiMax, SubDet  etMode = CALO  )

static

Definition at line 114 of file CaloTools.cc.

Referenced by l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::isoCalEgHwFootPrint(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), and l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent().

                                              {
  return calHwEtSum(iEta, iPhi, towers, localEtaMin, localEtaMax, localPhiMin, localPhiMax, kHFEnd, etMode);
}

int l1t::CaloTools::calHwEtSum ( int  iEta, int  iPhi, const std::vector< l1t::CaloTower > &  towers, int  localEtaMin, int  localEtaMax, int  localPhiMin, int  localPhiMax, int  iEtaAbsMax, SubDet  etMode = CALO  )

static

Definition at line 125 of file CaloTools.cc.

References funct::abs(), ECAL, HCAL, l1t::CaloTower::hwEtEm(), l1t::CaloTower::hwEtHad(), l1t::L1Candidate::hwPt(), l1t::CaloStage2Nav::offsetIEta(), l1t::CaloStage2Nav::offsetIPhi(), and hgcalTowerProducer_cfi::tower.

                                              {
  int hwEtSum = 0;
  for (int etaNr = localEtaMin; etaNr <= localEtaMax; etaNr++) {
    for (int phiNr = localPhiMin; phiNr <= localPhiMax; phiNr++) {
      int towerIEta = l1t::CaloStage2Nav::offsetIEta(iEta, etaNr);
      int towerIPhi = l1t::CaloStage2Nav::offsetIPhi(iPhi, phiNr);
      if (abs(towerIEta) <= iEtaAbsMax) {
        const l1t::CaloTower& tower = getTower(towers, towerIEta, towerIPhi);
        if (etMode == ECAL)
          hwEtSum += tower.hwEtEm();
        else if (etMode == HCAL)
          hwEtSum += tower.hwEtHad();
        else if (etMode == CALO)
          hwEtSum += tower.hwPt();
      }
    }
  }
  return hwEtSum;
}

size_t l1t::CaloTools::calNrTowers ( int  iEtaMin, int  iEtaMax, int  iPhiMin, int  iPhiMax, const std::vector< l1t::CaloTower > &  towers, int  minHwEt, int  maxHwEt, SubDet  etMode = CALO  )

static

Definition at line 153 of file CaloTools.cc.

References caloEta(), l1t::CaloStage2Nav::currIEta(), l1t::CaloStage2Nav::currIPhi(), l1t::CaloStage2Nav::east(), ECAL, getTower(), HCAL, l1t::CaloTower::hwEtEm(), l1t::CaloTower::hwEtHad(), l1t::L1Candidate::hwPt(), l1t::CaloStage2Nav::north(), l1t::CaloStage2Nav::resetIPhi(), and hgcalTowerProducer_cfi::tower.

Referenced by l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent().

                                                  {
  size_t nrTowers = 0;
  l1t::CaloStage2Nav nav(iEtaMin, iPhiMin);
  while (nav.currIEta() <= iEtaMax) {
    bool finishPhi = false;
    while (!finishPhi) {
      const l1t::CaloTower& tower =
          l1t::CaloTools::getTower(towers, CaloTools::caloEta(nav.currIEta()), nav.currIPhi());
      int towerHwEt = 0;
      if (etMode == ECAL)
        towerHwEt += tower.hwEtEm();
      else if (etMode == HCAL)
        towerHwEt += tower.hwEtHad();
      else if (etMode == CALO)
        towerHwEt += tower.hwPt();
      if (towerHwEt >= minHwEt && towerHwEt <= maxHwEt)
        nrTowers++;
      finishPhi = (nav.currIPhi() == iPhiMax);
      nav.north();
    }
    nav.east();
    nav.resetIPhi();
  }
  return nrTowers;
}

int l1t::CaloTools::caloEta ( int  ietaMP )

static

Definition at line 234 of file CaloTools.cc.

Referenced by L1TStage2InputPatternWriter::analyze(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::calibratedPt(), calNrTowers(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::donutPUEstimate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getChunkyRing(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getSumEtEtaMap(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), l1t::stage2::CaloTowerUnpacker::unpack(), l1t::stage2::MPUnpacker_0x10010033::unpack(), l1t::stage2::MPUnpacker_0x10010010::unpack(), l1t::stage2::MPUnpacker::unpack(), and l1t::stage2::MPUnpacker_0x1001000b::unpack().

                                   {
  if (mpEta >= kHFBegin)
    return mpEta + 1;
  else if (mpEta <= -1 * kHFBegin)
    return mpEta - 1;
  else
    return mpEta;
}

size_t l1t::CaloTools::caloTowerHash ( int  iEta, int  iPhi  )

static

Definition at line 81 of file CaloTools.cc.

References funct::abs().

Referenced by getTower(), insertTower(), l1t::stage2::CaloTowerPacker::pack(), and L1TStage2Layer2Producer::produce().

                                                     {
  if (!isValidIEtaIPhi(iEta, iPhi))
    return caloTowerHashMax();
  else {
    const int absIEta = abs(iEta);
    if (absIEta > kHFEnd)
      return kNrTowers;
    else if (absIEta <= kHBHEEnd) {  //HBHE
      int iEtaNoZero = iEta;
      if (iEta > 0)
        iEtaNoZero--;
      return (iEtaNoZero + kHBHEEnd) * kHBHENrPhi + iPhi - 1;
    } else {                          //HF
      int iEtaIndex = iEta + kHFEnd;  //iEta=-32 is 0
      if (iEta > 0)
        iEtaIndex = iEta - kHBHEEnd + (kHFEnd - kHBHEEnd) - 1;  //but iEta=29 is 4
      return iEtaIndex * kHFNrPhi + iPhi / kHFPhiSeg + kNrHBHETowers;
    }
  }
}

size_t l1t::CaloTools::caloTowerHashMax ( )

static

Definition at line 102 of file CaloTools.cc.

Referenced by L1TStage2Layer2Producer::produce().

{ return kNrTowers; }

l1t::EGamma l1t::CaloTools::egP4Demux ( l1t::EGamma &  eg )

static

Definition at line 322 of file CaloTools.cc.

References l1t::EGamma::footprintEt(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::EGamma::isoEt(), l1t::EGamma::nTT(), l1t::EGamma::rawEt(), l1t::EGamma::setTowerIPhi(), l1t::EGamma::shape(), l1t::EGamma::towerHoE(), l1t::EGamma::towerIEta(), and l1t::EGamma::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                               {
  l1t::EGamma tmpEG(p4Demux(&eg), eg.hwPt(), eg.hwEta(), eg.hwPhi(), eg.hwQual(), eg.hwIso());
  tmpEG.setTowerIPhi(eg.towerIPhi());
  tmpEG.setTowerIEta(eg.towerIEta());
  tmpEG.setRawEt(eg.rawEt());
  tmpEG.setIsoEt(eg.isoEt());
  tmpEG.setFootprintEt(eg.footprintEt());
  tmpEG.setNTT(eg.nTT());
  tmpEG.setShape(eg.shape());
  tmpEG.setTowerHoE(eg.towerHoE());

  return tmpEG;
}

l1t::EGamma l1t::CaloTools::egP4MP ( l1t::EGamma &  eg )

static

Definition at line 374 of file CaloTools.cc.

References l1t::EGamma::footprintEt(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::EGamma::isoEt(), l1t::EGamma::nTT(), l1t::EGamma::rawEt(), l1t::EGamma::setTowerIPhi(), l1t::EGamma::shape(), l1t::EGamma::towerHoE(), l1t::EGamma::towerIEta(), and l1t::EGamma::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                            {
  l1t::EGamma tmpEG(p4MP(&eg), eg.hwPt(), eg.hwEta(), eg.hwPhi(), eg.hwQual(), eg.hwIso());
  tmpEG.setTowerIPhi(eg.towerIPhi());
  tmpEG.setTowerIEta(eg.towerIEta());
  tmpEG.setRawEt(eg.rawEt());
  tmpEG.setIsoEt(eg.isoEt());
  tmpEG.setFootprintEt(eg.footprintEt());
  tmpEG.setNTT(eg.nTT());
  tmpEG.setShape(eg.shape());
  tmpEG.setTowerHoE(eg.towerHoE());

  return tmpEG;
}

l1t::EtSum l1t::CaloTools::etSumP4Demux ( l1t::EtSum &  etsum )

static

Definition at line 364 of file CaloTools.cc.

References l1t::EtSum::getType(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), and l1t::L1Candidate::hwQual().

Referenced by L1TStage2Layer2Producer::produce().

                                                   {
  return l1t::EtSum(p4Demux(&etsum), etsum.getType(), etsum.hwPt(), etsum.hwEta(), etsum.hwPhi(), etsum.hwQual());
}

l1t::EtSum l1t::CaloTools::etSumP4MP ( l1t::EtSum &  etsum )

static

Definition at line 416 of file CaloTools.cc.

References l1t::EtSum::getType(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), and l1t::L1Candidate::hwQual().

Referenced by L1TStage2Layer2Producer::produce().

                                                {
  return l1t::EtSum(p4MP(&etsum), etsum.getType(), etsum.hwPt(), etsum.hwEta(), etsum.hwPhi(), etsum.hwQual());
}

const l1t::CaloCluster & l1t::CaloTools::getCluster ( const std::vector< l1t::CaloCluster > &  clusters, int  iEta, int  iPhi  )

static

Definition at line 70 of file CaloTools.cc.

Referenced by l1t::Stage2Layer2ClusterAlgorithmFirmwareImp1::filtering().

                                                                                                            {
  for (size_t clusterNr = 0; clusterNr < clusters.size(); clusterNr++) {
    if (clusters[clusterNr].hwEta() == iEta && clusters[clusterNr].hwPhi() == iPhi)
      return clusters[clusterNr];
  }
  return nullCluster_;
}

const l1t::CaloTower & l1t::CaloTools::getTower ( const std::vector< l1t::CaloTower > &  towers, int  iEta, int  iPhi  )

static

Definition at line 43 of file CaloTools.cc.

References funct::abs(), caloTowerHash(), and kHFEnd.

Referenced by L1TStage2InputPatternWriter::analyze(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::calibratedPt(), calNrTowers(), l1t::Stage2Layer2ClusterAlgorithmFirmwareImp1::clustering(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::donutPUEstimate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getChunkyRing(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getSumEtEtaMap(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::is3x3Maximum(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::makeSecClusters(), l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2ClusterAlgorithmFirmwareImp1::refining().

                                                                                                    {
  if (abs(iEta) > CaloTools::kHFEnd)
    return nullTower_;

  size_t towerIndex = CaloTools::caloTowerHash(iEta, iPhi);
  if (towerIndex < towers.size()) {
    if (towers[towerIndex].hwEta() != iEta ||
        towers[towerIndex].hwPhi() !=
            iPhi) {  //it failed, this is bad, but we will not log the error due to policy and silently attempt to do a brute force search instead
      //std::cout <<"error, tower "<<towers[towerIndex].hwEta()<<" "<<towers[towerIndex].hwPhi()<<" does not match "<<iEta<<" "<<iPhi<<" index "<<towerIndex<<" nr towrs "<<towers.size()<<std::endl;
      for (size_t towerNr = 0; towerNr < towers.size(); towerNr++) {
        if (towers[towerNr].hwEta() == iEta && towers[towerNr].hwPhi() == iPhi)
          return towers[towerNr];
      }
    } else
      return towers[towerIndex];

  } else {  // in case the vector of towers do not contain all the towers (towerIndex can be > towers.size())
    for (size_t towerNr = 0; towerNr < towers.size(); towerNr++) {
      if (towers[towerNr].hwEta() == iEta && towers[towerNr].hwPhi() == iPhi)
        return towers[towerNr];
    }
  }

  return nullTower_;
}

unsigned int l1t::CaloTools::gloriousDivision ( uint32_t  aNumerator, uint32_t  aDenominator  )

static

Definition at line 419 of file CaloTools.cc.

References mps_fire::result.

Referenced by l1t::Stage2Layer2DemuxSumsAlgoFirmwareImp1::processEvent().

                                                                                      {
  static const uint64_t lLut[] = {
      0,     16777215, 4194304, 1864135, 1048576, 671089, 466034, 342392, 262144, 207126, 167772, 138655, 116508, 99273,
      85598, 74565,    65536,   58053,   51782,   46474,  41943,  38044,  34664,  31715,  29127,  26844,  24818,  23014,
      21400, 19949,    18641,   17458,   16384,   15406,  14513,  13696,  12945,  12255,  11619,  11030,  10486,  9980,
      9511,  9074,     8666,    8285,    7929,    7595,   7282,   6988,   6711,   6450,   6205,   5973,   5754,   5546,
      5350,  5164,     4987,    4820,    4660,    4509,   4365,   4227,   4096,   3971,   3852,   3737,   3628,   3524,
      3424,  3328,     3236,    3148,    3064,    2983,   2905,   2830,   2758,   2688,   2621,   2557,   2495,   2435,
      2378,  2322,     2268,    2217,    2166,    2118,   2071,   2026,   1982,   1940,   1899,   1859,   1820,   1783,
      1747,  1712,     1678,    1645,    1613,    1581,   1551,   1522,   1493,   1465,   1438,   1412,   1387,   1362,
      1337,  1314,     1291,    1269,    1247,    1226,   1205,   1185,   1165,   1146,   1127,   1109,   1091,   1074,
      1057,  1040,     1024,    1008,    993,     978,    963,    948,    934,    921,    907,    894,    881,    868,
      856,   844,      832,     820,     809,     798,    787,    776,    766,    756,    746,    736,    726,    717,
      707,   698,      689,     681,     672,     664,    655,    647,    639,    631,    624,    616,    609,    602,
      594,   587,      581,     574,     567,     561,    554,    548,    542,    536,    530,    524,    518,    512,
      506,   501,      496,     490,     485,     480,    475,    470,    465,    460,    455,    450,    446,    441,
      437,   432,      428,     424,     419,     415,    411,    407,    403,    399,    395,    392,    388,    384,
      380,   377,      373,     370,     366,     363,    360,    356,    353,    350,    347,    344,    340,    337,
      334,   331,      328,     326,     323,     320,    317,    314,    312,    309,    306,    304,    301,    299,
      296,   294,      291,     289,     286,     284,    282,    280,    277,    275,    273,    271,    268,    266,
      264,   262,      260,     258,     256,     254,    252,    250,    248,    246,    244,    243,    241,    239,
      237,   235,      234,     232,     230,     228,    227,    225,    223,    222,    220,    219,    217,    216,
      214,   212,      211,     209,     208,     207,    205,    204,    202,    201,    199,    198,    197,    195,
      194,   193,      191,     190,     189,     188,    186,    185,    184,    183,    182,    180,    179,    178,
      177,   176,      175,     173,     172,     171,    170,    169,    168,    167,    166,    165,    164,    163,
      162,   161,      160,     159,     158,     157,    156,    155,    154,    153,    152,    151,    150,    149,
      149,   148,      147,     146,     145,     144,    143,    143,    142,    141,    140,    139,    139,    138,
      137,   136,      135,     135,     134,     133,    132,    132,    131,    130,    129,    129,    128,    127,
      127,   126,      125,     125,     124,     123,    123,    122,    121,    121,    120,    119,    119,    118,
      117,   117,      116,     116,     115,     114,    114,    113,    113,    112,    111,    111,    110,    110,
      109,   109,      108,     108,     107,     106,    106,    105,    105,    104,    104,    103,    103,    102,
      102,   101,      101,     100,     100,     99,     99,     98,     98,     97,     97,     96,     96,     96,
      95,    95,       94,      94,      93,      93,     92,     92,     92,     91,     91,     90,     90,     89,
      89,    89,       88,      88,      87,      87,     87,     86,     86,     85,     85,     85,     84,     84,
      84,    83,       83,      82,      82,      82,     81,     81,     81,     80,     80,     80,     79,     79,
      79,    78,       78,      78,      77,      77,     77,     76,     76,     76,     75,     75,     75,     74,
      74,    74,       73,      73,      73,      73,     72,     72,     72,     71,     71,     71,     70,     70,
      70,    70,       69,      69,      69,      68,     68,     68,     68,     67,     67,     67,     67,     66,
      66,    66,       66,      65,      65,      65,     65,     64};

  // Firmware uses 18bit integers - make sure we are in the same range
  aNumerator &= 0x3FFFF;
  aDenominator &= 0x3FFFF;

  // Shift the denominator to optimise the polynomial expansion
  // I limit the shift to half the denominator size in the firmware to save on resources
  uint32_t lBitShift(0);
  for (; lBitShift != 9; ++lBitShift) {
    if (aDenominator & 0x20000)
      break;  // There is a 1 in the MSB
    aDenominator <<= 1;
  }

  // The magical polynomial expansion Voodoo
  uint64_t lInverseDenominator(((aDenominator & 0x3FE00) - (aDenominator & 0x001FF)) * (lLut[aDenominator >> 9]));

  // Save on DSPs by throwing away a bunch of LSBs
  lInverseDenominator >>= 17;

  // Multiply the numerator by the inverse denominator
  uint64_t lResult(aNumerator * lInverseDenominator);

  // Add two bits to the result, to make the Voodoo below work (saves us having an if-else on the shift direction)
  lResult <<= 2;

  // Restore the scale by taking into account the bitshift applied above.
  // We are now 18 bit left-shifted, so the 18 LSBs are effectively the fractional part...

  uint32_t aFractional = (lResult >>= (9 - lBitShift)) & 0x3FFFF;
  // ...and the top 18 bits are the integer part

  // uint32_t aInteger    = ( lResult >>= 18 ) & 0x3FFFF;

  unsigned int result = aFractional >> 10;

  return result;

  // Simples!
}

int l1t::CaloTools::gtEta ( int  ieta )

static

Definition at line 304 of file CaloTools.cc.

References PVValHelper::eta.

Referenced by l1t::Stage2Layer2DemuxJetAlgoFirmwareImp1::processEvent().

                                {
  double eta = towerEta(ieta);
  return round(eta / kGTEtaLSB);
}

int l1t::CaloTools::gtPhi ( int  ieta, int  iphi  )

static

Definition at line 309 of file CaloTools.cc.

References M_PI.

Referenced by l1t::Stage2Layer2DemuxJetAlgoFirmwareImp1::processEvent().

                                          {
  double phi = towerPhi(ieta, iphi);
  if (phi < 0)
    phi = phi + 2 * M_PI;
  return round(phi / kGTPhiLSB);
}

bool l1t::CaloTools::insertTower ( std::vector< l1t::CaloTower > &  towers, const l1t::CaloTower &  tower  )

static

Definition at line 32 of file CaloTools.cc.

References caloTowerHash(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), and hgcalTowerProducer_cfi::tower.

                                                                                         {
  size_t towerIndex = CaloTools::caloTowerHash(tower.hwEta(), tower.hwPhi());
  if (towers.size() > towerIndex) {
    towers.at(towerIndex) = tower;
    return true;
  } else
    return false;
}

bool l1t::CaloTools::isValidIEtaIPhi ( int  iEta, int  iPhi  )

static

Definition at line 104 of file CaloTools.cc.

References funct::abs().

                                                     {
  size_t absIEta = abs(iEta);
  if (iPhi <= 0 || iPhi > kNPhi)
    return false;
  if (absIEta == 0 || absIEta > kHFEnd)
    return false;
  //if(absIEta>kHBHEEnd && iPhi%kHFPhiSeg!=1) return false;
  return true;
}

l1t::Jet l1t::CaloTools::jetP4Demux ( l1t::Jet &  jet )

static

Definition at line 349 of file CaloTools.cc.

References l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::Jet::puDonutEt(), l1t::Jet::puEt(), l1t::Jet::rawEt(), l1t::Jet::seedEt(), l1t::Jet::setTowerIPhi(), l1t::Jet::towerIEta(), and l1t::Jet::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                           {
  l1t::Jet tmpJet(p4Demux(&jet), jet.hwPt(), jet.hwEta(), jet.hwPhi(), jet.hwQual());
  tmpJet.setTowerIPhi(jet.towerIPhi());
  tmpJet.setTowerIEta(jet.towerIEta());
  tmpJet.setRawEt(jet.rawEt());
  tmpJet.setSeedEt(jet.seedEt());
  tmpJet.setPUEt(jet.puEt());
  tmpJet.setPUDonutEt(0, jet.puDonutEt(0));
  tmpJet.setPUDonutEt(1, jet.puDonutEt(1));
  tmpJet.setPUDonutEt(2, jet.puDonutEt(2));
  tmpJet.setPUDonutEt(3, jet.puDonutEt(3));

  return tmpJet;
}

l1t::Jet l1t::CaloTools::jetP4MP ( l1t::Jet &  jet )

static

Definition at line 401 of file CaloTools.cc.

References l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::Jet::puDonutEt(), l1t::Jet::puEt(), l1t::Jet::rawEt(), l1t::Jet::seedEt(), l1t::Jet::setTowerIPhi(), l1t::Jet::towerIEta(), and l1t::Jet::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                        {
  l1t::Jet tmpJet(p4MP(&jet), jet.hwPt(), jet.hwEta(), jet.hwPhi(), jet.hwQual());
  tmpJet.setTowerIPhi(jet.towerIPhi());
  tmpJet.setTowerIEta(jet.towerIEta());
  tmpJet.setRawEt(jet.rawEt());
  tmpJet.setSeedEt(jet.seedEt());
  tmpJet.setPUEt(jet.puEt());
  tmpJet.setPUDonutEt(0, jet.puDonutEt(0));
  tmpJet.setPUDonutEt(1, jet.puDonutEt(1));
  tmpJet.setPUDonutEt(2, jet.puDonutEt(2));
  tmpJet.setPUDonutEt(3, jet.puDonutEt(3));

  return tmpJet;
}

int l1t::CaloTools::mpEta ( int  ieta )

static

Definition at line 224 of file CaloTools.cc.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibrate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::donutPUEstimate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getChunkyRing(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2DemuxJetAlgoFirmwareImp1::processEvent().

                                {
  if (ieta > kHFBegin)
    return ieta - 1;
  else if (ieta < -1 * kHFBegin)
    return ieta + 1;
  else
    return ieta;
}

math::PtEtaPhiMLorentzVector l1t::CaloTools::p4Demux ( l1t::L1Candidate *  cand )

static

Definition at line 317 of file CaloTools.cc.

References l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), and l1t::L1Candidate::hwPt().

Referenced by l1t::stage2::EtSumUnpacker::unpack(), l1t::stage2::TauUnpacker::unpack(), l1t::stage2::JetUnpacker::unpack(), l1t::stage2::EtSumUnpacker_0x10010057::unpack(), and l1t::stage2::EGammaUnpacker::unpack().

                                                                     {
  return math::PtEtaPhiMLorentzVector(
      cand->hwPt() * kGTEtLSB + 1.E-6, cand->hwEta() * kGTEtaLSB, cand->hwPhi() * kGTPhiLSB, 0.);
}

math::PtEtaPhiMLorentzVector l1t::CaloTools::p4MP ( l1t::L1Candidate *  cand )

static

Definition at line 369 of file CaloTools.cc.

References l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), and l1t::L1Candidate::hwPt().

Referenced by l1t::stage2::MPUnpacker::unpack(), l1t::stage2::MPUnpacker_0x10010033::unpack(), l1t::stage2::MPUnpacker_0x10010010::unpack(), and l1t::stage2::MPUnpacker_0x1001000b::unpack().

                                                                  {
  return math::PtEtaPhiMLorentzVector(
      cand->hwPt() * 0.5 + 1.E-6, towerEta(cand->hwEta()), towerPhi(cand->hwEta(), cand->hwPhi()), 0.);
}

int l1t::CaloTools::regionEta ( int  ieta )

static

Definition at line 244 of file CaloTools.cc.

References funct::abs(), and reco::ceil().

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibrate().

                                    {
  // outside HF
  if (abs(ieta) > kHFEnd)
    return (ieta < 0 ? 0 : 21);

  // inside HBHE
  if (abs(ieta) <= kHFBegin) {
    if (ieta < 0)
      return 11 - ceil(double(abs(ieta) / 4.));
    else
      return ceil(double(abs(ieta) / 4.)) + 10;
  }

  // in HF
  if (ieta < 0)
    return 4 - ceil(double(abs(ieta) - 29) / 4.);
  else
    return ceil(double(abs(ieta) - 29) / 4.) + 17;
}

l1t::Tau l1t::CaloTools::tauP4Demux ( l1t::Tau &  tau )

static

Definition at line 336 of file CaloTools.cc.

References l1t::Tau::hasEM(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::Tau::isMerged(), l1t::Tau::isoEt(), l1t::Tau::nTT(), l1t::Tau::rawEt(), l1t::Tau::setTowerIPhi(), l1t::Tau::towerIEta(), and l1t::Tau::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                           {
  l1t::Tau tmpTau(p4Demux(&tau), tau.hwPt(), tau.hwEta(), tau.hwPhi(), tau.hwQual(), tau.hwIso());
  tmpTau.setTowerIPhi(tau.towerIPhi());
  tmpTau.setTowerIEta(tau.towerIEta());
  tmpTau.setRawEt(tau.rawEt());
  tmpTau.setIsoEt(tau.isoEt());
  tmpTau.setNTT(tau.nTT());
  tmpTau.setHasEM(tau.hasEM());
  tmpTau.setIsMerged(tau.isMerged());

  return tmpTau;
}

l1t::Tau l1t::CaloTools::tauP4MP ( l1t::Tau &  tau )

static

Definition at line 388 of file CaloTools.cc.

References l1t::Tau::hasEM(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), l1t::Tau::isMerged(), l1t::Tau::isoEt(), l1t::Tau::nTT(), l1t::Tau::rawEt(), l1t::Tau::setTowerIPhi(), l1t::Tau::towerIEta(), and l1t::Tau::towerIPhi().

Referenced by L1TStage2Layer2Producer::produce().

                                        {
  l1t::Tau tmpTau(p4MP(&tau), tau.hwPt(), tau.hwEta(), tau.hwPhi(), tau.hwQual(), tau.hwIso());
  tmpTau.setTowerIPhi(tau.towerIPhi());
  tmpTau.setTowerIEta(tau.towerIEta());
  tmpTau.setRawEt(tau.rawEt());
  tmpTau.setIsoEt(tau.isoEt());
  tmpTau.setNTT(tau.nTT());
  tmpTau.setHasEM(tau.hasEM());
  tmpTau.setIsMerged(tau.isMerged());

  return tmpTau;
}

float l1t::CaloTools::towerEta ( int  ieta )

static

Definition at line 201 of file CaloTools.cc.

References PVValHelper::eta, and jetcorrextractor::sign().

Referenced by l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent(), L1TowerCalibrator::produce(), L1EGCrystalClusterEmulatorProducer::produce(), and L1TPFCaloProducer::readHcalDigis_().

                                     {
  std::pair<float, float> bounds = towerEtaBounds(ieta);
  float eta = (bounds.second + bounds.first) / 2.;
  float sign = ieta > 0 ? 1. : -1.;
  return sign * eta;
}

std::pair< float, float > l1t::CaloTools::towerEtaBounds ( int  ieta )

staticprivate

Definition at line 186 of file CaloTools.cc.

References funct::abs().

                                                           {
  if (ieta == 0)
    ieta = 1;
  if (ieta > kHFEnd)
    ieta = kHFEnd;
  if (ieta < (-1 * kHFEnd))
    ieta = -1 * kHFEnd;
  //const float towerEtas[33] = {0,0.087,0.174,0.261,0.348,0.435,0.522,0.609,0.696,0.783,0.870,0.957,1.044,1.131,1.218,1.305,1.392,1.479,1.566,1.653,1.740,1.830,1.930,2.043,2.172,2.322,2.5,2.650,3.000,3.5,4.0,4.5,5.0};
  const float towerEtas[42] = {0,     0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.696, 0.783, 0.870,
                               0.957, 1.044, 1.131, 1.218, 1.305, 1.392, 1.479, 1.566, 1.653, 1.740, 1.830,
                               1.930, 2.043, 2.172, 2.322, 2.5,   2.650, 2.853, 3.139, 3.314, 3.489, 3.664,
                               3.839, 4.013, 4.191, 4.363, 4.538, 4.716, 4.889, 5.191, 5.191};
  return std::make_pair(towerEtas[abs(ieta) - 1], towerEtas[abs(ieta)]);
}

float l1t::CaloTools::towerEtaSize ( int  ieta )

static

Definition at line 215 of file CaloTools.cc.

References findQualityFiles::size.

Referenced by l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), and l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent().

                                         {
  std::pair<float, float> bounds = towerEtaBounds(ieta);
  float size = (bounds.second - bounds.first);
  return size;
}

float l1t::CaloTools::towerPhi ( int  ieta, int  iphi  )

static

Definition at line 208 of file CaloTools.cc.

References M_PI.

Referenced by l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent(), L1TowerCalibrator::produce(), L1EGCrystalClusterEmulatorProducer::produce(), and L1TPFCaloProducer::readHcalDigis_().

                                               {
  float phi = (float(iphi) - 0.5) * towerPhiSize(ieta);
  if (phi > M_PI)
    phi = phi - (2 * M_PI);
  return phi;
}

float l1t::CaloTools::towerPhiSize ( int  ieta )

static

Definition at line 221 of file CaloTools.cc.

References M_PI.

Referenced by l1t::Stage2Layer2TauAlgorithmFirmwareImp1::merging(), and l1t::Stage2Layer2EGammaAlgorithmFirmwareImp1::processEvent().

{ return 2. * M_PI / kNPhi; }

Member Data Documentation

const int64_t l1t::CaloTools::cos_coeff

static

Initial value:

= {
    1023,  1019,  1007,  988,  961,  927,  886,  838,  784,  723,  658,  587,  512,  432,  350,  265,  178,   89,
    0,     -89,   -178,  -265, -350, -432, -512, -587, -658, -723, -784, -838, -886, -927, -961, -988, -1007, -1019,
    -1023, -1019, -1007, -988, -961, -927, -886, -838, -784, -723, -658, -587, -512, -432, -350, -265, -178,  -89,
    0,     89,    178,   265,  350,  432,  511,  587,  658,  723,  784,  838,  886,  927,  961,  988,  1007,  1019}

Definition at line 137 of file CaloTools.h.

Referenced by l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent().

const int l1t::CaloTools::emul_to_data_sum_index_map

static

Initial value:

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

Definition at line 141 of file CaloTools.h.

Referenced by L1TStage2CaloLayer2Comp::compareSums().

const float l1t::CaloTools::kGTEtaLSB = 0.0435

staticprivate

Definition at line 152 of file CaloTools.h.

const float l1t::CaloTools::kGTEtLSB = 0.5

staticprivate

Definition at line 154 of file CaloTools.h.

const float l1t::CaloTools::kGTPhiLSB = 0.0435

staticprivate

Definition at line 153 of file CaloTools.h.

const int l1t::CaloTools::kHBHEEnd = 28

static

Definition at line 38 of file CaloTools.h.

const int l1t::CaloTools::kHBHENrPhi = 72

static

Definition at line 43 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::donutPUEstimate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getChunkyRing(), l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent().

const int l1t::CaloTools::kHFBegin = 29

static

Definition at line 39 of file CaloTools.h.

Referenced by l1t::CaloStage2Nav::offsetIEta(), l1t::stage2::CaloTowerPacker::pack(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), and L1TowerCalibrator::produce().

const int l1t::CaloTools::kHFEnd = 41

static

Definition at line 40 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::donutPUEstimate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::getChunkyRing(), getTower(), l1t::CaloStage2Nav::offsetIEta(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent(), L1TowerCalibrator::produce(), l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::Stage2Layer2EtSumAlgorithmFirmwareImp1(), and l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::Stage2Layer2JetSumAlgorithmFirmwareImp1().

const int l1t::CaloTools::kHFNrPhi = 72 / kHFPhiSeg

static

Definition at line 42 of file CaloTools.h.

const int l1t::CaloTools::kHFPhiSeg = 1

static

Definition at line 41 of file CaloTools.h.

const int l1t::CaloTools::kNPhi = 72

static

Definition at line 44 of file CaloTools.h.

const int l1t::CaloTools::kNrHBHETowers = kHBHEEnd * kHBHENrPhi * 2

static

Definition at line 46 of file CaloTools.h.

const int l1t::CaloTools::kNrTowers = ((kHFEnd - kHFBegin + 1) * kHFNrPhi + kHBHEEnd * kHBHENrPhi) * 2

static

Definition at line 45 of file CaloTools.h.

const int l1t::CaloTools::kSatEcal = 510

static

Definition at line 53 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), and l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent().

const int l1t::CaloTools::kSatHcal = 509

static

Definition at line 52 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), and l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent().

const int l1t::CaloTools::kSatJet = 65535

static

Definition at line 57 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::calibrate(), l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), and l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent().

const int l1t::CaloTools::kSatTower = 511

static

Definition at line 54 of file CaloTools.h.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), and l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent().

const l1t::CaloCluster l1t::CaloTools::nullCluster_

staticprivate

Definition at line 150 of file CaloTools.h.

const l1t::CaloTower l1t::CaloTools::nullTower_

staticprivate

Definition at line 148 of file CaloTools.h.

const int64_t l1t::CaloTools::sin_coeff

static

Initial value:

= {
    0,     89,    178,   265,  350,  432,  512,  587,  658,  723,  784,  838,  886,  927,  961,  988,  1007,  1019,
    1023,  1019,  1007,  988,  961,  927,  886,  838,  784,  723,  658,  587,  512,  432,  350,  265,  178,   89,
    0,     -89,   -178,  -265, -350, -432, -512, -587, -658, -723, -784, -838, -886, -927, -961, -988, -1007, -1019,
    -1023, -1019, -1007, -988, -961, -927, -886, -838, -784, -723, -658, -587, -512, -432, -350, -265, -178,  -89}

Definition at line 138 of file CaloTools.h.

Referenced by l1t::Stage2Layer2EtSumAlgorithmFirmwareImp1::processEvent(), and l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent().