ATLAS Offline Software
Public Types | Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
ActsTrk::MeasurementToTrackParticleDecoration Class Reference

#include <MeasurementToTrackParticleDecoration.h>

Inheritance diagram for ActsTrk::MeasurementToTrackParticleDecoration:
Collaboration diagram for ActsTrk::MeasurementToTrackParticleDecoration:

Public Types

enum  Subdetector {
  INVALID_DETECTOR =-1, INNERMOST_PIXEL, PIXEL, STRIP,
  N_SUBDETECTORS
}
 
enum  Region { INVALID_REGION =-1, BARREL, ENDCAP }
 
enum  MeasurementType {
  INVALID_MEASUREMENT =-1, HIT, OUTLIER, HOLE,
  BIASED, UNBIASED
}
 

Public Member Functions

 MeasurementToTrackParticleDecoration (const std::string &name, ISvcLocator *pSvcLocator)
 
virtual ~MeasurementToTrackParticleDecoration ()=default
 
virtual StatusCode initialize () override
 
virtual StatusCode execute (const EventContext &ctx) const override
 
virtual StatusCode finalize () override
 
float getChi2Contribution (const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &state) const
 
std::pair< Acts::BoundVector, Acts::BoundMatrix > getUnbiasedTrackParameters (const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &state, bool useSmoothed=true) const
 
float evaluatePull (const float residual, const float measurementCovariance, const float trackParameterCovariance, const bool evaluateUnbiased) const
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual bool isClonable () const override
 Specify if the algorithm is clonable. More...
 
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
 
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. More...
 
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. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
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. More...
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

ToolHandle< IActsTrackingGeometryToolm_trackingGeometryTool {this, "TrackingGeometryTool", ""}
 
SG::ReadHandleKey< xAOD::TrackParticleContainerm_trackParticlesKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementRegionKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementDetectorKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementLayerKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_chi2HitPredictedKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_chi2HitFilteredKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementTypeKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementPhiWidthKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementEtaWidthKey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_residualLocXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_pullLocXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementLocXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_trackParameterLocXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementLocCovXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_trackParameterLocCovXkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_residualLocYkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_pullLocYkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementLocYkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_trackParameterLocYkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_measurementLocCovYkey
 
SG::WriteDecorHandleKey< xAOD::TrackParticleContainerm_trackParameterLocCovYkey
 
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 37 of file MeasurementToTrackParticleDecoration.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

◆ MeasurementType

Enumerator
INVALID_MEASUREMENT 
HIT 
OUTLIER 
HOLE 
BIASED 
UNBIASED 

Definition at line 52 of file MeasurementToTrackParticleDecoration.h.

52  {
54  };

◆ Region

Enumerator
INVALID_REGION 
BARREL 
ENDCAP 

Definition at line 49 of file MeasurementToTrackParticleDecoration.h.

49  {
51  };

◆ Subdetector

Enumerator
INVALID_DETECTOR 
INNERMOST_PIXEL 
PIXEL 
STRIP 
N_SUBDETECTORS 

Definition at line 46 of file MeasurementToTrackParticleDecoration.h.

Constructor & Destructor Documentation

◆ MeasurementToTrackParticleDecoration()

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

Definition at line 21 of file MeasurementToTrackParticleDecoration.cxx.

21  :
22  AthReentrantAlgorithm(name,pSvcLocator)
23  {}

◆ ~MeasurementToTrackParticleDecoration()

virtual ActsTrk::MeasurementToTrackParticleDecoration::~MeasurementToTrackParticleDecoration ( )
virtualdefault

Member Function Documentation

◆ cardinality()

unsigned int AthReentrantAlgorithm::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.

Override this to return 0 for reentrant algorithms.

Definition at line 55 of file AthReentrantAlgorithm.cxx.

56 {
57  return 0;
58 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType  
)
inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType  
)
inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType  
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ 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.

95 { return m_detStore; }

◆ evaluatePull()

float ActsTrk::MeasurementToTrackParticleDecoration::evaluatePull ( const float  residual,
const float  measurementCovariance,
const float  trackParameterCovariance,
const bool  evaluateUnbiased 
) const

Definition at line 398 of file MeasurementToTrackParticleDecoration.cxx.

401  {
402  float correlation = evaluateUnbiased ? 1. : -1.;
403  float residualCovariance = measurementCovariance + correlation*trackParameterCovariance;
404  if (residualCovariance<=0.) {
405  // If the total covariance is non-positive return 0
406  ATH_MSG_DEBUG("--- Total covariance for pull evaluation is non-positive! Returning pulls = 0!");
407  return 0.;
408  }
409  return 0.001 * residual/std::sqrt(residualCovariance);
410  }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( ) const
inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode ActsTrk::MeasurementToTrackParticleDecoration::execute ( const EventContext &  ctx) const
overridevirtual

Definition at line 54 of file MeasurementToTrackParticleDecoration.cxx.

