Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Classes | Functions
FEHelpers Namespace Reference

Classes

class  FillNeutralFlowElements
 

Functions

TLorentzVector getVertexCorrectedFourVec (const xAOD::FlowElement &fe, const xAOD::Vertex &vertexToCorrectTo)
 
TLorentzVector getVertexCorrectedFourVec (const xAOD::FlowElement &fe, const TVector3 &vertexToCorrectTo)
 
void vertexCorrectTheFourVector (const xAOD::FlowElement &fe, const TVector3 &vertexToCorrectTo, TLorentzVector &theFourVector)
 
bool getClusterMoment (const xAOD::FlowElement &fe, xAOD::CaloCluster::MomentType momentType, float &value)
 
std::vector< float > getEnergiesPerSampling (const xAOD::FlowElement &fe)
 
std::string getClusterMomentName (xAOD::CaloCluster::MomentType momentType)
 
xAOD::Type::ObjectType signalToXAODType (const xAOD::FlowElement &fe)
 

Function Documentation

◆ getClusterMoment()

bool FEHelpers::getClusterMoment ( const xAOD::FlowElement &  fe,
xAOD::CaloCluster::MomentType  momentType,
float &  value 
)

Definition at line 51 of file FEHelpers.cxx.

51  {
52 
53  // First try to retrieve directly from the cluster
54  if(fe.nOtherObjects() >= 1){
55  const xAOD::CaloCluster* cluster = dynamic_cast<const xAOD::CaloCluster*>(fe.otherObject(0));
56  if(cluster != nullptr){
57  double tmpValue(-99.);
58  if(cluster->retrieveMoment(momentType, tmpValue)){
59  value = static_cast<float>(tmpValue);
60  return true;
61  }
62  }
63  }
64 
65  // Cluster or moment from it was unavailable, try to retrieve from the FE itself instead
66  std::string momentName = getClusterMomentName(momentType);
67  if(!momentName.empty()){
68  const SG::AuxElement::ConstAccessor<float> acc(momentName);
69  if(acc.isAvailable(fe)){
70  value = acc(fe);
71  return true;
72  }
73  }
74 
75  return false;
76  }

◆ getClusterMomentName()

std::string FEHelpers::getClusterMomentName ( xAOD::CaloCluster::MomentType  momentType)

Definition at line 147 of file FEHelpers.cxx.

147  {
148 
149  switch(momentType){
150  case xAOD::CaloCluster::ENG_FRAC_CORE: return "ENG_FRAC_CORE";
151  case xAOD::CaloCluster::FIRST_ENG_DENS: return "FIRST_ENG_DENS";
152  case xAOD::CaloCluster::CENTER_LAMBDA: return "CENTER_LAMBDA";
153  case xAOD::CaloCluster::SECOND_R: return "SECOND_R";
154  case xAOD::CaloCluster::DELTA_ALPHA: return "DELTA_ALPHA";
155  case xAOD::CaloCluster::LATERAL: return "LATERAL";
156  case xAOD::CaloCluster::LONGITUDINAL: return "LONGITUDINAL";
157  case xAOD::CaloCluster::SECOND_LAMBDA: return "SECOND_LAMBDA";
158  case xAOD::CaloCluster::ISOLATION: return "ISOLATION";
159  case xAOD::CaloCluster::ENG_FRAC_MAX: return "ENG_FRAC_MAX";
160  case xAOD::CaloCluster::ENG_BAD_CELLS: return "ENG_BAD_CELLS";
161  case xAOD::CaloCluster::N_BAD_CELLS: return "N_BAD_CELLS";
162  case xAOD::CaloCluster::BADLARQ_FRAC: return "BADLARQ_FRAC";
163  case xAOD::CaloCluster::ENG_POS: return "ENG_POS";
164  case xAOD::CaloCluster::SIGNIFICANCE: return "SIGNIFICANCE";
165  case xAOD::CaloCluster::CELL_SIGNIFICANCE: return "CELL_SIGNIFICANCE";
166  case xAOD::CaloCluster::CELL_SIG_SAMPLING: return "CELL_SIG_SAMPLING";
167  case xAOD::CaloCluster::AVG_LAR_Q: return "AVG_LAR_Q";
168  case xAOD::CaloCluster::AVG_TILE_Q: return "AVG_TILE_Q";
169  case xAOD::CaloCluster::EM_PROBABILITY: return "EM_PROBABILITY";
170  case xAOD::CaloCluster::ENG_CALIB_TOT: return "ENG_CALIB_TOT";
171  case xAOD::CaloCluster::ENG_CALIB_FRAC_EM: return "ENG_CALIB_FRAC_EM";
172  case xAOD::CaloCluster::ENG_CALIB_FRAC_HAD: return "ENG_CALIB_FRAC_HAD";
173  case xAOD::CaloCluster::ENG_CALIB_FRAC_REST: return "ENG_CALIB_FRAC_REST";
174  default: return "";
175  }
176  }

