ATLAS Offline Software
Loading...
Searching...
No Matches
MergedElectronDetailsDecorator.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3*/
4
5
6#include "DerivationFrameworkHiggs/MergedElectronDetailsDecorator.h"
7#include "xAODTracking/TrackParticleContainer.h"
8#include "xAODTruth/TruthParticleContainer.h"
9#include "xAODTruth/TruthVertexContainer.h"
10#include "xAODCaloEvent/CaloClusterContainer.h"
11
12#include "FourMomUtils/P4Helpers.h"
13#include "xAODTruth/xAODTruthHelpers.h"
14#include "ElectronPhotonSelectorTools/ElectronSelectorHelpers.h"
15#include "AthContainers/Decorator.h"
16
17#include "GaudiKernel/EventContext.h"
18
19#include <vector>
20
21namespace DerivationFramework {
22
23 StatusCode MergedElectronDetailsDecorator::initialize() {
24
25 ATH_MSG_INFO("Initialize " );
26
27 if (m_emExtrapolationTool.retrieve().isFailure()) {
28 ATH_MSG_FATAL("Failed to retrieve tool: " << m_emExtrapolationTool);
29 return StatusCode::FAILURE;
30 }
31 ATH_MSG_INFO("Retrieved tool: " << m_emExtrapolationTool);
32
33 if (m_VertexFitter.retrieve().isFailure()) {
34 ATH_MSG_FATAL("Failed to retrieve tool: " << m_VertexFitter);
35 return StatusCode::FAILURE;
36 }
37 ATH_MSG_INFO("Retrieved tool: " << m_VertexFitter);
38
39 if (m_V0Tools.retrieve().isFailure()) {
40 ATH_MSG_FATAL("Failed to retrieve tool: " << m_V0Tools);
41 return StatusCode::FAILURE;
42 }
43 ATH_MSG_INFO("Retrieved tool: " << m_V0Tools);
44
45 ATH_CHECK(m_eventInfoKey.initialize());
46
47 ATH_CHECK(m_electronKey.initialize());
48 ATH_CHECK(m_caloMgrKey.initialize());
49 return StatusCode::SUCCESS;
50 }
51
52
53 StatusCode MergedElectronDetailsDecorator::addBranches(const EventContext& ctx) const{
54
55
56 // Retrieve the xAOD event info
57 SG::ReadHandle<xAOD::ElectronContainer> electrons (m_electronKey, ctx);
58
59 SG::ReadHandle<xAOD::EventInfo> eventInfo (m_eventInfoKey, ctx);
60
61 bool isMC = false;
62 if( eventInfo->eventType(xAOD::EventInfo::IS_SIMULATION) )
63 isMC = true;
64
65
66 for( const auto *const el : *electrons ){
67 fillTrackDetails(el,isMC);
68 fillVertexDetails(el);
69 fillClusterDetails(el);
70 }
71
72 return StatusCode::SUCCESS;
73
74 } // addBranches
75
76 void DerivationFramework::MergedElectronDetailsDecorator::fillTruthDetails( std::vector<float>& trkMatchTrk,
77 const xAOD::TrackParticle* tp,
78 const xAOD::CaloCluster* cluster) const{
79
80 trkMatchTrk.clear();
81 const xAOD::TruthParticle* truthPart = xAOD::TruthHelpers::getTruthParticle( *tp );
82
83 if( !truthPart || !truthPart->hasProdVtx() )
84 return;
85
86 const auto *prodVtx = truthPart->prodVtx();
87
88 Amg::Vector3D pos( prodVtx->x(), prodVtx->y(), prodVtx->z() );
89 Amg::Vector3D mom( truthPart->px(), truthPart->py(), truthPart->pz() );
90
91 Trk::CurvilinearParameters truthTP( pos, mom, truthPart->charge() );
92
93 const EventContext& ctx = Gaudi::Hive::currentContext();
94 float etaAtCalo, phiAtCalo;
95 if( m_emExtrapolationTool->getEtaPhiAtCalo ( ctx,
96 &truthTP,
97 &etaAtCalo,
98 &phiAtCalo) ){
99 trkMatchTrk.push_back( cluster->etaBE(2) - etaAtCalo );
100 trkMatchTrk.push_back( P4Helpers::deltaPhi(cluster->phiBE(2), phiAtCalo) );
101 }
102
103 return;
104
105 }
106
107 void DerivationFramework::MergedElectronDetailsDecorator::fillMatchDetails( std::vector<float>& trkMatchTrk,
108 const xAOD::TrackParticle* tp,
109 const xAOD::CaloCluster* cluster) const {
110
111 SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{ m_caloMgrKey };
112 const CaloDetDescrManager* caloDDMgr = *caloMgrHandle;
113
114 //Reset values
115 for(unsigned int i(0); i<trkMatchTrk.size(); ++i ){
116 trkMatchTrk[i] = -999;
117 }
118
119 std::array<double,4> eta = { -999.0, -999.0, -999.0, -999.0 };
120 std::array<double,4> phi = { -999.0, -999.0, -999.0, -999.0 };
121 std::array<double,4> deltaEta = { -999.0, -999.0, -999.0, -999.0 };
122 std::array<double,4> deltaPhi = { -999.0, -999.0, -999.0, -999.0 };
123
124 const EventContext& ctx = Gaudi::Hive::currentContext();
125 std::pair<std::vector<CaloSampling::CaloSample>,
126 std::vector<std::unique_ptr<Trk::Surface>>>
127 layersAndSurfaces =
128 m_emExtrapolationTool->getClusterLayerSurfaces(*cluster, *caloDDMgr);
129
130 if (m_emExtrapolationTool
131 ->getMatchAtCalo(ctx,
132 *cluster,
133 *tp,
134 layersAndSurfaces.first,
135 layersAndSurfaces.second,
136 eta,
137 phi,
138 deltaEta,
139 deltaPhi,
140 IEMExtrapolationTools::fromPerigee)
141 .isSuccess()) // Perigee
142 {
143 trkMatchTrk[0] = deltaEta[1];
144 trkMatchTrk[1] = deltaEta[2];
145 trkMatchTrk[2] = deltaPhi[1];
146 trkMatchTrk[3] = deltaPhi[2];
147 }
148
149 if (m_emExtrapolationTool
150 ->getMatchAtCalo(ctx,
151 *cluster,
152 *tp,
153 layersAndSurfaces.first,
154 layersAndSurfaces.second,
155 eta,
156 phi,
157 deltaEta,
158 deltaPhi,
159 IEMExtrapolationTools::fromLastMeasurement)
160 .isSuccess()) // Last Measurement
161 {
162 trkMatchTrk[4] = deltaEta[1];
163 trkMatchTrk[5] = deltaEta[2];
164 trkMatchTrk[6] = deltaPhi[1];
165 trkMatchTrk[7] = deltaPhi[2];
166 }
167
168 if(m_emExtrapolationTool->getMatchAtCalo ( ctx,
169 *cluster,
170 *tp,
171 layersAndSurfaces.first,
172 layersAndSurfaces.second,
173 eta,
174 phi,
175 deltaEta,
176 deltaPhi,
177 IEMExtrapolationTools::fromPerigeeRescaled).isSuccess()) //Last Measurement
178 {
179 trkMatchTrk[8] = deltaPhi[2];
180 }
181 }
182
183
184 int DerivationFramework::MergedElectronDetailsDecorator::nSiHits( const xAOD::TrackParticle * tp )
185 {
186 uint8_t dummy(-1);
187 int nPix = tp->summaryValue( dummy, xAOD::numberOfPixelHits )? dummy : 0;
188 int nPix_DS = tp->summaryValue( dummy, xAOD::numberOfPixelDeadSensors )? dummy : 0;
189 int nSCT = tp->summaryValue( dummy, xAOD::numberOfSCTHits )? dummy : 0;
190 int nSCT_DS = tp->summaryValue( dummy, xAOD::numberOfSCTDeadSensors )? dummy : 0;
191 return nPix + nPix_DS + nSCT + nSCT_DS;
192 }
193
194
195 void DerivationFramework::MergedElectronDetailsDecorator::fillTrackDetails(const xAOD::Electron* el,bool isMC) const
196 {
197 std::vector<float> trkMatchTrkP_dEta1(el->nTrackParticles(),-999);
198 std::vector<float> trkMatchTrkP_dEta2(el->nTrackParticles(),-999);
199 std::vector<float> trkMatchTrkP_dPhi1(el->nTrackParticles(),-999);
200 std::vector<float> trkMatchTrkP_dPhi2(el->nTrackParticles(),-999);
201 std::vector<float> trkMatchTrkLM_dEta1(el->nTrackParticles(),-999);
202 std::vector<float> trkMatchTrkLM_dEta2(el->nTrackParticles(),-999);
203 std::vector<float> trkMatchTrkLM_dPhi1(el->nTrackParticles(),-999);
204 std::vector<float> trkMatchTrkLM_dPhi2(el->nTrackParticles(),-999);
205 std::vector<float> trkMatchTrkR_dPhi2(el->nTrackParticles(),-999);
206
207 std::vector<float> trueMatch_dEta2(el->nTrackParticles(),-999);
208 std::vector<float> trueMatch_dPhi2(el->nTrackParticles(),-999);
209
210 const auto *caloCluster = el->caloCluster();
211 if( caloCluster && caloCluster->pt() > m_minET ){
212 std::vector<float> trkMatch(9,-999);
213 for( unsigned int i(0); i < el->nTrackParticles(); ++i ){
214 const auto *trackParticle = el->trackParticle( i );
215 if(trackParticle){
216 fillMatchDetails( trkMatch, trackParticle, caloCluster);
217 trkMatchTrkP_dEta1[i] = trkMatch[0];
218 trkMatchTrkP_dEta2[i] = trkMatch[1];
219 trkMatchTrkP_dPhi1[i] = trkMatch[2];
220 trkMatchTrkP_dPhi2[i] = trkMatch[3];
221 trkMatchTrkLM_dEta1[i] = trkMatch[4];
222 trkMatchTrkLM_dEta2[i] = trkMatch[5];
223 trkMatchTrkLM_dPhi1[i] = trkMatch[6];
224 trkMatchTrkLM_dPhi2[i] = trkMatch[7];
225 trkMatchTrkR_dPhi2[i] = trkMatch[8];
226
227 if(isMC){
228 fillTruthDetails( trkMatch, trackParticle, caloCluster );
229 if( trkMatch.size() == 2 ){
230 trueMatch_dEta2[i] = trkMatch[0];
231 trueMatch_dPhi2[i] = trkMatch[1];
232 }
233 }
234 }
235 }
236 }
237
238 static const SG::Decorator<std::vector<float> > TrackMatchingP_dEta1Dec("TrackMatchingP_dEta1");
239 static const SG::Decorator<std::vector<float> > TrackMatchingP_dEta2Dec("TrackMatchingP_dEta2");
240 static const SG::Decorator<std::vector<float> > TrackMatchingP_dPhi1Dec("TrackMatchingP_dPhi1");
241 static const SG::Decorator<std::vector<float> > TrackMatchingP_dPhi2Dec("TrackMatchingP_dPhi2");
242 static const SG::Decorator<std::vector<float> > TrackMatchingLM_dEta1Dec("TrackMatchingLM_dEta1");
243 static const SG::Decorator<std::vector<float> > TrackMatchingLM_dEta2Dec("TrackMatchingLM_dEta2");
244 static const SG::Decorator<std::vector<float> > TrackMatchingLM_dPhi1Dec("TrackMatchingLM_dPhi1");
245 static const SG::Decorator<std::vector<float> > TrackMatchingLM_dPhi2Dec("TrackMatchingLM_dPhi2");
246 static const SG::Decorator<std::vector<float> > TrackMatchingR_dPhi2Dec("TrackMatchingR_dPhi2");
247 static const SG::Decorator<std::vector<float> > TrackMatchingTrue_dEta2Dec("TrackMatchingTrue_dEta2");
248 static const SG::Decorator<std::vector<float> > TrackMatchingTrue_dPhi2Dec("TrackMatchingTrue_dPhi2");
249 TrackMatchingP_dEta1Dec(*el) = trkMatchTrkP_dEta1;
250 TrackMatchingP_dEta2Dec(*el) = trkMatchTrkP_dEta2;
251 TrackMatchingP_dPhi1Dec(*el) = trkMatchTrkP_dPhi1;
252 TrackMatchingP_dPhi2Dec(*el) = trkMatchTrkP_dPhi2;
253 TrackMatchingLM_dEta1Dec(*el) = trkMatchTrkLM_dEta1;
254 TrackMatchingLM_dEta2Dec(*el) = trkMatchTrkLM_dEta2;
255 TrackMatchingLM_dPhi1Dec(*el) = trkMatchTrkLM_dPhi1;
256 TrackMatchingLM_dPhi2Dec(*el) = trkMatchTrkLM_dPhi2;
257 TrackMatchingR_dPhi2Dec(*el) = trkMatchTrkR_dPhi2;
258 TrackMatchingTrue_dEta2Dec(*el) = trueMatch_dEta2;
259 TrackMatchingTrue_dPhi2Dec(*el) = trueMatch_dPhi2;
260 }
261
262 void DerivationFramework::MergedElectronDetailsDecorator::fillClusterDetails(const xAOD::Electron* el)
263 {
264 const auto *caloCluster = el->caloCluster();
265
266 std::vector<float> subCluster_E;
267 std::vector<float> subCluster_dEta;
268 std::vector<float> subCluster_dPhi;
269
270 static const SG::Accessor<std::vector<ElementLink<xAOD::CaloClusterContainer> > > clusterLinksAcc("constituentClusterLinks");
271 if(caloCluster && clusterLinksAcc.isAvailable(*caloCluster) ){
272 std::vector<ElementLink<xAOD::CaloClusterContainer> > clusterLinks = clusterLinksAcc(*caloCluster);
273 for( const auto& link : clusterLinks){
274 if( link.isValid() ){
275 subCluster_E.push_back( (*link)->e() );
276 subCluster_dEta.push_back( caloCluster->eta() - (*link)->eta() );
277 float dphi = P4Helpers::deltaPhi( caloCluster->phi(), (*link)->phi() );
278 subCluster_dPhi.push_back(dphi);
279 }
280 }
281 }
282 static const SG::Decorator<std::vector<float> > SubCluster_EDec("SubCluster_E");
283 static const SG::Decorator<std::vector<float> > SubCluster_dEtaDec("SubCluster_dEta");
284 static const SG::Decorator<std::vector<float> > SubCluster_dPhiDec("SubCluster_dPhi");
285 SubCluster_EDec(*el) = subCluster_E;
286 SubCluster_dEtaDec(*el) = subCluster_dEta;
287 SubCluster_dPhiDec(*el) = subCluster_dPhi;
288 }
289
290
291 void DerivationFramework::MergedElectronDetailsDecorator::fillVertexDetails(const xAOD::Electron* el) const
292 {
293 const EventContext& ctx = Gaudi::Hive::currentContext();
294 const auto *caloCluster = el->caloCluster();
295
296 float vtxR = -999;
297 float vtxZ = -999;
298 float vtxM = -999;
299 float vtxP = -999;
300 int vtxTrkParticleIndex1 = -999;
301 int vtxTrkParticleIndex2 = -999;
302 float vtxTrkParticle1_dEta2 = -999;
303 float vtxTrkParticle1_dPhi2 = -999;
304 float vtxTrkParticle2_dEta2 = -999;
305 float vtxTrkParticle2_dPhi2 = -999;
306 float vtxRerr = -999;
307 float vtxZerr = -999;
308 float vtxMerr = -999;
309 float vtxPerr = -999;
310 float vtxChi2 = -999;
311 int vtxNdof = -999;
312 float vtxdEta = -999;
313 float vtxdPhi = -999;
314 float vtxE = -999;
315 float vtxPhi = -999;
316 float vtxEta = -999;
317
318 if( caloCluster && caloCluster->pt() > m_minET ){
319 const xAOD::TrackParticle* trk1 = nullptr;
320 const xAOD::TrackParticle* trk2 = nullptr;
321 int trkIndex1 = -999;
322 int trkIndex2 = -999;
323 for( unsigned int i(0); i < el->nTrackParticles(); ++i ){
324 const auto *trackParticle = el->trackParticle( i );
325 if( nSiHits(trackParticle) >= 7 && ElectronSelectorHelpers::passBLayerRequirement(*trackParticle) ) {
326 if ( trk1 == nullptr ){
327 trk1 = trackParticle;
328 trkIndex1 = i;
329 } else if( trk2 == nullptr ){
330 if( trk1->charge() == trackParticle->charge() )
331 continue;
332 trk2 = trackParticle;
333 trkIndex2 = i;
334 break;
335 }
336 }
337 }
338 if( trk2 != nullptr){
339 // fit vertex
340
341 // Take the initial guess as the first tracks global position
342 const Trk::Perigee& perigeeParameters = trk1->perigeeParameters();
343 Amg::Vector3D startingPoint = perigeeParameters.position();
344
345 std::vector<const xAOD::TrackParticle*> trksToFit;
346 trksToFit.push_back( trk1 );
347 trksToFit.push_back( trk2 );
348 vtxTrkParticleIndex1 = trkIndex1;
349 vtxTrkParticleIndex2 = trkIndex2;
350
351 std::unique_ptr<xAOD::Vertex> myVertex( m_VertexFitter->fit( trksToFit, startingPoint ) );
352
353 if(myVertex){
354
355 vtxR = m_V0Tools->rxy( myVertex.get() );
356 vtxRerr = m_V0Tools->rxyError( myVertex.get() );
357 vtxZ = myVertex->z();
358 vtxZerr = myVertex->covariancePosition()(2,2) > 0 ? sqrt( myVertex->covariancePosition()(2,2) ): myVertex->covariancePosition()(2,2);
359
360 std::vector<double> masses = { 0.511e-3, 0.511e-3 };
361 vtxM = m_V0Tools->invariantMass( myVertex.get(), masses);
362 vtxMerr = m_V0Tools->invariantMassError( myVertex.get(), masses);
363
364 vtxP = m_V0Tools->pT( myVertex.get() );
365 vtxPerr = m_V0Tools->pTError( myVertex.get() );
366
367 vtxChi2 = myVertex->chiSquared();
368 vtxNdof = myVertex->numberDoF();
369
370 m_emExtrapolationTool->getEtaPhiAtCalo( ctx, myVertex.get(), &vtxdEta, &vtxdPhi);
371
372 vtxdPhi = P4Helpers::deltaPhi( vtxdPhi, caloCluster->phiBE(2) );
373 vtxdEta -= caloCluster->etaBE(2);
374
375 xAOD::TrackParticle::FourMom_t vertex4P = m_V0Tools->V04Momentum(myVertex.get(), ParticleConstants::electronMassInMeV);
376 vtxE = vertex4P.E();
377 vtxPhi = vertex4P.Phi();
378 vtxEta = vertex4P.Eta();
379 double momentumScaleFactor = caloCluster->e() / vtxE;
380
381 const auto& perigeeParameters = trk1->perigeeParameters();
382 Amg::Vector3D pos = perigeeParameters.position();
383 Amg::Vector3D mom = perigeeParameters.momentum();
384 mom *= momentumScaleFactor;
385 Trk::CurvilinearParameters scaledTrk( pos, mom, trk1->charge() );
386
387 float etaAtCalo, phiAtCalo;
388 if( m_emExtrapolationTool->getEtaPhiAtCalo (ctx,
389 &scaledTrk,
390 &etaAtCalo,
391 &phiAtCalo) ){
392 vtxTrkParticle1_dEta2 = caloCluster->etaBE(2) - etaAtCalo;
393 vtxTrkParticle1_dPhi2 = P4Helpers::deltaPhi(caloCluster->phiBE(2), phiAtCalo) ;
394 }
395
396 const auto& perigeeParameters2 = trk2->perigeeParameters();
397 pos = perigeeParameters2.position();
398 mom = perigeeParameters2.momentum();
399 mom *= momentumScaleFactor;
400 scaledTrk = Trk::CurvilinearParameters( pos, mom, trk2->charge() );
401
402 if( m_emExtrapolationTool->getEtaPhiAtCalo (ctx,
403 &scaledTrk,
404 &etaAtCalo,
405 &phiAtCalo) ){
406 vtxTrkParticle2_dEta2 = caloCluster->etaBE(2) - etaAtCalo;
407 vtxTrkParticle2_dPhi2 = P4Helpers::deltaPhi(caloCluster->phiBE(2), phiAtCalo) ;
408 }
409
410 } else {
411 vtxChi2 = 0;
412 vtxNdof = 0;
413 }
414 }
415 }
416
417 static const SG::Decorator<float> vtxRDec("vtxR");
418 static const SG::Decorator<float> vtxRerrDec("vtxRerr");
419 static const SG::Decorator<float> vtxZDec("vtxZ");
420 static const SG::Decorator<float> vtxZerrDec("vtxZerr");
421 static const SG::Decorator<float> vtxMDec("vtxM");
422 static const SG::Decorator<float> vtxMerrDec("vtxMerr");
423 static const SG::Decorator<float> vtxPtDec("vtxPt");
424 static const SG::Decorator<float> vtxPterrDec("vtxPterr");
425 static const SG::Decorator<float> vtxEDec("vtxE");
426 static const SG::Decorator<float> vtxEtaDec("vtxEta");
427 static const SG::Decorator<float> vtxPhiDec("vtxPhi");
428 static const SG::Decorator<int> vtxTrkParticleIndex1Dec("vtxTrkParticleIndex1");
429 static const SG::Decorator<int> vtxTrkParticleIndex2Dec("vtxTrkParticleIndex2");
430 static const SG::Decorator<float> vtxTrkParticle1_dPhi2Dec("vtxTrkParticle1_dPhi2");
431 static const SG::Decorator<float> vtxTrkParticle1_dEta2Dec("vtxTrkParticle1_dEta2");
432 static const SG::Decorator<float> vtxTrkParticle2_dPhi2Dec("vtxTrkParticle2_dPhi2");
433 static const SG::Decorator<float> vtxTrkParticle2_dEta2Dec("vtxTrkParticle2_dEta2");
434 static const SG::Decorator<float> vtxChi2Dec("vtxChi2");
435 static const SG::Decorator<int> vtxNdofDec("vtxNdof");
436 static const SG::Decorator<float> vtxdEtaDec("vtxdEta");
437 static const SG::Decorator<float> vtxdPhiDec("vtxdPhi");
438
439 vtxRDec(*el) = vtxR;
440 vtxRerrDec(*el) = vtxRerr;
441 vtxZDec(*el) = vtxZ;
442 vtxZerrDec(*el) = vtxZerr;
443 vtxMDec(*el) = vtxM;
444 vtxMerrDec(*el) = vtxMerr;
445 vtxPtDec(*el) = vtxP;
446 vtxPterrDec(*el)= vtxPerr;
447 vtxEDec(*el) = vtxE;
448 vtxEtaDec(*el) = vtxEta;
449 vtxPhiDec(*el) = vtxPhi;
450 vtxTrkParticleIndex1Dec(*el) = vtxTrkParticleIndex1;
451 vtxTrkParticleIndex2Dec(*el) = vtxTrkParticleIndex2;
452 vtxTrkParticle1_dPhi2Dec(*el) = vtxTrkParticle1_dPhi2;
453 vtxTrkParticle1_dEta2Dec(*el) = vtxTrkParticle1_dEta2;
454 vtxTrkParticle2_dPhi2Dec(*el) = vtxTrkParticle2_dPhi2;
455 vtxTrkParticle2_dEta2Dec(*el) = vtxTrkParticle2_dEta2;
456 vtxChi2Dec(*el) = vtxChi2;
457 vtxNdofDec(*el) = vtxNdof;
458 vtxdEtaDec(*el) = vtxdEta;
459 vtxdPhiDec(*el) = vtxdPhi;
460 }
461
462
463
464
465} // namespace
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
deltaPhi
Scalar deltaPhi(const MatrixBase< Derived > &vec) const
Definition AmgMatrixBasePlugin.h:112
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition AthMsgStreamMacros.h:34
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
Decorator.h
Helper class to provide type-safe access to aux data.
CaloDetDescrManager
This class provides the client interface for accessing the detector description information common to...
Definition CaloDetDescrManager.h:469
DerivationFramework::MergedElectronDetailsDecorator::m_caloMgrKey
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
Definition MergedElectronDetailsDecorator.h:64
DerivationFramework::MergedElectronDetailsDecorator::fillTrackDetails
void fillTrackDetails(const xAOD::Electron *el, bool isMC) const
Definition MergedElectronDetailsDecorator.cxx:195
DerivationFramework::MergedElectronDetailsDecorator::fillTruthDetails
void fillTruthDetails(std::vector< float > &trkMatchTrk, const xAOD::TrackParticle *tp, const xAOD::CaloCluster *cluster) const
Definition MergedElectronDetailsDecorator.cxx:76
DerivationFramework::MergedElectronDetailsDecorator::m_emExtrapolationTool
PublicToolHandle< IEMExtrapolationTools > m_emExtrapolationTool
Definition MergedElectronDetailsDecorator.h:50
DerivationFramework::MergedElectronDetailsDecorator::addBranches
virtual StatusCode addBranches(const EventContext &ctx) const override final
Definition MergedElectronDetailsDecorator.cxx:53
DerivationFramework::MergedElectronDetailsDecorator::m_electronKey
SG::ReadHandleKey< xAOD::ElectronContainer > m_electronKey
Definition MergedElectronDetailsDecorator.h:60
DerivationFramework::MergedElectronDetailsDecorator::fillMatchDetails
void fillMatchDetails(std::vector< float > &trkMatchTrk, const xAOD::TrackParticle *tp, const xAOD::CaloCluster *cluster) const
Definition MergedElectronDetailsDecorator.cxx:107
IEMExtrapolationTools::fromPerigee
@ fromPerigee
from the perigee of TrackParticle
Definition IEMExtrapolationTools.h:37
IEMExtrapolationTools::fromLastMeasurement
@ fromLastMeasurement
from the last measurement of TrackParticle
Definition IEMExtrapolationTools.h:35
IEMExtrapolationTools::fromPerigeeRescaled
@ fromPerigeeRescaled
from the perigee of TrackParticle recaled by Ecluster
Definition IEMExtrapolationTools.h:39
SG::Accessor
Helper class to provide type-safe access to aux data.
Definition Control/AthContainers/AthContainers/Accessor.h:68
SG::ConstAccessor::isAvailable
bool isAvailable(const ELT &e) const
Test to see if this variable exists in the store.
SG::Decorator
Helper class to provide type-safe access to aux data.
Definition Decorator.h:59
xAOD::CaloCluster_v1::phiBE
float phiBE(const unsigned layer) const
Get the phi in one layer of the EM Calo.
Definition CaloCluster_v1.cxx:634
xAOD::CaloCluster_v1::etaBE
float etaBE(const unsigned layer) const
Get the eta in one layer of the EM Calo.
Definition CaloCluster_v1.cxx:628
xAOD::EventInfo_v1::IS_SIMULATION
@ IS_SIMULATION
true: simulation, false: data
Definition EventInfo_v1.h:151
xAOD::TrackParticle_v1::perigeeParameters
const Trk::Perigee & perigeeParameters() const
Returns the Trk::MeasuredPerigee track parameters.
Definition TrackParticle_v1.cxx:415
xAOD::TrackParticle_v1::FourMom_t
IParticle::FourMom_t FourMom_t
Definition of the 4-momentum type.
Definition TrackParticle_v1.h:73
xAOD::TrackParticle_v1::charge
float charge() const
Returns the charge.
Definition TrackParticle_v1.cxx:143
xAOD::TruthParticle_v1::hasProdVtx
bool hasProdVtx() const
Check for a production vertex on this particle.
Definition TruthParticle_v1.cxx:69
xAOD::TruthParticle_v1::px
float px() const
The x component of the particle's momentum.
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition TruthParticle_v1.cxx:75
xAOD::TruthParticle_v1::py
float py() const
The y component of the particle's momentum.
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
xAOD::TruthParticle_v1::pz
float pz() const
The z component of the particle's momentum.
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition GeoPrimitives.h:47
DerivationFramework
THE reconstruction tool.
Definition ParticleSortingAlg.h:24
ElectronSelectorHelpers::passBLayerRequirement
bool passBLayerRequirement(const xAOD::TrackParticle &tp)
return true if effective number of BL hits + outliers is at least one
Definition ElectronSelectorHelpers.cxx:59
P4Helpers::deltaPhi
double deltaPhi(double phiA, double phiB)
delta Phi in range [-pi,pi[
Definition P4Helpers.h:34
ParticleConstants::PDG2011::electronMassInMeV
constexpr double electronMassInMeV
the mass of the electron (in MeV)
Definition ParticleConstants.h:26
Trk::Perigee
ParametersT< TrackParametersDim, Charged, PerigeeSurface > Perigee
Definition Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:33
Trk::CurvilinearParameters
CurvilinearParametersT< TrackParametersDim, Charged, PlaneSurface > CurvilinearParameters
Definition Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:29
xAOD::TruthHelpers::getTruthParticle
const xAOD::TruthParticle * getTruthParticle(const xAOD::IParticle &p)
Return the truthParticle associated to the given IParticle (if any)
Definition xAODTruthHelpers.cxx:25
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
xAOD::numberOfSCTDeadSensors
@ numberOfSCTDeadSensors
number of dead SCT sensors crossed [unit8_t].
Definition TrackingPrimitives.h:274
xAOD::numberOfSCTHits
@ numberOfSCTHits
number of hits in SCT [unit8_t].
Definition TrackingPrimitives.h:269
xAOD::numberOfPixelHits
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
Definition TrackingPrimitives.h:260
xAOD::numberOfPixelDeadSensors
@ numberOfPixelDeadSensors
number of dead pixel sensors crossed [unit8_t].
Definition TrackingPrimitives.h:267
xAOD::Electron
Electron_v1 Electron
Definition of the current "egamma version".
Definition Event/xAOD/xAODEgamma/xAODEgamma/Electron.h:17
Generated on for ATLAS Offline Software by doxygen 1.14.0