ATLAS Offline Software
Loading...
Searching...
No Matches
ActsTrk::TrackFindingBaseAlg Class Reference

#include <TrackFindingBaseAlg.h>

Inheritance diagram for ActsTrk::TrackFindingBaseAlg:

Classes

struct  MeasurementSelectorConfig
struct  DetectorContextHolder
struct  TrackFindingDefaultOptions
struct  BranchState
struct  CKF_pimpl

Public Member Functions

 TrackFindingBaseAlg (const std::string &name, ISvcLocator *pSvcLocator)
 ~TrackFindingBaseAlg ()
virtual StatusCode initialize () override
virtual StatusCode finalize () override
virtual StatusCode execute (const EventContext &ctx) const override
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
virtual bool filterPassed (const EventContext &ctx) const
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Types

enum  EStat : std::size_t {
  kNTotalSeeds , kNoTrackParam , kNUsedSeeds , kNoTrack ,
  kNDuplicateSeeds , kNNoEstimatedParams , kNOutputTracks , kNSeedRefitFailure ,
  kNRejectedRefinedSeeds , kNSelectedTracks , kNStoppedTracksMaxHoles , kMultipleBranches ,
  kNoSecond , kNStoppedTracksMinPt , kNStoppedTracksMaxEta , kNTotalSharedHits ,
  kNForcedSeedMeasurements , kNResolvedTracks , kNStat
}
using TrackFinderOptions = Acts::CombinatorialKalmanFilterOptions<detail::RecoTrackContainer>
using EventStats = std::vector<std::array<unsigned int, kNStat>>
using TrkProxy = detail::RecoTrackContainer::TrackProxy
using BranchStopperResult = Acts::CombinatorialKalmanFilterBranchStopperResult

Protected Member Functions

CKF_pimpltrackFinder ()
const CKF_pimpltrackFinder () const
StatusCode initializeMeasurementSelector ()
std::unique_ptr< ActsTrk::IMeasurementSelectorsetMeasurementSelector (const EventContext &ctx, const detail::TrackFindingMeasurements &measurements, TrackFinderOptions &options) const
 Setup and attach measurement selector to KF options.
TrackFindingDefaultOptions getDefaultOptions (const EventContext &ctx, const DetectorContextHolder &detContext, const detail::TrackFindingMeasurements &measurements, const Acts::PerigeeSurface *pSurface) const
 Get CKF options for first and second pass + pointer to MeasurementSelector.
template<class HandleArrayKeyType, class ContainerType>
StatusCode getContainersFromKeys (const EventContext &ctx, HandleArrayKeyType &handleKeyArray, std::vector< const ContainerType * > &outputContainers, std::size_t &sum) const
 Take the array of handle keys and for each key retrieve containers, then append them to the output vector.
const Acts::TrackSelector::Config & getCuts (double eta) const
 Retrieves track selector configuration for given eta value.
template<class MeasurementSource>
std::variant< std::unique_ptr< Acts::BoundTrackParameters >, TrackFindingBaseAlg::EStatdoRefit (const MeasurementSource &measurement, const Acts::BoundTrackParameters &initialParameters, const DetectorContextHolder &detContext, const bool paramsAtOutermostSurface) const
 Perform Kalman Filter fit and update given initialParameters.
std::vector< typename detail::RecoTrackContainer::TrackProxy > doTwoWayTrackFinding (const detail::RecoTrackStateContainerProxy &firstMeasurement, const TrkProxy &trackProxy, detail::RecoTrackContainer &tracksContainerTemp, const TrackFinderOptions &options) const
 Perform two-way track finding.
BranchStopperResult stopBranch (const detail::RecoTrackContainer::TrackProxy &track, const detail::RecoTrackContainer::TrackStateProxy &trackState, const Acts::TrackSelector::EtaBinnedConfig &trackSelectorCfg, const Acts::GeometryContext &tgContext, const detail::MeasurementIndex &measurementIndex, const std::size_t typeIndex, EventStats::value_type &event_stat_category_i) const
 Branch stopper.
void checkCounts (const detail::RecoTrackContainer::TrackProxy &track) const
std::array< bool, 3 > selectCounts (const detail::RecoTrackContainer::TrackProxy &track, double eta) const
bool selectCountsFinal (const detail::RecoTrackContainer::TrackProxy &track) const
const Acts::Logger & logger () const
 Private access to the logger.
void initStatTables ()
void copyStats (const EventStats &event_stat) const
void printStatTables () const
std::size_t nSeedCollections () const
std::size_t seedCollectionStride () const
std::size_t getStatCategory (std::size_t seed_collection, float eta) const
std::size_t computeStatSum (std::size_t seed_collection, EStat counter_i, const EventStats &stat) const
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.

Static Protected Member Functions

static xAOD::UncalibMeasType measurementType (const detail::RecoTrackContainer::TrackStateProxy &trackState)
static void addCountsAndProperties (detail::RecoTrackContainer &tracksContainer, bool add_counts)
static void initCounts (const detail::RecoTrackContainer::TrackProxy &track)
static void updateCounts (const detail::RecoTrackContainer::TrackProxy &track, Acts::ConstTrackStateTypeMap typeFlags, xAOD::UncalibMeasType detType)

Protected Attributes

struct ActsTrk::TrackFindingBaseAlg::MeasurementSelectorConfig m_measurementSelectorConfig
std::unique_ptr< CKF_pimplm_trackFinder
detail::xAODUncalibMeasSurfAcc m_unalibMeasSurfAcc {}
ToolHandle< GenericMonitoringToolm_monTool {this, "MonTool", "", "Monitoring tool"}
ToolHandle< ActsTrk::IExtrapolationToolm_extrapolationTool {this, "ExtrapolationTool", ""}
PublicToolHandle< ActsTrk::ITrackingGeometryToolm_trackingGeometryTool {this, "TrackingGeometryTool", ""}
ToolHandle< ActsTrk::TrackStatePrinterToolm_trackStatePrinter {this, "TrackStatePrinter", "", "optional track state printer"}
ToolHandle< ActsTrk::IFitterToolm_fitterTool {this, "FitterTool", "", "Fitter Tool for Seeds"}
ToolHandle< ActsTrk::IPixelOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_pixelCalibTool {this, "PixelCalibrator", "", "Opt. pixel measurement calibrator"}
ToolHandle< ActsTrk::IStripOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_stripCalibTool {this, "StripCalibrator", "", "Opt. strip measurement calibrator"}
ToolHandle< ActsTrk::IHGTDOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_hgtdCalibTool {this, "HGTDCalibrator", "", "Opt. HGTD measurement calibrator"}
SG::WriteHandleKey< ActsTrk::TrackContainerm_trackContainerKey {this, "ACTSTracksLocation", "", "Output track collection (ActsTrk variant)"}
ActsTrk::MutableTrackContainerHandlesHelper m_tracksBackendHandlesHelper {this}
Gaudi::Property< unsigned int > m_maxPropagationStep {this, "maxPropagationStep", 1000, "Maximum number of steps for one propagate call"}
Gaudi::Property< std::vector< double > > m_etaBins {this, "etaBins", {}, "bins in |eta| to specify variable selections"}
Gaudi::Property< std::vector< double > > m_chi2CutOff {this, "chi2CutOff", {}, "MeasurementSelector: maximum local chi2 contribution"}
Gaudi::Property< std::vector< double > > m_chi2OutlierCutOff {this, "chi2OutlierCutOff", {}, "MeasurementSelector: maximum local chi2 contribution for outlier"}
Gaudi::Property< std::vector< size_t > > m_numMeasurementsCutOff {this, "numMeasurementsCutOff", {}, "MeasurementSelector: maximum number of associated measurements on a single surface"}
Gaudi::Property< double > m_edgeHoleBorderWidth {this,"EdgeHoleBorderWidth", 0. ,"Width of the border of surfaces (mm) within which holes are not counted as holes."}
Gaudi::Property< std::vector< std::size_t > > m_ptMinMeasurements {this, "ptMinMeasurements", {}, "if specified for the given seed collection, applies ptMin cut in branch stopper once ptMinMinMeasurements have been encountered"}
Gaudi::Property< std::vector< std::size_t > > m_absEtaMaxMeasurements {this, "absEtaMaxMeasurements", {}, "if specified for the given seed collection, applies absEtaMax cut in branch stopper once absEtaMaxMeasurements have been encountered"}
Gaudi::Property< bool > m_doBranchStopper {this, "doBranchStopper", true, "use branch stopper"}
Gaudi::Property< bool > m_doTwoWay {this, "doTwoWay", true, "run CKF twice, first with forward propagation with smoothing, then with backward propagation"}
Gaudi::Property< bool > m_inflateCovarianceTwoWay {this, "inflateCovarianceTwoWay", false, "inflate covariance matrix at the beginning of two-way track finding"}
Gaudi::Property< double > m_twoWayinflateCovarianceFactor {this, "twoWayinflateCovarianceFactor", 1.0, "factor to multiply the initial covariance matrix at the beginning of two-way track finding"}
Gaudi::Property< double > m_branchStopperPtMinFactor {this, "branchStopperPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the branch stopper"}
Gaudi::Property< double > m_seedRefitPtMinFactor {this, "seedRefitPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the seed refit"}
Gaudi::Property< double > m_branchStopperAbsEtaMaxExtra {this, "branchStopperAbsEtaMaxExtra", 0.0, "increase absEtaMax cut when used in the branch stopper"}
Gaudi::Property< double > m_branchStopperMeasCutReduce {this, "branchStopperMeasCutReduce", 2, "how much to reduce the minMeas requirement for the branch stopper"}
Gaudi::Property< double > m_branchStopperAbsEtaMeasCut {this, "branchStopperAbsEtaMeasCut", 1.2, "the minimum |eta| to apply the reduction to the minMeas requirement for the branch stopper"}
Gaudi::Property< std::vector< double > > m_phiMin {this, "phiMin", {}, "TrackSelector: phiMin"}
Gaudi::Property< std::vector< double > > m_phiMax {this, "phiMax", {}, "TrackSelector: phiMax"}
Gaudi::Property< std::vector< double > > m_etaMin {this, "etaMin", {}, "TrackSelector: etaMin"}
Gaudi::Property< std::vector< double > > m_etaMax {this, "etaMax", {}, "TrackSelector: etaMax"}
Gaudi::Property< double > m_absEtaMin {this, "absEtaMin", 0.0, "TrackSelector: absEtaMin"}
Gaudi::Property< double > m_absEtaMax {this, "absEtaMax", std::numeric_limits<double>::max(), "TrackSelector: absEtaMax"}
Gaudi::Property< std::vector< double > > m_ptMin {this, "ptMin", {}, "TrackSelector: ptMin"}
Gaudi::Property< std::vector< double > > m_ptMax {this, "ptMax", {}, "TrackSelector: ptMax"}
Gaudi::Property< std::vector< double > > m_d0Min {this, "d0Min", {}, "TrackSelector: d0Min"}
Gaudi::Property< std::vector< double > > m_d0Max {this, "d0Max", {}, "TrackSelector: d0Max"}
Gaudi::Property< std::vector< double > > m_z0Min {this, "z0Min", {}, "TrackSelector: z0Min"}
Gaudi::Property< std::vector< double > > m_z0Max {this, "z0Max", {}, "TrackSelector: z0Max"}
Gaudi::Property< std::vector< std::size_t > > m_minMeasurements {this, "minMeasurements", {}, "TrackSelector: minMeasurements"}
Gaudi::Property< std::vector< std::size_t > > m_maxHoles {this, "maxHoles", {}, "TrackSelector: maxHoles"}
Gaudi::Property< std::vector< std::size_t > > m_maxOutliers {this, "maxOutliers", {}, "TrackSelector: maxOutliers"}
Gaudi::Property< std::vector< std::size_t > > m_maxSharedHits {this, "maxSharedHits", {}, "TrackSelector: maxSharedHits"}
Gaudi::Property< std::vector< double > > m_maxChi2 {this, "maxChi2", {}, "TrackSelector: maxChi2"}
Gaudi::Property< bool > m_addCounts {this, "addCounts", true, "keep separate pixel, strip and hgtd counts and apply the following cuts"}
Gaudi::Property< bool > m_checkCounts {this, "checkCounts", false, "check consistency among track state counts"}
Gaudi::Property< std::vector< std::size_t > > m_minPixelHits {this, "minPixelHits", {}, "minimum number of pixel hits"}
Gaudi::Property< std::vector< std::size_t > > m_minStripHits {this, "minStripHits", {}, "minimum number of strip hits"}
Gaudi::Property< std::vector< std::size_t > > m_minHgtdHits {this, "minHgtdHits", {}, "minimum number of hgtd hits"}
Gaudi::Property< std::vector< std::size_t > > m_maxPixelHoles {this, "maxPixelHoles", {}, "maximum number of pixel holes"}
Gaudi::Property< std::vector< std::size_t > > m_maxStripHoles {this, "maxStripHoles", {}, "maximum number of strip holes"}
Gaudi::Property< std::vector< std::size_t > > m_maxHgtdHoles {this, "maxHgtdHoles", {}, "maximum number of hgtd holes"}
Gaudi::Property< std::vector< std::size_t > > m_maxPixelOutliers {this, "maxPixelOutliers", {}, "maximum number of pixel outliers"}
Gaudi::Property< std::vector< std::size_t > > m_maxStripOutliers {this, "maxStripOutliers", {}, "maximum number of strip outliers"}
Gaudi::Property< std::vector< std::size_t > > m_maxHgtdOutliers {this, "maxHgtdOutliers", {}, "maximum number of hgtd outliers"}
Gaudi::Property< std::vector< std::uint32_t > > m_endOfWorldVolumeIds {this, "EndOfTheWorldVolumeIds", {}, ""}
Gaudi::Property< std::vector< float > > m_statEtaBins {this, "StatisticEtaBins", {-4, -2.6, -2, 0, 2., 2.6, 4}, "Gather statistics separately for these bins."}
Gaudi::Property< std::vector< std::string > > m_seedLabels {this, "SeedLabels", {}, "One label per seed key used in outputs"}
Gaudi::Property< bool > m_dumpAllStatEtaBins {this, "DumpEtaBinsForAll", false, "Dump eta bins of all statistics counter."}
std::unique_ptr< const Acts::Logger > m_logger
 logging instance
bool m_useAbsEtaForStat = false
bool m_showResolvedStats = false
std::mutex m_mutex ATLAS_THREAD_SAFE
std::vector< std::array< std::size_t, kNStat > > m_stat ATLAS_THREAD_SAFE {}

Static Protected Attributes

static constexpr BranchState s_branchState {}

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

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

Private Attributes

DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 46 of file TrackFindingBaseAlg.h.

Member Typedef Documentation

◆ BranchStopperResult

using ActsTrk::TrackFindingBaseAlg::BranchStopperResult = Acts::CombinatorialKalmanFilterBranchStopperResult
protected

Definition at line 282 of file TrackFindingBaseAlg.h.

◆ EventStats

using ActsTrk::TrackFindingBaseAlg::EventStats = std::vector<std::array<unsigned int, kNStat>>
protected

Definition at line 184 of file TrackFindingBaseAlg.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ TrackFinderOptions

using ActsTrk::TrackFindingBaseAlg::TrackFinderOptions = Acts::CombinatorialKalmanFilterOptions<detail::RecoTrackContainer>
protected

Definition at line 56 of file TrackFindingBaseAlg.h.

◆ TrkProxy

using ActsTrk::TrackFindingBaseAlg::TrkProxy = detail::RecoTrackContainer::TrackProxy
protected

Definition at line 257 of file TrackFindingBaseAlg.h.

Member Enumeration Documentation

◆ EStat

enum ActsTrk::TrackFindingBaseAlg::EStat : std::size_t
protected
Enumerator
kNTotalSeeds 
kNoTrackParam 
kNUsedSeeds 
kNoTrack 
kNDuplicateSeeds 
kNNoEstimatedParams 
kNOutputTracks 
kNSeedRefitFailure 
kNRejectedRefinedSeeds 
kNSelectedTracks 
kNStoppedTracksMaxHoles 
kMultipleBranches 
kNoSecond 
kNStoppedTracksMinPt 
kNStoppedTracksMaxEta 
kNTotalSharedHits 
kNForcedSeedMeasurements 
kNResolvedTracks 
kNStat 

Definition at line 161 of file TrackFindingBaseAlg.h.

162 {
166 kNoTrack,
175 kNoSecond,
181 kNStat
182 };

Constructor & Destructor Documentation

◆ TrackFindingBaseAlg()

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

Definition at line 24 of file TrackFindingBaseAlg.cxx.

24: AthReentrantAlgorithm(name, pSvcLocator) {}

◆ ~TrackFindingBaseAlg()

ActsTrk::TrackFindingBaseAlg::~TrackFindingBaseAlg ( )
default

Member Function Documentation

◆ addCountsAndProperties()

void ActsTrk::TrackFindingBaseAlg::addCountsAndProperties ( detail::RecoTrackContainer & tracksContainer,
bool add_counts )
staticprotected

Definition at line 411 of file TrackFindingBaseAlg.cxx.

412 {
413 if (addCounts) {
414 tracksContainer.addColumn<unsigned int>("nPixelHits");
415 tracksContainer.addColumn<unsigned int>("nStripHits");
416 tracksContainer.addColumn<unsigned int>("nHgtdHits");
417 tracksContainer.addColumn<unsigned int>("nPixelHoles");
418 tracksContainer.addColumn<unsigned int>("nStripHoles");
419 tracksContainer.addColumn<unsigned int>("nHgtdHoles");
420 tracksContainer.addColumn<unsigned int>("nPixelOutliers");
421 tracksContainer.addColumn<unsigned int>("nStripOutliers");
422 tracksContainer.addColumn<unsigned int>("nHgtdOutliers");
423 }
425 }
static void addFitterTypeProperty(track_container_t &tracksContainer)
add fitter column to the track container

◆ checkCounts()

void ActsTrk::TrackFindingBaseAlg::checkCounts ( const detail::RecoTrackContainer::TrackProxy & track) const
protected

Definition at line 470 of file TrackFindingBaseAlg.cxx.

470 {
471 // This check will fail if there are other types (HGTD, MS?) of hits, holes, or outliers.
472 // The check can be removed when it is no longer appropriate.
473 if (track.nMeasurements() != s_branchState.nPixelHits(track) + s_branchState.nStripHits(track) + s_branchState.nHgtdHits(track))
474 ATH_MSG_WARNING("mismatched hit count: total (" << track.nMeasurements()
475 << ") != pixel (" << s_branchState.nPixelHits(track)
476 << ") + strip (" << s_branchState.nStripHits(track)
477 << ") + hgtd (" << s_branchState.nHgtdHits(track)
478 << ")");
479 if (track.nHoles() != s_branchState.nPixelHoles(track) + s_branchState.nStripHoles(track) + s_branchState.nHgtdHoles(track))
480 ATH_MSG_WARNING("mismatched hole count: total (" << track.nHoles()
481 << ") < pixel (" << s_branchState.nPixelHoles(track)
482 << ") + strip (" << s_branchState.nStripHoles(track)
483 << ") + hgtd (" << s_branchState.nHgtdHoles(track)
484 << ")");
485 if (track.nOutliers() != s_branchState.nPixelOutliers(track) + s_branchState.nStripOutliers(track) + s_branchState.nHgtdOutliers(track))
486 ATH_MSG_WARNING("mismatched outlier count: total (" << track.nOutliers()
487 << ") != pixel (" << s_branchState.nPixelOutliers(track)
488 << ") + strip (" << s_branchState.nStripOutliers(track)
489 << ") + hgtd (" << s_branchState.nHgtdOutliers(track)
490 << ")");
491 };
#define ATH_MSG_WARNING(x)
static constexpr BranchState s_branchState

◆ computeStatSum()

std::size_t ActsTrk::TrackFindingBaseAlg::computeStatSum ( std::size_t seed_collection,
EStat counter_i,
const EventStats & stat ) const
protected

Definition at line 780 of file TrackFindingBaseAlg.cxx.

780 {
781 std::size_t out = 0u;
782 for (std::size_t category_i = seed_collection * seedCollectionStride();
783 category_i < (seed_collection + 1) * seedCollectionStride();
784 ++category_i)
785 {
786 assert(category_i < stat.size());
787 out += stat[category_i][counter_i];
788 }
789 return out;
790 }
std::size_t seedCollectionStride() const
@ u
Enums for curvilinear frames.
Definition ParamDefs.h:77

◆ copyStats()

void ActsTrk::TrackFindingBaseAlg::copyStats ( const EventStats & event_stat) const
protected

Definition at line 566 of file TrackFindingBaseAlg.cxx.

566 {
567 std::lock_guard<std::mutex> lock(m_mutex);
568 std::size_t category_i = 0;
569 for (const std::array<unsigned int, kNStat> &src_stat : event_stat)
570 {
571 std::array<std::size_t, kNStat> &dest_stat = m_stat[category_i++];
572 for (std::size_t i = 0; i < src_stat.size(); ++i)
573 {
574 assert(i < dest_stat.size());
575 dest_stat[i] += src_stat[i];
576 }
577 }
578 }
virtual void lock()=0
Interface to allow an object to lock itself when made const in SG.

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

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

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

◆ doRefit()

template<class MeasurementSource>
std::variant< std::unique_ptr< Acts::BoundTrackParameters >, TrackFindingBaseAlg::EStat > ActsTrk::TrackFindingBaseAlg::doRefit ( const MeasurementSource & measurement,
const Acts::BoundTrackParameters & initialParameters,
const DetectorContextHolder & detContext,
const bool paramsAtOutermostSurface ) const
protected

Perform Kalman Filter fit and update given initialParameters.

Template Parameters
MeasurementSourceType of measurement source: ActsTrk::Seed or (in future) ActsTrk::ProtoTrack
Parameters
measurementMeasurement source for KF
initialParametersParameters to use in KF
detContextStruct holding geometry, magnetic field and calibration contexts
paramsAtOutermostSurfaceFlag for searching in reverse direction
Returns
Unique pointer to updated parameters

◆ doTwoWayTrackFinding()

std::vector< typename detail::RecoTrackContainer::TrackProxy > ActsTrk::TrackFindingBaseAlg::doTwoWayTrackFinding ( const detail::RecoTrackStateContainerProxy & firstMeasurement,
const TrkProxy & trackProxy,
detail::RecoTrackContainer & tracksContainerTemp,
const TrackFinderOptions & options ) const
protected

Perform two-way track finding.

Parameters
addTrackFunction or lambda which adds newly found track to container
trackProxyTrack proxy object
tracksContainerTempTrack proxy container
optionsFit options
tgContextGeometry context
reverseSearchFlag for searching in reverse direction
seedTypeSeed type (only used for warnings/debug printouts)
iseedSeed number (only used for warnings/debug printouts)
itrackTrack number (only used for warnings/debug printouts)
Returns
Number of found tracks

Definition at line 243 of file TrackFindingBaseAlg.cxx.

