ATLAS Offline Software
Loading...
Searching...
No Matches
JiveXML::xAODTrackParticleRetriever Class Reference

Retrieves the TrackParticle objects (InDetTrackParticles etc.) More...

#include <xAODTrackParticleRetriever.h>

Inheritance diagram for JiveXML::xAODTrackParticleRetriever:
Collaboration diagram for JiveXML::xAODTrackParticleRetriever:

Public Member Functions

 xAODTrackParticleRetriever (const std::string &type, const std::string &name, const IInterface *parent)
 Standard Constructor.
virtual StatusCode retrieve (ToolHandle< IFormatTool > &FormatTool)
 For each TrackParticle collections retrieve basic parameters.
virtual StatusCode initialize ()
const DataMap getData (const xAOD::TrackParticleContainer *)
 Puts the variables into a DataMap.
virtual std::string dataTypeName () const
 Return the name of the data type that is generated by this retriever.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Static Public Member Functions

static const InterfaceID & interfaceID ()
 Return the interface identifier.

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

const std::string m_typeName = "Track"
 The data type that is generated by this retriever.
SG::ReadHandleKeyArray< xAOD::TrackParticleContainerm_keys {this,"TrackParticleCollections", {"InDetTrackParticles","InDetLargeD0TrackParticles","CombinedMuonTrackParticles","GSFTrackParticles"}, "Name of the Track Particle containers that will be written out to the JiveXML file"}
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Retrieves the TrackParticle objects (InDetTrackParticles etc.)

  • Properties
    • TrackParticleCollections
  • Retrieved Data
    • Usual four-vectors: phi, eta, et etc.
    • Associations for clusters and tracks via ElementLink: key/index scheme

Definition at line 29 of file xAODTrackParticleRetriever.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ xAODTrackParticleRetriever()

JiveXML::xAODTrackParticleRetriever::xAODTrackParticleRetriever ( const std::string & type,
const std::string & name,
const IInterface * parent )

Standard Constructor.

This is the standard AthAlgTool constructor.

Parameters
typeAlgTool type name
nameAlgTool instance name
parentAlgTools parent owning this tool

Definition at line 18 of file xAODTrackParticleRetriever.cxx.

18 :
19 AthAlgTool(type,name,parent){}
AthAlgTool()
Default constructor:

Member Function Documentation

◆ dataTypeName()

virtual std::string JiveXML::xAODTrackParticleRetriever::dataTypeName ( ) const
inlinevirtual

Return the name of the data type that is generated by this retriever.

Implements JiveXML::IDataRetriever.

Definition at line 42 of file xAODTrackParticleRetriever.h.

42{ return m_typeName; };
const std::string m_typeName
The data type that is generated by this retriever.

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ getData()

const DataMap JiveXML::xAODTrackParticleRetriever::getData ( const xAOD::TrackParticleContainer * TrackParticleCont)

Puts the variables into a DataMap.

Retrieve basic parameters, mainly four-vectors, for each collection.

Also association with clusters and tracks (ElementLink).

Definition at line 59 of file xAODTrackParticleRetriever.cxx.

