ATLAS Offline Software
Loading...
Searching...
No Matches
InDetPhysValTruthDecoratorAlg Class Reference

#include <InDetPhysValTruthDecoratorAlg.h>

Inheritance diagram for InDetPhysValTruthDecoratorAlg:

Public Member Functions

 InDetPhysValTruthDecoratorAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~InDetPhysValTruthDecoratorAlg ()
virtual StatusCode initialize ()
virtual StatusCode finalize ()
virtual StatusCode execute (const EventContext &ctx) const
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
virtual void setFilterPassed (bool state, const EventContext &ctx) const
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 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

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

enum  { kPixel , kSCT , kNClusterTypes }
enum  EDecorations {
  kDecorD0 , kDecorZ0 , kDecorPhi , kDecorTheta ,
  kDecorZ0st , kDecorQOverP , kDecorProdR , kDecorProdZ ,
  kDecorNSilHits , kNDecorators
}
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

bool decorateTruth (const xAOD::TruthParticle &particle, std::vector< std::pair< SG::WriteDecorHandle< xAOD::TruthParticleContainer, float >, bool > > &float_decor, const Amg::Vector3D &beamPos, const std::vector< std::array< uint16_t, kNClusterTypes > > &counts) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

PublicToolHandle< Trk::IExtrapolatorm_extrapolator {this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""}
SG::ReadDecorHandleKeyArray< xAOD::EventInfom_beamSpotDecoKey
PublicToolHandle< IAthSelectionToolm_truthSelectionTool {this,"TruthSelectionTool","",""}
std::mutex m_mutex
CutFlow m_cutFlow ATLAS_THREAD_SAFE
std::atomic< std::size_t > m_nMissingTruthParticles = 0u
std::atomic< bool > m_errorEmitted {false}
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthParticleName {this, "TruthParticleContainerName", "TruthParticles", ""}
 TruthParticle container's name needed to create decorators.
Gaudi::Property< std::string > m_prefix {this, "Prefix", "", "Decoration prefix to avoid clashes."}
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainerm_truthPixelClusterName {this, "PixelClusterContainerName", "PixelClusters", ""}
 TruthPixelClusterContainer and TruthSCTClusterContainer needed for truth silicon hit cut.
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainerm_truthSCTClusterName {this, "SCTClusterContainerName", "SCT_Clusters", ""}
SG::ReadDecorHandleKey< xAOD::TruthParticleContainerm_truthParticleIndexDecor {this, "TruthParticleIndexDecoration", "origTruthIndex", "decoration name for the original truth particle index."}
std::vector< std::pair< SG::WriteDecorHandleKey< xAOD::TruthParticleContainer >, SG::AuxElement::ConstAccessor< float > > > m_decor
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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

Definition at line 35 of file InDetPhysValTruthDecoratorAlg.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ anonymous enum

◆ EDecorations

Constructor & Destructor Documentation

◆ InDetPhysValTruthDecoratorAlg()

InDetPhysValTruthDecoratorAlg::InDetPhysValTruthDecoratorAlg ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 29 of file InDetPhysValTruthDecoratorAlg.cxx.

29 :
30 AthReentrantAlgorithm(name, pSvcLocator)
31{
32}

◆ ~InDetPhysValTruthDecoratorAlg()

InDetPhysValTruthDecoratorAlg::~InDetPhysValTruthDecoratorAlg ( )
virtual

Definition at line 34 of file InDetPhysValTruthDecoratorAlg.cxx.

34 {
35 // nop
36}

Member Function Documentation

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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>

◆ decorateTruth()

bool InDetPhysValTruthDecoratorAlg::decorateTruth ( const xAOD::TruthParticle & particle,
std::vector< std::pair< SG::WriteDecorHandle< xAOD::TruthParticleContainer, float >, bool > > & float_decor,
const Amg::Vector3D & beamPos,
const std::vector< std::array< uint16_t, kNClusterTypes > > & counts ) const
private

Definition at line 197 of file InDetPhysValTruthDecoratorAlg.cxx.

