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", m_truthParticleName, "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 194 of file InDetPhysValTruthDecoratorAlg.cxx.

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

81 {
82 SG::ReadHandle<xAOD::TruthParticleContainer> ptruth(m_truthParticleName, ctx);
83 if ((not ptruth.isValid())) {
84 return StatusCode::FAILURE;
85 }
86 std::size_t ptruth_size=ptruth->size();
87
88 std::vector<unsigned int> truthIndexMap;
89 if (!m_truthParticleIndexDecor.empty()) {
90 SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int> decor_index(m_truthParticleIndexDecor, ctx);
91 unsigned int max_size=0;
92 assert( ptruth_size < std::numeric_limits<unsigned int>::max());
93 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
94 max_size = std::max( max_size, decor_index(*truth_particle) );
95 }
96 if (max_size>=std::numeric_limits<unsigned int>::max()) {
97 ATH_MSG_ERROR("Truth index exceed max allowed range.");
98 return StatusCode::FAILURE;
99 }
100 ++max_size;
101 truthIndexMap.resize( max_size, std::numeric_limits<unsigned int>::max());
102 unsigned int new_index=0;
103 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
104 truthIndexMap.at( decor_index(*truth_particle) ) = new_index;
105 ++new_index;
106 }
107 }
108 else {
109 truthIndexMap.reserve(ptruth_size);
110 for (unsigned int i=0; i<ptruth_size; ++i) {
111 truthIndexMap.push_back(i);
112 }
113 }
114
115 std::vector< IDPVM::OptionalDecoration<xAOD::TruthParticleContainer,float> >
116 float_decor( IDPVM::createDecoratorsIfNeeded(*ptruth, m_decor, ctx, msgLvl(MSG::DEBUG)) );
117
119 std::vector< std::array<uint16_t, kNClusterTypes> > tp_clustercount;
120 tp_clustercount.resize(ptruth_size,std::array<uint16_t,kNClusterTypes>{});
121 unsigned int missing_truth_particle=0u;
122 //Loop over the pixel and sct clusters to fill the truth barcode - cluster count maps
123 SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> sctClusters(m_truthSCTClusterName, ctx);
124 SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> pixelClusters(m_truthPixelClusterName, ctx);
125 //only decorate the truth particles with truth silicon hits if both containers are available
126 if (sctClusters.isValid() && pixelClusters.isValid()) {
127 for (const auto *const sct : *sctClusters) {
128 const xAOD::TrackMeasurementValidation* sctCluster = sct;
129 static const SG::AuxElement::ConstAccessor< std::vector<unsigned int> > truthIndexAcc("truth_index");
130 if (truthIndexAcc.isAvailable(*sctCluster)) {
131 const std::vector<unsigned int> &truth_indices = truthIndexAcc(*sctCluster);
132 for (auto index : truth_indices) {
133 if (index != std::numeric_limits<unsigned int>::max()) {
134 if (index < truthIndexMap.size() && truthIndexMap[index] != std::numeric_limits<unsigned int>::max()) {
135 ++tp_clustercount.at(truthIndexMap[index])[kSCT];
136 }
137 else {
138 ++missing_truth_particle;
139 }
140 }
141 }
142 }
143 } // Loop over SCT clusters
144
145 for (const auto *const pix : *pixelClusters) {
146 const xAOD::TrackMeasurementValidation* pixCluster = pix;
147 static const SG::AuxElement::ConstAccessor< std::vector<unsigned int> > truthIndexAcc("truth_index");
148 if (truthIndexAcc.isAvailable(*pixCluster)) {
149 const std::vector<unsigned int> &truth_indices = truthIndexAcc(*pixCluster);
150 for (auto index : truth_indices) {
151 if (index != std::numeric_limits<unsigned int>::max()) {
152 if (index < truthIndexMap.size() && truthIndexMap[index] != std::numeric_limits<unsigned int>::max()) {
153 ++tp_clustercount.at(truthIndexMap[index])[kPixel];
154 }
155 else {
156 ++missing_truth_particle;
157 }
158 }
159 }
160 }
161 } // Loop over PIX clusters
162 }
163 m_nMissingTruthParticles += missing_truth_particle;
164
165 if (not float_decor.empty()) {
166 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosX(m_beamSpotDecoKey[0], ctx);
167 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosY(m_beamSpotDecoKey[1], ctx);
168 SG::ReadDecorHandle<xAOD::EventInfo, float> beamPosZ(m_beamSpotDecoKey[2], ctx);
169 Amg::Vector3D beamPos = Amg::Vector3D(beamPosX(0), beamPosY(0), beamPosZ(0));
170
171 if ( m_truthSelectionTool.get() ) {
172 CutFlow tmp_cut_flow(m_truthSelectionTool->nCuts());
173 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
174 auto passed = m_truthSelectionTool->accept(truth_particle);
175 tmp_cut_flow.update( passed.missingCuts() );
176 if (not passed) continue;
177 decorateTruth(*truth_particle, float_decor, beamPos, tp_clustercount);
178 }
179 std::lock_guard<std::mutex> lock(m_mutex);
180 m_cutFlow.merge(std::move(tmp_cut_flow));
181 }
182 else {
183 for (const xAOD::TruthParticle *truth_particle : *ptruth) {
184 decorateTruth(*truth_particle, float_decor, beamPos, tp_clustercount);
185 }
186 }
187 }
188 return StatusCode::SUCCESS;
189}
#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 69 of file InDetPhysValTruthDecoratorAlg.cxx.

69 {
70 if (not m_truthSelectionTool.name().empty()) {
71 std::lock_guard<std::mutex> lock(m_mutex);
72 ATH_MSG_DEBUG( "Truth selection cut flow : " << m_cutFlow.report(m_truthSelectionTool->names()) );
73 }
75 ATH_MSG_INFO( "Clusters which reference missing / thinned truth particles : " << m_nMissingTruthParticles );
76 }
77 return StatusCode::SUCCESS;
78}
#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());
51
52 std::vector<std::string> decor_names(kNDecorators);
53 decor_names[kDecorD0]="d0";
54 decor_names[kDecorZ0]="z0";
55 decor_names[kDecorPhi]="phi";
56 decor_names[kDecorTheta]="theta";
57 decor_names[kDecorZ0st]="z0st";
58 decor_names[kDecorQOverP]="qOverP";
59 decor_names[kDecorProdR]="prodR";
60 decor_names[kDecorProdZ]="prodZ";
61 decor_names[kDecorNSilHits]="nSilHits";
62
64 assert( m_decor.size() == kNDecorators);
65 return StatusCode::SUCCESS;
66}
#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", m_truthParticleName, "origTruthIndex", "decoration name for the original truth particle index."}
private

Definition at line 82 of file InDetPhysValTruthDecoratorAlg.h.

83{this, "TruthParticleIndexDecoration", m_truthParticleName, "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: