ATLAS Offline Software
Loading...
Searching...
No Matches
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 }
xAOD::CaloCluster_v1::retrieveMoment
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
Definition CaloCluster_v1.cxx:662
xAOD::FlowElement_v1::nOtherObjects
std::size_t nOtherObjects() const
Definition FlowElement_v1.cxx:192
xAOD::FlowElement_v1::otherObject
const xAOD::IParticle * otherObject(std::size_t i) const
Definition FlowElement_v1.cxx:196
FEHelpers::getClusterMomentName
std::string getClusterMomentName(xAOD::CaloCluster::MomentType momentType)
Definition FEHelpers.cxx:147
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19

◆ 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 }
xAOD::CaloCluster_v1::DELTA_ALPHA
@ DELTA_ALPHA
Angular shower axis deviation ( ) from IP-to-Center.
Definition CaloCluster_v1.h:133
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
xAOD::CaloCluster_v1::SECOND_LAMBDA
@ SECOND_LAMBDA
Second Moment in .
Definition CaloCluster_v1.h:127
xAOD::CaloCluster_v1::CELL_SIGNIFICANCE
@ CELL_SIGNIFICANCE
Cell significance = E/sig of the cell with the largest |E|/sig.
Definition CaloCluster_v1.h:162
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::EM_PROBABILITY
@ EM_PROBABILITY
Classification probability to be em-like.
Definition CaloCluster_v1.h:176
xAOD::CaloCluster_v1::LATERAL
@ LATERAL
Normalized lateral moment.
Definition CaloCluster_v1.h:140
xAOD::CaloCluster_v1::N_BAD_CELLS
@ N_BAD_CELLS
number of bad cells
Definition CaloCluster_v1.h:152
xAOD::CaloCluster_v1::LONGITUDINAL
@ LONGITUDINAL
Normalized longitudinal moment.
Definition CaloCluster_v1.h:141
xAOD::CaloCluster_v1::ENG_FRAC_MAX
@ ENG_FRAC_MAX
Energy fraction of hottest cell.
Definition CaloCluster_v1.h:143
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
xAOD::CaloCluster_v1::SECOND_R
@ SECOND_R
Second Moment in .
Definition CaloCluster_v1.h:126
xAOD::CaloCluster_v1::FIRST_ENG_DENS
@ FIRST_ENG_DENS
First Moment in E/V.
Definition CaloCluster_v1.h:146
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
xAOD::CaloCluster_v1::CENTER_LAMBDA
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
Definition CaloCluster_v1.h:139
xAOD::CaloCluster_v1::SIGNIFICANCE
@ SIGNIFICANCE
Cluster significance.
Definition CaloCluster_v1.h:160
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::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::ENG_BAD_CELLS
@ ENG_BAD_CELLS
Total em-scale energy of bad cells in this cluster.
Definition CaloCluster_v1.h:151
xAOD::CaloCluster_v1::ENG_CALIB_TOT
@ ENG_CALIB_TOT
Calibration Hit energy inside the cluster.
Definition CaloCluster_v1.h:198
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::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
xAOD::CaloCluster_v1::ISOLATION
@ ISOLATION
Energy weighted fraction of non-clustered perimeter cells.
Definition CaloCluster_v1.h:149

◆ 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 }
SG::AuxElement::ConstAccessor
SG::ConstAccessor< T, ALLOC > ConstAccessor
Definition AuxElement.h:569
result

◆ 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 }
xAOD::FlowElement_v1::pt
virtual double pt() const override
xAOD::FlowElement_v1::m
virtual double m() const override
The invariant mass of the particle.
xAOD::FlowElement_v1::phi
virtual double phi() const override
The azimuthal angle ( ) of the particle.
xAOD::FlowElement_v1::eta
virtual double eta() const override
The pseudorapidity ( ) of the particle.
FEHelpers::vertexCorrectTheFourVector
void vertexCorrectTheFourVector(const xAOD::FlowElement &fe, const TVector3 &vertexToCorrectTo, TLorentzVector &theFourVector)
Definition FEHelpers.cxx:27

◆ 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 }
xAOD::Vertex_v1::z
float z() const
Returns the z position.
xAOD::Vertex_v1::y
float y() const
Returns the y position.
xAOD::Vertex_v1::x
float x() const
Returns the x position.
FEHelpers::getVertexCorrectedFourVec
TLorentzVector getVertexCorrectedFourVec(const xAOD::FlowElement &fe, const xAOD::Vertex &vertexToCorrectTo)
Definition FEHelpers.cxx:13

◆ 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 }
xAOD::FlowElement_v1::signalType
signal_t signalType() const
xAOD::FlowElement_v1::CaloCluster
@ CaloCluster
Definition FlowElement_v1.h:42
xAOD::FlowElement_v1::PFlow
@ PFlow
Definition FlowElement_v1.h:45
xAOD::FlowElement_v1::Muon
@ Muon
Definition FlowElement_v1.h:44
xAOD::FlowElement_v1::UFO
@ UFO
Definition FlowElement_v1.h:53
xAOD::FlowElement_v1::Track
@ Track
Definition FlowElement_v1.h:43
xAOD::FlowElement_v1::TCC
@ TCC
Definition FlowElement_v1.h:50
xAOD::FlowElement_v1::signal_t
unsigned long signal_t
Definition FlowElement_v1.h:30
xAODType::TrackParticle
@ TrackParticle
The object is a charged track particle.
Definition ObjectType.h:43
xAODType::ParticleFlow
@ ParticleFlow
The object is a particle-flow object.
Definition ObjectType.h:41
xAODType::Other
@ Other
An object not falling into any of the other categories.
Definition ObjectType.h:34
xAODType::CaloCluster
@ CaloCluster
The object is a calorimeter cluster.
Definition ObjectType.h:39
xAODType::Muon
@ Muon
The object is a muon.
Definition ObjectType.h:48
xAODType::TrackCaloCluster
@ TrackCaloCluster
The object is a track-calo-cluster.
Definition ObjectType.h:51

◆ 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 }
Generated on for ATLAS Offline Software by doxygen 1.14.0