200 {
201 ATH_MSG_VERBOSE("Decorate truth with d0 etc");
202 if (particle.isNeutral()) {
203 return false;
204 }
205 const EventContext& ctx = Gaudi::Hive::currentContext();
206 const Amg::Vector3D momentum(particle.px(), particle.py(), particle.pz());
207 const int pid(particle.pdgId());
208 double charge = particle.charge();
209
210 if (std::isnan(charge)) {
211 ATH_MSG_DEBUG("charge not found on particle with pid " << pid);
212 return false;
213 }
214
215 // @TODO float?
216 float nSiHits = std::accumulate(counts.at(particle.index()).begin(), counts[particle.index()].end(), 0u);
217
218 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorNSilHits], nSiHits);
219
220 const xAOD::TruthVertex* ptruthVertex(nullptr);
221 try{
222 ptruthVertex = particle.prodVtx();
223 } catch (const std::exception& e) {
224 if (not m_errorEmitted) {
225 ATH_MSG_WARNING("A non existent production vertex was requested in calculating the track parameters d0 etc");
226 }
227 m_errorEmitted = true;
228 return false;
229 }
230 if (!ptruthVertex) {
231 ATH_MSG_DEBUG("A production vertex pointer was retrieved, but it is NULL");
232 return false;
233 }
234 const auto xPos = ptruthVertex->x();
235 const auto yPos = ptruthVertex->y();
236 const auto z_truth = ptruthVertex->z();
237 const Amg::Vector3D position(xPos, yPos, z_truth);
238 const float prodR_truth = std::sqrt(xPos * xPos + yPos * yPos);
239 // delete ptruthVertex;ptruthVertex=0;
240 const Trk::CurvilinearParameters cParameters(position, momentum, charge);
241
242 Trk::PerigeeSurface persf(beamPos);
243
244 std::unique_ptr<const Trk::TrackParameters> tP ( m_extrapolator->extrapolate(ctx,
245 cParameters,
246 persf, Trk::anyDirection, false) );
247 if (tP) {
248 float d0_truth = tP->parameters()[Trk::d0];
249 float theta_truth = tP->parameters()[Trk::theta];
250 float z0_truth = tP->parameters()[Trk::z0];
251 float phi_truth = tP->parameters()[Trk::phi];
252 float qOverP_truth = tP->parameters()[Trk::qOverP]; // P or Pt ??
253 float z0st_truth = z0_truth * std::sin(theta_truth);
254
255 // 'safeDecorator' used to prevent a crash in case of adding something which pre-exists.
256 // behaviour chosen is to reject quietly
257 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorD0],d0_truth);
258 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorZ0],z0_truth);
259 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorPhi],phi_truth);
260 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorTheta],theta_truth);
261 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorZ0st],z0st_truth);
262 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorQOverP],qOverP_truth);
263 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorProdR],prodR_truth);
264 IDPVM::decorateOrRejectQuietly(particle,float_decor[kDecorProdZ],z_truth);
265
266 return true;
267 } else {
268 ATH_MSG_DEBUG("The TrackParameters pointer for this TruthParticle is NULL");
269 return false;
270 }
271}
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
double charge(const T &p)
Definition AtlasPID.h:997
PublicToolHandle< Trk::IExtrapolator > m_extrapolator
Eigen::Matrix< double, 3, 1 > Vector3D
void decorateOrRejectQuietly(const T_Cont_Elm &particle, OptionalDecoration< T_Cont, T > &decorator, const T &value)
float nSiHits(const U &p)
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
@ anyDirection
CurvilinearParametersT< TrackParametersDim, Charged, PlaneSurface > CurvilinearParameters
@ theta
Definition ParamDefs.h:66
@ qOverP
perigee
Definition ParamDefs.h:67
@ phi
Definition ParamDefs.h:75
@ d0
Definition ParamDefs.h:63
@ z0
Definition ParamDefs.h:64
TruthVertex_v1 TruthVertex
Typedef to implementation.
Definition TruthVertex.h:15

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::evtStore ( )
inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode InDetPhysValTruthDecoratorAlg::execute ( const EventContext & ctx) const
virtual

truthbarcode-cluster maps to be pre-stored at event level

Definition at line 84 of file InDetPhysValTruthDecoratorAlg.cxx.