◆ getEnergiesPerSampling()

std::vector< float > FEHelpers::getEnergiesPerSampling ( const xAOD::FlowElement &  fe)

Definition at line 78 of file FEHelpers.cxx.

78  {
79 
80  const static SG::AuxElement::ConstAccessor< std::vector<float> > accEPerSampling("ePerSampling");
81  if(accEPerSampling.isAvailable(fe)) return accEPerSampling(fe);
82 
83  // This wasn't stored as a vector, try to retrieve the individual elements
84  const static SG::AuxElement::ConstAccessor<float> accPreSamplerB("LAYERENERGY_PreSamplerB");
85  const static SG::AuxElement::ConstAccessor<float> accEMB1("LAYERENERGY_EMB1");
86  const static SG::AuxElement::ConstAccessor<float> accEMB2("LAYERENERGY_EMB2");
87  const static SG::AuxElement::ConstAccessor<float> accEMB3("LAYERENERGY_EMB3");
88  const static SG::AuxElement::ConstAccessor<float> accPreSamplerE("LAYERENERGY_PreSamplerE");
89  const static SG::AuxElement::ConstAccessor<float> accEME1("LAYERENERGY_EME1");
90  const static SG::AuxElement::ConstAccessor<float> accEME2("LAYERENERGY_EME2");
91  const static SG::AuxElement::ConstAccessor<float> accEME3("LAYERENERGY_EME3");
92  const static SG::AuxElement::ConstAccessor<float> accHEC0("LAYERENERGY_HEC0");
93  const static SG::AuxElement::ConstAccessor<float> accHEC1("LAYERENERGY_HEC1");
94  const static SG::AuxElement::ConstAccessor<float> accHEC2("LAYERENERGY_HEC2");
95  const static SG::AuxElement::ConstAccessor<float> accHEC3("LAYERENERGY_HEC3");
96  const static SG::AuxElement::ConstAccessor<float> accTileBar0("LAYERENERGY_TileBar0");
97  const static SG::AuxElement::ConstAccessor<float> accTileBar1("LAYERENERGY_TileBar1");
98  const static SG::AuxElement::ConstAccessor<float> accTileBar2("LAYERENERGY_TileBar2");
99  const static SG::AuxElement::ConstAccessor<float> accTileGap1("LAYERENERGY_TileGap1");
100  const static SG::AuxElement::ConstAccessor<float> accTileGap2("LAYERENERGY_TileGap2");
101  const static SG::AuxElement::ConstAccessor<float> accTileGap3("LAYERENERGY_TileGap3");
102  const static SG::AuxElement::ConstAccessor<float> accTileExt0("LAYERENERGY_TileExt0");
103  const static SG::AuxElement::ConstAccessor<float> accTileExt1("LAYERENERGY_TileExt1");
104  const static SG::AuxElement::ConstAccessor<float> accTileExt2("LAYERENERGY_TileExt2");
105  const static SG::AuxElement::ConstAccessor<float> accFCAL0("LAYERENERGY_FCAL0");
106  const static SG::AuxElement::ConstAccessor<float> accFCAL1("LAYERENERGY_FCAL1");
107  const static SG::AuxElement::ConstAccessor<float> accFCAL2("LAYERENERGY_FCAL2");
108  const static SG::AuxElement::ConstAccessor<float> accMINIFCAL0("LAYERENERGY_MINIFCAL0");
109  const static SG::AuxElement::ConstAccessor<float> accMINIFCAL1("LAYERENERGY_MINIFCAL1");
110  const static SG::AuxElement::ConstAccessor<float> accMINIFCAL2("LAYERENERGY_MINIFCAL2");
111  const static SG::AuxElement::ConstAccessor<float> accMINIFCAL3("LAYERENERGY_MINIFCAL3");
112 
113  // Build the vector up from the individual samplings
114  std::vector<float> result(28, 0);
115  result[0] = accPreSamplerB(fe);
116  result[1] = accEMB1(fe);
117  result[2] = accEMB2(fe);
118  result[3] = accEMB3(fe);
119  result[4] = accPreSamplerE(fe);
120  result[5] = accEME1(fe);
121  result[6] = accEME2(fe);
122  result[7] = accEME3(fe);
123  result[8] = accHEC0(fe);
124  result[9] = accHEC1(fe);
125  result[10] = accHEC2(fe);
126  result[11] = accHEC3(fe);
127  result[12] = accTileBar0(fe);
128  result[13] = accTileBar1(fe);
129  result[14] = accTileBar2(fe);
130  result[15] = accTileGap1(fe);
131  result[16] = accTileGap2(fe);
132  result[17] = accTileGap3(fe);
133  result[18] = accTileExt0(fe);
134  result[19] = accTileExt1(fe);
135  result[20] = accTileExt2(fe);
136  result[21] = accFCAL0(fe);
137  result[22] = accFCAL1(fe);
138  result[23] = accFCAL2(fe);
139  result[24] = accMINIFCAL0(fe);
140  result[25] = accMINIFCAL1(fe);
141  result[26] = accMINIFCAL2(fe);
142  result[27] = accMINIFCAL3(fe);
143 
144  return result;
145  }

