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.
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
ToolHandle< GenericMonitoringToolm_monTool {this, "MonTool", "", "Monitoring tool"}
ServiceHandle< ActsTrk::ITrackingGeometrySvcm_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}
ToolHandle< ActsTrk::TrackStatePrinterToolm_trackStatePrinter {this, "TrackStatePrinter", "", "optional track state printer"}
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}
ContextUtility m_ctxProvider {this}
 Utility to fetch the geometry, magnetic field and calibration context in the event.
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_recordMaterialStates {this, "recordMaterialStates", false, "record track states on surfaces which have material but no measurement; keeping them off is faster, but they have to be on for the multi-component stepper and for anything downstream which needs the full track state sequence"}
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 42 of file TrackFindingBaseAlg.h.

Member Typedef Documentation

◆ BranchStopperResult

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

Definition at line 259 of file TrackFindingBaseAlg.h.

◆ EventStats

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

Definition at line 179 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 52 of file TrackFindingBaseAlg.h.

◆ TrkProxy

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

Definition at line 234 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 156 of file TrackFindingBaseAlg.h.

157 {
161 kNoTrack,
170 kNoSecond,
176 kNStat
177 };

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 410 of file TrackFindingBaseAlg.cxx.

411 {
412 if (addCounts) {
413 tracksContainer.addColumn<unsigned int>("nPixelHits");
414 tracksContainer.addColumn<unsigned int>("nStripHits");
415 tracksContainer.addColumn<unsigned int>("nHgtdHits");
416 tracksContainer.addColumn<unsigned int>("nPixelHoles");
417 tracksContainer.addColumn<unsigned int>("nStripHoles");
418 tracksContainer.addColumn<unsigned int>("nHgtdHoles");
419 tracksContainer.addColumn<unsigned int>("nPixelOutliers");
420 tracksContainer.addColumn<unsigned int>("nStripOutliers");
421 tracksContainer.addColumn<unsigned int>("nHgtdOutliers");
422 }
424 }
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 469 of file TrackFindingBaseAlg.cxx.

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

779 {
780 std::size_t out = 0u;
781 for (std::size_t category_i = seed_collection * seedCollectionStride();
782 category_i < (seed_collection + 1) * seedCollectionStride();
783 ++category_i)
784 {
785 assert(category_i < stat.size());
786 out += stat[category_i][counter_i];
787 }
788 return out;
789 }
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 565 of file TrackFindingBaseAlg.cxx.

565 {
566 std::lock_guard<std::mutex> lock(m_mutex);
567 std::size_t category_i = 0;
568 for (const std::array<unsigned int, kNStat> &src_stat : event_stat)
569 {
570 std::array<std::size_t, kNStat> &dest_stat = m_stat[category_i++];
571 for (std::size_t i = 0; i < src_stat.size(); ++i)
572 {
573 assert(i < dest_stat.size());
574 dest_stat[i] += src_stat[i];
575 }
576 }
577 }
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.

◆ 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 242 of file TrackFindingBaseAlg.cxx.

245 {
246 if (not m_doTwoWay) return {};
247
248 // Create initial parameters for the propagation
249 Acts::BoundTrackParameters secondInitialParameters = trackProxy.createParametersFromState(detail::RecoConstTrackStateContainerProxy{firstMeasurement});
250 if (!secondInitialParameters.referenceSurface().insideBounds(secondInitialParameters.localPosition())) { // #3751
251 return {};
252 }
253
254 // First, inflate the covariance matrix if configured
256 ATH_MSG_DEBUG("Inflating covariance matrix for second track finding with factor = " << m_twoWayinflateCovarianceFactor.value());
257 ATH_MSG_VERBOSE("Original parameters before inflation: \n" << secondInitialParameters);
258
259 auto inflatedCovariance = secondInitialParameters.covariance().value();
260 inflatedCovariance *= m_twoWayinflateCovarianceFactor;
261
262 const auto& origSurface = secondInitialParameters.referenceSurface();
263 auto surfacePtr = const_cast<Acts::Surface&>(origSurface).shared_from_this();
264
265 Acts::BoundTrackParameters newParams(
266 std::static_pointer_cast<const Acts::Surface>(std::move(surfacePtr)),
267 secondInitialParameters.parameters(),
268 std::make_optional(inflatedCovariance),
269 secondInitialParameters.particleHypothesis());
270 secondInitialParameters = std::move(newParams);
271
272 ATH_MSG_VERBOSE("Inflated covariance matrix : \n" << secondInitialParameters.covariance().value());
273 }
274
275 auto rootBranch = tracksContainerTemp.makeTrack();
276 rootBranch.copyFromWithoutStates(trackProxy); // #3534
277
278 // perform track finding
279 auto secondResult =
280 trackFinder().ckf.findTracks(secondInitialParameters, options, tracksContainerTemp, rootBranch);
281 if (not secondResult.ok()) {
282 return {};
283 }
284 return secondResult.value();
285 }
#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 170 of file TrackFindingBaseAlg.cxx.