246 {
247 if (not m_doTwoWay) return {};
248
249 // Create initial parameters for the propagation
250 Acts::BoundTrackParameters secondInitialParameters = trackProxy.createParametersFromState(detail::RecoConstTrackStateContainerProxy{firstMeasurement});
251 if (!secondInitialParameters.referenceSurface().insideBounds(secondInitialParameters.localPosition())) { // #3751
252 return {};
253 }
254
255 // First, inflate the covariance matrix if configured
257 ATH_MSG_DEBUG("Inflating covariance matrix for second track finding with factor = " << m_twoWayinflateCovarianceFactor.value());
258 ATH_MSG_VERBOSE("Original parameters before inflation: \n" << secondInitialParameters);
259
260 auto inflatedCovariance = secondInitialParameters.covariance().value();
261 inflatedCovariance *= m_twoWayinflateCovarianceFactor;
262
263 const auto& origSurface = secondInitialParameters.referenceSurface();
264 auto surfacePtr = const_cast<Acts::Surface&>(origSurface).shared_from_this();
265
266 Acts::BoundTrackParameters newParams(
267 std::static_pointer_cast<const Acts::Surface>(std::move(surfacePtr)),
268 secondInitialParameters.parameters(),
269 std::make_optional(inflatedCovariance),
270 secondInitialParameters.particleHypothesis());
271 secondInitialParameters = std::move(newParams);
272
273 ATH_MSG_VERBOSE("Inflated covariance matrix : \n" << secondInitialParameters.covariance().value());
274 }
275
276 auto rootBranch = tracksContainerTemp.makeTrack();
277 rootBranch.copyFromWithoutStates(trackProxy); // #3534
278
279 // perform track finding
280 auto secondResult =
281 trackFinder().ckf.findTracks(secondInitialParameters, options, tracksContainerTemp, rootBranch);
282 if (not secondResult.ok()) {
283 return {};
284 }
285 return secondResult.value();
286 }
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_DEBUG(x)
Gaudi::Property< bool > m_doTwoWay
Gaudi::Property< double > m_twoWayinflateCovarianceFactor
Gaudi::Property< bool > m_inflateCovarianceTwoWay
RecoTrackStateContainer::ConstTrackStateProxy RecoConstTrackStateContainerProxy
const xAOD::UncalibratedMeasurement * firstMeasurement(const xAOD::MuonSegment &segment, const bool skipOutlier=true)
Retrieves the first measurement associated with the segment.

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

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

Reimplemented in ActsTrk::HGTDTrackExtensionAlg, and ActsTrk::TrackFindingAlg.

Definition at line 173 of file TrackFindingBaseAlg.cxx.

173 {
174 ATH_MSG_FATAL("execute() method from the base class was called! Implement proper execute() method in the derived class!");
175
176 return StatusCode::FAILURE;
177 }
#define ATH_MSG_FATAL(x)

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

Return the list of extra output dependencies.

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

Definition at line 89 of file AthCommonAlgorithm.cxx.

54{
55 // If we didn't find any symlinks to add, just return the collection
56 // from the base class. Otherwise, return the extended collection.
57 if (!m_extendedExtraObjects.empty()) {
59 }
61}
Common base class for algorithms.

◆ filterPassed()

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

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

93 {
94 return execState( ctx ).filterPassed();
95 }
virtual bool filterPassed(const EventContext &ctx) const
Get filter decision:

◆ finalize()

StatusCode ActsTrk::TrackFindingBaseAlg::finalize ( )
overridevirtual

Reimplemented in ActsTrk::TrackFindingAlg.

Definition at line 179 of file TrackFindingBaseAlg.cxx.

179 {
181
182 return StatusCode::SUCCESS;
183 }

◆ getContainersFromKeys()

template<class HandleArrayKeyType, class ContainerType>
StatusCode ActsTrk::TrackFindingBaseAlg::getContainersFromKeys ( const EventContext & ctx,
HandleArrayKeyType & handleKeyArray,
std::vector< const ContainerType * > & outputContainers,
std::size_t & sum ) const
protected

Take the array of handle keys and for each key retrieve containers, then append them to the output vector.

Template Parameters
HandleArrayKeyTypeType of the list of handle keys
ContainerTypeType of the output container
Parameters
ctxEvent context
handleKeyArrayList of handle keys
outputContainersVector of output containers
sumNumber of total elements in all retrieved containers
Returns
Status code

◆ getCuts()

const Acts::TrackSelector::Config & ActsTrk::TrackFindingBaseAlg::getCuts ( double eta) const
protected

Retrieves track selector configuration for given eta value.

Parameters
etatrack candidate eta value

Definition at line 234 of file TrackFindingBaseAlg.cxx.

