Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Namespaces | Classes | Typedefs | Functions
ActsTrk::detail Namespace Reference

Namespaces

 FitterHelperFunctions
 
 HitCategory
 
 MakeDerivedVariant
 

Classes

struct  accepted_decoration_types
 
class  AnalogueClusteringToolImpl
 
class  AtlasMeasurementContainerList
 
struct  CKF_config
 
struct  Decoration
 
class  DuplicateSeedDetector
 
class  GenMeasurementRangeList
 
class  GenUncalibSourceLinkAccessor
 Accessor for the above source link container. More...
 
class  HitSummaryData
 Helper class to gather hit summary information for e.g. More...
 
class  MeasurementCalibratorBase
 
class  MeasurementIndex
 
struct  MeasurementRange
 
class  OnTrackCalibrator
 
class  RefittingCalibrator
 
class  SharedHitCounter
 
class  SumOfValues
 Helper class to gather statistics and compute the biased variance. More...
 
class  TrackFindingMeasurements
 
class  TrkMeasurementCalibrator
 

Typedefs

using SetterType = std::function< std::any(SG::IAuxStore *, ActsTrk::IndexType, SG::auxid_t)>
 
using GetterType = std::function< const std::any(const SG::IConstAuxStore *, ActsTrk::IndexType, SG::auxid_t)>
 
using CopierType = std::function< void(SG::IAuxStore *, ActsTrk::IndexType, SG::auxid_t, const std::any &)>
 
using MeasurementRangeList = GenMeasurementRangeList< AtlasMeasurementContainerList >
 
using UncalibSourceLinkAccessor = GenUncalibSourceLinkAccessor< MeasurementRangeList >
 
using RecoTrackContainer = Acts::TrackContainer< Acts::VectorTrackContainer, Acts::VectorMultiTrajectory >
 
using RecoTrackContainerProxy = RecoTrackContainer::TrackProxy
 
using RecoTrackStateContainer = Acts::VectorMultiTrajectory
 
using RecoTrackStateContainerProxy = RecoTrackStateContainer::TrackStateProxy
 
using RecoConstTrackStateContainerProxy = RecoTrackStateContainer::ConstTrackStateProxy
 
using Stepper = Acts::SympyStepper
 Adapted from Acts Examples/Algorithms/TrackFinding/src/TrackFindingAlgorithmFunction.cpp. More...
 
using Navigator = Acts::Navigator
 
using Propagator = Acts::Propagator< Stepper, Navigator >
 
using CKF = Acts::CombinatorialKalmanFilter< Propagator, RecoTrackContainer >
 
using Extrapolator = Propagator
 

Functions

template<typename T >
const std::any constDecorationGetter (const SG::IConstAuxStore *container, ActsTrk::IndexType idx, SG::auxid_t decorationId)
 
template<typename T >
const std::any decorationGetter (const SG::IAuxStore *container, ActsTrk::IndexType idx, SG::auxid_t decorationId)
 
template<typename T >
std::any decorationSetter (SG::IAuxStore *container, ActsTrk::IndexType idx, SG::auxid_t decorationId)
 
template<typename T >
void decorationCopier (SG::IAuxStore *dst, ActsTrk::IndexType dst_idx, SG::auxid_t decorationId, const std::any &src_ptr)
 
std::vector< Decoration > restoreDecorations (const SG::IConstAuxStore *container, const std::set< std::string > &staticVariables)
 
template<typename T >
bool build (const std::type_info *typeInfo, const std::string &name, std::vector< Decoration > &decorations)
 
std::tuple< std::vector< int >, std::vector< std::vector< int > >, std::vector< std::vector< float > > > getSDOInformation (const std::vector< Identifier > &rdoList, const InDetSimDataCollection &sdoCollection)
 
std::tuple< std::vector< float >, std::vector< float >, std::vector< int >, std::vector< int >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float > > getSiHitInformation (const InDetDD::SiDetectorElement &element, const std::vector< SiHit > &matchingHits)
 
std::unique_ptr< ActsTrk::IMeasurementSelector > getMeasurementSelector (const ActsTrk::IOnBoundStateCalibratorTool *onTrackCalibratorTool, const ActsTrk::detail::MeasurementRangeList &measurementRanges, const std::vector< float > &etaBinsf, const std::vector< std::pair< float, float > > &chi2CutOffOutlier, const std::vector< size_t > &numMeasurementsCutOff)
 
double computeDtDs (const Acts::ParticleHypothesis &hypothesis, double qop)
 Helper function to compute dt/ds Helper function to compute the derivative of the time as function of the path length. More...
 
Acts::FreeToPathMatrix computeFreeToPathDerivatives (const Acts::Vector3 &direction, double qop, const Acts::Vector3 &bfield, const Acts::ParticleHypothesis &particle_hypothesis)
 Compute the path length derivatives for the free/bound to curvilinear paramter transform. More...
 