59 {
60
61 ATH_MSG_DEBUG("in getData()");
62
64
65 DataVect d0; d0.reserve(TrackParticleCont->size());
66 DataVect z0; z0.reserve(TrackParticleCont->size());
67 DataVect pt; pt.reserve(TrackParticleCont->size());
68 DataVect phi0; phi0.reserve(TrackParticleCont->size());
69 DataVect cotTheta; cotTheta.reserve(TrackParticleCont->size());
70 DataVect label; label.reserve(TrackParticleCont->size());
71 DataVect nBLayerHits; nBLayerHits.reserve(TrackParticleCont->size());
72 DataVect nPixHits; nPixHits.reserve(TrackParticleCont->size());
73 DataVect nSCTHits; nSCTHits.reserve(TrackParticleCont->size());
74 DataVect nTRTHits; nTRTHits.reserve(TrackParticleCont->size());
75
76 xAOD::TrackParticleContainer::const_iterator TrackParticleItr = TrackParticleCont->begin();
77 xAOD::TrackParticleContainer::const_iterator TrackParticleItrE = TrackParticleCont->end();
78
79 int counter = 0;
80 float charge = 0.;
81 float myQOverP = 0.;
82 double countHits = 0.;
83 std::string labelStr = "unknownHits";
84
85 for (; TrackParticleItr != TrackParticleItrE; ++TrackParticleItr) {
86
87 ATH_MSG_VERBOSE(" TrackParticle #" << counter++ << " : d0 = " << (*TrackParticleItr)->d0() << ", z0 = "
88 << (*TrackParticleItr)->z0() << ", pt[GeV] = " << (*TrackParticleItr)->pt()*(1./CLHEP::GeV)
89 << ", phi = " << (*TrackParticleItr)->phi()
90 << ", qOverP = " << (*TrackParticleItr)->qOverP()
91 << ", abs(qOverP) = " << fabs((*TrackParticleItr)->qOverP()));
92
93 // Event/xAOD/xAODTrackingCnv/trunk/src/TrackParticleCnvAlg.cxx#L190 Info from Nick Styles
94
96 uint8_t numberOfBLayerHits_tmp=0;
97 if ( (*TrackParticleItr)->summaryValue(numberOfBLayerHits_tmp,xAOD::numberOfBLayerHits) ){
98 numberOfBLayerHits += numberOfBLayerHits_tmp;
99 }
101 uint8_t numberOfPixelHits_tmp = 0;
102 if ( (*TrackParticleItr)->summaryValue(numberOfPixelHits_tmp,xAOD::numberOfPixelHits)){
103 numberOfPixelHits += numberOfPixelHits_tmp;
104 }
106 uint8_t numberOfTRTHits_tmp = 0;
107 if ( (*TrackParticleItr)->summaryValue(numberOfTRTHits_tmp,xAOD::numberOfTRTHits)){
108 numberOfTRTHits += numberOfTRTHits_tmp;
109 }
111 uint8_t numberOfSCTHits_tmp = 0;
112 if ( (*TrackParticleItr)->summaryValue(numberOfSCTHits_tmp,xAOD::numberOfSCTHits)){
113 numberOfSCTHits += numberOfSCTHits_tmp;
114 }
115 labelStr = "_PixelHits"+DataType( (double)numberOfPixelHits ).toString()
116 + "_SCTHits"+DataType( (double)numberOfSCTHits ).toString()
117 + "_BLayerHits"+DataType( (double)numberOfBLayerHits ).toString()
118 + "_TRTHits"+DataType( (double)numberOfTRTHits ).toString() ;
119
120 countHits = (double)numberOfBLayerHits + (double)numberOfPixelHits
121 + (double)numberOfSCTHits + (double)numberOfTRTHits;
122
123 ATH_MSG_VERBOSE(" TrackParticle #" << counter
124 << " BLayer hits: " << (double)numberOfBLayerHits
125 << ", Pixel hits: " << (double)numberOfPixelHits
126 << ", SCT hits: " << (double)numberOfSCTHits
127 << ", TRT hits: " << (double)numberOfTRTHits
128 << ", Total hits: " << countHits
129 << "; Label: " << labelStr);
130
131 nBLayerHits.emplace_back( DataType( (double)numberOfBLayerHits ));
132 nPixHits.emplace_back( DataType( (double)numberOfPixelHits ));
133 nSCTHits.emplace_back( DataType( (double)numberOfSCTHits ));
134 nTRTHits.emplace_back( DataType( (double)numberOfTRTHits ));
135
136 label.emplace_back( DataType( labelStr ).toString() );
137 d0.emplace_back(DataType((*TrackParticleItr)->d0()/CLHEP::cm));
138 z0.emplace_back(DataType((*TrackParticleItr)->z0()/CLHEP::cm));
139 phi0.emplace_back(DataType((*TrackParticleItr)->phi()));
140
141 // pt is always positive, get charge from qOverP
142 myQOverP = (*TrackParticleItr)->qOverP() ;
143 if (fabs(myQOverP) != myQOverP){
144 charge = -1.;
145 }else{
146 charge = +1.;
147 }
148 pt.emplace_back(DataType(( charge*(*TrackParticleItr)->pt() )/CLHEP::GeV));
149
150 if ( (*TrackParticleItr)->theta() == 0.) {
151 cotTheta.emplace_back(DataType(9999.));
152 } else {
153 cotTheta.emplace_back(DataType(1./tan((*TrackParticleItr)->theta())));
154 }
155 } // end TrackParticleIterator
156
157 // four-vectors
158 DataMap["d0"] = d0;
159 DataMap["z0"] = z0;
160 DataMap["pt"] = pt;
161 DataMap["phi0"] = phi0;
162 DataMap["cotTheta"] = cotTheta;
163 DataMap["label"] = label;
164 DataMap["nBLayerHits"] = nBLayerHits;
165 DataMap["nPixHits"] = nPixHits;
166 DataMap["nSCTHits"] = nSCTHits;
167 DataMap["nTRTHits"] = nTRTHits;
168
169 ATH_MSG_DEBUG(dataTypeName() << " retrieved with " << d0.size() << " entries");
170
171 //All collections retrieved okay
172 return DataMap;
173
174 }
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_DEBUG(x)
double charge(const T &p)
Definition AtlasPID.h:997
OFFLINE_FRAGMENTS_NAMESPACE::PointerType DataType
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
size_type size() const noexcept
Returns the number of elements in the collection.
virtual std::string dataTypeName() const
Return the name of the data type that is generated by this retriever.
std::string label(const std::string &format, int i)
Definition label.h:19
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
float nSCTHits(const U &p)
float nTRTHits(const U &p)
std::map< std::string, DataVect > DataMap
Definition DataType.h:59
std::vector< DataType > DataVect
Defines a map with a key and a vector of DataType objects e.g.
Definition DataType.h:58
@ numberOfTRTHits
number of TRT hits [unit8_t].
@ numberOfBLayerHits
these are the hits in the first pixel layer, i.e.
@ numberOfSCTHits
number of hits in SCT [unit8_t].
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].