234 {
235 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
236 // return the last bin for |eta|>=4 or nan
237 return (!(std::abs(eta) < trackSelectorCfg.absEtaEdges.back())) ? trackSelectorCfg.cutSets.back()
238 : (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? trackSelectorCfg.cutSets.front()
239 : trackSelectorCfg.getCuts(eta);
240 };
Scalar eta() const
pseudorapidity method

◆ getDefaultOptions()

TrackFindingBaseAlg::TrackFindingDefaultOptions ActsTrk::TrackFindingBaseAlg::getDefaultOptions ( const EventContext & ctx,
const DetectorContextHolder & detContext,
const detail::TrackFindingMeasurements & measurements,
const Acts::PerigeeSurface * pSurface ) const
protected

Get CKF options for first and second pass + pointer to MeasurementSelector.

Parameters
detContextObject holding detector-related context
measurements<easurements container used in MeasurementSelector
pSurfaceRaw pointer to perigee surface

Definition at line 206 of file TrackFindingBaseAlg.cxx.

210 {
211 Acts::PropagatorPlainOptions plainOptions{detContext.geometry, detContext.magField};
212 plainOptions.maxSteps = m_maxPropagationStep;
213 plainOptions.direction = Acts::Direction::Forward();
214 plainOptions.endOfWorldVolumeIds = m_endOfWorldVolumeIds;
215
216 // Set the CombinatorialKalmanFilter options
217 TrackFinderOptions options(detContext.geometry, detContext.magField, detContext.calib,
218 trackFinder().ckfExtensions, plainOptions, pSurface);
219
220 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = setMeasurementSelector(ctx, measurements, options);
221
222 Acts::PropagatorPlainOptions plainSecondOptions{detContext.geometry, detContext.magField};
223 plainSecondOptions.maxSteps = m_maxPropagationStep;
224 plainSecondOptions.direction = plainOptions.direction.invert();
225
226 TrackFinderOptions secondOptions(detContext.geometry, detContext.magField, detContext.calib,
227 options.extensions, plainSecondOptions, pSurface);
228 secondOptions.targetSurface = pSurface;
229 secondOptions.skipPrePropagationUpdate = true;
230
231 return {std::move(options), std::move(secondOptions), std::move(measurementSelector)};
232 };
Gaudi::Property< unsigned int > m_maxPropagationStep
Gaudi::Property< std::vector< std::uint32_t > > m_endOfWorldVolumeIds
Acts::CombinatorialKalmanFilterOptions< detail::RecoTrackContainer > TrackFinderOptions
std::unique_ptr< ActsTrk::IMeasurementSelector > setMeasurementSelector(const EventContext &ctx, const detail::TrackFindingMeasurements &measurements, TrackFinderOptions &options) const
Setup and attach measurement selector to KF options.

◆ getStatCategory()

std::size_t ActsTrk::TrackFindingBaseAlg::getStatCategory ( std::size_t seed_collection,
float eta ) const
protected

Definition at line 771 of file TrackFindingBaseAlg.cxx.

771 {
772 std::vector<float>::const_iterator bin_iter = std::upper_bound(m_statEtaBins.begin(),
773 m_statEtaBins.end(),
774 m_useAbsEtaForStat ? std::abs(eta) : eta);
775 std::size_t category_i = seed_collection * seedCollectionStride() + static_cast<std::size_t>(bin_iter - m_statEtaBins.begin());
776 assert(category_i < m_stat.size());
777 return category_i;
778 }
Gaudi::Property< std::vector< float > > m_statEtaBins

◆ initCounts()

void ActsTrk::TrackFindingBaseAlg::initCounts ( const detail::RecoTrackContainer::TrackProxy & track)
staticprotected

Definition at line 427 of file TrackFindingBaseAlg.cxx.

428 {
429 s_branchState.nPixelHits(track) = 0;
430 s_branchState.nStripHits(track) = 0;
431 s_branchState.nHgtdHits(track) = 0;
432 s_branchState.nPixelHoles(track) = 0;
433 s_branchState.nStripHoles(track) = 0;
434 s_branchState.nHgtdHoles(track) = 0;
435 s_branchState.nPixelOutliers(track) = 0;
436 s_branchState.nStripOutliers(track) = 0;
437 s_branchState.nHgtdOutliers(track) = 0;
438 }

◆ initialize()

StatusCode ActsTrk::TrackFindingBaseAlg::initialize ( )
overridevirtual

Reimplemented in ActsTrk::HGTDTrackExtensionAlg, and ActsTrk::TrackFindingAlg.

Definition at line 28 of file TrackFindingBaseAlg.cxx.

28 {
29 ATH_MSG_DEBUG("Properties Summary:");
41 ATH_MSG_DEBUG(" " << m_phiMin);
42 ATH_MSG_DEBUG(" " << m_phiMax);
43 ATH_MSG_DEBUG(" " << m_etaMin);
44 ATH_MSG_DEBUG(" " << m_etaMax);
47 ATH_MSG_DEBUG(" " << m_ptMin);
48 ATH_MSG_DEBUG(" " << m_ptMax);
49 ATH_MSG_DEBUG(" " << m_d0Min);
50 ATH_MSG_DEBUG(" " << m_d0Max);
51 ATH_MSG_DEBUG(" " << m_z0Min);
52 ATH_MSG_DEBUG(" " << m_z0Max);
70
71 m_logger = makeActsAthenaLogger(this, "Acts");
72
73 // Read and Write handles
74 ATH_CHECK(m_trackContainerKey.initialize());
76
77 ATH_CHECK(m_monTool.retrieve(EnableTool{not m_monTool.empty()}));
80 ATH_CHECK(m_trackStatePrinter.retrieve(EnableTool{not m_trackStatePrinter.empty()}));
81 ATH_CHECK(m_fitterTool.retrieve());
82 ATH_CHECK(m_pixelCalibTool.retrieve(EnableTool{not m_pixelCalibTool.empty()}));
83 ATH_CHECK(m_stripCalibTool.retrieve(EnableTool{not m_stripCalibTool.empty()}));
84 ATH_CHECK(m_hgtdCalibTool.retrieve(EnableTool{not m_hgtdCalibTool.empty()}));
85
86 auto magneticField = std::make_unique<ATLASMagneticFieldWrapper>();
87 auto trackingGeometry = m_trackingGeometryTool->trackingGeometry();
88
89 detail::Stepper stepper(std::move(magneticField));
90 detail::Navigator::Config config{trackingGeometry};
91 config.resolvePassive = false;
92 config.resolveMaterial = true;
93 config.resolveSensitive = true;
94 detail::Navigator navigator(config, logger().cloneWithSuffix("Navigator"));
95 detail::Propagator propagator(std::move(stepper), std::move(navigator), logger().cloneWithSuffix("Prop"));
96
97 // Using the CKF propagator as extrapolator
98 detail::Extrapolator extrapolator = propagator;
99
100 // m_etaBins (from flags.Tracking.ActiveConfig.etaBins) includes a dummy first and last bin, which we ignore
101 std::vector<double> absEtaEdges;
102 if (m_etaBins.size() <= 2)
103 {
104 absEtaEdges.reserve(2ul);
105 absEtaEdges.push_back(0.0);
106 absEtaEdges.push_back(std::numeric_limits<double>::infinity());
107 }
108 else
109 {
110 absEtaEdges.reserve(m_etaBins.size());
111 absEtaEdges.push_back(m_absEtaMin);
112 absEtaEdges.insert(absEtaEdges.end(), m_etaBins.begin() + 1, m_etaBins.end() - 1);
113 absEtaEdges.push_back(m_absEtaMax);
114 }
115
116 auto setCut = [](auto &cfgVal, const auto &cuts, size_t ind) -> void
117 {
118 if (cuts.empty())
119 return;
120 cfgVal = (ind < cuts.size()) ? cuts[ind] : cuts[cuts.size() - 1];
121 };
122
123 Acts::TrackSelector::EtaBinnedConfig trackSelectorCfg{std::move(absEtaEdges)};
124 if (m_etaBins.size() <= 2)
125 {
126 assert(trackSelectorCfg.cutSets.size() == 1);
127 trackSelectorCfg.cutSets[0].absEtaMin = m_absEtaMin;
128 trackSelectorCfg.cutSets[0].absEtaMax = m_absEtaMax;
129 }
130 size_t cutIndex = 0;
131 for (auto &cfg : trackSelectorCfg.cutSets)
132 {
133 setCut(cfg.phiMin, m_phiMin, cutIndex);
134 setCut(cfg.phiMax, m_phiMax, cutIndex);
135 setCut(cfg.etaMin, m_etaMin, cutIndex);
136 setCut(cfg.etaMax, m_etaMax, cutIndex);
137 setCut(cfg.ptMin, m_ptMin, cutIndex);
138 setCut(cfg.ptMax, m_ptMax, cutIndex);
139 setCut(cfg.loc0Min, m_d0Min, cutIndex);
140 setCut(cfg.loc0Max, m_d0Max, cutIndex);
141 setCut(cfg.loc1Min, m_z0Min, cutIndex);
142 setCut(cfg.loc1Max, m_z0Max, cutIndex);
143 setCut(cfg.minMeasurements, m_minMeasurements, cutIndex);
144 setCut(cfg.maxHoles, m_maxHoles, cutIndex);
145 setCut(cfg.maxOutliers, m_maxOutliers, cutIndex);
146 setCut(cfg.maxSharedHits, m_maxSharedHits, cutIndex);
147 setCut(cfg.maxChi2, m_maxChi2, cutIndex);
148 ++cutIndex;
149 }
150
151 ATH_MSG_DEBUG(trackSelectorCfg);
152
153 // initializer measurement selector and connect it to the delegates of the track finder optins
155
156 detail::CKF_config ckfConfig{
157 std::move(extrapolator),
158 detail::CKF{std::move(propagator), logger().cloneWithSuffix("CKF")},
159 {},
160 Acts::TrackSelector{trackSelectorCfg}};
161
162 m_trackFinder = std::make_unique<CKF_pimpl>(std::move(ckfConfig));
163
165
166 m_unalibMeasSurfAcc = detail::xAODUncalibMeasSurfAcc {m_trackingGeometryTool.get()};
167
169
170 return StatusCode::SUCCESS;
171 }
#define ATH_CHECK
Evaluate an expression and check for errors.
std::unique_ptr< const Acts::Logger > makeActsAthenaLogger(IMessageSvc *svc, const std::string &name, int level, std::optional< std::string > parent_name)
Gaudi::Property< std::vector< std::size_t > > m_maxStripHoles
Gaudi::Property< double > m_branchStopperAbsEtaMeasCut
Gaudi::Property< double > m_branchStopperPtMinFactor
Gaudi::Property< std::vector< std::size_t > > m_maxOutliers
std::unique_ptr< const Acts::Logger > m_logger
logging instance
Gaudi::Property< std::vector< double > > m_etaBins
Gaudi::Property< std::vector< size_t > > m_numMeasurementsCutOff
Gaudi::Property< std::vector< std::size_t > > m_maxSharedHits
ToolHandle< ActsTrk::TrackStatePrinterTool > m_trackStatePrinter
Gaudi::Property< std::vector< double > > m_etaMax
Gaudi::Property< std::vector< double > > m_ptMin
Gaudi::Property< std::vector< std::size_t > > m_minPixelHits
Gaudi::Property< std::vector< std::size_t > > m_maxStripOutliers
SG::WriteHandleKey< ActsTrk::TrackContainer > m_trackContainerKey
Gaudi::Property< std::vector< double > > m_ptMax
Gaudi::Property< std::vector< double > > m_chi2CutOff
Gaudi::Property< double > m_seedRefitPtMinFactor
Gaudi::Property< std::vector< std::size_t > > m_maxPixelHoles
Gaudi::Property< std::vector< double > > m_phiMin
Gaudi::Property< double > m_branchStopperMeasCutReduce
ToolHandle< ActsTrk::IPixelOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_pixelCalibTool
Gaudi::Property< std::vector< std::size_t > > m_minStripHits
ActsTrk::MutableTrackContainerHandlesHelper m_tracksBackendHandlesHelper
Gaudi::Property< double > m_branchStopperAbsEtaMaxExtra
ToolHandle< GenericMonitoringTool > m_monTool
Gaudi::Property< std::vector< std::size_t > > m_absEtaMaxMeasurements
Gaudi::Property< std::vector< double > > m_etaMin
Gaudi::Property< std::vector< double > > m_d0Min
PublicToolHandle< ActsTrk::ITrackingGeometryTool > m_trackingGeometryTool
Gaudi::Property< std::vector< double > > m_z0Max
Gaudi::Property< double > m_absEtaMax
ToolHandle< ActsTrk::IExtrapolationTool > m_extrapolationTool
Gaudi::Property< std::vector< double > > m_z0Min
Gaudi::Property< bool > m_doBranchStopper
Gaudi::Property< std::vector< std::size_t > > m_minMeasurements
Gaudi::Property< std::vector< double > > m_maxChi2
const Acts::Logger & logger() const
Private access to the logger.
Gaudi::Property< double > m_absEtaMin
Gaudi::Property< std::vector< std::size_t > > m_maxPixelOutliers
Gaudi::Property< bool > m_checkCounts
Gaudi::Property< std::vector< std::size_t > > m_maxHoles
ToolHandle< ActsTrk::IStripOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_stripCalibTool
std::unique_ptr< CKF_pimpl > m_trackFinder
ToolHandle< ActsTrk::IHGTDOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_hgtdCalibTool
detail::xAODUncalibMeasSurfAcc m_unalibMeasSurfAcc
Gaudi::Property< bool > m_addCounts
Gaudi::Property< std::vector< double > > m_phiMax
Gaudi::Property< std::vector< std::size_t > > m_ptMinMeasurements
Gaudi::Property< std::vector< double > > m_d0Max
ToolHandle< ActsTrk::IFitterTool > m_fitterTool
Gaudi::Property< std::vector< double > > m_chi2OutlierCutOff
Acts::Result< void > gainMatrixUpdate(const Acts::GeometryContext &gctx, typename trajectory_t::TrackStateProxy trackState, const Acts::Logger &logger)
Acts::SympyStepper Stepper
Adapted from Acts Examples/Algorithms/TrackFinding/src/TrackFindingAlgorithmFunction....
Acts::CombinatorialKalmanFilter< Propagator, RecoTrackContainer > CKF
std::string prefixFromTrackContainerName(const std::string &tracks)
Parse TrackContainer name to get the prefix for backends The name has to contain XYZTracks,...
Acts::CombinatorialKalmanFilterExtensions< RecoTrackContainer > ckfExtensions

◆ initializeMeasurementSelector()

StatusCode ActsTrk::TrackFindingBaseAlg::initializeMeasurementSelector ( )
protected

Definition at line 519 of file TrackFindingBaseAlg.cxx.

519 {
520 std::vector<std::pair<float, float> > &chi2CutOffOutlier = m_measurementSelectorConfig.m_chi2CutOffOutlier;
521 chi2CutOffOutlier .reserve( m_chi2CutOff.size() );
522 if (!m_chi2OutlierCutOff.empty()) {
523 if (m_chi2CutOff.size() != m_chi2OutlierCutOff.size()) {
524 ATH_MSG_ERROR("Outlier chi2 cut off provided but number of elements does not agree with"
525 " chi2 cut off for measurements which however is required: "
526 << m_chi2CutOff.size() << " != " << m_chi2OutlierCutOff.size());
527 return StatusCode::FAILURE;
528 }
529 }
530 unsigned int idx=0;
531 for (const auto &elm : m_chi2CutOff) {
532 chi2CutOffOutlier.push_back( std::make_pair(static_cast<float>(elm),
533 idx < m_chi2OutlierCutOff.size()
534 ? static_cast<float>(m_chi2OutlierCutOff[idx])
535 : std::numeric_limits<float>::max()) );
536 ++idx;
537 }
538 if (m_etaBins.size() > 2) {
539 std::vector<float> &etaBinsf = m_measurementSelectorConfig.m_etaBins;
540 etaBinsf.assign(m_etaBins.begin() + 1, m_etaBins.end() - 1);
541 }
542
543 return /*m_measurementSelector ?*/ StatusCode::SUCCESS /*: StatusCode::FAILURE*/;
544 }
#define ATH_MSG_ERROR(x)
struct ActsTrk::TrackFindingBaseAlg::MeasurementSelectorConfig m_measurementSelectorConfig

◆ initStatTables()

void ActsTrk::TrackFindingBaseAlg::initStatTables ( )
protected

Definition at line 548 of file TrackFindingBaseAlg.cxx.

548 {
549 if (!m_statEtaBins.empty())
550 {
552 float last_eta = m_statEtaBins[0];
553 for (float eta : m_statEtaBins)
554 {
555 if (eta < last_eta)
556 {
557 ATH_MSG_FATAL("Eta bins for statistics counter not in ascending order.");
558 }
559 last_eta = eta;
560 }
561 }
562 m_stat.resize(nSeedCollections() * seedCollectionStride());
563 }

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

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const
inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

68 {
69 return true;
70 }

◆ logger()

const Acts::Logger & ActsTrk::TrackFindingBaseAlg::logger ( ) const
inlineprotected

Private access to the logger.

Definition at line 329 of file TrackFindingBaseAlg.h.

330 {
331 return *m_logger;
332 }

◆ measurementType()

xAOD::UncalibMeasType ActsTrk::TrackFindingBaseAlg::measurementType ( const detail::RecoTrackContainer::TrackStateProxy & trackState)
staticprotected

Definition at line 288 of file TrackFindingBaseAlg.cxx.

288 {
289 if (trackState.hasReferenceSurface()) {
290 if (const auto *actsDetElem = dynamic_cast<const ISurfacePlacement*>(trackState.referenceSurface().surfacePlacement())) {
291 switch (actsDetElem->detectorType()) {
298 default:
299 break;
300 }
301 }
302 }
304 }
@ Pixel
Inner detector legacy.

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ nSeedCollections()

std::size_t ActsTrk::TrackFindingBaseAlg::nSeedCollections ( ) const
inlineprotected

Definition at line 342 of file TrackFindingBaseAlg.h.

342 {
343 return m_seedLabels.size();
344 }
Gaudi::Property< std::vector< std::string > > m_seedLabels

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

◆ printStatTables()

void ActsTrk::TrackFindingBaseAlg::printStatTables ( ) const
protected

Definition at line 581 of file TrackFindingBaseAlg.cxx.

581 {
582 if (msgLvl(MSG::INFO))
583 {
584 std::vector<std::string> stat_labels =
586 {
587 std::make_pair(kNTotalSeeds, "Input seeds"),
588 std::make_pair(kNoTrackParam, "No track parameters"),
589 std::make_pair(kNUsedSeeds, "Used seeds"),
590 std::make_pair(kNoTrack, "Cannot find track"),
591 std::make_pair(kNDuplicateSeeds, "Duplicate seeds"),
592 std::make_pair(kNNoEstimatedParams, "Initial param estimation failed"),
593 std::make_pair(kNRejectedRefinedSeeds, "Rejected refined parameters"),
594 std::make_pair(kNSeedRefitFailure, "Seed refit Kalman fit failure"),
595 std::make_pair(kNOutputTracks, "CKF tracks"),
596 std::make_pair(kNSelectedTracks, "selected tracks"),
597 std::make_pair(kNResolvedTracks, "resolved tracks"),
598 std::make_pair(kNStoppedTracksMaxHoles, "Stopped tracks reaching max holes"),
599 std::make_pair(kMultipleBranches, "Seeds with more than one branch"),
600 std::make_pair(kNoSecond, "Tracks failing second CKF"),
601 std::make_pair(kNStoppedTracksMinPt, "Stopped tracks below pT cut"),
602 std::make_pair(kNStoppedTracksMaxEta, "Stopped tracks above max eta"),
603 std::make_pair(kNTotalSharedHits, "Total shared hits"),
604 std::make_pair(kNForcedSeedMeasurements, "Total forced measurements")
605 });
606 assert(stat_labels.size() == kNStat);
607 std::vector<std::string> categories;
608 categories.reserve(m_seedLabels.size() + 1);
609 categories.insert(categories.end(), m_seedLabels.begin(), m_seedLabels.end());
610 categories.push_back("ALL");
611
612 std::vector<std::string> eta_labels;
613 eta_labels.reserve(m_statEtaBins.size() + 2);
614 for (std::size_t eta_bin_i = 0; eta_bin_i < m_statEtaBins.size() + 2; ++eta_bin_i)
615 {
616 eta_labels.push_back(TableUtils::makeEtaBinLabel(m_statEtaBins, eta_bin_i, m_useAbsEtaForStat));
617 }
618
619 // vector used as 3D array stat[ eta_bin ][ stat_i ][ seed_type]
620 // stat_i = [0, kNStat)
621 // eta_bin = [0, m_statEtaBins.size()+2 ); eta_bin == m_statEtaBinsSize()+1 means sum of all etaBins
622 // seed_type = [0, nSeedCollections()+1) seed_type == nSeedCollections() means sum of all seed collections
623 std::vector<std::size_t> stat =
625 m_statEtaBins.size() + 1,
626 m_stat);
627
628 // the extra columns and rows for the projections are addeded internally:
629 std::size_t stat_stride =
631 m_statEtaBins.size() + 1,
632 kNStat);
633 std::size_t eta_stride =
635 m_statEtaBins.size() + 1,
636 kNStat);
637 std::stringstream table_out;
638
639 if (m_dumpAllStatEtaBins.value())
640 {
641 // dump for each counter a table with one row per eta bin
642 std::size_t max_label_width = TableUtils::maxLabelWidth(stat_labels) + TableUtils::maxLabelWidth(eta_labels);
643 for (std::size_t stat_i = 0; stat_i < kNStat; ++stat_i)
644 {
645 std::size_t dest_idx_offset = stat_i * stat_stride;
646 table_out << makeTable(stat, dest_idx_offset, eta_stride,
647 eta_labels,
648 categories)
649 .columnWidth(10)
650 // only dump the footer for the last eta bin i.e. total
651 .dumpHeader(stat_i == 0)
652 .dumpFooter(stat_i + 1 == kNStat)
653 .separateLastRow(true) // separate the sum of all eta bins
654 .minLabelWidth(max_label_width)
655 .labelPrefix(stat_labels.at(stat_i));
656 }
657 }
658 else
659 {
660 // dump one table with one row per counter showing the total eta range
661 for (std::size_t eta_bin_i = (m_dumpAllStatEtaBins.value() ? 0 : m_statEtaBins.size() + 1);
662 eta_bin_i < m_statEtaBins.size() + 2;
663 ++eta_bin_i)
664 {
665 std::size_t dest_idx_offset = eta_bin_i * eta_stride;
666 table_out << makeTable(stat, dest_idx_offset, stat_stride,
667 stat_labels,
668 categories,
669 eta_labels.at(eta_bin_i))
670 .columnWidth(10)
671 // only dump the footer for the last eta bin i.e. total
672 .dumpFooter(!m_dumpAllStatEtaBins.value() || eta_bin_i == m_statEtaBins.size() + 1);
673 }
674 }
675 ATH_MSG_INFO("statistics:\n"
676 << table_out.str());
677 table_out.str("");
678
679 // define retios first element numerator, second element denominator
680 // each element contains a vector of counter and a multiplier e.g. +- 1
681 // ratios are computed as (sum_i stat[stat_i] * multiplier_i ) / (sum_j stat[stat_j] * multiplier_j )
682 auto [ratio_labels, ratio_def] =
684 std::vector<TableUtils::SummandDefinition>{
688 // no track counted as used but want to include it as failed
690 }, // failed seeds i.e. seeds which are not duplicates but did not produce a track
691 std::vector<TableUtils::SummandDefinition>{TableUtils::defineSummand(kNTotalSeeds, 1)}),
693 TableUtils::defineSimpleRatio("Rejected refined params / seeds", kNRejectedRefinedSeeds, kNTotalSeeds),
695 TableUtils::defineSimpleRatio("selected tracks / used seeds", kNSelectedTracks, kNUsedSeeds),
697 TableUtils::defineSimpleRatio("resolved tracks / used seeds", kNResolvedTracks, kNUsedSeeds),
698 TableUtils::defineSimpleRatio("branched tracks / used seeds", kMultipleBranches, kNUsedSeeds),
699 TableUtils::defineSimpleRatio("no 2nd CKF / CKF tracks", kNoSecond, kNOutputTracks),
701 TableUtils::defineSimpleRatio("forced measurements / used seeds", kNForcedSeedMeasurements, kNUsedSeeds)});
702
703 std::vector<float> ratio = TableUtils::computeRatios(ratio_def,
704 nSeedCollections() + 1,
705 m_statEtaBins.size() + 2,
706 stat);
707
708 // the extra columns and rows for the projections are _not_ added internally
709 std::size_t ratio_stride = TableUtils::ratioStride(nSeedCollections() + 1,
710 m_statEtaBins.size() + 2,
711 ratio_def);
712 std::size_t ratio_eta_stride = TableUtils::subCategoryStride(nSeedCollections() + 1,
713 m_statEtaBins.size() + 2,
714 ratio_def);
715
716 std::size_t max_label_width = TableUtils::maxLabelWidth(ratio_labels) + TableUtils::maxLabelWidth(eta_labels);
717 if (m_dumpAllStatEtaBins.value())
718 {
719 // show for each ratio a table with one row per eta bin
720 for (std::size_t ratio_i = 0; ratio_i < ratio_labels.size(); ++ratio_i)
721 {
722 table_out << makeTable(ratio,
723 ratio_i * ratio_stride,
724 ratio_eta_stride,
725 eta_labels,
726 categories)
727 .columnWidth(10)
728 // only dump the footer for the last eta bin i.e. total
729 .dumpHeader(ratio_i == 0)
730 .dumpFooter(ratio_i + 1 == ratio_labels.size())
731 .separateLastRow(true) // separate the sum of las
732 .minLabelWidth(max_label_width)
733 .labelPrefix(ratio_labels.at(ratio_i));
734 }
735 }
736 else
737 {
738 // dump one table with one row per ratio showing the total eta range
739 table_out << makeTable(ratio,
740 (m_statEtaBins.size() + 1) * ratio_eta_stride + 0 * ratio_stride,
741 ratio_stride,
742 ratio_labels,
743 categories)
744 .columnWidth(10)
745 // only dump the footer for the last eta bin i.e. total
746 .minLabelWidth(max_label_width)
747 .dumpFooter(false);
748
749 // also dump a table for final tracks over used seeds (ratio_i==6 or 4) showing one row per eta bin
750 eta_labels.erase(eta_labels.end() - 1); // drop last line of table which shows again all eta bins summed.
751 std::size_t ratio_i = m_showResolvedStats ? 6 : 4;
752 table_out << makeTable(ratio,
753 ratio_i * ratio_stride,
754 ratio_eta_stride,
755 eta_labels,
756 categories)
757 .columnWidth(10)
758 .dumpHeader(false)
759 // only dump the footer for the last eta bin i.e. total
760 .dumpFooter(!m_dumpAllStatEtaBins.value() || ratio_i + 1 == ratio_labels.size())
761 .separateLastRow(false)
762 .minLabelWidth(max_label_width)
763 .labelPrefix(ratio_labels.at(ratio_i));
764 }
765
766 ATH_MSG_INFO("Ratios:\n"
767 << table_out.str());
768 }
769 }
#define ATH_MSG_INFO(x)
TableUtils::StatTable< T > makeTable(const std::array< T, N > &counter, const std::array< std::string, N > &label)
Definition TableUtils.h:543
Gaudi::Property< bool > m_dumpAllStatEtaBins
bool msgLvl(const MSG::Level lvl) const
std::tuple< std::vector< std::string >, std::vector< RatioDefinition > > splitRatioDefinitionsAndLabels(std::initializer_list< std::tuple< std::string, RatioDefinition > > a_ratio_list)
Definition TableUtils.h:468
SummandDefinition defineSummand(T counter_idx, int multiplier)
Definition TableUtils.h:422
std::size_t maxLabelWidth(const T_Collection &col)
Definition TableUtils.h:310
RatioDefinition defineSimpleRatio(T numerator, T denominator)
Definition TableUtils.h:404
std::vector< std::string > makeLabelVector(T_index n_entries, std::initializer_list< std::pair< T_index, T_string > > a_list)
Definition TableUtils.h:295
constexpr std::size_t subCategoryStride(const std::size_t categories, const std::size_t sub_categories, const std::size_t n_counter)
Definition TableUtils.h:324
constexpr std::size_t counterStride(const std::size_t categories, const std::size_t sub_categories, const std::size_t n_counter)
Definition TableUtils.h:329
std::string makeEtaBinLabel(const std::vector< float > &eta_bins, std::size_t eta_bin_i, bool abs_eta=false)
Definition TableUtils.h:534
std::tuple< std::string, RatioDefinition > makeRatioDefinition(std::string &&name, std::vector< SummandDefinition > &&numerator, std::vector< SummandDefinition > &&denominator)
Definition TableUtils.h:458
std::vector< float > computeRatios(const std::vector< RatioDefinition > &ratio_def, const std::size_t categories, const std::size_t sub_categories, const std::vector< std::size_t > &counter)
constexpr std::size_t ratioStride(const std::size_t categories, const std::size_t sub_categories, const std::vector< RatioDefinition > &ratio_def)
Definition TableUtils.h:493
std::vector< T_Output > createCounterArrayWithProjections(const std::size_t categories, const std::size_t sub_categories, const std::vector< std::array< T_Input, N > > &input_counts)
Definition TableUtils.h:342

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ seedCollectionStride()

