CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Static Public Member Functions | Static Public Attributes | Private Member Functions | Static Private Member Functions | Static Private Attributes
l1t::CaloTools Class Reference

#include <CaloTools.h>

Public Types

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

Static Public Member Functions

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::CaloClustergetCluster (const std::vector< l1t::CaloCluster > &clusters, int iEta, int iPhi)
 
static const l1t::CaloTowergetTower (const std::vector< l1t::CaloTower > &towers, int iEta, int iPhi)
 
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 int32_t cos_coeff [72] = {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, 512, 587, 658, 723, 784, 838, 886, 927, 961, 988, 1007, 1019}
 
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 int32_t sin_coeff [72] = {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}
 

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

Enumerator
ECAL 
HCAL 
CALO 

Definition at line 50 of file CaloTools.h.

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

Constructor & Destructor Documentation

l1t::CaloTools::CaloTools ( )
inlineprivate

Definition at line 33 of file CaloTools.h.

33 {}
l1t::CaloTools::~CaloTools ( )
inlineprivate

Definition at line 34 of file CaloTools.h.

34 {}

Member Function Documentation

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 100 of file CaloTools.cc.

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

103 {
104 
105  return calHwEtSum(iEta,iPhi,towers,localEtaMin,localEtaMax,localPhiMin,localPhiMax,kHFEnd,etMode);
106 }
static int calHwEtSum(int iEta, int iPhi, const std::vector< l1t::CaloTower > &towers, int localEtaMin, int localEtaMax, int localPhiMin, int localPhiMax, SubDet etMode=CALO)
Definition: CaloTools.cc:100
static const int kHFEnd
Definition: CaloTools.h:41
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 108 of file CaloTools.cc.

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

111 {
112  int hwEtSum=0;
113  for(int etaNr=localEtaMin;etaNr<=localEtaMax;etaNr++){
114  for(int phiNr=localPhiMin;phiNr<=localPhiMax;phiNr++){
115 
116  int towerIEta = l1t::CaloStage2Nav::offsetIEta(iEta,etaNr);
117  int towerIPhi = l1t::CaloStage2Nav::offsetIPhi(iPhi,phiNr);
118  if(abs(towerIEta)<=iEtaAbsMax){
119  const l1t::CaloTower& tower = getTower(towers,towerIEta,towerIPhi);
120  if(etMode==ECAL) hwEtSum+=tower.hwEtEm();
121  else if(etMode==HCAL) hwEtSum+=tower.hwEtHad();
122  else if(etMode==CALO) hwEtSum+=tower.hwPt();
123  }
124  }
125  }
126  return hwEtSum;
127 }
static int offsetIEta(int iEta, int offset)
Definition: CaloStage2Nav.h:43
int hwEtEm() const
Definition: CaloTower.cc:64
static int offsetIPhi(int iPhi, int offset)
Definition: CaloStage2Nav.h:32
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
int hwPt() const
Definition: L1Candidate.cc:69
static const l1t::CaloTower & getTower(const std::vector< l1t::CaloTower > &towers, int iEta, int iPhi)
Definition: CaloTools.cc:30
int hwEtHad() const
Definition: CaloTower.cc:69
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 130 of file CaloTools.cc.

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

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

131 {
132  size_t nrTowers=0;
133  l1t::CaloStage2Nav nav(iEtaMin,iPhiMin);
134  while(nav.currIEta()<=iEtaMax){
135  bool finishPhi = false;
136  while(!finishPhi){
137  const l1t::CaloTower& tower = l1t::CaloTools::getTower(towers,nav.currIEta(),nav.currIPhi());
138  int towerHwEt =0;
139  if(etMode==ECAL) towerHwEt+=tower.hwEtEm();
140  else if(etMode==HCAL) towerHwEt+=tower.hwEtHad();
141  else if(etMode==CALO) towerHwEt+=tower.hwPt();
142  if(towerHwEt>=minHwEt && towerHwEt<=maxHwEt) nrTowers++;
143  finishPhi = (nav.currIPhi() == iPhiMax);
144  nav.north();
145  }
146  nav.east();
147  nav.resetIPhi();
148  }
149  return nrTowers;
150 }
int hwEtEm() const
Definition: CaloTower.cc:64
int hwPt() const
Definition: L1Candidate.cc:69
static const l1t::CaloTower & getTower(const std::vector< l1t::CaloTower > &towers, int iEta, int iPhi)
Definition: CaloTools.cc:30
int hwEtHad() const
Definition: CaloTower.cc:69
int l1t::CaloTools::caloEta ( int  ietaMP)
static