84 {
85 SG::ReadHandle<xAOD::TruthParticleContainer> ptruth(m_truthParticleName, ctx);
86 if ((not ptruth.isValid())) {
87 return StatusCode::FAILURE;
88 }
89 std::size_t ptruth_size=ptruth->size();
90
91 std::vector<unsigned int> truthIndexMap;
92 if (!m_truthParticleIndexDecor.empty()) {
93 SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int> decor_index(m_truthParticleIndexDecor, ctx);
94 unsigned int max_size=0;
95 assert( ptruth_size < std::numeric_limits<unsigned int>::max());
96 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
97 max_size = std::max( max_size, decor_index(*truth_particle) );
98 }
99 if (max_size>=std::numeric_limits<unsigned int>::max()) {
100 ATH_MSG_ERROR("Truth index exceed max allowed range.");
101 return StatusCode::FAILURE;
102 }
103 ++max_size;
104 truthIndexMap.resize( max_size, std::numeric_limits<unsigned int>::max());
105 unsigned int new_index=0;
106 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
107 truthIndexMap.at( decor_index(*truth_particle) ) = new_index;
108 ++new_index;
109 }
110 }
111 else {
112 truthIndexMap.reserve(ptruth_size);
113 for (unsigned int i=0; i<ptruth_size; ++i) {
114 truthIndexMap.push_back(i);
115 }
116 }
117
118 std::vector< IDPVM::OptionalDecoration<xAOD::TruthParticleContainer,float> >
119 float_decor( IDPVM::createDecoratorsIfNeeded(*ptruth, m_decor, ctx, msgLvl(MSG::DEBUG)) );
120
122 std::vector< std::array<uint16_t, kNClusterTypes> > tp_clustercount;
123 tp_clustercount.resize(ptruth_size,std::array<uint16_t,kNClusterTypes>{});
124 unsigned int missing_truth_particle=0u;
125 //Loop over the pixel and sct clusters to fill the truth barcode - cluster count maps
126 SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> sctClusters(m_truthSCTClusterName, ctx);
127 SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> pixelClusters(m_truthPixelClusterName, ctx);
128 //only decorate the truth particles with truth silicon hits if both containers are available
129 if (sctClusters.isValid() && pixelClusters.isValid()) {
130 for (const auto *const sct : *sctClusters) {
131 const xAOD::TrackMeasurementValidation* sctCluster = sct;
132 static const SG::AuxElement::ConstAccessor< std::vector<unsigned int> > truthIndexAcc("truth_index");
133 if (truthIndexAcc.isAvailable(*sctCluster)) {
134 const std::vector<unsigned int> &truth_indices = truthIndexAcc(*sctCluster);
135 for (auto index : truth_indices) {
136 if (index != std::numeric_limits<unsigned int>::max()) {
137 if (index < truthIndexMap.size() && truthIndexMap[index] != std::numeric_limits<unsigned int>::max()) {
138 ++tp_clustercount.at(truthIndexMap[index])[kSCT];
139 }
140 else {
141 ++missing_truth_particle;
142 }
143 }
144 }
145 }
146 } // Loop over SCT clusters
147
148 for (const auto *const pix : *pixelClusters) {
149 const xAOD::TrackMeasurementValidation* pixCluster = pix;
150 static const SG::AuxElement::ConstAccessor< std::vector<unsigned int> > truthIndexAcc("truth_index");
151 if (truthIndexAcc.isAvailable(*pixCluster)) {
152 const std::vector<unsigned int> &truth_indices = truthIndexAcc(*pixCluster);
153 for (auto index : truth_indices) {
154 if (index != std::numeric_limits<unsigned int>::max()) {
155 if (index < truthIndexMap.size() && truthIndexMap[index] != std::numeric_limits<unsigned int>::max()) {
156 ++tp_clustercount.at(truthIndexMap[index])[kPixel];
157 }
158 else {
159 ++missing_truth_particle;
160 }
161 }
162 }
163 }
164 } // Loop over PIX clusters
165 }
166 m_nMissingTruthParticles += missing_truth_particle;
167
168 if (not float_decor.empty()) {
169 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosX(m_beamSpotDecoKey[0], ctx);
170 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosY(m_beamSpotDecoKey[1], ctx);
171 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosZ(m_beamSpotDecoKey[2], ctx);
172 Amg::Vector3D beamPos = Amg::Vector3D(beamPosX(0), beamPosY(0), beamPosZ(0));
173
174 if ( m_truthSelectionTool.get() ) {
175 CutFlow tmp_cut_flow(m_truthSelectionTool->nCuts());
176 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
177 auto passed = m_truthSelectionTool->accept(truth_particle);
178 tmp_cut_flow.update( passed.missingCuts() );
179 if (not passed) continue;
180 decorateTruth(*truth_particle, float_decor, beamPos, tp_clustercount);
181 }
182 std::lock_guard<std::mutex> lock(m_mutex);
183 m_cutFlow.merge(std::move(tmp_cut_flow));
184 }
185 else {
186 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
187 decorateTruth(*truth_particle, float_decor, beamPos, tp_clustercount);
188 }
189 }
190 }
191 return StatusCode::SUCCESS;
192}
#define ATH_MSG_ERROR(x)
bool passed(DecisionID id, const DecisionIDContainer &)
checks if required decision ID is in the set of IDs in the container
bool msgLvl(const MSG::Level lvl) const
PublicToolHandle< IAthSelectionTool > m_truthSelectionTool
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthParticleIndexDecor
std::vector< std::pair< SG::WriteDecorHandleKey< xAOD::TruthParticleContainer >, SG::AuxElement::ConstAccessor< float > > > m_decor
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainer > m_truthPixelClusterName
TruthPixelClusterContainer and TruthSCTClusterContainer needed for truth silicon hit cut.
std::atomic< std::size_t > m_nMissingTruthParticles
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainer > m_truthSCTClusterName
SG::ReadDecorHandleKeyArray< xAOD::EventInfo > m_beamSpotDecoKey
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthParticleName
TruthParticle container's name needed to create decorators.
bool decorateTruth(const xAOD::TruthParticle &particle, std::vector< std::pair< SG::WriteDecorHandle< xAOD::TruthParticleContainer, float >, bool > > &float_decor, const Amg::Vector3D &beamPos, const std::vector< std::array< uint16_t, kNClusterTypes > > &counts) const
SG::ConstAccessor< T, ALLOC > ConstAccessor
Definition AuxElement.h:569
std::vector< OptionalDecoration< T_Cont, T > > createDecoratorsIfNeeded(const T_Cont &container, const std::vector< WriteKeyAccessorPair< T_Cont, T > > &keys, const EventContext &ctx, bool verbose=false)
@ u
Enums for curvilinear frames.
Definition ParamDefs.h:77
TrackMeasurementValidation_v1 TrackMeasurementValidation
Reference the current persistent version:
TruthParticle_v1 TruthParticle
Typedef to implementation.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ finalize()