std::optional< Acts::BoundMatrix > convertActsBoundCovToCurvilinearParam (const Acts::GeometryContext &tgContext, const Acts::BoundTrackParameters &param, const Acts::Vector3 &magnFieldVect, const Acts::ParticleHypothesis &particle_hypothesis)
 Convert the covariance of the given Acts track parameters into curvilinear parameterisation. More...
 
std::array< unsigned int, 4 > expectedLayerPattern (const EventContext &ctx, const IActsExtrapolationTool &extrapolator, Acts::BoundTrackParameters perigee_parameters, double pathLimit)
 Extrapolate from the perigee outwards and gather information which detector layers should have hits. More...
 
void gatherTrackSummaryData (const ActsTrk::TrackContainer &tracksContainer, const typename ActsTrk::TrackContainer::ConstTrackProxy &track, const std::array< const InDetDD::SiDetectorElementCollection *, to_underlying(xAOD::UncalibMeasType::nTypes)> &siDetEleColl, const std::array< unsigned short, to_underlying(xAOD::UncalibMeasType::nTypes)> &measurement_to_summary_type, SumOfValues &chi2_stat_out, HitSummaryData &hit_info_out, std::vector< ActsTrk::TrackStateBackend::ConstTrackStateProxy::IndexType > &param_state_idx_out, std::array< std::array< uint8_t, to_underlying(HitCategory::N)>, to_underlying(xAOD::UncalibMeasType::nTypes)> &special_hit_counts_out)
 Helper to gather track summary information from the track states of the specified track. More...
 
template<typename T_EnumClass >
constexpr std::underlying_type< T_EnumClass >::type to_underlying (T_EnumClass an_enum)
 Helper to convert class enum into an integer. More...
 

Typedef Documentation

◆ CKF

using ActsTrk::detail::CKF = typedef Acts::CombinatorialKalmanFilter<Propagator, RecoTrackContainer>

Definition at line 33 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ CopierType

using ActsTrk::detail::CopierType = typedef std::function<void(SG::IAuxStore*, ActsTrk::IndexType, SG::auxid_t, const std::any&)>

Definition at line 20 of file Decoration.h.

◆ Extrapolator

using ActsTrk::detail::Extrapolator = typedef Propagator

Definition at line 34 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ GetterType

using ActsTrk::detail::GetterType = typedef std::function<const std::any( const SG::IConstAuxStore*, ActsTrk::IndexType, SG::auxid_t)>

Definition at line 18 of file Decoration.h.

◆ MeasurementRangeList

using ActsTrk::detail::MeasurementRangeList = typedef GenMeasurementRangeList< AtlasMeasurementContainerList >

Definition at line 212 of file AtlasUncalibSourceLinkAccessor.h.

◆ Navigator

using ActsTrk::detail::Navigator = typedef Acts::Navigator

Definition at line 31 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ Propagator

using ActsTrk::detail::Propagator = typedef Acts::Propagator<Stepper, Navigator>

Definition at line 32 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ RecoConstTrackStateContainerProxy

using ActsTrk::detail::RecoConstTrackStateContainerProxy = typedef RecoTrackStateContainer::ConstTrackStateProxy

Definition at line 26 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ RecoTrackContainer

using ActsTrk::detail::RecoTrackContainer = typedef Acts::TrackContainer<Acts::VectorTrackContainer, Acts::VectorMultiTrajectory>

Definition at line 21 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ RecoTrackContainerProxy

using ActsTrk::detail::RecoTrackContainerProxy = typedef RecoTrackContainer::TrackProxy

Definition at line 23 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ RecoTrackStateContainer

using ActsTrk::detail::RecoTrackStateContainer = typedef Acts::VectorMultiTrajectory

Definition at line 24 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ RecoTrackStateContainerProxy

using ActsTrk::detail::RecoTrackStateContainerProxy = typedef RecoTrackStateContainer::TrackStateProxy

Definition at line 25 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ SetterType

using ActsTrk::detail::SetterType = typedef std::function<std::any(SG::IAuxStore*, ActsTrk::IndexType, SG::auxid_t)>

Definition at line 16 of file Decoration.h.

◆ Stepper

using ActsTrk::detail::Stepper = typedef Acts::SympyStepper

Adapted from Acts Examples/Algorithms/TrackFinding/src/TrackFindingAlgorithmFunction.cpp.

Definition at line 30 of file Tracking/Acts/ActsTrackReconstruction/src/detail/Definitions.h.

◆ UncalibSourceLinkAccessor

using ActsTrk::detail::UncalibSourceLinkAccessor = typedef GenUncalibSourceLinkAccessor< MeasurementRangeList >

Definition at line 213 of file AtlasUncalibSourceLinkAccessor.h.

Function Documentation

◆ build()

template<typename T >
bool ActsTrk::detail::build ( const std::type_info *  typeInfo,
const std::string &  name,
std::vector< Decoration > &  decorations 
)

Definition at line 11 of file Decoration.cxx.