55  {
56  ATH_MSG_DEBUG("Executing MeasurementToTrackParticleDecoration...");
57 
58  auto tgContext = m_trackingGeometryTool->getGeometryContext(ctx).context();
59 
66 
69 
76 
83 
85  ATH_CHECK(trackParticlesHandle.isValid());
86  const xAOD::TrackParticleContainer *track_particles = trackParticlesHandle.cptr();
87 
88  static const SG::AuxElement::ConstAccessor<ElementLink<ActsTrk::TrackContainer> > actsTrackLink("actsTrack");
89 
90  for (const xAOD::TrackParticle *track_particle : *track_particles) {
91  ElementLink<ActsTrk::TrackContainer> link_to_track = actsTrackLink(*track_particle);
92  ATH_CHECK(link_to_track.isValid());
93 
94  // to ensure that the code does not suggest something stupid (i.e. creating an unnecessary copy)
96  std::optional<ActsTrk::TrackContainer::ConstTrackProxy> >::value);
97  std::optional<ActsTrk::TrackContainer::ConstTrackProxy> optional_track = *link_to_track;
98 
99  if ( not optional_track.has_value() ) {
100  ATH_MSG_WARNING("Invalid track link for particle " << track_particle->index() << ". Skipping track..");
101  continue;
102  }
103 
104  ATH_MSG_DEBUG("Track link found for track particle with index " << track_particle->index());
105  ActsTrk::TrackContainer::ConstTrackProxy track = optional_track.value();
106 
107  std::vector<int>& regions{measurementRegionHandle(*track_particle)};
108  regions.reserve(track.nMeasurements());
109  std::vector<int>& detectors{measurementDetectorHandle(*track_particle)};
110  detectors.reserve(track.nMeasurements());
111  std::vector<int>& layers{measurementLayerHandle(*track_particle)};
112  layers.reserve(track.nMeasurements());
113  std::vector<int>& types{measurementTypeHandle(*track_particle)};
114  types.reserve(track.nMeasurements());
115  std::vector<float>& predchi2s{chi2HitPredictedHandle(*track_particle)};
116  predchi2s.reserve(track.nMeasurements());
117  std::vector<float>& filtchi2s{chi2HitFilteredHandle(*track_particle)};
118  filtchi2s.reserve(track.nMeasurements());
119 
120  std::vector<int>& sizesPhi{measurementPhiWidthHandle(*track_particle)};
121  sizesPhi.reserve(track.nMeasurements());
122  std::vector<int>& sizesEta{measurementEtaWidthHandle(*track_particle)};
123  sizesEta.reserve(track.nMeasurements());
124  std::vector<float>& residualsLocX{residualLocXhandle(*track_particle)};
125  residualsLocX.reserve(track.nMeasurements());
126  std::vector<float>& pullsLocX{pullLocXhandle(*track_particle)};
127  pullsLocX.reserve(track.nMeasurements());
128  std::vector<float>& measurementsLocX{measurementLocXhandle(*track_particle)};
129  measurementsLocX.reserve(track.nMeasurements());
130  std::vector<float>& trackParametersLocX{trackParameterLocXhandle(*track_particle)};
131  trackParametersLocX.reserve(track.nMeasurements());
132  std::vector<float>& measurementsLocCovX{measurementLocCovXhandle(*track_particle)};
133  measurementsLocCovX.reserve(track.nMeasurements());
134  std::vector<float>& trackParametersLocCovX{trackParameterLocCovXhandle(*track_particle)};
135  trackParametersLocCovX.reserve(track.nMeasurements());
136  std::vector<float>& residualsLocY{residualLocYhandle(*track_particle)};
137  residualsLocY.reserve(track.nMeasurements());
138  std::vector<float>& pullsLocY{pullLocYhandle(*track_particle)};
139  pullsLocY.reserve(track.nMeasurements());
140  std::vector<float>& measurementsLocY{measurementLocYhandle(*track_particle)};
141  measurementsLocY.reserve(track.nMeasurements());
142  std::vector<float>& trackParametersLocY{trackParameterLocYhandle(*track_particle)};
143  trackParametersLocY.reserve(track.nMeasurements());
144  std::vector<float>& measurementsLocCovY{measurementLocCovYhandle(*track_particle)};
145  measurementsLocCovY.reserve(track.nMeasurements());
146  std::vector<float>& trackParametersLocCovY{trackParameterLocCovYhandle(*track_particle)};
147  trackParametersLocCovY.reserve(track.nMeasurements());
148 
149  for (const auto state : track.trackStatesReversed()) {
150 
151  auto flag = state.typeFlags();
152  // consider holes and measurements (also outliers)
153  bool anyHit = flag.test(Acts::TrackStateFlag::HoleFlag) or flag.test(Acts::TrackStateFlag::MeasurementFlag);
154  if (not anyHit) {
155  ATH_MSG_DEBUG("--- This is not a hit measurement, skipping...");
156  continue;
157  }
158 
159  // starting with invalid values and setting them where needed.
160 
161  int detector = -999;
162  int region = -999;
163  int layer = -999;
164  int type = -999;
165  float chi2_hit_predicted = -999.;
166  float chi2_hit_filtered = -999.;
167  int sizePhi = -999;
168  int sizeEta = -999;
169  float residualLocX = -999.;
170  float pullLocX = -999.;
171  float measurementLocX = -999.;
172  float trackParameterLocX = -999.;
173  float measurementLocCovX = -999.;
174  float trackParameterLocCovX = -999;
175  float residualLocY = -999.;
176  float pullLocY = -999.;
177  float measurementLocY = -999.;
178  float trackParameterLocY = -999.;
179  float measurementLocCovY = -999.;
180  float trackParameterLocCovY = -999.;
181  bool isAnnulusBound = false;
182 
183  // Get the measurement type
184  if (flag.test(Acts::TrackStateFlag::HoleFlag)) {
185  type = MeasurementType::HOLE;
186  ATH_MSG_DEBUG("--- This is a hole");
187  } else if (flag.test(Acts::TrackStateFlag::OutlierFlag)) {
188  type = MeasurementType::OUTLIER;
189  ATH_MSG_DEBUG("--- This is an outlier");
190  } else {
191  type = MeasurementType::HIT;
192  ATH_MSG_DEBUG("--- This is a hit");
193  }
194 
195  // Check the location of the state
196  if (state.hasReferenceSurface() and state.referenceSurface().associatedDetectorElement()) {
197  const ActsDetectorElement * detectorElement = dynamic_cast<const ActsDetectorElement *>(state.referenceSurface().associatedDetectorElement());
198  if (!detectorElement) {
199  ATH_MSG_WARNING("--- TrackState reference surface returned an invalid associated detector element");
200  continue;
201  }
202  const InDetDD::SiDetectorElement * siliconDetectorElement = dynamic_cast<const InDetDD::SiDetectorElement *>(detectorElement->upstreamDetectorElement());
203  if (!siliconDetectorElement) {
204  ATH_MSG_WARNING("--- TrackState associated detector element is not silicon");
205  continue;
206  }
207 
208  const Acts::AnnulusBounds* annulusBounds = dynamic_cast<const Acts::AnnulusBounds*>(&(state.referenceSurface().bounds()));
209  isAnnulusBound = annulusBounds ? true : false;
210 
211  if (siliconDetectorElement) {
212  Identifier detectorIdentifier = siliconDetectorElement->identify();
213  if (siliconDetectorElement->isPixel()) {
214  const PixelID* pixel_id = static_cast<const PixelID *>(siliconDetectorElement->getIdHelper());
215  int layer_disk = pixel_id->layer_disk(detectorIdentifier);
216  layer = layer_disk;
217  if (pixel_id->barrel_ec(detectorIdentifier) == 0) {
218  if (layer_disk == 0) {
219  detector = Subdetector::INNERMOST_PIXEL;
220  } else detector = Subdetector::PIXEL;
221  region = Region::BARREL;
222  } else {
224  region = Region::ENDCAP;
225  }
226  } else if (siliconDetectorElement->isSCT()) {
227  const SCT_ID* sct_id = static_cast<const SCT_ID *>(siliconDetectorElement->getIdHelper());
229  region = sct_id->barrel_ec(detectorIdentifier) == 0 ?
231  layer = sct_id->layer_disk(detectorIdentifier);;
232  } else ATH_MSG_WARNING("--- Unknown detector type - It is not pixel nor strip detecor element!");
233  } else ATH_MSG_WARNING("--- Missing silicon detector element!");
234  } else ATH_MSG_WARNING("--- Missing reference surface or associated detector element!");
235 
236 
237 
238  // If I have a measurement (hit or outlier) then proceed with computing the residuals / pulls
239 
240  if (type == MeasurementType::OUTLIER || type == MeasurementType::HIT) {
241 
242  //Get the Chi2 computation
243 
244  if (type == MeasurementType::HIT) {
245  chi2_hit_filtered = state.chi2();
246  }
247 
248  if (state.hasUncalibratedSourceLink()) {
249  chi2_hit_predicted = getChi2Contribution(state);
250  }
251 
252  // Skip all states without smoothed parameters or without projector
253  if (!state.hasSmoothed() || !state.hasProjector())
254  continue;
255 
256  // Calling effective Calibrated has some runtime overhead
257  const auto &calibratedParameters = state.effectiveCalibrated();
258  const auto &calibratedCovariance = state.effectiveCalibratedCovariance();
259 
260  // We evaluate the unbiased parameters for:
261  // - measurements added to the fit. For outliers, the measurement is not part of the fit, hence track parameters are already unbiased
262  // - if the filtered parameters and the projector exist.
263  bool evaluateUnbiased = (!flag.test(Acts::TrackStateFlag::OutlierFlag));
264 
265  if (evaluateUnbiased) {
266  ATH_MSG_DEBUG("--- Good for unbiased parameters evaluation!");
267  type = MeasurementType::UNBIASED;
268  // if unbiased, access the associated uncalibrated measurement and store the size
269  if (state.hasUncalibratedSourceLink()) {
270  ATLASUncalibSourceLink sourceLink = state.getUncalibratedSourceLink().template get<ATLASUncalibSourceLink>();
271  const xAOD::UncalibratedMeasurement &uncalibratedMeasurement = getUncalibratedMeasurement(sourceLink);
272  const xAOD::UncalibMeasType measurementType = uncalibratedMeasurement.type();
273  if (measurementType == xAOD::UncalibMeasType::PixelClusterType) {
274  auto pixelCluster = static_cast<const xAOD::PixelCluster *>(&uncalibratedMeasurement);
275  sizePhi = pixelCluster->channelsInPhi();
276  sizeEta = pixelCluster->channelsInEta();
277  } else if (measurementType == xAOD::UncalibMeasType::StripClusterType) {
278  auto stripCluster = static_cast<const xAOD::StripCluster *>(&uncalibratedMeasurement);
279  sizePhi = stripCluster->channelsInPhi();
280  } else {
281  ATH_MSG_DEBUG("xAOD::UncalibratedMeasurement is neither xAOD::PixelCluster nor xAOD::StripCluster");
282  }
283  }
284  }
285 
286  const auto& [unbiasedParameters, unbiasedCovariance] =
287  evaluateUnbiased ? Acts::calculateUnbiasedParametersCovariance(state) : std::make_pair(state.parameters(), state.covariance());
288 
289  measurementLocX = calibratedParameters[Acts::eBoundLoc0];
290  measurementLocCovX = calibratedCovariance(Acts::eBoundLoc0, Acts::eBoundLoc0);
291 
292  if (!isAnnulusBound) {
293 
294  trackParameterLocX = unbiasedParameters[Acts::eBoundLoc0];
295  residualLocX = (measurementLocX - trackParameterLocX) / 1_um; //in um
296  trackParameterLocCovX = unbiasedCovariance(Acts::eBoundLoc0, Acts::eBoundLoc0);
297  }
298  else {
299  // TODO:: use directly phi instead of r*phi in the future
300 
301  float locR = unbiasedParameters[Acts::eBoundLoc0];
302  float covR = unbiasedCovariance(Acts::eBoundLoc0,Acts::eBoundLoc0);
303  float locphi = unbiasedParameters[Acts::eBoundLoc1];
304  float covphi = unbiasedCovariance(Acts::eBoundLoc1,Acts::eBoundLoc1);
305  float covRphi = unbiasedCovariance(Acts::eBoundLoc0,Acts::eBoundLoc1);
306 
307  trackParameterLocX = locphi;
308  residualLocX = locR * (measurementLocX - trackParameterLocX) / 1_um;
309  // Compute the error on the local rphi
310  trackParameterLocCovX = locR*locR*covphi + locphi*locphi*covR + 2*locphi*locR*covRphi;
311 
312  // Rescale the error of the measurement to Rphi.
313  measurementLocCovX = locR*locR * measurementLocCovX;
314  }
315 
316  pullLocX = evaluatePull(residualLocX, measurementLocCovX,
317  trackParameterLocCovX, evaluateUnbiased);
318 
319  if (state.calibratedSize() == 2) {
320  measurementLocY = calibratedParameters[Acts::eBoundLoc1];
321  trackParameterLocY = unbiasedParameters[Acts::eBoundLoc1];
322  residualLocY = (measurementLocY - trackParameterLocY) / 1_um;
323  measurementLocCovY = calibratedCovariance(Acts::eBoundLoc1, Acts::eBoundLoc1);
324  trackParameterLocCovY = unbiasedCovariance(Acts::eBoundLoc1, Acts::eBoundLoc1);
325  pullLocY = evaluatePull(residualLocY, measurementLocCovY,
326  trackParameterLocCovY, evaluateUnbiased);
327  }
328 
329  } // hit or outliers
330 
331  else if (type == MeasurementType::HOLE) {
332 
333  // Get the predicted position on sensor
334  auto pred = state.predicted();
335  trackParameterLocX = pred[Acts::eBoundLoc0];
336  trackParameterLocY = pred[Acts::eBoundLoc1];
337  }
338 
339  // Always fill with this information
340 
341  regions.push_back(region);
342  detectors.push_back(detector);
343  layers.push_back(layer);
344  types.push_back(type);
345  predchi2s.push_back(chi2_hit_predicted);
346  filtchi2s.push_back(chi2_hit_filtered);
347  sizesPhi.push_back(sizePhi);
348  sizesEta.push_back(sizeEta);
349  residualsLocX.push_back(residualLocX);
350  pullsLocX.push_back(pullLocX);
351  measurementsLocX.push_back(measurementLocX);
352  trackParametersLocX.push_back(trackParameterLocX);
353  measurementsLocCovX.push_back(measurementLocCovX);
354  trackParametersLocCovX.push_back(trackParameterLocCovX);
355  residualsLocY.push_back(residualLocY);
356  pullsLocY.push_back(pullLocY);
357  measurementsLocY.push_back(measurementLocY);
358  trackParametersLocY.push_back(trackParameterLocY);
359  measurementsLocCovY.push_back(measurementLocCovY);
360  trackParametersLocCovY.push_back(trackParameterLocCovY);
361 
362 
363  } // loop on states
364 
365  } // loop on tracks
366 
367  return StatusCode::SUCCESS;
368  }