StatusCode InDetPhysValTruthDecoratorAlg::finalize ( )
virtual

Definition at line 72 of file InDetPhysValTruthDecoratorAlg.cxx.

72 {
73 if (not m_truthSelectionTool.name().empty()) {
74 std::lock_guard<std::mutex> lock(m_mutex);
75 ATH_MSG_DEBUG( "Truth selection cut flow : " << m_cutFlow.report(m_truthSelectionTool->names()) );
76 }
78 ATH_MSG_INFO( "Clusters which reference missing / thinned truth particles : " << m_nMissingTruthParticles );
79 }
80 return StatusCode::SUCCESS;
81}
#define ATH_MSG_INFO(x)

◆ initialize()

StatusCode InDetPhysValTruthDecoratorAlg::initialize ( )
virtual

Definition at line 39 of file InDetPhysValTruthDecoratorAlg.cxx.

39 {
40 ATH_CHECK(m_extrapolator.retrieve());
41 ATH_CHECK(m_beamSpotDecoKey.initialize());
42 ATH_CHECK( m_truthPixelClusterName.initialize() );
43 ATH_CHECK( m_truthSCTClusterName.initialize() );
44 ATH_CHECK( m_truthSelectionTool.retrieve( EnableTool { not m_truthSelectionTool.name().empty() } ) );
45 if (not m_truthSelectionTool.name().empty() ) {
46 m_cutFlow = CutFlow(m_truthSelectionTool->nCuts() );
47 }
48
49 ATH_CHECK( m_truthParticleName.initialize());
50 if (!m_truthParticleIndexDecor.key().empty()) {
52 }
54
55 std::vector<std::string> decor_names(kNDecorators);
56 decor_names[kDecorD0]="d0";
57 decor_names[kDecorZ0]="z0";
58 decor_names[kDecorPhi]="phi";
59 decor_names[kDecorTheta]="theta";
60 decor_names[kDecorZ0st]="z0st";
61 decor_names[kDecorQOverP]="qOverP";
62 decor_names[kDecorProdR]="prodR";
63 decor_names[kDecorProdZ]="prodZ";
64 decor_names[kDecorNSilHits]="nSilHits";
65
67 assert( m_decor.size() == kNDecorators);
68 return StatusCode::SUCCESS;
69}
#define ATH_CHECK
Evaluate an expression and check for errors.
Gaudi::Property< std::string > m_prefix
void createDecoratorKeysAndAccessor(T_Parent &parent, const SG::ReadHandleKey< T_Cont > &container_key, const std::string &prefix, const std::vector< std::string > &decor_names, std::vector< WriteKeyAccessorPair< T_Cont, T > > &decor_out)

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

◆ msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

77{
78 return BaseAlg::sysExecute (ctx);
79}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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 }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

CutFlow m_cutFlow InDetPhysValTruthDecoratorAlg::ATLAS_THREAD_SAFE
mutableprivate

Definition at line 63 of file InDetPhysValTruthDecoratorAlg.h.