11  {
12  if (*typeInfo == typeid(T)) {
13  decorations.emplace_back(
14  decoration<T>(name, ActsTrk::detail::constDecorationGetter<T>,
15  ActsTrk::detail::decorationCopier<T>));
16  return true;
17  }
18  return false;
19 }

◆ computeDtDs()

double ActsTrk::detail::computeDtDs ( const Acts::ParticleHypothesis &  hypothesis,
double  qop 
)
inline

Helper function to compute dt/ds Helper function to compute the derivative of the time as function of the path length.

Definition at line 16 of file CurvilinearCovarianceHelper.h.

16  {
17  return std::hypot(1.,hypothesis.mass() / hypothesis.extractMomentum(qop));
18  }

◆ computeFreeToPathDerivatives()

Acts::FreeToPathMatrix ActsTrk::detail::computeFreeToPathDerivatives ( const Acts::Vector3 &  direction,
double  qop,
const Acts::Vector3 &  bfield,
const Acts::ParticleHypothesis &  particle_hypothesis 
)
inline

Compute the path length derivatives for the free/bound to curvilinear paramter transform.

Parameters
directionthe direction of trajectory at the location in question
qopq/p of the particle at the location in question in Acts units
bfieldthe magnetic field at the location in question in Acts units
particle_hypothesisthe particle hypothesis e.g. Acts::ParticleHypothesis::pion()
Returns
path length derivatives ( dr(...)/ds, dt/ds, dr(...)/ds2, d qop/ds [== 0] )

Definition at line 27 of file CurvilinearCovarianceHelper.h.

31  {
32  Acts::FreeToPathMatrix path_length_deriv;
33  static_assert(path_length_deriv.cols() == 8); // ensure that all elements are initialized
34  path_length_deriv.segment<3>(Acts::eFreePos0) = direction;
35  path_length_deriv(0,Acts::eFreeTime) = computeDtDs(particle_hypothesis,qop);
36  path_length_deriv.segment<3>(Acts::eFreeDir0) = (qop * direction.cross(bfield)).transpose();
37  path_length_deriv(0,Acts::eFreeQOverP) = 0.;
38  return path_length_deriv;
39  }

◆ constDecorationGetter()

template<typename T >
const std::any ActsTrk::detail::constDecorationGetter ( const SG::IConstAuxStore *  container,
ActsTrk::IndexType  idx,
SG::auxid_t  decorationId 
)

Definition at line 49 of file Decoration.h.

51  {
52  const void* data = container->getData(decorationId);
53  return &(static_cast<const T*>(data)[idx]);
54 }

◆ convertActsBoundCovToCurvilinearParam()

std::optional<Acts::BoundMatrix> ActsTrk::detail::convertActsBoundCovToCurvilinearParam ( const Acts::GeometryContext &  tgContext,
const Acts::BoundTrackParameters &  param,
const Acts::Vector3 &  magnFieldVect,
const Acts::ParticleHypothesis &  particle_hypothesis 
)
inline

Convert the covariance of the given Acts track parameters into curvilinear parameterisation.

Parameters
tgContextthe current geometry context
paramthe Acts track parameters
magnFieldVectthe vector of the magnetic field at the position of the track parameters in Acts units
particle_hypothesisthe particle hypothesis e.g. Acts::ParticleHypothesis::pion() Usage: 
  
     using namespace Acts::UnitLiterals;
     MagField::AtlasFieldCache fieldCache;
     ...
     Amg::Vector3D magnFieldVect;
     fieldCache.getField(param.position(tgContext).data(), magnFieldVect.data());
     magnFieldVect *= 1000_T;
     auto cov = convertActsBoundCovToCurvilinearParam(tgContext, param, magnFieldVect);
     if (cov.has_value()) { ... }

Definition at line 61 of file CurvilinearCovarianceHelper.h.

65  {
66  if (param.covariance().has_value()) {
67  Acts::FreeVector freeParams = Acts::transformBoundToFreeParameters(
68  param.referenceSurface(), tgContext, param.parameters());
69  Acts::Vector3 position = freeParams.segment<3>(Acts::eFreePos0);
70  Acts::Vector3 direction = freeParams.segment<3>(Acts::eFreeDir0);
71  Acts::BoundMatrix b2c;
72  Acts::detail::boundToCurvilinearTransportJacobian(direction, // magnFieldVect.normalized(),
73  param.referenceSurface().boundToFreeJacobian(tgContext, position, direction),
74  Acts::FreeMatrix::Identity(),
75  computeFreeToPathDerivatives(direction,
76  param.parameters()[Acts::eBoundQOverP],
77  magnFieldVect,
78  particle_hypothesis),
79  b2c);
80 
81  return b2c * param.covariance().value() * b2c.transpose();
82  }
83  else {
84  return {};
85  }
86  }

◆ decorationCopier()