Definition at line 201 of file CaloTools.cc.

Referenced by l1t::Stage2Layer2JetAlgorithmFirmwareImp1::create(), l1t::stage2::CaloTowerUnpacker::unpack(), l1t::stage2::MPUnpacker_0x1001000b::unpack(), l1t::stage2::MPUnpacker_0x10010010::unpack(), and l1t::stage2::MPUnpacker::unpack().

201  {
202 
203  if (mpEta>=kHFBegin) return mpEta+1;
204  else if (mpEta<=-1*kHFBegin) return mpEta-1;
205  else return mpEta;
206 
207 }
static int mpEta(int ieta)
Definition: CaloTools.cc:191
static const int kHFBegin
Definition: CaloTools.h:40
size_t l1t::CaloTools::caloTowerHash ( int  iEta,
int  iPhi 
)
static

Definition at line 64 of file CaloTools.cc.

References funct::abs().

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

65 {
66 
67  if(!isValidIEtaIPhi(iEta,iPhi)) return caloTowerHashMax();
68  else{
69  const int absIEta = abs(iEta);
70  if(absIEta>kHFEnd) return kNrTowers;
71  else if(absIEta<=kHBHEEnd){ //HBHE
72  int iEtaNoZero=iEta;
73  if(iEta>0) iEtaNoZero--;
74  return (iEtaNoZero+kHBHEEnd)*kHBHENrPhi+iPhi-1;
75  }else{ //HF
76  int iEtaIndex = iEta+kHFEnd; //iEta=-32 is 0
77  if(iEta>0) iEtaIndex= iEta-kHBHEEnd+(kHFEnd-kHBHEEnd)-1; //but iEta=29 is 4
78  return iEtaIndex*kHFNrPhi+iPhi/kHFPhiSeg + kNrHBHETowers;
79  }
80  }
81 }
static const int kHFPhiSeg
Definition: CaloTools.h:42
static bool isValidIEtaIPhi(int iEta, int iPhi)
Definition: CaloTools.cc:90
static const int kHBHEEnd
Definition: CaloTools.h:39
static const int kNrHBHETowers
Definition: CaloTools.h:47
static const int kHFNrPhi
Definition: CaloTools.h:43
static const int kHFEnd
Definition: CaloTools.h:41
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
static const int kNrTowers
Definition: CaloTools.h:46
static const int kHBHENrPhi
Definition: CaloTools.h:44
static size_t caloTowerHashMax()
Definition: CaloTools.cc:84
size_t l1t::CaloTools::caloTowerHashMax ( )
static

Definition at line 84 of file CaloTools.cc.

Referenced by L1TStage2Layer2Producer::produce().

85 {
86  return kNrTowers;
87 }
static const int kNrTowers
Definition: CaloTools.h:46
l1t::EGamma l1t::CaloTools::egP4Demux ( l1t::EGamma eg)
static

Definition at line 266 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

266  {
267 
268  return l1t::EGamma( p4Demux(&eg),
269  eg.hwPt(),
270  eg.hwEta(),
271  eg.hwPhi(),
272  eg.hwQual(),
273  eg.hwIso() );
274 
275 }
int hwPhi() const
Definition: L1Candidate.cc:79
int hwIso() const
Definition: L1Candidate.cc:84
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
static math::PtEtaPhiMLorentzVector p4Demux(l1t::L1Candidate *)
Definition: CaloTools.cc:256
l1t::EGamma l1t::CaloTools::egP4MP ( l1t::EGamma eg)
static

Definition at line 324 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