std::size_t ActsTrk::TrackFindingBaseAlg::seedCollectionStride ( ) const
inlineprotected

Definition at line 345 of file TrackFindingBaseAlg.h.

345 {
346 return m_statEtaBins.size() + 1;
347 }

◆ selectCounts()

std::array< bool, 3 > ActsTrk::TrackFindingBaseAlg::selectCounts ( const detail::RecoTrackContainer::TrackProxy & track,
double eta ) const
protected

Definition at line 493 of file TrackFindingBaseAlg.cxx.

493 {
494 bool enoughMeasurements = true, tooManyHoles = false, tooManyOutliers = false;
495 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
496 std::size_t etaBin = (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? 0
497 : (std::abs(eta) >= trackSelectorCfg.absEtaEdges.back()) ? trackSelectorCfg.absEtaEdges.size() - 1
498 : trackSelectorCfg.binIndex(eta);
499 auto cutMin = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
500 return !cutSet.empty() && (val < (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
501 };
502 auto cutMax = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
503 return !cutSet.empty() && (val > (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
504 };
505
506 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nPixelHits(track), m_minPixelHits);
507 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nStripHits(track), m_minStripHits);
508 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nHgtdHits(track), m_minHgtdHits);
509 tooManyHoles = tooManyHoles || cutMax(s_branchState.nPixelHoles(track), m_maxPixelHoles);
510 tooManyHoles = tooManyHoles || cutMax(s_branchState.nStripHoles(track), m_maxStripHoles);
511 tooManyHoles = tooManyHoles || cutMax(s_branchState.nHgtdHoles(track), m_maxHgtdHoles);
512 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nPixelOutliers(track), m_maxPixelOutliers);
513 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nStripOutliers(track), m_maxStripOutliers);
514 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nHgtdOutliers(track), m_maxHgtdOutliers);
515
516 return {enoughMeasurements, tooManyHoles, tooManyOutliers};
517 }
size_t size() const
Number of registered mappings.
Gaudi::Property< std::vector< std::size_t > > m_maxHgtdHoles
Gaudi::Property< std::vector< std::size_t > > m_maxHgtdOutliers
Gaudi::Property< std::vector< std::size_t > > m_minHgtdHits
return m_collEvts back().back().max_entries
unsigned int binIndex(const T &val, const std::vector< T > &partitions)
general utility function to return bin index given a value and the upper endpoints of each bin
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap etaBin