◆ initialize()

StatusCode JiveXML::xAODTrackParticleRetriever::initialize ( )
virtual

Definition at line 22 of file xAODTrackParticleRetriever.cxx.

22 {
23 ATH_CHECK(m_keys.initialize());
24 return StatusCode::SUCCESS;
25 }
#define ATH_CHECK
Evaluate an expression and check for errors.
SG::ReadHandleKeyArray< xAOD::TrackParticleContainer > m_keys

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

const InterfaceID & JiveXML::IDataRetriever::interfaceID ( )
inlinestaticinherited

Return the interface identifier.

Definition at line 40 of file IDataRetriever.h.

40{ return IID_IDataRetriever; }
static const InterfaceID IID_IDataRetriever("JiveXML::IDataRetriever", 3, 0)

◆ msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ retrieve()

StatusCode JiveXML::xAODTrackParticleRetriever::retrieve ( ToolHandle< IFormatTool > & FormatTool)
virtual

For each TrackParticle collections retrieve basic parameters.

Parameters
FormatToolthe tool that will create formated output from the DataMap

Implements JiveXML::IDataRetriever.

Definition at line 31 of file xAODTrackParticleRetriever.cxx.

31 {
32
33 ATH_MSG_DEBUG("In retrieve()");
34
35
36 // Loop through the keys and retrieve the corresponding data
37 for (const auto& key : m_keys) {
38 SG::ReadHandle<xAOD::TrackParticleContainer> cont(key);
39 if (cont.isValid()) {
40 DataMap data = getData(&(*cont));
41 if (FormatTool->AddToEvent(dataTypeName(), key.key() + "_xAOD", &data).isFailure()) {
42 ATH_MSG_WARNING("Failed to add collection " << key.key());
43 } else {
44 ATH_MSG_DEBUG(" (" << key.key() << ") retrieved");
45 }
46 } else {
47 ATH_MSG_WARNING("Collection " << key.key() << " not found in SG");
48 }
49 }
50
51 return StatusCode::SUCCESS;
52 }
#define ATH_MSG_WARNING(x)
char data[hepevt_bytes_allocation_ATLAS]
Definition HepEvt.cxx:11
const DataMap getData(const xAOD::TrackParticleContainer *)
Puts the variables into a DataMap.

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_keys

SG::ReadHandleKeyArray<xAOD::TrackParticleContainer> JiveXML::xAODTrackParticleRetriever::m_keys {this,"TrackParticleCollections", {"InDetTrackParticles","InDetLargeD0TrackParticles","CombinedMuonTrackParticles","GSFTrackParticles"}, "Name of the Track Particle containers that will be written out to the JiveXML file"}
private

Definition at line 49 of file xAODTrackParticleRetriever.h.

49{this,"TrackParticleCollections", {"InDetTrackParticles","InDetLargeD0TrackParticles","CombinedMuonTrackParticles","GSFTrackParticles"}, "Name of the Track Particle containers that will be written out to the JiveXML file"};

◆ m_typeName

const std::string JiveXML::xAODTrackParticleRetriever::m_typeName = "Track"
private

The data type that is generated by this retriever.

Definition at line 47 of file xAODTrackParticleRetriever.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.


The documentation for this class was generated from the following files: