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 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
 

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

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

99 {
100 
101  return calHwEtSum(iEta,iPhi,towers,localEtaMin,localEtaMax,localPhiMin,localPhiMax,kHFEnd,etMode);
102 }
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:96
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 104 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().

107 {
108  int hwEtSum=0;
109  for(int etaNr=localEtaMin;etaNr<=localEtaMax;etaNr++){
110  for(int phiNr=localPhiMin;phiNr<=localPhiMax;phiNr++){
111 
112  int towerIEta = l1t::CaloStage2Nav::offsetIEta(iEta,etaNr);
113  int towerIPhi = l1t::CaloStage2Nav::offsetIPhi(iPhi,phiNr);
114  if(abs(towerIEta)<=iEtaAbsMax){
115  const l1t::CaloTower& tower = getTower(towers,towerIEta,towerIPhi);
116  if(etMode==ECAL) hwEtSum+=tower.hwEtEm();
117  else if(etMode==HCAL) hwEtSum+=tower.hwEtHad();
118  else if(etMode==CALO) hwEtSum+=tower.hwPt();
119  }
120  }
121  }
122  return hwEtSum;
123 }
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:26
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 126 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().

127 {
128  size_t nrTowers=0;
129  l1t::CaloStage2Nav nav(iEtaMin,iPhiMin);
130  while(nav.currIEta()<=iEtaMax){
131  bool finishPhi = false;
132  while(!finishPhi){
133  const l1t::CaloTower& tower = l1t::CaloTools::getTower(towers,nav.currIEta(),nav.currIPhi());
134  int towerHwEt =0;
135  if(etMode==ECAL) towerHwEt+=tower.hwEtEm();
136  else if(etMode==HCAL) towerHwEt+=tower.hwEtHad();
137  else if(etMode==CALO) towerHwEt+=tower.hwPt();
138  if(towerHwEt>=minHwEt && towerHwEt<=maxHwEt) nrTowers++;
139  finishPhi = (nav.currIPhi() == iPhiMax);
140  nav.north();
141  }
142  nav.east();
143  nav.resetIPhi();
144  }
145  return nrTowers;
146 }
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:26
int hwEtHad() const
Definition: CaloTower.cc:69
int l1t::CaloTools::caloEta ( int  ietaMP)
static

Definition at line 197 of file CaloTools.cc.

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

197  {
198 
199  if (mpEta>kHFBegin) return mpEta+1;
200  else if (mpEta<-1*kHFBegin) return mpEta-1;
201  else return mpEta;
202 
203 }
static int mpEta(int ieta)
Definition: CaloTools.cc:187
static const int kHFBegin
Definition: CaloTools.h:40
size_t l1t::CaloTools::caloTowerHash ( int  iEta,
int  iPhi 
)
static

Definition at line 60 of file CaloTools.cc.

References funct::abs().

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

61 {
62 
63  if(!isValidIEtaIPhi(iEta,iPhi)) return caloTowerHashMax();
64  else{
65  const int absIEta = abs(iEta);
66  if(absIEta>kHFEnd) return kNrTowers;
67  else if(absIEta<=kHBHEEnd){ //HBHE
68  int iEtaNoZero=iEta;
69  if(iEta>0) iEtaNoZero--;
70  return (iEtaNoZero+kHBHEEnd)*kHBHENrPhi+iPhi-1;
71  }else{ //HF
72  int iEtaIndex = iEta+kHFEnd; //iEta=-32 is 0
73  if(iEta>0) iEtaIndex= iEta-kHBHEEnd+(kHFEnd-kHBHEEnd)-1; //but iEta=29 is 4
74  return iEtaIndex*kHFNrPhi+iPhi/kHFPhiSeg + kNrHBHETowers;
75  }
76  }
77 }
static const int kHFPhiSeg
Definition: CaloTools.h:42
static bool isValidIEtaIPhi(int iEta, int iPhi)
Definition: CaloTools.cc:86
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:80
size_t l1t::CaloTools::caloTowerHashMax ( )
static

Definition at line 80 of file CaloTools.cc.