◆ 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 & AthReentrantAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 79 of file AthReentrantAlgorithm.cxx.

80 {
81  // If we didn't find any symlinks to add, just return the collection
82  // from the base class. Otherwise, return the extended collection.
83  if (!m_extendedExtraObjects.empty()) {
85  }
86  return Algorithm::extraOutputDeps();
87 }

◆ filterPassed()

virtual bool AthReentrantAlgorithm::filterPassed ( const EventContext &  ctx) const
inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

135  {
136  return execState( ctx ).filterPassed();
137  }

◆ finalize()

virtual StatusCode ActsTrk::MeasurementToTrackParticleDecoration::finalize ( )
inlineoverridevirtual

Definition at line 44 of file MeasurementToTrackParticleDecoration.h.

44 { return StatusCode::SUCCESS;};

◆ getChi2Contribution()

float ActsTrk::MeasurementToTrackParticleDecoration::getChi2Contribution ( const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &  state) const

Definition at line 370 of file MeasurementToTrackParticleDecoration.cxx.

370  {
371 
372  auto pred = state.predicted();
373  auto predC = state.predictedCovariance();
374 
375  return Acts::visit_measurement(
376  state.calibratedSize(),
377  [&]<std::size_t measdim>(std::integral_constant<std::size_t, measdim>) {
378  Acts::FixedBoundSubspaceHelper<measdim> subspaceHelper =
379  state.template projectorSubspaceHelper<measdim>();
380 
381  // TODO use subspace helper for projection instead
382  auto H = subspaceHelper.projector();
383 
384  const auto calibrated = state.template calibrated<measdim>();
385  const auto calibratedCov = state.template calibratedCovariance<measdim>();
386 
387  auto residual = (H * pred - calibrated).eval();
388  auto rescov = (H * predC * H.transpose() + calibratedCov).eval();
389 
390  return ((residual.transpose() * rescov.inverse() * residual).eval())(0,0);
391  });
392 
393 
394 
395  }