template<typename T >
void ActsTrk::detail::decorationCopier ( SG::IAuxStore *  dst,
ActsTrk::IndexType  dst_idx,
SG::auxid_t  decorationId,
const std::any &  src_ptr 
)

Definition at line 76 of file Decoration.h.

77  {
78  *std::any_cast<T*>(decorationSetter<T>(dst, dst_idx, decorationId)) =
79  *std::any_cast<const T*>(src_ptr);
80 }

◆ decorationGetter()

template<typename T >
const std::any ActsTrk::detail::decorationGetter ( const SG::IAuxStore *  container,
ActsTrk::IndexType  idx,
SG::auxid_t  decorationId 
)

Definition at line 57 of file Decoration.h.

59  {
60  const void* data = container->getData(decorationId);
61  return &(static_cast<T*>(data)[idx]);
62 }

◆ decorationSetter()

template<typename T >
std::any ActsTrk::detail::decorationSetter ( SG::IAuxStore *  container,
ActsTrk::IndexType  idx,
SG::auxid_t  decorationId 
)

Definition at line 67 of file Decoration.h.

68  {
69  assert (idx < container->size());
70  // The size requested for the decoration must match the size of the container.
71  void* data = container->getData(decorationId, container->size(), container->size());
72  return &(static_cast<T*>(data)[idx]);
73 }

◆ expectedLayerPattern()

std::array< unsigned int, 4 > ActsTrk::detail::expectedLayerPattern ( const EventContext &  ctx,
const IActsExtrapolationTool &  extrapolator,
Acts::BoundTrackParameters  perigee_parameters,
double  pathLimit 
)

Extrapolate from the perigee outwards and gather information which detector layers should have hits.

Parameters
ctxthe current athena EvetContext @paran extrapolator referece to the Acts extrapolation tool
perigee_parametersthe Acts defining track parameters (parameters at the perigee).
pathLimitthe maximum path length to extrapolate to.
Returns
array of layer pattern words: two pairs for pixel and strips containing the barrel and end caps bit pattern with one bit per detector layer. Will extrapolate from the perigee to the given cylinder volume surface and gather the layers which are crossed and on which hits are expected.

Definition at line 14 of file ExpectedHitUtils.cxx.

17  {
18  ActsPropagationOutput result = extrapolator.propagationSteps(ctx,
19  perigee_parameters,
20  Acts::Direction::Forward(),
21  pathLimit);
22  std::array<unsigned int,4> expected_layer_pattern {0u,0u,0u,0u};
23  for (const Acts::detail::Step &step : result.first ) {
24  // @TODO boundary check ?, does layer number match layer numbering in athena ?
25  // @TODO filter out dead modules
26  if (step.geoID.sensitive()) {
27  if (step.surface && step.surface->associatedDetectorElement()) {
28  const ActsDetectorElement *
29  actsDetEl = dynamic_cast<const ActsDetectorElement *>(step.surface->associatedDetectorElement());
30  if (actsDetEl) {
31  const InDetDD::SiDetectorElement *
32  detEl = dynamic_cast<const InDetDD::SiDetectorElement *>(actsDetEl->upstreamDetectorElement());
33  if (detEl) {
34  if (detEl->isPixel()) {
35  InDetDD::DetectorType type = detEl->design().type();
36  unsigned int pattern_idx = type==InDetDD::PixelBarrel ? 0 : 1;
37  const PixelID* pixel_id = static_cast<const PixelID *>(detEl->getIdHelper());
38  expected_layer_pattern[pattern_idx] |= (1<< pixel_id->layer_disk(detEl->identify()) );
39  }
40  else if (detEl->isSCT()) {
41  // @TODO SCT/strip layer number should be shifted by total number of pixel layers ...
42  const SCT_ID* sct_id = static_cast<const SCT_ID *>(detEl->getIdHelper());
43  unsigned int pattern_idx = detEl->isBarrel() ? 2 : 3;
44  expected_layer_pattern[pattern_idx] |= (1<< sct_id->layer_disk(detEl->identify()) );
45  }
46  }
47  }
48  }
49  }
50  }
51  return expected_layer_pattern;
52  }

◆ gatherTrackSummaryData()

void ActsTrk::detail::gatherTrackSummaryData ( const ActsTrk::TrackContainer &  tracksContainer,
const typename ActsTrk::TrackContainer::ConstTrackProxy &  track,
const std::array< const InDetDD::SiDetectorElementCollection *, to_underlying(xAOD::UncalibMeasType::nTypes)> &  siDetEleColl,
const std::array< unsigned short, to_underlying(xAOD::UncalibMeasType::nTypes)> &  measurement_to_summary_type,
SumOfValues &  chi2_stat_out,
HitSummaryData &  hit_info_out,
std::vector< ActsTrk::TrackStateBackend::ConstTrackStateProxy::IndexType > &  param_state_idx_out,
std::array< std::array< uint8_t, to_underlying(HitCategory::N)>, to_underlying(xAOD::UncalibMeasType::nTypes)> &  special_hit_counts_out 
)