170 {
171 ATH_MSG_FATAL("execute() method from the base class was called! Implement proper execute() method in the derived class!");
172
173 return StatusCode::FAILURE;
174 }
#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 176 of file TrackFindingBaseAlg.cxx.

176 {
178
179 return StatusCode::SUCCESS;
180 }

◆ 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 233 of file TrackFindingBaseAlg.cxx.

233 {
234 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
235 // return the last bin for |eta|>=4 or nan
236 return (!(std::abs(eta) < trackSelectorCfg.absEtaEdges.back())) ? trackSelectorCfg.cutSets.back()
237 : (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? trackSelectorCfg.cutSets.front()
238 : trackSelectorCfg.getCuts(eta);
239 };
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 203 of file TrackFindingBaseAlg.cxx.

207 {
208 Acts::PropagatorPlainOptions plainOptions{detContext.geometry, detContext.magField};
209 plainOptions.maxSteps = m_maxPropagationStep;
210 plainOptions.direction = Acts::Direction::Forward();
211 plainOptions.endOfWorldVolumeIds = m_endOfWorldVolumeIds;
212
213 // Set the CombinatorialKalmanFilter options
214 TrackFinderOptions options(detContext.geometry, detContext.magField, detContext.calib,
215 trackFinder().ckfExtensions, plainOptions, pSurface);
216 options.recordMaterialStates = m_recordMaterialStates;
217
218 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = setMeasurementSelector(ctx, measurements, options);
219
220 Acts::PropagatorPlainOptions plainSecondOptions{detContext.geometry, detContext.magField};
221 plainSecondOptions.maxSteps = m_maxPropagationStep;
222 plainSecondOptions.direction = plainOptions.direction.invert();
223
224 TrackFinderOptions secondOptions(detContext.geometry, detContext.magField, detContext.calib,
225 options.extensions, plainSecondOptions, pSurface);
226 secondOptions.targetSurface = pSurface;
227 secondOptions.skipPrePropagationUpdate = true;
228 secondOptions.recordMaterialStates = m_recordMaterialStates;
229
230 return {std::move(options), std::move(secondOptions), std::move(measurementSelector)};
231 };
Gaudi::Property< unsigned int > m_maxPropagationStep
Gaudi::Property< std::vector< std::uint32_t > > m_endOfWorldVolumeIds
Gaudi::Property< bool > m_recordMaterialStates
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 770 of file TrackFindingBaseAlg.cxx.

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

◆ initCounts()

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

Definition at line 426 of file TrackFindingBaseAlg.cxx.

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

◆ 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()}));
79 ATH_CHECK(m_ctxProvider.initialize());
80 ATH_CHECK(m_trackStatePrinter.retrieve(EnableTool{not m_trackStatePrinter.empty()}));
81 ATH_CHECK(m_pixelCalibTool.retrieve(EnableTool{not m_pixelCalibTool.empty()}));
82 ATH_CHECK(m_stripCalibTool.retrieve(EnableTool{not m_stripCalibTool.empty()}));
83 ATH_CHECK(m_hgtdCalibTool.retrieve(EnableTool{not m_hgtdCalibTool.empty()}));
84
85 auto magneticField = std::make_unique<ATLASMagneticFieldWrapper>();
86 auto trackingGeometry = m_trackingGeometrySvc->trackingGeometry();
87
88 detail::Stepper stepper(std::move(magneticField));
89 detail::Navigator::Config config{std::move(trackingGeometry)};
90 config.resolvePassive = false;
91 config.resolveMaterial = true;
92 config.resolveSensitive = true;
93 detail::Navigator navigator(std::move(config), logger().cloneWithSuffix("Navigator"));
94 detail::Propagator propagator(std::move(stepper), std::move(navigator), logger().cloneWithSuffix("Prop"));
95
96 // Using the CKF propagator as extrapolator
97 detail::Extrapolator extrapolator = propagator;
98
99 // m_etaBins (from flags.Tracking.ActiveConfig.etaBins) includes a dummy first and last bin, which we ignore
100 std::vector<double> absEtaEdges;
101 if (m_etaBins.size() <= 2)
102 {
103 absEtaEdges.reserve(2ul);
104 absEtaEdges.push_back(0.0);
105 absEtaEdges.push_back(std::numeric_limits<double>::infinity());
106 }
107 else
108 {
109 absEtaEdges.reserve(m_etaBins.size());
110 absEtaEdges.push_back(m_absEtaMin);
111 absEtaEdges.insert(absEtaEdges.end(), m_etaBins.begin() + 1, m_etaBins.end() - 1);
112 absEtaEdges.push_back(m_absEtaMax);
113 }
114
115 auto setCut = [](auto &cfgVal, const auto &cuts, size_t ind) -> void
116 {
117 if (cuts.empty())
118 return;
119 cfgVal = (ind < cuts.size()) ? cuts[ind] : cuts[cuts.size() - 1];
120 };
121
122 Acts::TrackSelector::EtaBinnedConfig trackSelectorCfg{std::move(absEtaEdges)};
123 if (m_etaBins.size() <= 2)
124 {
125 assert(trackSelectorCfg.cutSets.size() == 1);
126 trackSelectorCfg.cutSets[0].absEtaMin = m_absEtaMin;
127 trackSelectorCfg.cutSets[0].absEtaMax = m_absEtaMax;
128 }
129 size_t cutIndex = 0;
130 for (auto &cfg : trackSelectorCfg.cutSets)
131 {
132 setCut(cfg.phiMin, m_phiMin, cutIndex);
133 setCut(cfg.phiMax, m_phiMax, cutIndex);
134 setCut(cfg.etaMin, m_etaMin, cutIndex);
135 setCut(cfg.etaMax, m_etaMax, cutIndex);
136 setCut(cfg.ptMin, m_ptMin, cutIndex);
137 setCut(cfg.ptMax, m_ptMax, cutIndex);
138 setCut(cfg.loc0Min, m_d0Min, cutIndex);
139 setCut(cfg.loc0Max, m_d0Max, cutIndex);
140 setCut(cfg.loc1Min, m_z0Min, cutIndex);
141 setCut(cfg.loc1Max, m_z0Max, cutIndex);
142 setCut(cfg.minMeasurements, m_minMeasurements, cutIndex);
143 setCut(cfg.maxHoles, m_maxHoles, cutIndex);
144 setCut(cfg.maxOutliers, m_maxOutliers, cutIndex);
145 setCut(cfg.maxSharedHits, m_maxSharedHits, cutIndex);
146 setCut(cfg.maxChi2, m_maxChi2, cutIndex);
147 ++cutIndex;
148 }
149
150 ATH_MSG_DEBUG(trackSelectorCfg);
151
152 // initializer measurement selector and connect it to the delegates of the track finder optins
154
155 detail::CKF_config ckfConfig{
156 std::move(extrapolator),
157 detail::CKF{std::move(propagator), logger().cloneWithSuffix("CKF")},
158 {},
159 Acts::TrackSelector{trackSelectorCfg}};
160
161 m_trackFinder = std::make_unique<CKF_pimpl>(std::move(ckfConfig));
162
164
166
167 return StatusCode::SUCCESS;
168 }
#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)
Definition Logger.cxx:64
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
Gaudi::Property< std::vector< double > > m_z0Max
Gaudi::Property< double > m_absEtaMax
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
ServiceHandle< ActsTrk::ITrackingGeometrySvc > m_trackingGeometrySvc
ContextUtility m_ctxProvider
Utility to fetch the geometry, magnetic field and calibration context in the event.
ToolHandle< ActsTrk::IHGTDOnTrackCalibratorTool< detail::RecoTrackStateContainer > > m_hgtdCalibTool
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
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 518 of file TrackFindingBaseAlg.cxx.

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