◆ selectCountsFinal()

bool ActsTrk::TrackFindingBaseAlg::selectCountsFinal ( const detail::RecoTrackContainer::TrackProxy & track) const
protected

Definition at line 792 of file TrackFindingBaseAlg.cxx.

792 {
793 if (not m_addCounts) return true;
794 double eta = -std::log(std::tan(0.5 * track.theta()));
795 auto [enoughMeasurementsPS, tooManyHolesPS, tooManyOutliersPS] = selectCounts(track, eta);
796 return enoughMeasurementsPS && !tooManyHolesPS && !tooManyOutliersPS;
797 }
std::array< bool, 3 > selectCounts(const detail::RecoTrackContainer::TrackProxy &track, double eta) const

◆ setFilterPassed()

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

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

99 {
101 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const
Set filter decision:

◆ setMeasurementSelector()

std::unique_ptr< ActsTrk::IMeasurementSelector > ActsTrk::TrackFindingBaseAlg::setMeasurementSelector ( const EventContext & ctx,
const detail::TrackFindingMeasurements & measurements,
TrackFinderOptions & options ) const
nodiscardprotected

Setup and attach measurement selector to KF options.

Parameters
ctxthe current event context
measurementsmeasurements container used in MeasurementSelector
optionsKalman filter options
Returns
unique_ptr to MeasurementSelector

Definition at line 185 of file TrackFindingBaseAlg.cxx.

188 {
189
190 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = ActsTrk::detail::getMeasurementSelector(
191 ctx,
192 m_pixelCalibTool.isEnabled() ? &(*m_pixelCalibTool) : nullptr,
193 m_stripCalibTool.isEnabled() ? &(*m_stripCalibTool) : nullptr,
194 m_hgtdCalibTool.isEnabled() ? &(*m_hgtdCalibTool) : nullptr,
195 measurements.measurementRanges(),
197 m_measurementSelectorConfig.m_chi2CutOffOutlier,
199 m_edgeHoleBorderWidth.value());
200
201 measurementSelector->connect(&options.extensions.createTrackStates);
202
203 return measurementSelector;
204 }
Gaudi::Property< double > m_edgeHoleBorderWidth
std::unique_ptr< ActsTrk::IMeasurementSelector > getMeasurementSelector(const EventContext &ctx, const ActsTrk::IPixelOnTrackCalibratorTool< detail::RecoTrackStateContainer > *pixelOnTrackCalibratorTool, const ActsTrk::IStripOnTrackCalibratorTool< detail::RecoTrackStateContainer > *stripOnTrackCalibratorTool, const ActsTrk::IHGTDOnTrackCalibratorTool< detail::RecoTrackStateContainer > *hgtdOnTrackCalibratorTool, 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 edge_hole_border_width)

◆ stopBranch()

TrackFindingBaseAlg::BranchStopperResult ActsTrk::TrackFindingBaseAlg::stopBranch ( const detail::RecoTrackContainer::TrackProxy & track,
const detail::RecoTrackContainer::TrackStateProxy & trackState,
const Acts::TrackSelector::EtaBinnedConfig & trackSelectorCfg,
const Acts::GeometryContext & tgContext,
const detail::MeasurementIndex & measurementIndex,
const std::size_t typeIndex,
EventStats::value_type & event_stat_category_i ) const
protected

Branch stopper.

Parameters
trackTrack proxy object
trackStateTrack state proxy object
trackSelectorCfgTrack selector configuration
tgContextGeometry context
measurementIndexMeasurement index
typeIndexType index
event_stat_category_iEvent statistics for current category
Returns
BranchStopperResult

Definition at line 306 of file TrackFindingBaseAlg.cxx.

313 {
314 if (m_addCounts) {
315 updateCounts(track, trackState.typeFlags(),
316 measurementType(trackState));
317 if (m_checkCounts) {
318 checkCounts(track);
319 }
320 }
321
322 if (m_trackStatePrinter.isSet()) {
323 m_trackStatePrinter->printTrackState(tgContext, trackState,
324 measurementIndex, true);
325 }
326
327 if (!m_doBranchStopper) {
328 return BranchStopperResult::Continue;
329 }
330
331 const auto &parameters = trackState.hasFiltered() ? trackState.filtered()
332 : trackState.predicted();
333 double eta = -std::log(std::tan(0.5 * parameters[Acts::eBoundTheta]));
334 const auto &cutSet = getCuts(eta);
335
336 if (typeIndex < m_ptMinMeasurements.size() &&
337 !(track.nMeasurements() < m_ptMinMeasurements[typeIndex])) {
338 double pT = std::sin(parameters[Acts::eBoundTheta]) /
339 parameters[Acts::eBoundQOverP];
340 if (std::abs(pT) < cutSet.ptMin * m_branchStopperPtMinFactor) {
341 ++event_stat_category_i[kNStoppedTracksMinPt];
342 ATH_MSG_DEBUG("CkfBranchStopper: drop branch with q*pT="
343 << pT << " after " << track.nMeasurements()
344 << " measurements");
345 return BranchStopperResult::StopAndDrop;
346 }
347 }
348
349 if (typeIndex < m_absEtaMaxMeasurements.size() &&
350 !(track.nMeasurements() < m_absEtaMaxMeasurements[typeIndex]) &&
351 !(std::abs(eta) < trackSelectorCfg.absEtaEdges.back() +
353 ++event_stat_category_i[kNStoppedTracksMaxEta];
354 ATH_MSG_DEBUG("CkfBranchStopper: drop branch with eta="
355 << eta << " after " << track.nMeasurements()
356 << " measurements");
357 return BranchStopperResult::StopAndDrop;
358 }
359
360
361 // In the pixel endcap regions relax the requirement for minMeasurements before cutting the branch off
362 auto minMeasurementsBranchStop = std::abs(eta) > m_branchStopperAbsEtaMeasCut ? cutSet.minMeasurements - m_branchStopperMeasCutReduce : cutSet.minMeasurements;
363 bool enoughMeasurements = (track.nMeasurements() >= minMeasurementsBranchStop);
364 bool tooManyHoles = (track.nHoles() > cutSet.maxHoles);
365 bool tooManyOutliers = (track.nOutliers() > cutSet.maxOutliers);
366
367 if (m_addCounts) {
368 auto [enoughMeasurementsPS, tooManyHolesPS, tooManyOutliersPS] =
369 selectCounts(track, eta);
370 enoughMeasurements = enoughMeasurements && enoughMeasurementsPS;
371 tooManyHoles = tooManyHoles || tooManyHolesPS;
372 tooManyOutliers = tooManyOutliers || tooManyOutliersPS;
373 }
374
375 if (!(tooManyHoles || tooManyOutliers)) {
376 return BranchStopperResult::Continue;
377 }
378
379 if (!enoughMeasurements) {
380 ++event_stat_category_i[kNStoppedTracksMaxHoles];
381 }
382
383 if (m_addCounts) {
384 ATH_MSG_DEBUG("CkfBranchStopper: stop and "
385 << (enoughMeasurements ? "keep" : "drop")
386 << " branch with nHoles=" << track.nHoles() << " ("
387 << s_branchState.nPixelHoles(track) << " pixel+"
388 << s_branchState.nStripHoles(track) << " strip+"
389 << s_branchState.nHgtdHoles(track)
390 << " hgtd), nOutliers=" << track.nOutliers() << " ("
391 << s_branchState.nPixelOutliers(track) << "+"
392 << s_branchState.nStripOutliers(track) << "+"
393 << s_branchState.nHgtdOutliers(track)
394 << "), nMeasurements=" << track.nMeasurements() << " ("
395 << s_branchState.nPixelHits(track) << "+"
396 << s_branchState.nStripHits(track) << "+"
397 << s_branchState.nHgtdHits(track) << ")");
398 } else {
399 ATH_MSG_DEBUG("CkfBranchStopper: stop and "
400 << (enoughMeasurements ? "keep" : "drop")
401 << " branch with nHoles=" << track.nHoles()
402 << ", nOutliers=" << track.nOutliers()
403 << ", nMeasurements=" << track.nMeasurements());
404 }
405
406 return enoughMeasurements ? BranchStopperResult::StopAndKeep
407 : BranchStopperResult::StopAndDrop;
408 }
static xAOD::UncalibMeasType measurementType(const detail::RecoTrackContainer::TrackStateProxy &trackState)
void checkCounts(const detail::RecoTrackContainer::TrackProxy &track) const
const Acts::TrackSelector::Config & getCuts(double eta) const
Retrieves track selector configuration for given eta value.
static void updateCounts(const detail::RecoTrackContainer::TrackProxy &track, Acts::ConstTrackStateTypeMap typeFlags, xAOD::UncalibMeasType detType)

◆ sysExecute()

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

Execute an algorithm.

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

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

41{
42 return BaseAlg::sysExecute (ctx);
43}

◆ sysInitialize()

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

Override sysInitialize.

Override sysInitialize from the base class.

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

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

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

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

71 {
73
74 if (sc.isFailure()) {
75 return sc;
76 }
77
78 ServiceHandle<ICondSvc> cs("CondSvc",name());
79 for (auto h : outputHandles()) {
80 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
81 // do this inside the loop so we don't create the CondSvc until needed
82 if ( cs.retrieve().isFailure() ) {
83 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
85 }
86 if (cs->regHandle(this,*h).isFailure()) {
88 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
89 << " with CondSvc");
90 }
91 }
92 }
93 return sc;
94}
virtual StatusCode sysInitialize() override
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

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

◆ trackFinder() [1/2]

TrackFindingBaseAlg::CKF_pimpl & ActsTrk::TrackFindingBaseAlg::trackFinder ( )
protected