Helper to gather track summary information from the track states of the specified track.

Parameters
tracksContainerthe Acts track container
tracka track of the given acts track container for which the summary information is to be gathered
siDetEleCollarray of SiDetectorElement collections per measurement type.
measurement_to_summary_typea LUT to map measurement types to the corresponding summary type for the measurement counts
chi2_stat_outoutput of the per track state chi-squared sums and squared sums to compute the per state chi2 variance.
hit_info_outoutput of the gathered measurement statistics per detector region, layer, ... .
param_state_idx_outoutput vector to be filled with the state index of all track states which are not holes.
special_hit_counts_outarrays to count holes (and @TODO dead sensors) per measurement type.

Definition at line 14 of file HitSummaryDataUtils.cxx.

25  {
26  chi2_stat_out.reset();
27 
28  hit_info_out.reset();
29  param_state_idx_out.clear();
30 
31  const auto lastMeasurementIndex = track.tipIndex();
32  tracksContainer.trackStateContainer().visitBackwards(
33  lastMeasurementIndex,
34  [&siDetEleColl,
35  &measurement_to_summary_type,
36  &chi2_stat_out,
37  &hit_info_out,
38  &param_state_idx_out,
39  &special_hit_counts_out
40  ](const typename ActsTrk::TrackStateBackend::ConstTrackStateProxy &state) -> void
41  {
42 
43  auto flag = state.typeFlags();
44  if (flag.test(Acts::TrackStateFlag::HoleFlag)) {
45  xAOD::UncalibMeasType det_type = xAOD::UncalibMeasType::Other;
46  if (state.hasReferenceSurface()) {
47  if (state.referenceSurface().associatedDetectorElement()) {
48  const ActsDetectorElement *
49  actsDetEl = dynamic_cast<const ActsDetectorElement *>(state.referenceSurface().associatedDetectorElement());
50  if (actsDetEl) {
51  const InDetDD::SiDetectorElement *
52  detEl = dynamic_cast<const InDetDD::SiDetectorElement *>(actsDetEl->upstreamDetectorElement());
53  if (detEl) {
54  if (detEl->isPixel()) {
55  det_type = xAOD::UncalibMeasType::PixelClusterType;
56  }
57  else if (detEl->isSCT()) {
58  det_type = xAOD::UncalibMeasType::StripClusterType;
59  }
60  }
61  }
62  }
63  }
64 
65  Acts::BoundaryTolerance tolerance = Acts::BoundaryTolerance::AbsoluteBound(
66  // local0 tolerance
67  (det_type == xAOD::UncalibMeasType::PixelClusterType
68  || state.referenceSurface().bounds().type() != Acts::SurfaceBounds::eAnnulus)
69  ? 0.0 : std::numeric_limits<double>::infinity(),
70  // local1 tolerance
71  (det_type == xAOD::UncalibMeasType::PixelClusterType
72  || state.referenceSurface().bounds().type() == Acts::SurfaceBounds::eAnnulus)
73  ? 0.0 : std::numeric_limits<double>::infinity()
74  );
75  Acts::Vector2 localPos(state.predicted()[Acts::eBoundLoc0],state.predicted()[Acts::eBoundLoc1]);
76  if (state.referenceSurface().insideBounds(localPos,tolerance)) {
77  // @TODO check whether detector element is dead..
78  // if (dead) {
79  // ++specialHitCounts.at(to_underlying(det_type)).at(HitCategory::DeadSensor);
80  // } else {
81  ++special_hit_counts_out.at(to_underlying(det_type)).at(HitCategory::Hole);
82  }
83 
84  }
85  else if (flag.test(Acts::TrackStateFlag::MeasurementFlag)) {
86  // do not consider material states
87  param_state_idx_out.push_back(state.index());
88  }
89  // @TODO dead elements
90 
91  if (state.hasUncalibratedSourceLink()) {
92  auto sl = state.getUncalibratedSourceLink().template get<ActsTrk::ATLASUncalibSourceLink>();
93  assert( sl != nullptr );
94  const xAOD::UncalibratedMeasurement &uncalibMeas = getUncalibratedMeasurement(sl);
95  if (measurement_to_summary_type.at(to_underlying(uncalibMeas.type())) < xAOD::numberOfTrackSummaryTypes ) {
96  if (static_cast<unsigned int>(to_underlying(uncalibMeas.type())) < siDetEleColl.size()) {
97  HitSummaryData::EHitSelection hit_selection = (flag.test(Acts::TrackStateFlag::OutlierFlag)
98  ? HitSummaryData::Outlier
99  : HitSummaryData::Hit);
100  if (flag.test(Acts::TrackStateFlag::SharedHitFlag)) {
101  hit_selection = HitSummaryData::EHitSelection(hit_selection | HitSummaryData::SharedHit);
102  }
103  hit_info_out.addHit(siDetEleColl[to_underlying(uncalibMeas.type())],
104  uncalibMeas.identifierHash(),
105  hit_selection);
106  }
107  }
108  if (state.calibratedSize()>0 && !flag.test(Acts::TrackStateFlag::OutlierFlag)) {
109  // from Tracking/TrkTools/TrkTrackSummaryTool/src/TrackSummaryTool.cxx
110  // processTrackState
111  double chi2add = std::min(state.chi2(),1e5f) / state.calibratedSize();
112  chi2_stat_out.add(chi2add );
113  }
114  }
115 
116  });
117  hit_info_out.computeSummaries();
118  }

◆ getMeasurementSelector()

std::unique_ptr< ActsTrk::IMeasurementSelector > ActsTrk::detail::getMeasurementSelector ( const ActsTrk::IOnBoundStateCalibratorTool *  onTrackCalibratorTool,
const ActsTrk::detail::MeasurementRangeList &  measurementRanges,
const std::vector< float > &  etaBinsf,
const std::vector< std::pair< float, float > > &  chi2CutOffOutlier,
const std::vector< size_t > &  numMeasurementsCutOff 
)

Definition at line 385 of file AtlasMeasurementSelector.cxx.

389  {
390 
391  // set calibrators per measurement container type (order does not matter);
392  ActsTrk::MeasurementCalibrator2 atl_measurement_calibrator(onTrackCalibratorTool);
393  using AtlMeasurementSelectorCuts = AtlasMeasurementSelectorCuts;
394 
395  using AtlMeasurementSelector = AtlasActsMeasurmentSelector<RecoTrackContainer>;
396  using AtlMeasurementSelectorConfig = AtlMeasurementSelector::TheAtlasMeasurementSelector::Config;
397 
398  std::unique_ptr<ActsTrk::IMeasurementSelector>
399  selector(new AtlMeasurementSelector(
400  std::move(atl_measurement_calibrator),
401  measurementRanges,
402  AtlMeasurementSelectorConfig{ {Acts::GeometryIdentifier(),
403  AtlMeasurementSelectorCuts{ etaBinsf,
404  chi2CutOffOutlier,
405  numMeasurementsCutOff} }}));
406  return selector;
407 
408 }

◆ getSDOInformation()

std::tuple< std::vector< int >, std::vector< std::vector< int > >, std::vector< std::vector< float > > > ActsTrk::detail::getSDOInformation ( const std::vector< Identifier > &  rdoList,
const InDetSimDataCollection &  sdoCollection 
)

Definition at line 13 of file Tracking/Acts/ActsObjectDecoration/src/detail/Utilities.cxx.

15  {
16  std::vector<int> sdo_word {};
17  std::vector< std::vector< int > > sdo_depositsBarcode {};
18  std::vector< std::vector< float > > sdo_depositsEnergy {};
19 
20  for (const Identifier hitIdentifier : rdoList) {
21  auto pos = sdoCollection.find(hitIdentifier);
22  if( pos == sdoCollection.end() ) continue;
23 
24  sdo_word.push_back( pos->second.word() ) ;
25 
26  std::vector<int> sdoDepBC(pos->second.getdeposits().size(), HepMC::INVALID_PARTICLE_ID);
27  std::vector<float> sdoDepEnergy(pos->second.getdeposits().size());
28 
29  unsigned int nDepos {0};
30  for (const auto& deposit: pos->second.getdeposits()) {
31  if (deposit.first) sdoDepBC[nDepos] = HepMC::barcode(deposit.first);
32  sdoDepEnergy[nDepos] = deposit.second;
33  ++nDepos;
34  }
35 
36  sdo_depositsBarcode.push_back( std::move(sdoDepBC) );
37  sdo_depositsEnergy.push_back( std::move(sdoDepEnergy) );
38  }
39 
40  return std::make_tuple(std::move(sdo_word),
41  std::move(sdo_depositsBarcode),
42  std::move(sdo_depositsEnergy));
43  }

◆ getSiHitInformation()

std::tuple< std::vector< float >, std::vector< float >, std::vector< int >, std::vector< int >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float >, std::vector< float > > ActsTrk::detail::getSiHitInformation ( const InDetDD::SiDetectorElement &  element,
const std::vector< SiHit > &  matchingHits 
)

Definition at line 57 of file Tracking/Acts/ActsObjectDecoration/src/detail/Utilities.cxx.