◆ initStatTables()

void ActsTrk::TrackFindingBaseAlg::initStatTables ( )
protected

Definition at line 547 of file TrackFindingBaseAlg.cxx.

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

◆ 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 306 of file TrackFindingBaseAlg.h.

307 {
308 return *m_logger;
309 }

◆ measurementType()

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

Definition at line 287 of file TrackFindingBaseAlg.cxx.

287 {
288 if (trackState.hasReferenceSurface()) {
289 if (const auto *actsDetElem = dynamic_cast<const ISurfacePlacement*>(trackState.referenceSurface().surfacePlacement())) {
290 switch (actsDetElem->detectorType()) {
297 default:
298 break;
299 }
300 }
301 }
303 }
@ 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 319 of file TrackFindingBaseAlg.h.

319 {
320 return m_seedLabels.size();
321 }
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 580 of file TrackFindingBaseAlg.cxx.

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

◆ 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 322 of file TrackFindingBaseAlg.h.

322 {
323 return m_statEtaBins.size() + 1;
324 }

◆ selectCounts()

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

Definition at line 492 of file TrackFindingBaseAlg.cxx.

492 {
493 bool enoughMeasurements = true, tooManyHoles = false, tooManyOutliers = false;
494 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
495 std::size_t etaBin = (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? 0
496 : (std::abs(eta) >= trackSelectorCfg.absEtaEdges.back()) ? trackSelectorCfg.absEtaEdges.size() - 1
497 : trackSelectorCfg.binIndex(eta);
498 auto cutMin = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
499 return !cutSet.empty() && (val < (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
500 };
501 auto cutMax = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
502 return !cutSet.empty() && (val > (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
503 };
504
505 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nPixelHits(track), m_minPixelHits);
506 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nStripHits(track), m_minStripHits);
507 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nHgtdHits(track), m_minHgtdHits);
508 tooManyHoles = tooManyHoles || cutMax(s_branchState.nPixelHoles(track), m_maxPixelHoles);
509 tooManyHoles = tooManyHoles || cutMax(s_branchState.nStripHoles(track), m_maxStripHoles);
510 tooManyHoles = tooManyHoles || cutMax(s_branchState.nHgtdHoles(track), m_maxHgtdHoles);
511 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nPixelOutliers(track), m_maxPixelOutliers);
512 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nStripOutliers(track), m_maxStripOutliers);
513 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nHgtdOutliers(track), m_maxHgtdOutliers);
514
515 return {enoughMeasurements, tooManyHoles, tooManyOutliers};
516 }
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 791 of file TrackFindingBaseAlg.cxx.

791 {
792 if (not m_addCounts) return true;
793 double eta = -std::log(std::tan(0.5 * track.theta()));
794 auto [enoughMeasurementsPS, tooManyHolesPS, tooManyOutliersPS] = selectCounts(track, eta);
795 return enoughMeasurementsPS && !tooManyHolesPS && !tooManyOutliersPS;
796 }
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 182 of file TrackFindingBaseAlg.cxx.

185 {
186
187 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = ActsTrk::detail::getMeasurementSelector(
188 ctx,
189 m_pixelCalibTool.isEnabled() ? &(*m_pixelCalibTool) : nullptr,
190 m_stripCalibTool.isEnabled() ? &(*m_stripCalibTool) : nullptr,
191 m_hgtdCalibTool.isEnabled() ? &(*m_hgtdCalibTool) : nullptr,
192 measurements.measurementRanges(),
194 m_measurementSelectorConfig.m_chi2CutOffOutlier,
196 m_edgeHoleBorderWidth.value());
197
198 measurementSelector->connect(&options.extensions.createTrackStates);
199
200 return measurementSelector;
201 }
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 305 of file TrackFindingBaseAlg.cxx.

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

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

◆ 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 331 of file TrackFindingBaseAlg.h.

331{};

◆ ATLAS_THREAD_SAFE [2/2]

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

Definition at line 330 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 123 of file TrackFindingBaseAlg.h.