◆ getUnbiasedTrackParameters()

std::pair<Acts::BoundVector, Acts::BoundMatrix> ActsTrk::MeasurementToTrackParticleDecoration::getUnbiasedTrackParameters ( const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &  state,
bool  useSmoothed = true 
) const

◆ initialize()

StatusCode ActsTrk::MeasurementToTrackParticleDecoration::initialize ( )
overridevirtual

◆ 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()

bool AthReentrantAlgorithm::isClonable ( ) const
overridevirtualinherited

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ 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();
383  PBASE::renounce (h);
384  }

◆ 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  {
365  handlesArray.renounce();
366  }

◆ setFilterPassed()

virtual void AthReentrantAlgorithm::setFilterPassed ( bool  state,
const EventContext &  ctx 
) const
inlinevirtualinherited

Definition at line 139 of file AthReentrantAlgorithm.h.

139  {
140  execState( ctx ).setFilterPassed( state );
141  }

◆ sysExecute()

StatusCode AthReentrantAlgorithm::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 67 of file AthReentrantAlgorithm.cxx.

68 {
69  return Gaudi::Algorithm::sysExecute (ctx);
70 }

◆ sysInitialize()

StatusCode AthReentrantAlgorithm::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 InputMakerBase, and HypoBase.

Definition at line 96 of file AthReentrantAlgorithm.cxx.