Definition at line 21 of file TrackFindingBaseAlg.cxx.

21{ return *m_trackFinder; }

◆ trackFinder() [2/2]

const TrackFindingBaseAlg::CKF_pimpl & ActsTrk::TrackFindingBaseAlg::trackFinder ( ) const
protected

Definition at line 22 of file TrackFindingBaseAlg.cxx.

22{ return *m_trackFinder; }

◆ updateCounts()

void ActsTrk::TrackFindingBaseAlg::updateCounts ( const detail::RecoTrackContainer::TrackProxy & track,
Acts::ConstTrackStateTypeMap typeFlags,
xAOD::UncalibMeasType detType )
staticprotected

Definition at line 440 of file TrackFindingBaseAlg.cxx.

442 {
444 if (typeFlags.isHole()) {
445 s_branchState.nPixelHoles(track)++;
446 } else if (typeFlags.isOutlier()) {
447 s_branchState.nPixelOutliers(track)++;
448 } else if (typeFlags.isMeasurement()) {
449 s_branchState.nPixelHits(track)++;
450 }
451 } else if (detType == xAOD::UncalibMeasType::StripClusterType) {
452 if (typeFlags.isHole()) {
453 s_branchState.nStripHoles(track)++;
454 } else if (typeFlags.isOutlier()) {
455 s_branchState.nStripOutliers(track)++;
456 } else if (typeFlags.isMeasurement()) {
457 s_branchState.nStripHits(track)++;
458 }
459 } else if (detType == xAOD::UncalibMeasType::HGTDClusterType) {
460 if (typeFlags.isHole()) {
461 s_branchState.nHgtdHoles(track)++;
462 } else if (typeFlags.isOutlier()) {
463 s_branchState.nHgtdOutliers(track)++;
464 } else if (typeFlags.isMeasurement()) {
465 s_branchState.nHgtdHits(track)++;
466 }
467 }
468 }

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

std::vector<std::array<std::size_t, kNStat> > m_stat ActsTrk::TrackFindingBaseAlg::ATLAS_THREAD_SAFE {}
mutableprotected

Definition at line 354 of file TrackFindingBaseAlg.h.

354{};

◆ ATLAS_THREAD_SAFE [2/2]

std::mutex m_mutex ActsTrk::TrackFindingBaseAlg::ATLAS_THREAD_SAFE
mutableprotected

Definition at line 353 of file TrackFindingBaseAlg.h.

◆ m_absEtaMax

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_absEtaMax {this, "absEtaMax", std::numeric_limits<double>::max(), "TrackSelector: absEtaMax"}
protected

Definition at line 128 of file TrackFindingBaseAlg.h.

128{this, "absEtaMax", std::numeric_limits<double>::max(), "TrackSelector: absEtaMax"};

◆ m_absEtaMaxMeasurements

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_absEtaMaxMeasurements {this, "absEtaMaxMeasurements", {}, "if specified for the given seed collection, applies absEtaMax cut in branch stopper once absEtaMaxMeasurements have been encountered"}
protected

Definition at line 110 of file TrackFindingBaseAlg.h.

110{this, "absEtaMaxMeasurements", {}, "if specified for the given seed collection, applies absEtaMax cut in branch stopper once absEtaMaxMeasurements have been encountered"};

◆ m_absEtaMin

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_absEtaMin {this, "absEtaMin", 0.0, "TrackSelector: absEtaMin"}
protected

Definition at line 127 of file TrackFindingBaseAlg.h.

127{this, "absEtaMin", 0.0, "TrackSelector: absEtaMin"};

◆ m_addCounts

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_addCounts {this, "addCounts", true, "keep separate pixel, strip and hgtd counts and apply the following cuts"}
protected

Definition at line 142 of file TrackFindingBaseAlg.h.

142{this, "addCounts", true, "keep separate pixel, strip and hgtd counts and apply the following cuts"};

◆ m_branchStopperAbsEtaMaxExtra

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_branchStopperAbsEtaMaxExtra {this, "branchStopperAbsEtaMaxExtra", 0.0, "increase absEtaMax cut when used in the branch stopper"}
protected

Definition at line 117 of file TrackFindingBaseAlg.h.

117{this, "branchStopperAbsEtaMaxExtra", 0.0, "increase absEtaMax cut when used in the branch stopper"};

◆ m_branchStopperAbsEtaMeasCut

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_branchStopperAbsEtaMeasCut {this, "branchStopperAbsEtaMeasCut", 1.2, "the minimum |eta| to apply the reduction to the minMeas requirement for the branch stopper"}
protected

Definition at line 119 of file TrackFindingBaseAlg.h.

119{this, "branchStopperAbsEtaMeasCut", 1.2, "the minimum |eta| to apply the reduction to the minMeas requirement for the branch stopper"};

◆ m_branchStopperMeasCutReduce

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_branchStopperMeasCutReduce {this, "branchStopperMeasCutReduce", 2, "how much to reduce the minMeas requirement for the branch stopper"}
protected

Definition at line 118 of file TrackFindingBaseAlg.h.

118{this, "branchStopperMeasCutReduce", 2, "how much to reduce the minMeas requirement for the branch stopper"};

◆ m_branchStopperPtMinFactor

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_branchStopperPtMinFactor {this, "branchStopperPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the branch stopper"}
protected

Definition at line 115 of file TrackFindingBaseAlg.h.

115{this, "branchStopperPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the branch stopper"};

◆ m_checkCounts

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_checkCounts {this, "checkCounts", false, "check consistency among track state counts"}
protected

Definition at line 143 of file TrackFindingBaseAlg.h.

143{this, "checkCounts", false, "check consistency among track state counts"};

◆ m_chi2CutOff

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_chi2CutOff {this, "chi2CutOff", {}, "MeasurementSelector: maximum local chi2 contribution"}
protected

Definition at line 105 of file TrackFindingBaseAlg.h.

105{this, "chi2CutOff", {}, "MeasurementSelector: maximum local chi2 contribution"};

◆ m_chi2OutlierCutOff

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_chi2OutlierCutOff {this, "chi2OutlierCutOff", {}, "MeasurementSelector: maximum local chi2 contribution for outlier"}
protected

Definition at line 106 of file TrackFindingBaseAlg.h.

106{this, "chi2OutlierCutOff", {}, "MeasurementSelector: maximum local chi2 contribution for outlier"};

◆ m_d0Max

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_d0Max {this, "d0Max", {}, "TrackSelector: d0Max"}
protected

Definition at line 132 of file TrackFindingBaseAlg.h.

132{this, "d0Max", {}, "TrackSelector: d0Max"};

◆ m_d0Min

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_d0Min {this, "d0Min", {}, "TrackSelector: d0Min"}
protected

Definition at line 131 of file TrackFindingBaseAlg.h.

131{this, "d0Min", {}, "TrackSelector: d0Min"};

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

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_doBranchStopper {this, "doBranchStopper", true, "use branch stopper"}
protected

Definition at line 111 of file TrackFindingBaseAlg.h.

111{this, "doBranchStopper", true, "use branch stopper"};

◆ m_doTwoWay

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_doTwoWay {this, "doTwoWay", true, "run CKF twice, first with forward propagation with smoothing, then with backward propagation"}
protected

Definition at line 112 of file TrackFindingBaseAlg.h.

112{this, "doTwoWay", true, "run CKF twice, first with forward propagation with smoothing, then with backward propagation"};

◆ m_dumpAllStatEtaBins

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_dumpAllStatEtaBins {this, "DumpEtaBinsForAll", false, "Dump eta bins of all statistics counter."}
protected

Definition at line 159 of file TrackFindingBaseAlg.h.

159{this, "DumpEtaBinsForAll", false, "Dump eta bins of all statistics counter."};

◆ m_edgeHoleBorderWidth

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_edgeHoleBorderWidth {this,"EdgeHoleBorderWidth", 0. ,"Width of the border of surfaces (mm) within which holes are not counted as holes."}
protected

Definition at line 108 of file TrackFindingBaseAlg.h.

108{this,"EdgeHoleBorderWidth", 0. ,"Width of the border of surfaces (mm) within which holes are not counted as holes."};

◆ m_endOfWorldVolumeIds

Gaudi::Property<std::vector<std::uint32_t> > ActsTrk::TrackFindingBaseAlg::m_endOfWorldVolumeIds {this, "EndOfTheWorldVolumeIds", {}, ""}
protected

Definition at line 154 of file TrackFindingBaseAlg.h.

154{this, "EndOfTheWorldVolumeIds", {}, ""};

◆ m_etaBins

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_etaBins {this, "etaBins", {}, "bins in |eta| to specify variable selections"}
protected

Definition at line 103 of file TrackFindingBaseAlg.h.

103{this, "etaBins", {}, "bins in |eta| to specify variable selections"};

◆ m_etaMax

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_etaMax {this, "etaMax", {}, "TrackSelector: etaMax"}
protected

Definition at line 126 of file TrackFindingBaseAlg.h.

126{this, "etaMax", {}, "TrackSelector: etaMax"};

◆ m_etaMin

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_etaMin {this, "etaMin", {}, "TrackSelector: etaMin"}
protected

Definition at line 125 of file TrackFindingBaseAlg.h.

125{this, "etaMin", {}, "TrackSelector: etaMin"};

◆ 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 AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_extrapolationTool

ToolHandle<ActsTrk::IExtrapolationTool> ActsTrk::TrackFindingBaseAlg::m_extrapolationTool {this, "ExtrapolationTool", ""}
protected

Definition at line 90 of file TrackFindingBaseAlg.h.

90{this, "ExtrapolationTool", ""};

◆ m_fitterTool

ToolHandle<ActsTrk::IFitterTool> ActsTrk::TrackFindingBaseAlg::m_fitterTool {this, "FitterTool", "", "Fitter Tool for Seeds"}
protected

Definition at line 93 of file TrackFindingBaseAlg.h.

93{this, "FitterTool", "", "Fitter Tool for Seeds"};

◆ m_hgtdCalibTool

ToolHandle<ActsTrk::IHGTDOnTrackCalibratorTool<detail::RecoTrackStateContainer> > ActsTrk::TrackFindingBaseAlg::m_hgtdCalibTool {this, "HGTDCalibrator", "", "Opt. HGTD measurement calibrator"}
protected

Definition at line 96 of file TrackFindingBaseAlg.h.

96{this, "HGTDCalibrator", "", "Opt. HGTD measurement calibrator"};

◆ m_inflateCovarianceTwoWay

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_inflateCovarianceTwoWay {this, "inflateCovarianceTwoWay", false, "inflate covariance matrix at the beginning of two-way track finding"}
protected

Definition at line 113 of file TrackFindingBaseAlg.h.

113{this, "inflateCovarianceTwoWay", false, "inflate covariance matrix at the beginning of two-way track finding"};

◆ m_logger

std::unique_ptr<const Acts::Logger> ActsTrk::TrackFindingBaseAlg::m_logger
protected

logging instance

Definition at line 335 of file TrackFindingBaseAlg.h.

◆ m_maxChi2

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_maxChi2 {this, "maxChi2", {}, "TrackSelector: maxChi2"}
protected

Definition at line 140 of file TrackFindingBaseAlg.h.

140{this, "maxChi2", {}, "TrackSelector: maxChi2"};

◆ m_maxHgtdHoles

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxHgtdHoles {this, "maxHgtdHoles", {}, "maximum number of hgtd holes"}
protected

Definition at line 149 of file TrackFindingBaseAlg.h.

149{this, "maxHgtdHoles", {}, "maximum number of hgtd holes"};