324  {
325 
326  return l1t::EGamma( p4MP(&eg),
327  eg.hwPt(),
328  eg.hwEta(),
329  eg.hwPhi(),
330  eg.hwQual(),
331  eg.hwIso() );
332 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:315
int hwPhi() const
Definition: L1Candidate.cc:79
int hwIso() const
Definition: L1Candidate.cc:84
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
l1t::EtSum l1t::CaloTools::etSumP4Demux ( l1t::EtSum etsum)
static

Definition at line 301 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().

301  {
302 
303  return l1t::EtSum( p4Demux(&etsum),
304  etsum.getType(),
305  etsum.hwPt(),
306  etsum.hwEta(),
307  etsum.hwPhi(),
308  etsum.hwQual() );
309 
310 }
int hwPhi() const
Definition: L1Candidate.cc:79
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
static math::PtEtaPhiMLorentzVector p4Demux(l1t::L1Candidate *)
Definition: CaloTools.cc:256
EtSumType getType() const
Definition: EtSum.cc:37
l1t::EtSum l1t::CaloTools::etSumP4MP ( l1t::EtSum etsum)
static

Definition at line 356 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().

356  {
357 
358  return l1t::EtSum( p4MP(&etsum),
359  etsum.getType(),
360  etsum.hwPt(),
361  etsum.hwEta(),
362  etsum.hwPhi(),
363  etsum.hwQual() );
364 
365 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:315
int hwPhi() const
Definition: L1Candidate.cc:79
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
EtSumType getType() const
Definition: EtSum.cc:37
const l1t::CaloCluster & l1t::CaloTools::getCluster ( const std::vector< l1t::CaloCluster > &  clusters,
int  iEta,
int  iPhi 
)
static

Definition at line 51 of file CaloTools.cc.

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

52 {
53  for(size_t clusterNr=0;clusterNr<clusters.size();clusterNr++){
54  if(clusters[clusterNr].hwEta()==iEta && clusters[clusterNr].hwPhi()==iPhi) return clusters[clusterNr];
55  }
56  return nullCluster_;
57 }
static const l1t::CaloCluster nullCluster_
Definition: CaloTools.h:113
const l1t::CaloTower & l1t::CaloTools::getTower ( const std::vector< l1t::CaloTower > &  towers,
int  iEta,
int  iPhi 
)
static

Definition at line 30 of file CaloTools.cc.

References caloTowerHash().

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

31 {
32  size_t towerIndex = CaloTools::caloTowerHash(iEta, iPhi);
33  if(towerIndex<towers.size()){
34  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
35  // std::cout <<"error, tower "<<towers[towerIndex].hwEta()<<" "<<towers[towerIndex].hwPhi()<<" does not match "<<iEta<<" "<<iPhi<<" index "<<towerIndex<<" nr towrs "<<towers.size()<<std::endl;
36  for(size_t towerNr=0;towerNr<towers.size();towerNr++){
37  if(towers[towerNr].hwEta()==iEta && towers[towerNr].hwPhi()==iPhi) return towers[towerNr];
38  }
39  }else return towers[towerIndex];
40 
41  }
42  else{// in case the vector of towers do not contain all the towers (towerIndex can be > towers.size())
43  for(size_t towerNr=0;towerNr<towers.size();towerNr++){
44  if(towers[towerNr].hwEta()==iEta && towers[towerNr].hwPhi()==iPhi) return towers[towerNr];
45  }
46  }
47 
48  return nullTower_;
49 }
static size_t caloTowerHash(int iEta, int iPhi)
Definition: CaloTools.cc:64
static const l1t::CaloTower nullTower_
Definition: CaloTools.h:112
int l1t::CaloTools::gtEta ( int  ieta)
static

Definition at line 236 of file CaloTools.cc.

References eta.

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

236  {
237 
238  double eta = towerEta(ieta);
239  return round ( eta / kGTEtaLSB );
240 
241 }
static float towerEta(int ieta)
Definition: CaloTools.cc:162
static const float kGTEtaLSB
Definition: CaloTools.h:115
int l1t::CaloTools::gtPhi ( int  ieta,
int  iphi 
)
static

Definition at line 243 of file CaloTools.cc.

References M_PI, and phi().

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

243  {
244 
245  double phi = towerPhi(ieta, iphi);
246  if (phi<0) phi = phi + 2*M_PI;
247  return round ( phi / kGTPhiLSB );
248 
249 }
static const float kGTPhiLSB
Definition: CaloTools.h:116
static float towerPhi(int ieta, int iphi)
Definition: CaloTools.cc:170
#define M_PI
bool l1t::CaloTools::insertTower ( std::vector< l1t::CaloTower > &  towers,
const l1t::CaloTower tower 
)
static

Definition at line 19 of file CaloTools.cc.

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

19  {
20  size_t towerIndex = CaloTools::caloTowerHash(tower.hwEta(), tower.hwPhi());
21  if (towers.size() > towerIndex) {
22  towers.at(towerIndex) = tower;
23  return true;
24  }
25  else return false;
26 }
int hwPhi() const
Definition: L1Candidate.cc:79
static size_t caloTowerHash(int iEta, int iPhi)
Definition: CaloTools.cc:64
int hwEta() const
Definition: L1Candidate.cc:74
bool l1t::CaloTools::isValidIEtaIPhi ( int  iEta,
int  iPhi 
)
static

Definition at line 90 of file CaloTools.cc.

References funct::abs().

91 {
92  size_t absIEta = abs(iEta);
93  if(iPhi<=0 || iPhi>kNPhi) return false;
94  if(absIEta==0 || absIEta>kHFEnd) return false;
95  //if(absIEta>kHBHEEnd && iPhi%kHFPhiSeg!=1) return false;
96  return true;
97 
98 }
static const int kNPhi
Definition: CaloTools.h:45
static const int kHFEnd
Definition: CaloTools.h:41
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
l1t::Jet l1t::CaloTools::jetP4Demux ( l1t::Jet jet)
static

Definition at line 290 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

290  {
291 
292  return l1t::Jet( p4Demux(&jet),
293  jet.hwPt(),
294  jet.hwEta(),
295  jet.hwPhi(),
296  jet.hwQual() );
297 
298 }
int hwPhi() const
Definition: L1Candidate.cc:79
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
static math::PtEtaPhiMLorentzVector p4Demux(l1t::L1Candidate *)
Definition: CaloTools.cc:256
l1t::Jet l1t::CaloTools::jetP4MP ( l1t::Jet jet)
static

Definition at line 346 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

346  {
347 
348  return l1t::Jet( p4MP(&jet),
349  jet.hwPt(),
350  jet.hwEta(),
351  jet.hwPhi(),
352  jet.hwQual() );
353 
354 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:315
int hwPhi() const
Definition: L1Candidate.cc:79
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
int l1t::CaloTools::mpEta ( int  ieta)
static

Definition at line 191 of file CaloTools.cc.

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

191  {
192 
193  if (ieta>kHFBegin) return ieta-1;
194  else if (ieta<-1*kHFBegin) return ieta+1;
195  else return ieta;
196 
197 }
static const int kHFBegin
Definition: CaloTools.h:40
math::PtEtaPhiMLorentzVector l1t::CaloTools::p4Demux ( l1t::L1Candidate cand)
static

Definition at line 256 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(), and l1t::stage2::EGammaUnpacker::unpack().

256  {
257 
258  return math::PtEtaPhiMLorentzVector( cand->hwPt() * kGTEtLSB + 1.E-6,
259  cand->hwEta() * kGTEtaLSB,
260  cand->hwPhi() * kGTPhiLSB,
261  0. ) ;
262 
263 }
static const float kGTPhiLSB
Definition: CaloTools.h:116
static const float kGTEtaLSB
Definition: CaloTools.h:115
int hwPhi() const
Definition: L1Candidate.cc:79
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
int hwEta() const
Definition: L1Candidate.cc:74
int hwPt() const
Definition: L1Candidate.cc:69
static const float kGTEtLSB
Definition: CaloTools.h:117
math::PtEtaPhiMLorentzVector l1t::CaloTools::p4MP ( l1t::L1Candidate cand)
static

Definition at line 315 of file CaloTools.cc.

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

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

315  {
316 
317  return math::PtEtaPhiMLorentzVector( cand->hwPt() * 0.5 + 1.E-6,
318  towerEta(cand->hwEta()),
319  towerPhi(cand->hwEta(), cand->hwPhi()),
320  0. ) ;
321 
322 }
static float towerEta(int ieta)
Definition: CaloTools.cc:162
static float towerPhi(int ieta, int iphi)
Definition: CaloTools.cc:170
int hwPhi() const
Definition: L1Candidate.cc:79
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
int hwEta() const
Definition: L1Candidate.cc:74
int hwPt() const
Definition: L1Candidate.cc:69
int l1t::CaloTools::regionEta ( int  ieta)
static

Definition at line 211 of file CaloTools.cc.

References funct::abs().

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

212 {
213 
214  // outside HF
215  if (abs(ieta) > kHFEnd)
216  return (ieta<0 ? 0 : 21);
217 
218  // inside HBHE
219  if (abs(ieta) <= kHFBegin)
220  {
221  if (ieta<0)
222  return 11 - ceil( double (abs(ieta) /4.) );
223  else
224  return ceil( double (abs(ieta) /4.) ) + 10;
225  }
226 
227  // in HF
228  if (ieta<0)
229  return 4 - ceil( double (abs(ieta)-29) /4. );
230  else
231  return ceil( double (abs(ieta)-29) /4. ) + 17;
232 
233 }
static const int kHFBegin
Definition: CaloTools.h:40
static const int kHFEnd
Definition: CaloTools.h:41
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
l1t::Tau l1t::CaloTools::tauP4Demux ( l1t::Tau tau)
static

Definition at line 278 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

278  {
279 
280  return l1t::Tau( p4Demux(&tau),
281  tau.hwPt(),
282  tau.hwEta(),
283  tau.hwPhi(),
284  tau.hwQual(),
285  tau.hwIso() );
286 
287 }
int hwPhi() const
Definition: L1Candidate.cc:79
int hwIso() const
Definition: L1Candidate.cc:84
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
static math::PtEtaPhiMLorentzVector p4Demux(l1t::L1Candidate *)
Definition: CaloTools.cc:256
l1t::Tau l1t::CaloTools::tauP4MP ( l1t::Tau tau)
static

Definition at line 335 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

335  {
336 
337  return l1t::Tau( p4MP(&tau),
338  tau.hwPt(),
339  tau.hwEta(),
340  tau.hwPhi(),
341  tau.hwQual(),
342  tau.hwIso() );
343 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:315
int hwPhi() const
Definition: L1Candidate.cc:79
int hwIso() const
Definition: L1Candidate.cc:84
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
int hwPt() const
Definition: L1Candidate.cc:69
float l1t::CaloTools::towerEta ( int  ieta)
static

Definition at line 162 of file CaloTools.cc.

References eta, and jetcorrextractor::sign().

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

163 {
164  std::pair<float,float> bounds = towerEtaBounds(ieta);
165  float eta = (bounds.second+bounds.first)/2.;
166  float sign = ieta>0 ? 1. : -1.;
167  return sign*eta;
168 }
double sign(double x)
static std::pair< float, float > towerEtaBounds(int ieta)
Definition: CaloTools.cc:152
std::pair< float, float > l1t::CaloTools::towerEtaBounds ( int  ieta)
staticprivate

Definition at line 152 of file CaloTools.cc.

References funct::abs().

153 {
154  if(ieta==0) ieta = 1;
155  if(ieta>kHFEnd) ieta = kHFEnd;
156  if(ieta<(-1*kHFEnd)) ieta = -1*kHFEnd;
157  //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};
158  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};
159  return std::make_pair( towerEtas[abs(ieta)-1],towerEtas[abs(ieta)] );
160 }
static const int kHFEnd
Definition: CaloTools.h:41
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
float l1t::CaloTools::towerEtaSize ( int  ieta)
static

Definition at line 177 of file CaloTools.cc.

References findQualityFiles::size.

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

178 {
179  std::pair<float,float> bounds = towerEtaBounds(ieta);
180  float size = (bounds.second-bounds.first);
181  return size;
182 }
static std::pair< float, float > towerEtaBounds(int ieta)
Definition: CaloTools.cc:152
tuple size
Write out results.
float l1t::CaloTools::towerPhi ( int  ieta,
int  iphi 
)
static

Definition at line 170 of file CaloTools.cc.

References M_PI, and phi().

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

171 {
172  float phi = (float(iphi)-0.5)*towerPhiSize(ieta);
173  if (phi > M_PI) phi = phi - (2*M_PI);
174  return phi;
175 }
#define M_PI
static float towerPhiSize(int ieta)
Definition: CaloTools.cc:184
float l1t::CaloTools::towerPhiSize ( int  ieta)
static

Definition at line 184 of file CaloTools.cc.

References M_PI.

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

185 {
186  return 2.*M_PI/kNPhi;
187 }
static const int kNPhi
Definition: CaloTools.h:45
#define M_PI

Member Data Documentation

const int l1t::CaloTools::cos_coeff = {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, 512, 587, 658, 723, 784, 838, 886, 927, 961, 988, 1007, 1019}
static
const float l1t::CaloTools::kGTEtaLSB = 0.0435
staticprivate

Definition at line 115 of file CaloTools.h.

const float l1t::CaloTools::kGTEtLSB = 0.5
staticprivate

Definition at line 117 of file CaloTools.h.

const float l1t::CaloTools::kGTPhiLSB = 0.0435
staticprivate

Definition at line 116 of file CaloTools.h.

const int l1t::CaloTools::kHBHEEnd =28
static

Definition at line 39 of file CaloTools.h.

const int l1t::CaloTools::kHBHENrPhi =72
static
const int l1t::CaloTools::kHFBegin =29
static
const int l1t::CaloTools::kHFEnd =41
static
const int l1t::CaloTools::kHFNrPhi =72/kHFPhiSeg
static

Definition at line 43 of file CaloTools.h.

const int l1t::CaloTools::kHFPhiSeg =1
static

Definition at line 42 of file CaloTools.h.

const int l1t::CaloTools::kNPhi =72
static

Definition at line 45 of file CaloTools.h.

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

Definition at line 47 of file CaloTools.h.

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

Definition at line 46 of file CaloTools.h.

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

Definition at line 113 of file CaloTools.h.

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

Definition at line 112 of file CaloTools.h.

const int l1t::CaloTools::sin_coeff = {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}
static