96  {
98 
99  if (sc.isFailure()) {
100  return sc;
101  }
102 
103  ServiceHandle<ICondSvc> cs("CondSvc",name());
104  for (auto h : outputHandles()) {
105  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
106  // do this inside the loop so we don't create the CondSvc until needed
107  if ( cs.retrieve().isFailure() ) {
108  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
109  return StatusCode::SUCCESS;
110  }
111  if (cs->regHandle(this,*h).isFailure()) {
112  sc = StatusCode::FAILURE;
113  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
114  << " with CondSvc");
115  }
116  }
117  }
118  return sc;
119 }

◆ 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) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_chi2HitFilteredKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_chi2HitFilteredKey
private
Initial value:
{
this, "Chi2HitFilteredKey", m_trackParticlesKey,"chi2_hit_filtered",
"Filtered Chi2 contribution for each hit"}

Definition at line 83 of file MeasurementToTrackParticleDecoration.h.

◆ m_chi2HitPredictedKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_chi2HitPredictedKey
private
Initial value:
{
this, "Chi2HitPredictedKey", m_trackParticlesKey,"chi2_hit_predicted",
"Predicted Chi2 contribution for each hit"}

Definition at line 80 of file MeasurementToTrackParticleDecoration.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_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 AthReentrantAlgorithm::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 153 of file AthReentrantAlgorithm.h.

◆ m_measurementDetectorKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementDetectorKey
private
Initial value:
{
this, "MeasurementDetectorKey", m_trackParticlesKey,"measurement_det",
"Decorate track particle with measurement detector id (innermost pix, pix, strip)"}

Definition at line 74 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementEtaWidthKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementEtaWidthKey
private
Initial value:
{
this, "MeasurementEtaWidthKey", m_trackParticlesKey,"hitResiduals_etaWidth",
"Decorate track particle with measurement cluster size (in eta)"}