123{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 104 of file TrackFindingBaseAlg.h.

104{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 122 of file TrackFindingBaseAlg.h.

122{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 137 of file TrackFindingBaseAlg.h.

137{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 112 of file TrackFindingBaseAlg.h.

112{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 114 of file TrackFindingBaseAlg.h.

114{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 113 of file TrackFindingBaseAlg.h.

113{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 110 of file TrackFindingBaseAlg.h.

110{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 138 of file TrackFindingBaseAlg.h.

138{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 99 of file TrackFindingBaseAlg.h.

99{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 100 of file TrackFindingBaseAlg.h.

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

◆ m_ctxProvider

ContextUtility ActsTrk::TrackFindingBaseAlg::m_ctxProvider {this}
protected

Utility to fetch the geometry, magnetic field and calibration context in the event.

Definition at line 94 of file TrackFindingBaseAlg.h.

94{this};

◆ m_d0Max

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

Definition at line 127 of file TrackFindingBaseAlg.h.

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

◆ m_d0Min

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

Definition at line 126 of file TrackFindingBaseAlg.h.

126{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 105 of file TrackFindingBaseAlg.h.

105{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 106 of file TrackFindingBaseAlg.h.

106{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 154 of file TrackFindingBaseAlg.h.

154{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 102 of file TrackFindingBaseAlg.h.

102{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 149 of file TrackFindingBaseAlg.h.

149{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 97 of file TrackFindingBaseAlg.h.

97{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 121 of file TrackFindingBaseAlg.h.

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

◆ m_etaMin

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

Definition at line 120 of file TrackFindingBaseAlg.h.

120{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_hgtdCalibTool

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

Definition at line 88 of file TrackFindingBaseAlg.h.

88{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 108 of file TrackFindingBaseAlg.h.

108{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 312 of file TrackFindingBaseAlg.h.

◆ m_maxChi2

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

Definition at line 135 of file TrackFindingBaseAlg.h.

135{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 144 of file TrackFindingBaseAlg.h.

144{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 147 of file TrackFindingBaseAlg.h.

147{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 132 of file TrackFindingBaseAlg.h.

132{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 133 of file TrackFindingBaseAlg.h.

133{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 142 of file TrackFindingBaseAlg.h.

142{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 145 of file TrackFindingBaseAlg.h.

145{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 96 of file TrackFindingBaseAlg.h.

96{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 134 of file TrackFindingBaseAlg.h.

134{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 143 of file TrackFindingBaseAlg.h.

143{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 146 of file TrackFindingBaseAlg.h.

146{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 141 of file TrackFindingBaseAlg.h.

141{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 131 of file TrackFindingBaseAlg.h.

131{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 139 of file TrackFindingBaseAlg.h.

139{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 140 of file TrackFindingBaseAlg.h.

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

◆ m_monTool

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

Definition at line 83 of file TrackFindingBaseAlg.h.

83{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 101 of file TrackFindingBaseAlg.h.

101{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 119 of file TrackFindingBaseAlg.h.

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

◆ m_phiMin

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

Definition at line 118 of file TrackFindingBaseAlg.h.

118{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 86 of file TrackFindingBaseAlg.h.

86{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 125 of file TrackFindingBaseAlg.h.

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

◆ m_ptMin

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

Definition at line 124 of file TrackFindingBaseAlg.h.

124{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 103 of file TrackFindingBaseAlg.h.

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

◆ m_recordMaterialStates

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_recordMaterialStates {this, "recordMaterialStates", false, "record track states on surfaces which have material but no measurement; keeping them off is faster, but they have to be on for the multi-component stepper and for anything downstream which needs the full track state sequence"}
protected

Definition at line 107 of file TrackFindingBaseAlg.h.

107{this, "recordMaterialStates", false, "record track states on surfaces which have material but no measurement; keeping them off is faster, but they have to be on for the multi-component stepper and for anything downstream which needs the full track state sequence"};

◆ 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 153 of file TrackFindingBaseAlg.h.

153{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 111 of file TrackFindingBaseAlg.h.

111{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 329 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 152 of file TrackFindingBaseAlg.h.

152{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 87 of file TrackFindingBaseAlg.h.

87{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 90 of file TrackFindingBaseAlg.h.

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

◆ m_trackFinder

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

Definition at line 79 of file TrackFindingBaseAlg.h.

◆ m_trackingGeometrySvc

ServiceHandle<ActsTrk::ITrackingGeometrySvc> ActsTrk::TrackFindingBaseAlg::m_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}
protected

Definition at line 84 of file TrackFindingBaseAlg.h.

84{this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"};

◆ m_tracksBackendHandlesHelper

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

Definition at line 91 of file TrackFindingBaseAlg.h.

91{this};

◆ m_trackStatePrinter

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

Definition at line 85 of file TrackFindingBaseAlg.h.

85{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 109 of file TrackFindingBaseAlg.h.

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

◆ m_useAbsEtaForStat

bool ActsTrk::TrackFindingBaseAlg::m_useAbsEtaForStat = false
protected

Definition at line 328 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 129 of file TrackFindingBaseAlg.h.

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

◆ m_z0Min

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

Definition at line 128 of file TrackFindingBaseAlg.h.

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

◆ s_branchState

BranchState ActsTrk::TrackFindingBaseAlg::s_branchState {}
staticconstexprprotected

Definition at line 293 of file TrackFindingBaseAlg.h.

293{};

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