◆ getVertexCorrectedFourVec() [1/2]

TLorentzVector FEHelpers::getVertexCorrectedFourVec ( const xAOD::FlowElement &  fe,
const TVector3 &  vertexToCorrectTo 
)

Definition at line 18 of file FEHelpers.cxx.

18  {
19 
20  TLorentzVector theNewVector(0.0,0.0,0.0,0.0);
21  theNewVector.SetPtEtaPhiM(fe.pt(), fe.eta(), fe.phi(), fe.m());
22 
23  vertexCorrectTheFourVector(fe, vertexToCorrectTo, theNewVector);
24  return theNewVector;
25  }

◆ getVertexCorrectedFourVec() [2/2]

TLorentzVector FEHelpers::getVertexCorrectedFourVec ( const xAOD::FlowElement &  fe,
const xAOD::Vertex &  vertexToCorrectTo 
)

Definition at line 13 of file FEHelpers.cxx.

13  {
14  TVector3 theVertexVector(vertexToCorrectTo.x(), vertexToCorrectTo.y(), vertexToCorrectTo.z());
15  return getVertexCorrectedFourVec(fe, theVertexVector);
16  }

◆ signalToXAODType()

xAOD::Type::ObjectType FEHelpers::signalToXAODType ( const xAOD::FlowElement &  fe)

Definition at line 178 of file FEHelpers.cxx.

178  {
179  xAOD::FlowElement::signal_t st = fe.signalType();
180  if(st & xAOD::FlowElement::PFlow) return xAOD::Type::ParticleFlow;
181  if(st & xAOD::FlowElement::CaloCluster) return xAOD::Type::CaloCluster;
182  if(st & xAOD::FlowElement::Track) return xAOD::Type::TrackParticle;
183  if(st & xAOD::FlowElement::Muon) return xAOD::Type::Muon;
184  if(st & xAOD::FlowElement::TCC) return xAOD::Type::TrackCaloCluster;
185  if(st & xAOD::FlowElement::UFO) return xAOD::Type::TrackCaloCluster; // UFO doesn't have its own xAOD type
186  return xAOD::Type::Other;
187  }

◆ vertexCorrectTheFourVector()

void FEHelpers::vertexCorrectTheFourVector ( const xAOD::FlowElement &  fe,
const TVector3 &  vertexToCorrectTo,
TLorentzVector &  theFourVector 
)

Definition at line 27 of file FEHelpers.cxx.