Referenced by L1TStage2Layer2Producer::produce().

81 {
82  return kNrTowers;
83 }
static const int kNrTowers
Definition: CaloTools.h:46
l1t::EGamma l1t::CaloTools::egP4Demux ( l1t::EGamma eg)
static

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

262  {
263 
264  return l1t::EGamma( p4Demux(&eg),
265  eg.hwPt(),
266  eg.hwEta(),
267  eg.hwPhi(),
268  eg.hwQual(),
269  eg.hwIso() );
270 
271 }
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:252
l1t::EGamma l1t::CaloTools::egP4MP ( l1t::EGamma eg)
static

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

320  {
321 
322  return l1t::EGamma( p4MP(&eg),
323  eg.hwPt(),
324  eg.hwEta(),
325  eg.hwPhi(),
326  eg.hwQual(),
327  eg.hwIso() );
328 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:311
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 297 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().

297  {
298 
299  return l1t::EtSum( p4Demux(&etsum),
300  etsum.getType(),
301  etsum.hwPt(),
302  etsum.hwEta(),
303  etsum.hwPhi(),
304  etsum.hwQual() );
305 
306 }
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:252
EtSumType getType() const
Definition: EtSum.cc:37
l1t::EtSum l1t::CaloTools::etSumP4MP ( l1t::EtSum etsum)
static

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

352  {
353 
354  return l1t::EtSum( p4MP(&etsum),
355  etsum.getType(),
356  etsum.hwPt(),
357  etsum.hwEta(),
358  etsum.hwPhi(),
359  etsum.hwQual() );
360 
361 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:311
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 47 of file CaloTools.cc.

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

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

Definition at line 26 of file CaloTools.cc.

References caloTowerHash().

Referenced by L1TStage2InputPatternWriter::analyze(), 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().

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

Definition at line 232 of file CaloTools.cc.

References eta.

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

232  {
233 
234  double eta = towerEta(ieta);
235  return round ( eta / kGTEtaLSB );
236 
237 }
static float towerEta(int ieta)
Definition: CaloTools.cc:158
static const float kGTEtaLSB
Definition: CaloTools.h:113
int l1t::CaloTools::gtPhi ( int  ieta,
int  iphi 
)
static

Definition at line 239 of file CaloTools.cc.

References M_PI, and phi().

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

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

Definition at line 15 of file CaloTools.cc.

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

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

Definition at line 86 of file CaloTools.cc.

References funct::abs().

87 {
88  size_t absIEta = abs(iEta);
89  if(iPhi<=0 || iPhi>kNPhi) return false;
90  if(absIEta==0 || absIEta>kHFEnd) return false;
91  //if(absIEta>kHBHEEnd && iPhi%kHFPhiSeg!=1) return false;
92  return true;
93 
94 }
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 286 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

286  {
287 
288  return l1t::Jet( p4Demux(&jet),
289  jet.hwPt(),
290  jet.hwEta(),
291  jet.hwPhi(),
292  jet.hwQual() );
293 
294 }
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:252
l1t::Jet l1t::CaloTools::jetP4MP ( l1t::Jet jet)
static

Definition at line 342 of file CaloTools.cc.

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

Referenced by L1TStage2Layer2Producer::produce().

342  {
343 
344  return l1t::Jet( p4MP(&jet),
345  jet.hwPt(),
346  jet.hwEta(),
347  jet.hwPhi(),
348  jet.hwQual() );
349 
350 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:311
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 187 of file CaloTools.cc.

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

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

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

252  {
253 
254  return math::PtEtaPhiMLorentzVector( cand->hwPt() * kGTEtLSB + 1.E-6,
255  cand->hwEta() * kGTEtaLSB,
256  cand->hwPhi() * kGTPhiLSB,
257  0. ) ;
258 
259 }
static const float kGTPhiLSB
Definition: CaloTools.h:114
static const float kGTEtaLSB
Definition: CaloTools.h:113
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:115
math::PtEtaPhiMLorentzVector l1t::CaloTools::p4MP ( l1t::L1Candidate cand)
static

Definition at line 311 of file CaloTools.cc.

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

Referenced by l1t::stage2::MPUnpacker::unpack().

311  {
312 
313  return math::PtEtaPhiMLorentzVector( cand->hwPt() * 0.5 + 1.E-6,
314  towerEta(cand->hwEta()),
315  towerPhi(cand->hwEta(), cand->hwPhi()),
316  0. ) ;
317 
318 }
static float towerEta(int ieta)
Definition: CaloTools.cc:158
static float towerPhi(int ieta, int iphi)
Definition: CaloTools.cc:166
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 207 of file CaloTools.cc.

References funct::abs().

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

208 {
209 
210  // outside HF
211  if (abs(ieta) > kHFEnd)
212  return (ieta<0 ? 0 : 21);
213 
214  // inside HBHE
215  if (abs(ieta) <= kHFBegin)
216  {
217  if (ieta<0)
218  return 11 - ceil( double (abs(ieta) /4.) );
219  else
220  return ceil( double (abs(ieta) /4.) ) + 10;
221  }
222 
223  // in HF
224  if (ieta<0)
225  return 4 - ceil( double (abs(ieta)-29) /4. );
226  else
227  return ceil( double (abs(ieta)-29) /4. ) + 17;
228 
229 }
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 274 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().

274  {
275 
276  return l1t::Tau( p4Demux(&tau),
277  tau.hwPt(),
278  tau.hwEta(),
279  tau.hwPhi(),
280  tau.hwQual(),
281  tau.hwIso() );
282 
283 }
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:252
l1t::Tau l1t::CaloTools::tauP4MP ( l1t::Tau tau)
static

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

331  {
332 
333  return l1t::Tau( p4MP(&tau),
334  tau.hwPt(),
335  tau.hwEta(),
336  tau.hwPhi(),
337  tau.hwQual(),
338  tau.hwIso() );
339 }
static math::PtEtaPhiMLorentzVector p4MP(l1t::L1Candidate *)
Definition: CaloTools.cc:311
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 158 of file CaloTools.cc.

References eta, and jetcorrextractor::sign().

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

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

Definition at line 148 of file CaloTools.cc.

References funct::abs().

149 {
150  if(ieta==0) ieta = 1;
151  if(ieta>kHFEnd) ieta = kHFEnd;
152  if(ieta<(-1*kHFEnd)) ieta = -1*kHFEnd;
153  //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};
154  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};
155  return std::make_pair( towerEtas[abs(ieta)-1],towerEtas[abs(ieta)] );
156 }
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 173 of file CaloTools.cc.

References findQualityFiles::size.

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

174 {
175  std::pair<float,float> bounds = towerEtaBounds(ieta);
176  float size = (bounds.second-bounds.first);
177  return size;
178 }
static std::pair< float, float > towerEtaBounds(int ieta)
Definition: CaloTools.cc:148
tuple size
Write out results.
float l1t::CaloTools::towerPhi ( int  ieta,
int  iphi 
)
static

Definition at line 166 of file CaloTools.cc.

References M_PI, and phi().

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

167 {
168  float phi = (float(iphi)-0.5)*towerPhiSize(ieta);
169  if (phi > M_PI) phi = phi - (2*M_PI);
170  return phi;
171 }
#define M_PI
static float towerPhiSize(int ieta)
Definition: CaloTools.cc:180
float l1t::CaloTools::towerPhiSize ( int  ieta)
static

Definition at line 180 of file CaloTools.cc.

References M_PI.

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

181 {
182  return 2.*M_PI/kNPhi;
183 }
static const int kNPhi
Definition: CaloTools.h:45
#define M_PI

Member Data Documentation

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

Definition at line 113 of file CaloTools.h.

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

Definition at line 115 of file CaloTools.h.

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

Definition at line 114 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

Definition at line 40 of file CaloTools.h.

Referenced by l1t::stage2::CaloTowerPacker::pack().

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 111 of file CaloTools.h.

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

Definition at line 110 of file CaloTools.h.