◆ m_beamSpotDecoKey

SG::ReadDecorHandleKeyArray<xAOD::EventInfo> InDetPhysValTruthDecoratorAlg::m_beamSpotDecoKey
private
Initial value:
{this, "BeamSpotDecoKeys",
{"EventInfo.beamPosX", "EventInfo.beamPosY", "EventInfo.beamPosZ"},
"Beamspot position decoration keys"}

Definition at line 54 of file InDetPhysValTruthDecoratorAlg.h.

55 {this, "BeamSpotDecoKeys",
56 {"EventInfo.beamPosX", "EventInfo.beamPosY", "EventInfo.beamPosZ"},
57 "Beamspot position decoration keys"};

◆ m_decor

std::vector< std::pair<SG::WriteDecorHandleKey<xAOD::TruthParticleContainer>,SG::AuxElement::ConstAccessor<float> > > InDetPhysValTruthDecoratorAlg::m_decor
private

Definition at line 98 of file InDetPhysValTruthDecoratorAlg.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_errorEmitted

std::atomic<bool> InDetPhysValTruthDecoratorAlg::m_errorEmitted {false}
mutableprivate

Definition at line 66 of file InDetPhysValTruthDecoratorAlg.h.

66{false};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_extrapolator

PublicToolHandle<Trk::IExtrapolator> InDetPhysValTruthDecoratorAlg::m_extrapolator {this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""}
private

Definition at line 52 of file InDetPhysValTruthDecoratorAlg.h.

53{this,"Extrapolator","Trk::Extrapolator/AtlasExtrapolator",""};

◆ m_mutex

std::mutex InDetPhysValTruthDecoratorAlg::m_mutex
mutableprivate

Definition at line 62 of file InDetPhysValTruthDecoratorAlg.h.

◆ m_nMissingTruthParticles

std::atomic<std::size_t> InDetPhysValTruthDecoratorAlg::m_nMissingTruthParticles = 0u
mutableprivate

Definition at line 65 of file InDetPhysValTruthDecoratorAlg.h.

◆ m_prefix

Gaudi::Property<std::string> InDetPhysValTruthDecoratorAlg::m_prefix {this, "Prefix", "", "Decoration prefix to avoid clashes."}
private

Definition at line 72 of file InDetPhysValTruthDecoratorAlg.h.

73{this, "Prefix", "", "Decoration prefix to avoid clashes."};

◆ m_truthParticleIndexDecor

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> InDetPhysValTruthDecoratorAlg::m_truthParticleIndexDecor {this, "TruthParticleIndexDecoration", "origTruthIndex", "decoration name for the original truth particle index."}
private

Definition at line 82 of file InDetPhysValTruthDecoratorAlg.h.

83{this, "TruthParticleIndexDecoration", "origTruthIndex", "decoration name for the original truth particle index."};

◆ m_truthParticleName

SG::ReadHandleKey<xAOD::TruthParticleContainer> InDetPhysValTruthDecoratorAlg::m_truthParticleName {this, "TruthParticleContainerName", "TruthParticles", ""}
private

TruthParticle container's name needed to create decorators.

Definition at line 69 of file InDetPhysValTruthDecoratorAlg.h.

70{this, "TruthParticleContainerName", "TruthParticles", ""};

◆ m_truthPixelClusterName

SG::ReadHandleKey<xAOD::TrackMeasurementValidationContainer> InDetPhysValTruthDecoratorAlg::m_truthPixelClusterName {this, "PixelClusterContainerName", "PixelClusters", ""}
private

TruthPixelClusterContainer and TruthSCTClusterContainer needed for truth silicon hit cut.

Definition at line 76 of file InDetPhysValTruthDecoratorAlg.h.

77{this, "PixelClusterContainerName", "PixelClusters", ""};

◆ m_truthSCTClusterName

SG::ReadHandleKey<xAOD::TrackMeasurementValidationContainer> InDetPhysValTruthDecoratorAlg::m_truthSCTClusterName {this, "SCTClusterContainerName", "SCT_Clusters", ""}
private

Definition at line 79 of file InDetPhysValTruthDecoratorAlg.h.

80{this, "SCTClusterContainerName", "SCT_Clusters", ""};

◆ m_truthSelectionTool

PublicToolHandle<IAthSelectionTool> InDetPhysValTruthDecoratorAlg::m_truthSelectionTool {this,"TruthSelectionTool","",""}
private

Definition at line 59 of file InDetPhysValTruthDecoratorAlg.h.

60{this,"TruthSelectionTool","",""};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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