27  {
28 
29  float clusterEta = theFourVector.Eta();
30 
31  const static SG::AuxElement::ConstAccessor<float> accCenterMag("CENTER_MAG");
32  float centerMag = accCenterMag(fe);
33 
34  float radius = centerMag/cosh(clusterEta);
35 
36  float EtaVertexCorrection = 0.0, PhiVertexCorrection = 0.0;
37  float clusterPhi = theFourVector.Phi();
38 
39  if (radius > 1.0 && centerMag > 1e-3){
40  EtaVertexCorrection = (-vertexToCorrectTo.Z()/cosh(clusterEta) + (vertexToCorrectTo.X()*cos(clusterPhi) + vertexToCorrectTo.Y()*sin(clusterPhi))*tanh(clusterEta))/radius;
41  PhiVertexCorrection = (vertexToCorrectTo.X()*sin(clusterPhi) - vertexToCorrectTo.Y()*cos(clusterPhi))/radius;
42  }
43 
44  float etaVertexCorrected = clusterEta + EtaVertexCorrection;
45  float p = std::sqrt(theFourVector.E()*theFourVector.E()-theFourVector.M()*theFourVector.M());
46  float ptVertexCorrected = p/cosh(etaVertexCorrected);
47  theFourVector.SetPtEtaPhiM(ptVertexCorrected, etaVertexCorrected, clusterPhi + PhiVertexCorrection, theFourVector.M());
48 
49  }
xAOD::CaloCluster_v1::SECOND_R
@ SECOND_R
Second Moment in .
Definition: CaloCluster_v1.h:126
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
FEHelpers::getVertexCorrectedFourVec
TLorentzVector getVertexCorrectedFourVec(const xAOD::FlowElement &fe, const TVector3 &vertexToCorrectTo)
Definition: FEHelpers.cxx:18
xAOD::Vertex_v1::x
float x() const
Returns the x position.
xAOD::FlowElement_v1::Track
@ Track
Definition: FlowElement_v1.h:43
xAOD::CaloCluster_v1::ENG_FRAC_CORE
@ ENG_FRAC_CORE
Energy fraction of the sum of the hottest cells in each sampling.
Definition: CaloCluster_v1.h:145
get_generator_info.result
result
Definition: get_generator_info.py:21
xAOD::CaloCluster_v1::EM_PROBABILITY
@ EM_PROBABILITY
Classification probability to be em-like.
Definition: CaloCluster_v1.h:176
xAOD::FlowElement_v1::m
virtual double m() const override
The invariant mass of the particle.
xAOD::FlowElement_v1::UFO
@ UFO
Definition: FlowElement_v1.h:53
xAOD::FlowElement_v1::CaloCluster
@ CaloCluster
Definition: FlowElement_v1.h:42
xAOD::FlowElement_v1::nOtherObjects
std::size_t nOtherObjects() const
Definition: FlowElement_v1.cxx:192
xAOD::CaloCluster_v1::AVG_LAR_Q
@ AVG_LAR_Q
Sum(E_cell_LAr^2 Q_cell_LAr)/Sum(E_cell_LAr^2)
Definition: CaloCluster_v1.h:166
FEHelpers::getClusterMomentName
std::string getClusterMomentName(xAOD::CaloCluster::MomentType momentType)
Definition: FEHelpers.cxx:147
xAOD::CaloCluster_v1::CENTER_LAMBDA
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
Definition: CaloCluster_v1.h:139
athena.value
value
Definition: athena.py:124
xAOD::FlowElement_v1::phi
virtual double phi() const override
The azimuthal angle ( ) of the particle.
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
xAOD::FlowElement_v1::pt
virtual double pt() const override
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
xAOD::FlowElement_v1::TCC
@ TCC
Definition: FlowElement_v1.h:50
xAOD::CaloCluster_v1::ENG_FRAC_MAX
@ ENG_FRAC_MAX
Energy fraction of hottest cell.
Definition: CaloCluster_v1.h:143
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19
xAOD::FlowElement_v1::signal_t
unsigned long signal_t
Definition: FlowElement_v1.h:30
xAOD::CaloCluster_v1::SECOND_LAMBDA
@ SECOND_LAMBDA
Second Moment in .
Definition: CaloCluster_v1.h:127
xAOD::TrackCaloCluster
TrackCaloCluster_v1 TrackCaloCluster
Reference the current persistent version:
Definition: TrackCaloCluster.h:12
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition: Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::FlowElement_v1::PFlow
@ PFlow
Definition: FlowElement_v1.h:45
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:62
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
xAOD::FlowElement_v1::Muon
@ Muon
Definition: FlowElement_v1.h:44
columnar::ContainerId::cluster
@ cluster
xAOD::CaloCluster_v1::CELL_SIG_SAMPLING
@ CELL_SIG_SAMPLING
CaloSample of the cell with the largest |E|/sig.
Definition: CaloCluster_v1.h:164
xAOD::CaloCluster_v1::ENG_CALIB_TOT
@ ENG_CALIB_TOT
Calibration Hit energy inside the cluster.
Definition: CaloCluster_v1.h:198
FEHelpers::vertexCorrectTheFourVector
void vertexCorrectTheFourVector(const xAOD::FlowElement &fe, const TVector3 &vertexToCorrectTo, TLorentzVector &theFourVector)
Definition: FEHelpers.cxx:27
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_EM
@ ENG_CALIB_FRAC_EM
Calibration Hit energy inside the cluster caused by e/gamma/pi0.
Definition: CaloCluster_v1.h:251
xAOD::CaloCluster_v1::ISOLATION
@ ISOLATION
Energy weighted fraction of non-clustered perimeter cells.
Definition: CaloCluster_v1.h:149
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_REST
@ ENG_CALIB_FRAC_REST
Calibration Hit energy inside the cluster caused by other particles.
Definition: CaloCluster_v1.h:256
xAOD::FlowElement_v1::signalType
signal_t signalType() const
AthenaPoolTestRead.acc
acc
Definition: AthenaPoolTestRead.py:16
xAODType::ParticleFlow
@ ParticleFlow
The object is a particle-flow object.
Definition: ObjectType.h:41
xAOD::Vertex_v1::z
float z() const
Returns the z position.
xAOD::CaloCluster_v1::FIRST_ENG_DENS
@ FIRST_ENG_DENS
First Moment in E/V.
Definition: CaloCluster_v1.h:146
xAOD::CaloCluster_v1::ENG_BAD_CELLS
@ ENG_BAD_CELLS
Total em-scale energy of bad cells in this cluster.
Definition: CaloCluster_v1.h:151
xAOD::CaloCluster_v1::AVG_TILE_Q
@ AVG_TILE_Q
Sum(E_cell_Tile^2 Q_cell_Tile)/Sum(E_cell_Tile^2)
Definition: CaloCluster_v1.h:168
ParticleGun_SamplingFraction.radius
radius
Definition: ParticleGun_SamplingFraction.py:96
Muon
struct TBPatternUnitContext Muon
xAOD::FlowElement_v1::eta
virtual double eta() const override
The pseudorapidity ( ) of the particle.
xAOD::CaloCluster_v1::SIGNIFICANCE
@ SIGNIFICANCE
Cluster significance.
Definition: CaloCluster_v1.h:160
xAOD::FlowElement_v1::otherObject
const xAOD::IParticle * otherObject(std::size_t i) const
Definition: FlowElement_v1.cxx:196
xAOD::Vertex_v1::y
float y() const
Returns the y position.
xAOD::CaloCluster_v1::LATERAL
@ LATERAL
Normalized lateral moment.
Definition: CaloCluster_v1.h:140
xAOD::CaloCluster_v1::DELTA_ALPHA
@ DELTA_ALPHA
Angular shower axis deviation ( ) from IP-to-Center.
Definition: CaloCluster_v1.h:133
xAOD::CaloCluster_v1::ENG_CALIB_FRAC_HAD
@ ENG_CALIB_FRAC_HAD
Calibration Hit energy inside the cluster caused by charged pi+ and pi-.
Definition: CaloCluster_v1.h:254
xAOD::CaloCluster_v1::CELL_SIGNIFICANCE
@ CELL_SIGNIFICANCE
Cell significance = E/sig of the cell with the largest |E|/sig.
Definition: CaloCluster_v1.h:162
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
xAOD::CaloCluster_v1::LONGITUDINAL
@ LONGITUDINAL
Normalized longitudinal moment.
Definition: CaloCluster_v1.h:141
xAOD::CaloCluster_v1::N_BAD_CELLS
@ N_BAD_CELLS
number of bad cells
Definition: CaloCluster_v1.h:152
xAOD::CaloCluster_v1::ENG_POS
@ ENG_POS
Total positive Energy of this cluster.
Definition: CaloCluster_v1.h:159
xAOD::CaloCluster_v1::BADLARQ_FRAC
@ BADLARQ_FRAC
Energy fraction of LAr cells with quality larger than a given cut.
Definition: CaloCluster_v1.h:158
Generated on Thu Apr 17 2025 21:53:31 for ATLAS Offline Software by   doxygen 1.8.18