Definition at line 92 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementLayerKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLayerKey
private
Initial value:
{
this, "MeasurementLayerKey", m_trackParticlesKey,"measurement_iLayer",
"Decorate track particle with measurement layer"}

Definition at line 77 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementLocCovXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocCovXkey
private
Initial value:
{
this, "MeasurementLocCovXkey", m_trackParticlesKey,"measurementLocCovX",
"Decorate track particle with local x measurement covariance"}

Definition at line 109 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementLocCovYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocCovYkey
private
Initial value:
{
this, "MeasurementLocCovYkey", m_trackParticlesKey,"measurementLocCovY",
"Decorate track particle with local y measurement covariance"}

Definition at line 128 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementLocXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocXkey
private
Initial value:
{
this, "MeasurementLocXkey", m_trackParticlesKey,"measurementLocX",
"Decorate track particle with measurement local x"}

Definition at line 103 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementLocYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocYkey
private
Initial value:
{
this, "MeasurementLocYkey", m_trackParticlesKey,"measurementLocY",
"Decorate track particle with measurement local y"}

Definition at line 122 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementPhiWidthKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementPhiWidthKey
private
Initial value:
{
this, "MeasurementPhiWidthKey", m_trackParticlesKey,"hitResiduals_phiWidth",
"Decorate track particle with measurement cluster size (in r-phi)"}

Definition at line 89 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementRegionKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementRegionKey
private
Initial value:
{
this, "MeasurementRegionKey", m_trackParticlesKey,"measurement_region",
"Decorate track particle with region of the measurement (barrel, ec)"}

Definition at line 71 of file MeasurementToTrackParticleDecoration.h.

◆ m_measurementTypeKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_measurementTypeKey
private
Initial value:
{
this, "MeasurementTypeKey", m_trackParticlesKey,"measurement_type",
"Decorate track particle with type of track state (outlier,hole, biased/unbiased)"}

Definition at line 86 of file MeasurementToTrackParticleDecoration.h.

◆ m_pullLocXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_pullLocXkey
private
Initial value:
{
this, "PullLocXkey", m_trackParticlesKey,"hitResiduals_pullLocX",
"Decorate track particle with unbiased pull in local x"}

Definition at line 100 of file MeasurementToTrackParticleDecoration.h.

◆ m_pullLocYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_pullLocYkey
private
Initial value:
{
this, "PullLocYkey", m_trackParticlesKey,"hitResiduals_pullLocY",
"Decorate track particle with unbiased pull in local y"}

Definition at line 119 of file MeasurementToTrackParticleDecoration.h.

◆ m_residualLocXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_residualLocXkey
private
Initial value:
{
this, "ResidualLocXkey", m_trackParticlesKey,"hitResiduals_residualLocX",
"Decorate track particle with unbiased residual in local x"}

Definition at line 97 of file MeasurementToTrackParticleDecoration.h.

◆ m_residualLocYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_residualLocYkey
private
Initial value:
{
this, "ResidualLocYkey", m_trackParticlesKey,"hitResiduals_residualLocY",
"Decorate track particle with unbiased residual in local y"}

Definition at line 116 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackingGeometryTool

ToolHandle<IActsTrackingGeometryTool> ActsTrk::MeasurementToTrackParticleDecoration::m_trackingGeometryTool {this, "TrackingGeometryTool", ""}
private

Definition at line 66 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackParameterLocCovXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocCovXkey
private
Initial value:
{
this, "TrackParameterLocCovXkey", m_trackParticlesKey,"trackParameterLocCovX",
"Decorate track particle with unbiased local x prediction covariance"}

Definition at line 112 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackParameterLocCovYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocCovYkey
private
Initial value:
{
this, "TrackParameterLocCovYkey", m_trackParticlesKey,"trackParameterLocCovY",
"Decorate track particle with unbiased local y prediction covariance"}

Definition at line 131 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackParameterLocXkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocXkey
private
Initial value:
{
this, "TrackParameterLocXkey", m_trackParticlesKey,"trackParamLocX",
"Decorate track particle with unbiased prediction in local x"}

Definition at line 106 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackParameterLocYkey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocYkey
private
Initial value:
{
this, "TrackParameterLocYkey", m_trackParticlesKey,"trackParamLocY",
"Decorate track particle with unbiased prediction in local y"}

Definition at line 125 of file MeasurementToTrackParticleDecoration.h.

◆ m_trackParticlesKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> ActsTrk::MeasurementToTrackParticleDecoration::m_trackParticlesKey
private
Initial value:
{
this, "TrackParticleKey", "InDetTrackParticles", "Input track particle collection"}

Definition at line 68 of file MeasurementToTrackParticleDecoration.h.