◆ m_maxHgtdOutliers

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxHgtdOutliers {this, "maxHgtdOutliers", {}, "maximum number of hgtd outliers"}
protected

Definition at line 152 of file TrackFindingBaseAlg.h.

152{this, "maxHgtdOutliers", {}, "maximum number of hgtd outliers"};

◆ m_maxHoles

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxHoles {this, "maxHoles", {}, "TrackSelector: maxHoles"}
protected

Definition at line 137 of file TrackFindingBaseAlg.h.

137{this, "maxHoles", {}, "TrackSelector: maxHoles"};

◆ m_maxOutliers

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxOutliers {this, "maxOutliers", {}, "TrackSelector: maxOutliers"}
protected

Definition at line 138 of file TrackFindingBaseAlg.h.

138{this, "maxOutliers", {}, "TrackSelector: maxOutliers"};

◆ m_maxPixelHoles

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxPixelHoles {this, "maxPixelHoles", {}, "maximum number of pixel holes"}
protected

Definition at line 147 of file TrackFindingBaseAlg.h.

147{this, "maxPixelHoles", {}, "maximum number of pixel holes"};

◆ m_maxPixelOutliers

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxPixelOutliers {this, "maxPixelOutliers", {}, "maximum number of pixel outliers"}
protected

Definition at line 150 of file TrackFindingBaseAlg.h.

150{this, "maxPixelOutliers", {}, "maximum number of pixel outliers"};

◆ m_maxPropagationStep

Gaudi::Property<unsigned int> ActsTrk::TrackFindingBaseAlg::m_maxPropagationStep {this, "maxPropagationStep", 1000, "Maximum number of steps for one propagate call"}
protected

Definition at line 102 of file TrackFindingBaseAlg.h.

102{this, "maxPropagationStep", 1000, "Maximum number of steps for one propagate call"};

◆ m_maxSharedHits

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxSharedHits {this, "maxSharedHits", {}, "TrackSelector: maxSharedHits"}
protected

Definition at line 139 of file TrackFindingBaseAlg.h.

139{this, "maxSharedHits", {}, "TrackSelector: maxSharedHits"};

◆ m_maxStripHoles

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxStripHoles {this, "maxStripHoles", {}, "maximum number of strip holes"}
protected

Definition at line 148 of file TrackFindingBaseAlg.h.

148{this, "maxStripHoles", {}, "maximum number of strip holes"};

◆ m_maxStripOutliers

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxStripOutliers {this, "maxStripOutliers", {}, "maximum number of strip outliers"}
protected

Definition at line 151 of file TrackFindingBaseAlg.h.

151{this, "maxStripOutliers", {}, "maximum number of strip outliers"};

◆ m_measurementSelectorConfig

struct ActsTrk::TrackFindingBaseAlg::MeasurementSelectorConfig ActsTrk::TrackFindingBaseAlg::m_measurementSelectorConfig
protected

◆ m_minHgtdHits

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_minHgtdHits {this, "minHgtdHits", {}, "minimum number of hgtd hits"}
protected

Definition at line 146 of file TrackFindingBaseAlg.h.

146{this, "minHgtdHits", {}, "minimum number of hgtd hits"};

◆ m_minMeasurements

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_minMeasurements {this, "minMeasurements", {}, "TrackSelector: minMeasurements"}
protected

Definition at line 136 of file TrackFindingBaseAlg.h.

136{this, "minMeasurements", {}, "TrackSelector: minMeasurements"};

◆ m_minPixelHits

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_minPixelHits {this, "minPixelHits", {}, "minimum number of pixel hits"}
protected

Definition at line 144 of file TrackFindingBaseAlg.h.

144{this, "minPixelHits", {}, "minimum number of pixel hits"};

◆ m_minStripHits

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_minStripHits {this, "minStripHits", {}, "minimum number of strip hits"}
protected

Definition at line 145 of file TrackFindingBaseAlg.h.

145{this, "minStripHits", {}, "minimum number of strip hits"};

◆ m_monTool

ToolHandle<GenericMonitoringTool> ActsTrk::TrackFindingBaseAlg::m_monTool {this, "MonTool", "", "Monitoring tool"}
protected

Definition at line 89 of file TrackFindingBaseAlg.h.

89{this, "MonTool", "", "Monitoring tool"};

◆ m_numMeasurementsCutOff

Gaudi::Property<std::vector<size_t> > ActsTrk::TrackFindingBaseAlg::m_numMeasurementsCutOff {this, "numMeasurementsCutOff", {}, "MeasurementSelector: maximum number of associated measurements on a single surface"}
protected

Definition at line 107 of file TrackFindingBaseAlg.h.

107{this, "numMeasurementsCutOff", {}, "MeasurementSelector: maximum number of associated measurements on a single surface"};

◆ m_phiMax

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_phiMax {this, "phiMax", {}, "TrackSelector: phiMax"}
protected

Definition at line 124 of file TrackFindingBaseAlg.h.

124{this, "phiMax", {}, "TrackSelector: phiMax"};

◆ m_phiMin

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_phiMin {this, "phiMin", {}, "TrackSelector: phiMin"}
protected

Definition at line 123 of file TrackFindingBaseAlg.h.

123{this, "phiMin", {}, "TrackSelector: phiMin"};

◆ m_pixelCalibTool

ToolHandle<ActsTrk::IPixelOnTrackCalibratorTool<detail::RecoTrackStateContainer> > ActsTrk::TrackFindingBaseAlg::m_pixelCalibTool {this, "PixelCalibrator", "", "Opt. pixel measurement calibrator"}
protected

Definition at line 94 of file TrackFindingBaseAlg.h.

94{this, "PixelCalibrator", "", "Opt. pixel measurement calibrator"};

◆ m_ptMax

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_ptMax {this, "ptMax", {}, "TrackSelector: ptMax"}
protected

Definition at line 130 of file TrackFindingBaseAlg.h.

130{this, "ptMax", {}, "TrackSelector: ptMax"};

◆ m_ptMin

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_ptMin {this, "ptMin", {}, "TrackSelector: ptMin"}
protected

Definition at line 129 of file TrackFindingBaseAlg.h.

129{this, "ptMin", {}, "TrackSelector: ptMin"};

◆ m_ptMinMeasurements

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_ptMinMeasurements {this, "ptMinMeasurements", {}, "if specified for the given seed collection, applies ptMin cut in branch stopper once ptMinMinMeasurements have been encountered"}
protected

Definition at line 109 of file TrackFindingBaseAlg.h.

109{this, "ptMinMeasurements", {}, "if specified for the given seed collection, applies ptMin cut in branch stopper once ptMinMinMeasurements have been encountered"};

◆ m_seedLabels

Gaudi::Property<std::vector<std::string> > ActsTrk::TrackFindingBaseAlg::m_seedLabels {this, "SeedLabels", {}, "One label per seed key used in outputs"}
protected

Definition at line 158 of file TrackFindingBaseAlg.h.

158{this, "SeedLabels", {}, "One label per seed key used in outputs"};

◆ m_seedRefitPtMinFactor

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_seedRefitPtMinFactor {this, "seedRefitPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the seed refit"}
protected

Definition at line 116 of file TrackFindingBaseAlg.h.

116{this, "seedRefitPtMinFactor", 1.0, "factor to multiply ptMin cut when used in the seed refit"};

◆ m_showResolvedStats

bool ActsTrk::TrackFindingBaseAlg::m_showResolvedStats = false
protected

Definition at line 352 of file TrackFindingBaseAlg.h.

◆ m_statEtaBins

Gaudi::Property<std::vector<float> > ActsTrk::TrackFindingBaseAlg::m_statEtaBins {this, "StatisticEtaBins", {-4, -2.6, -2, 0, 2., 2.6, 4}, "Gather statistics separately for these bins."}
protected

Definition at line 157 of file TrackFindingBaseAlg.h.

157{this, "StatisticEtaBins", {-4, -2.6, -2, 0, 2., 2.6, 4}, "Gather statistics separately for these bins."};

◆ m_stripCalibTool

ToolHandle<ActsTrk::IStripOnTrackCalibratorTool<detail::RecoTrackStateContainer> > ActsTrk::TrackFindingBaseAlg::m_stripCalibTool {this, "StripCalibrator", "", "Opt. strip measurement calibrator"}
protected

Definition at line 95 of file TrackFindingBaseAlg.h.

95{this, "StripCalibrator", "", "Opt. strip measurement calibrator"};

◆ m_trackContainerKey

SG::WriteHandleKey<ActsTrk::TrackContainer> ActsTrk::TrackFindingBaseAlg::m_trackContainerKey {this, "ACTSTracksLocation", "", "Output track collection (ActsTrk variant)"}
protected

Definition at line 98 of file TrackFindingBaseAlg.h.

98{this, "ACTSTracksLocation", "", "Output track collection (ActsTrk variant)"};

◆ m_trackFinder

std::unique_ptr<CKF_pimpl> ActsTrk::TrackFindingBaseAlg::m_trackFinder
protected

Definition at line 83 of file TrackFindingBaseAlg.h.

◆ m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> ActsTrk::TrackFindingBaseAlg::m_trackingGeometryTool {this, "TrackingGeometryTool", ""}
protected

Definition at line 91 of file TrackFindingBaseAlg.h.

91{this, "TrackingGeometryTool", ""};

◆ m_tracksBackendHandlesHelper

ActsTrk::MutableTrackContainerHandlesHelper ActsTrk::TrackFindingBaseAlg::m_tracksBackendHandlesHelper {this}
protected

Definition at line 99 of file TrackFindingBaseAlg.h.

99{this};

◆ m_trackStatePrinter

ToolHandle<ActsTrk::TrackStatePrinterTool> ActsTrk::TrackFindingBaseAlg::m_trackStatePrinter {this, "TrackStatePrinter", "", "optional track state printer"}
protected

Definition at line 92 of file TrackFindingBaseAlg.h.

92{this, "TrackStatePrinter", "", "optional track state printer"};

◆ m_twoWayinflateCovarianceFactor

Gaudi::Property<double> ActsTrk::TrackFindingBaseAlg::m_twoWayinflateCovarianceFactor {this, "twoWayinflateCovarianceFactor", 1.0, "factor to multiply the initial covariance matrix at the beginning of two-way track finding"}
protected

Definition at line 114 of file TrackFindingBaseAlg.h.

114{this, "twoWayinflateCovarianceFactor", 1.0, "factor to multiply the initial covariance matrix at the beginning of two-way track finding"};

◆ m_unalibMeasSurfAcc

detail::xAODUncalibMeasSurfAcc ActsTrk::TrackFindingBaseAlg::m_unalibMeasSurfAcc {}
protected

Definition at line 85 of file TrackFindingBaseAlg.h.

85{};

◆ m_useAbsEtaForStat

bool ActsTrk::TrackFindingBaseAlg::m_useAbsEtaForStat = false
protected

Definition at line 351 of file TrackFindingBaseAlg.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.

◆ m_z0Max

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_z0Max {this, "z0Max", {}, "TrackSelector: z0Max"}
protected

Definition at line 134 of file TrackFindingBaseAlg.h.

134{this, "z0Max", {}, "TrackSelector: z0Max"};

◆ m_z0Min

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_z0Min {this, "z0Min", {}, "TrackSelector: z0Min"}
protected

Definition at line 133 of file TrackFindingBaseAlg.h.

133{this, "z0Min", {}, "TrackSelector: z0Min"};

◆ s_branchState

BranchState ActsTrk::TrackFindingBaseAlg::s_branchState {}
staticconstexprprotected

Definition at line 316 of file TrackFindingBaseAlg.h.

316{};

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