59  {
60  int numHits = matchingHits.size();
61 
62  std::vector<float> sihit_energyDeposit(numHits, 0);
63  std::vector<float> sihit_meanTime(numHits, 0);
64  std::vector<int> sihit_barcode(numHits, 0);
65  std::vector<int> sihit_pdgid(numHits, 0);
66 
67  std::vector<float> sihit_startPosX(numHits, 0);
68  std::vector<float> sihit_startPosY(numHits, 0);
69  std::vector<float> sihit_startPosZ(numHits, 0);
70 
71  std::vector<float> sihit_endPosX(numHits, 0);
72  std::vector<float> sihit_endPosY(numHits, 0);
73  std::vector<float> sihit_endPosZ(numHits, 0);
74 
75  int hitNumber {0};
76  for ( const SiHit& sihit : matchingHits ) {
77  sihit_energyDeposit[hitNumber] = sihit.energyLoss() ;
78  sihit_meanTime[hitNumber] = sihit.meanTime() ;
79 
80  const HepMcParticleLink& HMPL = sihit.particleLink();
81  sihit_barcode[hitNumber] = HepMC::barcode(HMPL) ;
82  if( HMPL.isValid() ){
83  sihit_pdgid[hitNumber] = HMPL->pdg_id();
84  }
85 
86  // Convert Simulation frame into reco frame
87  const HepGeom::Point3D<double>& startPos=sihit.localStartPosition();
88 
89  Amg::Vector2D pos = element.hitLocalToLocal( startPos.z(), startPos.y() );
90  sihit_startPosX[hitNumber] = pos[0];
91  sihit_startPosY[hitNumber] = pos[1];
92  sihit_startPosZ[hitNumber] = startPos.x();
93 
94  const HepGeom::Point3D<double>& endPos=sihit.localEndPosition();
95  pos = element.hitLocalToLocal( endPos.z(), endPos.y() );
96  sihit_endPosX[hitNumber] = pos[0];
97  sihit_endPosY[hitNumber] = pos[1];
98  sihit_endPosZ[hitNumber] = endPos.x();
99  ++hitNumber;
100  }
101 
102  return std::make_tuple(std::move(sihit_energyDeposit), std::move(sihit_meanTime), std::move(sihit_barcode), std::move(sihit_pdgid),
103  std::move(sihit_startPosX), std::move(sihit_startPosY), std::move(sihit_startPosZ),
104  std::move(sihit_endPosX), std::move(sihit_endPosY), std::move(sihit_endPosZ));
105  }

◆ restoreDecorations()

std::vector< Decoration > ActsTrk::detail::restoreDecorations ( const SG::IConstAuxStore *  container,
const std::set< std::string > &  staticVariables 
)
  • container - source container to look for decorations
  • staticVariables - set of names of predefined variables for this container

Definition at line 21 of file Decoration.cxx.

23  {
24  std::vector<Decoration> decorations;
25  for (auto id : container->getAuxIDs()) {
26  const std::string name = SG::AuxTypeRegistry::instance().getName(id);
27  const std::type_info* typeInfo =
28  SG::AuxTypeRegistry::instance().getType(id);
29  if (staticVariables.count(name) == 1) {
30  continue;
31  }
32 
33  // try making decoration accessor of matching type
34  // there is a fixed set of supported types (as there is a fixed set
35  // available in MutableMTJ) setters are not needed so replaced by a
36  if ( build<float>(typeInfo, name, decorations) ||
37  build<double>(typeInfo, name, decorations) ||
38  build<short>(typeInfo, name, decorations) ||
39  build<int>(typeInfo, name, decorations) ||
40  build<uint8_t>(typeInfo, name, decorations) ||
41  build<uint16_t>(typeInfo, name, decorations) ||
42  build<uint32_t>(typeInfo, name, decorations) ||
43  build<uint64_t>(typeInfo, name, decorations) ||
44  build<int8_t>(typeInfo, name, decorations) ||
45  build<int16_t>(typeInfo, name, decorations) ||
46  build<int32_t>(typeInfo, name, decorations) ||
47  build<int64_t>(typeInfo, name, decorations)
48  ) {
49  continue;
50  }
51  throw std::runtime_error("Can't restore decoration of " + name +
52  " because it is of an unsupported type");
53  }
54  return decorations;
55 }

◆ to_underlying()

template<typename T_EnumClass >
constexpr std::underlying_type<T_EnumClass>::type ActsTrk::detail::to_underlying ( T_EnumClass  an_enum)
constexpr

Helper to convert class enum into an integer.

Definition at line 24 of file HitSummaryDataUtils.h.