◆ 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:
ActsTrk::MeasurementToTrackParticleDecoration::m_chi2HitFilteredKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_chi2HitFilteredKey
Definition: MeasurementToTrackParticleDecoration.h:83
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocCovYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementLocCovYkey
Definition: MeasurementToTrackParticleDecoration.h:128
SG::ReadHandle::cptr
const_pointer_type cptr()
Dereference the pointer.
ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_trackParameterLocXkey
Definition: MeasurementToTrackParticleDecoration.h:106
ActsTrk::MeasurementToTrackParticleDecoration::INVALID_MEASUREMENT
@ INVALID_MEASUREMENT
Definition: MeasurementToTrackParticleDecoration.h:53
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocCovXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementLocCovXkey
Definition: MeasurementToTrackParticleDecoration.h:109
ClusterSeg::residual
@ residual
Definition: ClusterNtuple.h:20
PixelID::barrel_ec
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: PixelID.h:619
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
xAOD::UncalibMeasType::StripClusterType
@ StripClusterType
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
TRTCalib_cfilter.detector
detector
Definition: TRTCalib_cfilter.py:241
ActsTrk::MeasurementToTrackParticleDecoration::INNERMOST_PIXEL
@ INNERMOST_PIXEL
Definition: MeasurementToTrackParticleDecoration.h:47
ActsTrk::MeasurementToTrackParticleDecoration::m_chi2HitPredictedKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_chi2HitPredictedKey
Definition: MeasurementToTrackParticleDecoration.h:80
AthCommonDataStore::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 > renounce(T &h)
Definition: AthCommonDataStore.h:380
module_driven_slicing.layers
layers
Definition: module_driven_slicing.py:114
ActsDetectorElement::upstreamDetectorElement
const GeoVDetectorElement * upstreamDetectorElement() const
Returns the underllying GeoModel detectorelement that this one is based on.
Definition: ActsDetectorElement.cxx:288
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
ActsTrk::MeasurementToTrackParticleDecoration::evaluatePull
float evaluatePull(const float residual, const float measurementCovariance, const float trackParameterCovariance, const bool evaluateUnbiased) const
Definition: MeasurementToTrackParticleDecoration.cxx:398
athena.value
value
Definition: athena.py:124
ActsTrk::MeasurementToTrackParticleDecoration::INVALID_DETECTOR
@ INVALID_DETECTOR
Definition: MeasurementToTrackParticleDecoration.h:47
RunPrintSiDetElements.detectors
detectors
Definition: RunPrintSiDetElements.py:59
ActsTrk::MeasurementToTrackParticleDecoration::BIASED
@ BIASED
Definition: MeasurementToTrackParticleDecoration.h:53
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:55
SCT_ID::barrel_ec
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: SCT_ID.h:728
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementLocXkey
Definition: MeasurementToTrackParticleDecoration.h:103
Trk::locR
@ locR
Definition: ParamDefs.h:44
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ActsTrk::MeasurementToTrackParticleDecoration::PIXEL
@ PIXEL
Definition: MeasurementToTrackParticleDecoration.h:47
ActsTrk::MeasurementToTrackParticleDecoration::OUTLIER
@ OUTLIER
Definition: MeasurementToTrackParticleDecoration.h:53
ActsTrk::MeasurementToTrackParticleDecoration::m_trackingGeometryTool
ToolHandle< IActsTrackingGeometryTool > m_trackingGeometryTool
Definition: MeasurementToTrackParticleDecoration.h:66
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementDetectorKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementDetectorKey
Definition: MeasurementToTrackParticleDecoration.h:74
InDetDD::SolidStateDetectorElementBase::getIdHelper
const AtlasDetectorID * getIdHelper() const
Returns the id helper (inline)
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
Trk::stripCluster
@ stripCluster
Definition: MeasurementType.h:23
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
SG::makeHandle
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Definition: ReadCondHandle.h:270
ActsTrk::MeasurementToTrackParticleDecoration::m_pullLocXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_pullLocXkey
Definition: MeasurementToTrackParticleDecoration.h:100
ENDCAP
@ ENDCAP
Definition: TRTRadiatorParameters.h:10
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
xAOD::UncalibratedMeasurement_v1
Definition: UncalibratedMeasurement_v1.h:13
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthReentrantAlgorithm::AthReentrantAlgorithm
AthReentrantAlgorithm()
Default constructor:
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
xAOD::UncalibratedMeasurement_v1::type
virtual xAOD::UncalibMeasType type() const =0
Returns the type of the measurement type as a simple enumeration.
ActsTrk::ATLASUncalibSourceLink
const xAOD::UncalibratedMeasurement * ATLASUncalibSourceLink
Definition: ATLASSourceLink.h:24
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
master.flag
bool flag
Definition: master.py:29
ActsTrk::MeasurementToTrackParticleDecoration::HIT
@ HIT
Definition: MeasurementToTrackParticleDecoration.h:53
ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocCovXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_trackParameterLocCovXkey
Definition: MeasurementToTrackParticleDecoration.h:112
xAOD::StripCluster_v1
Definition: StripCluster_v1.h:17
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
ActsDetectorElement
Definition: ActsDetectorElement.h:42
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
ActsTrk::MeasurementToTrackParticleDecoration::BARREL
@ BARREL
Definition: MeasurementToTrackParticleDecoration.h:50
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementTypeKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementTypeKey
Definition: MeasurementToTrackParticleDecoration.h:86
SG::ReadHandle::isValid
virtual bool isValid() override final
Can the handle be successfully dereferenced?
AthReentrantAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Definition: AthReentrantAlgorithm.h:153
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
InDetDD::SiDetectorElement::isPixel
bool isPixel() const
PixelID::layer_disk
int layer_disk(const Identifier &id) const
Definition: PixelID.h:626
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementRegionKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementRegionKey
Definition: MeasurementToTrackParticleDecoration.h:71
python.root_lsr_rank.types
types
Definition: root_lsr_rank.py:35
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocCovYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_trackParameterLocCovYkey
Definition: MeasurementToTrackParticleDecoration.h:131
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
ActsTrk::MeasurementToTrackParticleDecoration::m_residualLocXkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_residualLocXkey
Definition: MeasurementToTrackParticleDecoration.h:97
AthReentrantAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthReentrantAlgorithm.cxx:96
ActsTrk::getUncalibratedMeasurement
const xAOD::UncalibratedMeasurement & getUncalibratedMeasurement(const ATLASUncalibSourceLink &source_link)
Definition: ATLASSourceLink.h:26
SCT_ID::layer_disk
int layer_disk(const Identifier &id) const
Definition: SCT_ID.h:734
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
ActsTrk::MeasurementToTrackParticleDecoration::m_pullLocYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_pullLocYkey
Definition: MeasurementToTrackParticleDecoration.h:119
Trk::pixelCluster
@ pixelCluster
Definition: MeasurementType.h:22
a
TList * a
Definition: liststreamerinfos.cxx:10
xAOD::PixelCluster_v1
Definition: PixelCluster_v1.h:17
SG::WriteDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
h
SCT_ID
Definition: SCT_ID.h:68
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
LVL1TGCTrigger::STRIP
@ STRIP
Definition: TGCNumbering.h:52
ActsTrk::MeasurementToTrackParticleDecoration::UNBIASED
@ UNBIASED
Definition: MeasurementToTrackParticleDecoration.h:53
ActsTrk::MeasurementToTrackParticleDecoration::m_trackParameterLocYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_trackParameterLocYkey
Definition: MeasurementToTrackParticleDecoration.h:125
xAOD::UncalibMeasType
UncalibMeasType
Define the type of the uncalibrated measurement.
Definition: MeasurementDefs.h:24
ActsTrk::MeasurementToTrackParticleDecoration::m_residualLocYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_residualLocYkey
Definition: MeasurementToTrackParticleDecoration.h:116
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
BARREL
@ BARREL
Definition: TRTRadiatorParameters.h:10
physics_parameters.parameters
parameters
Definition: physics_parameters.py:144
ActsTrk::MeasurementToTrackParticleDecoration::getChi2Contribution
float getChi2Contribution(const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &state) const
Definition: MeasurementToTrackParticleDecoration.cxx:370
ActsTrk::MeasurementToTrackParticleDecoration::INVALID_REGION
@ INVALID_REGION
Definition: MeasurementToTrackParticleDecoration.h:50
InDetDD::SiDetectorElement::isSCT
bool isSCT() const
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
ActsTrk::MeasurementToTrackParticleDecoration::STRIP
@ STRIP
Definition: MeasurementToTrackParticleDecoration.h:47
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementPhiWidthKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementPhiWidthKey
Definition: MeasurementToTrackParticleDecoration.h:89
xAOD::track
@ track
Definition: TrackingPrimitives.h:512
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLayerKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementLayerKey
Definition: MeasurementToTrackParticleDecoration.h:77
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementLocYkey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementLocYkey
Definition: MeasurementToTrackParticleDecoration.h:122
ActsTrk::MeasurementToTrackParticleDecoration::m_trackParticlesKey
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticlesKey
Definition: MeasurementToTrackParticleDecoration.h:68
ActsTrk::MeasurementToTrackParticleDecoration::m_measurementEtaWidthKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_measurementEtaWidthKey
Definition: MeasurementToTrackParticleDecoration.h:92
PixelID
Definition: PixelID.h:67
ActsTrk::MeasurementToTrackParticleDecoration::ENDCAP
@ ENDCAP
Definition: MeasurementToTrackParticleDecoration.h:50
ActsTrk::MeasurementToTrackParticleDecoration::HOLE
@ HOLE
Definition: MeasurementToTrackParticleDecoration.h:53
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
ActsTrk::MeasurementToTrackParticleDecoration::N_SUBDETECTORS
@ N_SUBDETECTORS
Definition: MeasurementToTrackParticleDecoration.h:47
InDetDD::SolidStateDetectorElementBase::identify
virtual Identifier identify() const override final
identifier of this detector element (inline)
xAOD::UncalibMeasType::PixelClusterType
@ PixelClusterType
fitman.k
k
Definition: fitman.py:528
ServiceHandle< ICondSvc >
PIXEL
@ PIXEL
Definition: RegSelEnums.h:24
Identifier
Definition: IdentifierFieldParser.cxx:14