24  {
25  return static_cast<typename std::underlying_type<T_EnumClass>::type>(an_enum);
26  }
data
char data[hepevt_bytes_allocation_ATLAS]
Definition: HepEvt.cxx:11
SG::IConstAuxStore::getData
virtual const void * getData(SG::auxid_t auxid) const =0
Return the data vector for one aux data item.
get_generator_info.result
result
Definition: get_generator_info.py:21
SG::AuxTypeRegistry::instance
static AuxTypeRegistry & instance()
Return the singleton registry instance.
Definition: AuxTypeRegistry.cxx:639
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
SG::AuxTypeRegistry::getName
std::string getName(SG::auxid_t auxid) const
Return the name of an aux data item.
Definition: AuxTypeRegistry.cxx:881
ActsTrk::detail::computeFreeToPathDerivatives
Acts::FreeToPathMatrix computeFreeToPathDerivatives(const Acts::Vector3 &direction, double qop, const Acts::Vector3 &bfield, const Acts::ParticleHypothesis &particle_hypothesis)
Compute the path length derivatives for the free/bound to curvilinear paramter transform.
Definition: CurvilinearCovarianceHelper.h:27
HepMC::INVALID_PARTICLE_ID
constexpr int INVALID_PARTICLE_ID
Definition: MagicNumbers.h:57
ActsDetectorElement::upstreamDetectorElement
const GeoVDetectorElement * upstreamDetectorElement() const
Returns the underllying GeoModel detectorelement that this one is based on.
Definition: ActsDetectorElement.cxx:284
SG::IConstAuxStore::size
virtual size_t size() const =0
Return the number of elements in the store.
Trk::u
@ u
Enums for curvilinear frames.
Definition: ParamDefs.h:77
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
InDetDD::SolidStateDetectorElementBase::getIdHelper
const AtlasDetectorID * getIdHelper() const
Returns the id helper (inline)
SiHit
Definition: SiHit.h:19
python.setupRTTAlg.size
int size
Definition: setupRTTAlg.py:39
SG::AuxTypeRegistry::getType
const std::type_info * getType(SG::auxid_t auxid) const
Return the type of an aux data item.
Definition: AuxTypeRegistry.cxx:908
InDetDD::DetectorType
DetectorType
Definition: DetectorDesign.h:45
AtlasMeasurementSelectorCuts
Definition: MeasurementSelector.h:289
HepMC::barcode
int barcode(const T *p)
Definition: Barcode.h:16
master.flag
bool flag
Definition: master.py:29
InDetDD::DetectorDesign::type
virtual DetectorType type() const
Type of element.
Definition: DetectorDesign.cxx:101
ActsDetectorElement
Definition: ActsDetectorElement.h:42
ActsTrk::MeasurementCalibrator2
Definition: MeasurementCalibrator2.h:101
InDetDD::SolidStateDetectorElementBase::hitLocalToLocal
Amg::Vector2D hitLocalToLocal(double xEta, double xPhi) const
Simulation/Hit local frame to reconstruction local frame.
Definition: SolidStateDetectorElementBase.cxx:95
InDetDD::SiDetectorElement::isPixel
bool isPixel() const
tolerance
Definition: suep_shower.h:17
PixelID::layer_disk
int layer_disk(const Identifier &id) const
Definition: PixelID.h:626
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
InDetDD::PixelBarrel
@ PixelBarrel
Definition: DetectorDesign.h:46
SCT_ID::layer_disk
int layer_disk(const Identifier &id) const
Definition: SCT_ID.h:734
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
InDetDD::SiDetectorElement::isBarrel
bool isBarrel() const
python.Dumpers.typename
def typename(t)
Definition: Dumpers.py:194
ActsTrk::detail::computeDtDs
double computeDtDs(const Acts::ParticleHypothesis &hypothesis, double qop)
Helper function to compute dt/ds Helper function to compute the derivative of the time as function of...
Definition: CurvilinearCovarianceHelper.h:16
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:18
python.selector.AtlRunQuerySelectorLhcOlc.selector
selector
Definition: AtlRunQuerySelectorLhcOlc.py:611
SCT_ID
Definition: SCT_ID.h:68
SG::IAuxStore::getData
virtual void * getData(auxid_t auxid, size_t size, size_t capacity)=0
Return the data vector for one aux data item.
if
if(febId1==febId2)
Definition: LArRodBlockPhysicsV0.cxx:567
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
LArCellBinning.step
step
Definition: LArCellBinning.py:158
IActsExtrapolationTool::propagationSteps
virtual ActsPropagationOutput propagationSteps(const EventContext &ctx, const Acts::BoundTrackParameters &startParameters, Acts::Direction navDir=Acts::Direction::Forward(), double pathLimit=std::numeric_limits< double >::max()) const =0
InDetDD::SiDetectorElement::isSCT
bool isSCT() const
xAOD::track
@ track
Definition: TrackingPrimitives.h:513
InDetDD::SiDetectorElement::design
virtual const SiDetectorDesign & design() const override final
access to the local description (inline):
PixelID
Definition: PixelID.h:67
InDetDD::SolidStateDetectorElementBase::identify
virtual Identifier identify() const override final
identifier of this detector element (inline)
SG::IConstAuxStore::getAuxIDs
virtual const SG::auxid_set_t & getAuxIDs() const =0
Return a set of identifiers for existing data items in this store.
TSU::T
unsigned long long T
Definition: L1TopoDataTypes.h:35
ActsPropagationOutput
std::pair< std::vector< Acts::detail::Step >, ActsRecordedMaterial > ActsPropagationOutput
Finally the output of the propagation test.
Definition: IActsExtrapolationTool.h:25
xAOD::UncalibMeasType::PixelClusterType
@ PixelClusterType
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Mon Apr 21 2025 21:52:12 for ATLAS Offline Software by   doxygen 1.8.18