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_pimpl & trackFinder ()
const CKF_pimpl & trackFinder () const
StatusCode initializeMeasurementSelector ()
std::unique_ptr< ActsTrk::IMeasurementSelector > setMeasurementSelector (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_pimpl > m_trackFinder
ToolHandle< GenericMonitoringTool > m_monTool {this, "MonTool", "", "Monitoring tool"}
ServiceHandle< ActsTrk::ITrackingGeometrySvc > m_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}
ToolHandle< ActsTrk::TrackStatePrinterTool > m_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::TrackContainer > m_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_chi2CutOffSigma
Gaudi::Property< std::vector< double > > m_chi2OutlierCutOffSigma
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_doSecondPassBranchStopper {this, "doSecondPassBranchStopper", false, "run the branch stopper in the second pass, keeping instead of dropping"}
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_maxHolesAndOutliers {this, "maxHolesAndOutliers", {}, "TrackSelector: maximum number of holes and outliers together"}
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::uint32_t > m_envelopeConstraint
 Stop the propagation as soon as the envelope is left (Gen3 only).
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< StoreGateSvc > StoreGateSvc_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 269 of file TrackFindingBaseAlg.h.

◆ EventStats

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

Definition at line 190 of file TrackFindingBaseAlg.h.

◆ StoreGateSvc_t

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

Definition at line 376 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 245 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 167 of file TrackFindingBaseAlg.h.

168 {
172 kNoTrack,
181 kNoSecond,
187 kNStat
188 };

Constructor & Destructor Documentation

◆ TrackFindingBaseAlg()

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

Definition at line 28 of file TrackFindingBaseAlg.cxx.

28: 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 436 of file TrackFindingBaseAlg.cxx.

437 {
438 if (addCounts) {
439 tracksContainer.addColumn<unsigned int>("nPixelHits");
440 tracksContainer.addColumn<unsigned int>("nStripHits");
441 tracksContainer.addColumn<unsigned int>("nHgtdHits");
442 tracksContainer.addColumn<unsigned int>("nPixelHoles");
443 tracksContainer.addColumn<unsigned int>("nStripHoles");
444 tracksContainer.addColumn<unsigned int>("nHgtdHoles");
445 tracksContainer.addColumn<unsigned int>("nPixelOutliers");
446 tracksContainer.addColumn<unsigned int>("nStripOutliers");
447 tracksContainer.addColumn<unsigned int>("nHgtdOutliers");
448 }
450 }
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 495 of file TrackFindingBaseAlg.cxx.

495 {
496 // This check will fail if there are other types (HGTD, MS?) of hits, holes, or outliers.
497 // The check can be removed when it is no longer appropriate.
498 if (track.nMeasurements() != s_branchState.nPixelHits(track) + s_branchState.nStripHits(track) + s_branchState.nHgtdHits(track))
499 ATH_MSG_WARNING("mismatched hit count: total (" << track.nMeasurements()
500 << ") != pixel (" << s_branchState.nPixelHits(track)
501 << ") + strip (" << s_branchState.nStripHits(track)
502 << ") + hgtd (" << s_branchState.nHgtdHits(track)
503 << ")");
504 if (track.nHoles() != s_branchState.nPixelHoles(track) + s_branchState.nStripHoles(track) + s_branchState.nHgtdHoles(track))
505 ATH_MSG_WARNING("mismatched hole count: total (" << track.nHoles()
506 << ") < pixel (" << s_branchState.nPixelHoles(track)
507 << ") + strip (" << s_branchState.nStripHoles(track)
508 << ") + hgtd (" << s_branchState.nHgtdHoles(track)
509 << ")");
510 if (track.nOutliers() != s_branchState.nPixelOutliers(track) + s_branchState.nStripOutliers(track) + s_branchState.nHgtdOutliers(track))
511 ATH_MSG_WARNING("mismatched outlier count: total (" << track.nOutliers()
512 << ") != pixel (" << s_branchState.nPixelOutliers(track)
513 << ") + strip (" << s_branchState.nStripOutliers(track)
514 << ") + hgtd (" << s_branchState.nHgtdOutliers(track)
515 << ")");
516 };
#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 825 of file TrackFindingBaseAlg.cxx.

825 {
826 std::size_t out = 0u;
827 for (std::size_t category_i = seed_collection * seedCollectionStride();
828 category_i < (seed_collection + 1) * seedCollectionStride();
829 ++category_i)
830 {
831 assert(category_i < stat.size());
832 out += stat[category_i][counter_i];
833 }
834 return out;
835 }
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 611 of file TrackFindingBaseAlg.cxx.

611 {
612 std::lock_guard<std::mutex> lock(m_mutex);
613 std::size_t category_i = 0;
614 for (const std::array<unsigned int, kNStat> &src_stat : event_stat)
615 {
616 std::array<std::size_t, kNStat> &dest_stat = m_stat[category_i++];
617 for (std::size_t i = 0; i < src_stat.size(); ++i)
618 {
619 assert(i < dest_stat.size());
620 dest_stat[i] += src_stat[i];
621 }
622 }
623 }
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
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 267 of file TrackFindingBaseAlg.cxx.

270 {
271 if (not m_doTwoWay) return {};
272
273 // Create initial parameters for the propagation
274 Acts::BoundTrackParameters secondInitialParameters = trackProxy.createParametersFromState(detail::RecoConstTrackStateContainerProxy{firstMeasurement});
275 if (!secondInitialParameters.referenceSurface().insideBounds(secondInitialParameters.localPosition())) { // #3751
276 return {};
277 }
278
279 // First, inflate the covariance matrix if configured
281 ATH_MSG_DEBUG("Inflating covariance matrix for second track finding with factor = " << m_twoWayinflateCovarianceFactor.value());
282 ATH_MSG_VERBOSE("Original parameters before inflation: \n" << secondInitialParameters);
283
284 auto inflatedCovariance = secondInitialParameters.covariance().value();
285 inflatedCovariance *= m_twoWayinflateCovarianceFactor;
286
287 const auto& origSurface = secondInitialParameters.referenceSurface();
288 auto surfacePtr = const_cast<Acts::Surface&>(origSurface).shared_from_this();
289
290 Acts::BoundTrackParameters newParams(
291 std::static_pointer_cast<const Acts::Surface>(std::move(surfacePtr)),
292 secondInitialParameters.parameters(),
293 std::make_optional(inflatedCovariance),
294 secondInitialParameters.particleHypothesis());
295 secondInitialParameters = std::move(newParams);
296
297 ATH_MSG_VERBOSE("Inflated covariance matrix : \n" << secondInitialParameters.covariance().value());
298 }
299
300 auto rootBranch = tracksContainerTemp.makeTrack();
301 rootBranch.copyFromWithoutStates(trackProxy); // #3534
302
303 // perform track finding
304 auto secondResult =
305 trackFinder().ckf.findTracks(secondInitialParameters, options, tracksContainerTemp, rootBranch);
306 if (not secondResult.ok()) {
307 return {};
308 }
309 return secondResult.value();
310 }
#define ATH_MSG_DEBUG(x,...)
#define ATH_MSG_VERBOSE(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 194 of file TrackFindingBaseAlg.cxx.

194 {
195 ATH_MSG_FATAL("execute() method from the base class was called! Implement proper execute() method in the derived class!");
196
197 return StatusCode::FAILURE;
198 }
#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 200 of file TrackFindingBaseAlg.cxx.

200 {
202
203 return StatusCode::SUCCESS;
204 }

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

258 {
259 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
260 // return the last bin for |eta|>=4 or nan
261 return (!(std::abs(eta) < trackSelectorCfg.absEtaEdges.back())) ? trackSelectorCfg.cutSets.back()
262 : (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? trackSelectorCfg.cutSets.front()
263 : trackSelectorCfg.getCuts(eta);
264 };
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 227 of file TrackFindingBaseAlg.cxx.

231 {
232 Acts::PropagatorPlainOptions plainOptions{detContext.geometry, detContext.magField};
233 plainOptions.maxSteps = m_maxPropagationStep;
234 plainOptions.direction = Acts::Direction::Forward();
235 plainOptions.endOfWorldVolumeIds = m_endOfWorldVolumeIds;
236
237 // Set the CombinatorialKalmanFilter options
238 TrackFinderOptions options(detContext.geometry, detContext.magField, detContext.calib,
239 trackFinder().ckfExtensions, plainOptions);
240 options.recordMaterialStates = m_recordMaterialStates;
241
242 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = setMeasurementSelector(ctx, measurements, options);
243
244 Acts::PropagatorPlainOptions plainSecondOptions{detContext.geometry, detContext.magField};
245 plainSecondOptions.maxSteps = m_maxPropagationStep;
246 plainSecondOptions.direction = plainOptions.direction.invert();
247 plainSecondOptions.endOfWorldVolumeIds = m_endOfWorldVolumeIds;
248
249 TrackFinderOptions secondOptions(detContext.geometry, detContext.magField, detContext.calib,
250 options.extensions, plainSecondOptions);
251 secondOptions.targetSurface = pSurface;
252 secondOptions.skipPrePropagationUpdate = true;
253 secondOptions.recordMaterialStates = m_recordMaterialStates;
254
255 return {std::move(options), std::move(secondOptions), std::move(measurementSelector)};
256 };
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 816 of file TrackFindingBaseAlg.cxx.

816 {
817 std::vector<float>::const_iterator bin_iter = std::upper_bound(m_statEtaBins.begin(),
818 m_statEtaBins.end(),
819 m_useAbsEtaForStat ? std::abs(eta) : eta);
820 std::size_t category_i = seed_collection * seedCollectionStride() + static_cast<std::size_t>(bin_iter - m_statEtaBins.begin());
821 assert(category_i < m_stat.size());
822 return category_i;
823 }
Gaudi::Property< std::vector< float > > m_statEtaBins

◆ initCounts()

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

Definition at line 452 of file TrackFindingBaseAlg.cxx.

453 {
454 s_branchState.nPixelHits(track) = 0;
455 s_branchState.nStripHits(track) = 0;
456 s_branchState.nHgtdHits(track) = 0;
457 s_branchState.nPixelHoles(track) = 0;
458 s_branchState.nStripHoles(track) = 0;
459 s_branchState.nHgtdHoles(track) = 0;
460 s_branchState.nPixelOutliers(track) = 0;
461 s_branchState.nStripOutliers(track) = 0;
462 s_branchState.nHgtdOutliers(track) = 0;
463 }

◆ initialize()

StatusCode ActsTrk::TrackFindingBaseAlg::initialize ( )
overridevirtual

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

Definition at line 32 of file TrackFindingBaseAlg.cxx.

32 {
33 ATH_MSG_DEBUG("Properties Summary:");
45 ATH_MSG_DEBUG(" " << m_phiMin);
46 ATH_MSG_DEBUG(" " << m_phiMax);
47 ATH_MSG_DEBUG(" " << m_etaMin);
48 ATH_MSG_DEBUG(" " << m_etaMax);
51 ATH_MSG_DEBUG(" " << m_ptMin);
52 ATH_MSG_DEBUG(" " << m_ptMax);
53 ATH_MSG_DEBUG(" " << m_d0Min);
54 ATH_MSG_DEBUG(" " << m_d0Max);
55 ATH_MSG_DEBUG(" " << m_z0Min);
56 ATH_MSG_DEBUG(" " << m_z0Max);
75
76 m_logger = makeActsAthenaLogger(this, "Acts");
77
78 // Read and Write handles
79 ATH_CHECK(m_trackContainerKey.initialize());
81
82 ATH_CHECK(m_monTool.retrieve(EnableTool{not m_monTool.empty()}));
84 ATH_CHECK(m_ctxProvider.initialize());
85 ATH_CHECK(m_trackStatePrinter.retrieve(EnableTool{not m_trackStatePrinter.empty()}));
86 ATH_CHECK(m_pixelCalibTool.retrieve(EnableTool{not m_pixelCalibTool.empty()}));
87 ATH_CHECK(m_stripCalibTool.retrieve(EnableTool{not m_stripCalibTool.empty()}));
88 ATH_CHECK(m_hgtdCalibTool.retrieve(EnableTool{not m_hgtdCalibTool.empty()}));
89
90 auto magneticField = std::make_unique<ATLASMagneticFieldWrapper>();
91 auto trackingGeometry = m_trackingGeometrySvc->trackingGeometry();
92
93 if(trackingGeometry->geometryVersion() == Acts::TrackingGeometry::GeometryVersion::Gen3) {
94 m_endOfWorldVolumeIds.value().clear();
95 switch (m_envelopeConstraint.value()) {
96 using enum ActsTrk::SystemEnvelope;
97 case Acts::toUnderlying(ITkExit):
98 case Acts::toUnderlying(CaloExit): {
99 const Acts::TrackingVolume* volume = m_trackingGeometrySvc
100 ->getEnvelope(static_cast<SystemEnvelope>(m_envelopeConstraint.value()));
101 ATH_CHECK(volume != nullptr);
102 m_endOfWorldVolumeIds.value().push_back(volume->motherVolume()->geometryId().volume());
103 break;
104 } default: {
105 ATH_MSG_ERROR("Invalid configuration "<<m_envelopeConstraint);
106 return StatusCode::FAILURE;
107 }
108 }
109 }
110
111 detail::Stepper stepper(std::move(magneticField));
112 detail::Navigator::Config config{std::move(trackingGeometry)};
113 config.resolvePassive = false;
114 config.resolveMaterial = true;
115 config.resolveSensitive = true;
116 detail::Navigator navigator(std::move(config), logger().cloneWithSuffix("Navigator"));
117 detail::Propagator propagator(std::move(stepper), std::move(navigator), logger().cloneWithSuffix("Prop"));
118
119 // Using the CKF propagator as extrapolator
120 detail::Extrapolator extrapolator = propagator;
121
122 // m_etaBins (from flags.Tracking.ActiveConfig.etaBins) includes a dummy first and last bin, which we ignore
123 std::vector<double> absEtaEdges;
124 if (m_etaBins.size() <= 2)
125 {
126 absEtaEdges.reserve(2ul);
127 absEtaEdges.push_back(0.0);
128 absEtaEdges.push_back(std::numeric_limits<double>::infinity());
129 }
130 else
131 {
132 absEtaEdges.reserve(m_etaBins.size());
133 absEtaEdges.push_back(m_absEtaMin);
134 absEtaEdges.insert(absEtaEdges.end(), m_etaBins.begin() + 1, m_etaBins.end() - 1);
135 absEtaEdges.push_back(m_absEtaMax);
136 }
137
138 auto setCut = [](auto &cfgVal, const auto &cuts, size_t ind) -> void
139 {
140 if (cuts.empty())
141 return;
142 cfgVal = (ind < cuts.size()) ? cuts[ind] : cuts[cuts.size() - 1];
143 };
144
145 Acts::TrackSelector::EtaBinnedConfig trackSelectorCfg{std::move(absEtaEdges)};
146 if (m_etaBins.size() <= 2)
147 {
148 assert(trackSelectorCfg.cutSets.size() == 1);
149 trackSelectorCfg.cutSets[0].absEtaMin = m_absEtaMin;
150 trackSelectorCfg.cutSets[0].absEtaMax = m_absEtaMax;
151 }
152 size_t cutIndex = 0;
153 for (auto &cfg : trackSelectorCfg.cutSets)
154 {
155 setCut(cfg.phiMin, m_phiMin, cutIndex);
156 setCut(cfg.phiMax, m_phiMax, cutIndex);
157 setCut(cfg.etaMin, m_etaMin, cutIndex);
158 setCut(cfg.etaMax, m_etaMax, cutIndex);
159 setCut(cfg.ptMin, m_ptMin, cutIndex);
160 setCut(cfg.ptMax, m_ptMax, cutIndex);
161 setCut(cfg.loc0Min, m_d0Min, cutIndex);
162 setCut(cfg.loc0Max, m_d0Max, cutIndex);
163 setCut(cfg.loc1Min, m_z0Min, cutIndex);
164 setCut(cfg.loc1Max, m_z0Max, cutIndex);
165 setCut(cfg.minMeasurements, m_minMeasurements, cutIndex);
166 setCut(cfg.maxHoles, m_maxHoles, cutIndex);
167 setCut(cfg.maxOutliers, m_maxOutliers, cutIndex);
168 setCut(cfg.maxHolesAndOutliers, m_maxHolesAndOutliers, cutIndex);
169 setCut(cfg.maxSharedHits, m_maxSharedHits, cutIndex);
170 setCut(cfg.maxChi2, m_maxChi2, cutIndex);
171 ++cutIndex;
172 }
173
174 ATH_MSG_DEBUG(trackSelectorCfg);
175
176 // initializer measurement selector and connect it to the delegates of the track finder optins
177 ATH_CHECK(initializeMeasurementSelector());
178
179 detail::CKF_config ckfConfig{
180 std::move(extrapolator),
181 detail::CKF{std::move(propagator), logger().cloneWithSuffix("CKF")},
182 {},
183 Acts::TrackSelector{trackSelectorCfg}};
184
185 m_trackFinder = std::make_unique<CKF_pimpl>(std::move(ckfConfig));
186
188
189 initStatTables();
190
191 return StatusCode::SUCCESS;
192 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x,...)
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_chi2CutOffSigma
Gaudi::Property< std::vector< double > > m_ptMax
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
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
Gaudi::Property< std::vector< std::size_t > > m_maxHolesAndOutliers
Gaudi::Property< std::vector< double > > m_chi2OutlierCutOffSigma
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::uint32_t > m_envelopeConstraint
Stop the propagation as soon as the envelope is left (Gen3 only).
static Root::TMsgLogger logger("iLumiCalc")
Acts::Result< void > gainMatrixUpdate(const Acts::GeometryContext &gctx, typename trajectory_t::TrackStateProxy trackState, const Acts::Logger &logger)
std::string prefixFromTrackContainerName(const std::string &tracks)
Parse TrackContainer name to get the prefix for backends The name has to contain XYZTracks,...
SystemEnvelope
Define an enumeration to retrieve the envelope tracking volume from.
m_data push_back(elt)

◆ initializeMeasurementSelector()

StatusCode ActsTrk::TrackFindingBaseAlg::initializeMeasurementSelector ( )
protected

Definition at line 544 of file TrackFindingBaseAlg.cxx.

544 {
545 // the cuts for each eta bin are (measurement, outlier) pairs for each measurement dimension
546 auto &chi2CutOffOutlier = m_measurementSelectorConfig.m_chi2CutOffOutlier;
547 const std::vector<double> &cutOff = m_chi2CutOffSigma.value();
548 const std::vector<double> &outlierCutOff = m_chi2OutlierCutOffSigma.value();
549 chi2CutOffOutlier.reserve( cutOff.size() );
550 if (!outlierCutOff.empty()) {
551 if (cutOff.size() != outlierCutOff.size()) {
552 ATH_MSG_ERROR("Outlier chi2 cut off provided but number of elements does not agree with"
553 " chi2 cut off for measurements which however is required: "
554 << cutOff.size() << " != " << outlierCutOff.size());
555 return StatusCode::FAILURE;
556 }
557 }
558 // A cut at n sigma of a 1D Gaussian (two-sided) has the p-value erfc(n/sqrt(2)). The chi2 cut for
559 // a measurement of dimension DIM is the chi2 value with this p-value for DIM degrees of freedom.
560 // For DIM=1 this is n^2.
561 auto chi2Cut = [](double cut, std::size_t dim) -> float {
562 const double pValue = std::erfc(cut / std::numbers::sqrt2);
563 return static_cast<float>(boost::math::quantile(
564 boost::math::complement(boost::math::chi_squared_distribution<double>(dim), pValue)));
565 };
566 for (std::size_t idx = 0; idx < cutOff.size(); ++idx) {
567 std::array<std::pair<float, float>, 4> cuts{};
568 for (std::size_t dim = 1; dim < cuts.size(); ++dim) {
569 cuts[dim] = std::make_pair(chi2Cut(cutOff[idx], dim),
570 idx < outlierCutOff.size()
571 ? chi2Cut(outlierCutOff[idx], dim)
572 : std::numeric_limits<float>::max());
573 }
574 cuts[0] = cuts[1];
575 chi2CutOffOutlier.push_back(cuts);
576 ATH_MSG_DEBUG("Measurement selector eta bin " << idx << ": chi2 cuts for "
577 << cutOff[idx] << " / "
578 << (idx < outlierCutOff.size() ? outlierCutOff[idx] : std::numeric_limits<double>::infinity())
579 << " sigma: 1D " << cuts[1].first << " / " << cuts[1].second
580 << ", 2D " << cuts[2].first << " / " << cuts[2].second
581 << ", 3D " << cuts[3].first << " / " << cuts[3].second);
582 }
583 if (m_etaBins.size() > 2) {
584 std::vector<float> &etaBinsf = m_measurementSelectorConfig.m_etaBins;
585 etaBinsf.assign(m_etaBins.begin() + 1, m_etaBins.end() - 1);
586 }
587
588 return /*m_measurementSelector ?*/ StatusCode::SUCCESS /*: StatusCode::FAILURE*/;
589 }
struct ActsTrk::TrackFindingBaseAlg::MeasurementSelectorConfig m_measurementSelectorConfig

◆ initStatTables()

void ActsTrk::TrackFindingBaseAlg::initStatTables ( )
protected

Definition at line 593 of file TrackFindingBaseAlg.cxx.

593 {
594 if (!m_statEtaBins.empty())
595 {
597 float last_eta = m_statEtaBins[0];
598 for (float eta : m_statEtaBins)
599 {
600 if (eta < last_eta)
601 {
602 ATH_MSG_FATAL("Eta bins for statistics counter not in ascending order.");
603 }
604 last_eta = eta;
605 }
606 }
607 m_stat.resize(nSeedCollections() * seedCollectionStride());
608 }

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

317 {
318 return *m_logger;
319 }

◆ measurementType()

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

Definition at line 312 of file TrackFindingBaseAlg.cxx.

312 {
313 if (trackState.hasReferenceSurface()) {
314 if (const auto *actsDetElem = dynamic_cast<const ISurfacePlacement*>(trackState.referenceSurface().surfacePlacement())) {
315 switch (actsDetElem->detectorType()) {
322 default:
323 break;
324 }
325 }
326 }
328 }
@ 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 329 of file TrackFindingBaseAlg.h.

329 {
330 return m_seedLabels.size();
331 }
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 626 of file TrackFindingBaseAlg.cxx.

626 {
627 if (msgLvl(MSG::INFO))
628 {
629 std::vector<std::string> stat_labels =
631 {
632 std::make_pair(kNTotalSeeds, "Input seeds"),
633 std::make_pair(kNoTrackParam, "No track parameters"),
634 std::make_pair(kNUsedSeeds, "Used seeds"),
635 std::make_pair(kNoTrack, "Cannot find track"),
636 std::make_pair(kNDuplicateSeeds, "Duplicate seeds"),
637 std::make_pair(kNNoEstimatedParams, "Initial param estimation failed"),
638 std::make_pair(kNRejectedRefinedSeeds, "Rejected refined parameters"),
639 std::make_pair(kNSeedRefitFailure, "Seed refit Kalman fit failure"),
640 std::make_pair(kNOutputTracks, "CKF tracks"),
641 std::make_pair(kNSelectedTracks, "selected tracks"),
642 std::make_pair(kNResolvedTracks, "resolved tracks"),
643 std::make_pair(kNStoppedTracksMaxHoles, "Stopped tracks reaching max holes"),
644 std::make_pair(kMultipleBranches, "Seeds with more than one branch"),
645 std::make_pair(kNoSecond, "Tracks failing second CKF"),
646 std::make_pair(kNStoppedTracksMinPt, "Stopped tracks below pT cut"),
647 std::make_pair(kNStoppedTracksMaxEta, "Stopped tracks above max eta"),
648 std::make_pair(kNTotalSharedHits, "Total shared hits"),
649 std::make_pair(kNForcedSeedMeasurements, "Total forced measurements")
650 });
651 assert(stat_labels.size() == kNStat);
652 std::vector<std::string> categories;
653 categories.reserve(m_seedLabels.size() + 1);
654 categories.insert(categories.end(), m_seedLabels.begin(), m_seedLabels.end());
655 categories.push_back("ALL");
656
657 std::vector<std::string> eta_labels;
658 eta_labels.reserve(m_statEtaBins.size() + 2);
659 for (std::size_t eta_bin_i = 0; eta_bin_i < m_statEtaBins.size() + 2; ++eta_bin_i)
660 {
661 eta_labels.push_back(TableUtils::makeEtaBinLabel(m_statEtaBins, eta_bin_i, m_useAbsEtaForStat));
662 }
663
664 // vector used as 3D array stat[ eta_bin ][ stat_i ][ seed_type]
665 // stat_i = [0, kNStat)
666 // eta_bin = [0, m_statEtaBins.size()+2 ); eta_bin == m_statEtaBinsSize()+1 means sum of all etaBins
667 // seed_type = [0, nSeedCollections()+1) seed_type == nSeedCollections() means sum of all seed collections
668 std::vector<std::size_t> stat =
670 m_statEtaBins.size() + 1,
671 m_stat);
672
673 // the extra columns and rows for the projections are addeded internally:
674 std::size_t stat_stride =
676 m_statEtaBins.size() + 1,
677 kNStat);
678 std::size_t eta_stride =
680 m_statEtaBins.size() + 1,
681 kNStat);
682 std::stringstream table_out;
683
684 if (m_dumpAllStatEtaBins.value())
685 {
686 // dump for each counter a table with one row per eta bin
687 std::size_t max_label_width = TableUtils::maxLabelWidth(stat_labels) + TableUtils::maxLabelWidth(eta_labels);
688 for (std::size_t stat_i = 0; stat_i < kNStat; ++stat_i)
689 {
690 std::size_t dest_idx_offset = stat_i * stat_stride;
691 table_out << makeTable(stat, dest_idx_offset, eta_stride,
692 eta_labels,
693 categories)
694 .columnWidth(10)
695 // only dump the footer for the last eta bin i.e. total
696 .dumpHeader(stat_i == 0)
697 .dumpFooter(stat_i + 1 == kNStat)
698 .separateLastRow(true) // separate the sum of all eta bins
699 .minLabelWidth(max_label_width)
700 .labelPrefix(stat_labels.at(stat_i));
701 }
702 }
703 else
704 {
705 // dump one table with one row per counter showing the total eta range
706 for (std::size_t eta_bin_i = (m_dumpAllStatEtaBins.value() ? 0 : m_statEtaBins.size() + 1);
707 eta_bin_i < m_statEtaBins.size() + 2;
708 ++eta_bin_i)
709 {
710 std::size_t dest_idx_offset = eta_bin_i * eta_stride;
711 table_out << makeTable(stat, dest_idx_offset, stat_stride,
712 stat_labels,
713 categories,
714 eta_labels.at(eta_bin_i))
715 .columnWidth(10)
716 // only dump the footer for the last eta bin i.e. total
717 .dumpFooter(!m_dumpAllStatEtaBins.value() || eta_bin_i == m_statEtaBins.size() + 1);
718 }
719 }
720 ATH_MSG_INFO("statistics:\n"
721 << table_out.str());
722 table_out.str("");
723
724 // define retios first element numerator, second element denominator
725 // each element contains a vector of counter and a multiplier e.g. +- 1
726 // ratios are computed as (sum_i stat[stat_i] * multiplier_i ) / (sum_j stat[stat_j] * multiplier_j )
727 auto [ratio_labels, ratio_def] =
729 std::vector<TableUtils::SummandDefinition>{
733 // no track counted as used but want to include it as failed
735 }, // failed seeds i.e. seeds which are not duplicates but did not produce a track
736 std::vector<TableUtils::SummandDefinition>{TableUtils::defineSummand(kNTotalSeeds, 1)}),
738 TableUtils::defineSimpleRatio("Rejected refined params / seeds", kNRejectedRefinedSeeds, kNTotalSeeds),
740 TableUtils::defineSimpleRatio("selected tracks / used seeds", kNSelectedTracks, kNUsedSeeds),
742 TableUtils::defineSimpleRatio("resolved tracks / used seeds", kNResolvedTracks, kNUsedSeeds),
743 TableUtils::defineSimpleRatio("branched tracks / used seeds", kMultipleBranches, kNUsedSeeds),
744 TableUtils::defineSimpleRatio("no 2nd CKF / CKF tracks", kNoSecond, kNOutputTracks),
746 TableUtils::defineSimpleRatio("forced measurements / used seeds", kNForcedSeedMeasurements, kNUsedSeeds)});
747
748 std::vector<float> ratio = TableUtils::computeRatios(ratio_def,
749 nSeedCollections() + 1,
750 m_statEtaBins.size() + 2,
751 stat);
752
753 // the extra columns and rows for the projections are _not_ added internally
754 std::size_t ratio_stride = TableUtils::ratioStride(nSeedCollections() + 1,
755 m_statEtaBins.size() + 2,
756 ratio_def);
757 std::size_t ratio_eta_stride = TableUtils::subCategoryStride(nSeedCollections() + 1,
758 m_statEtaBins.size() + 2,
759 ratio_def);
760
761 std::size_t max_label_width = TableUtils::maxLabelWidth(ratio_labels) + TableUtils::maxLabelWidth(eta_labels);
762 if (m_dumpAllStatEtaBins.value())
763 {
764 // show for each ratio a table with one row per eta bin
765 for (std::size_t ratio_i = 0; ratio_i < ratio_labels.size(); ++ratio_i)
766 {
767 table_out << makeTable(ratio,
768 ratio_i * ratio_stride,
769 ratio_eta_stride,
770 eta_labels,
771 categories)
772 .columnWidth(10)
773 // only dump the footer for the last eta bin i.e. total
774 .dumpHeader(ratio_i == 0)
775 .dumpFooter(ratio_i + 1 == ratio_labels.size())
776 .separateLastRow(true) // separate the sum of las
777 .minLabelWidth(max_label_width)
778 .labelPrefix(ratio_labels.at(ratio_i));
779 }
780 }
781 else
782 {
783 // dump one table with one row per ratio showing the total eta range
784 table_out << makeTable(ratio,
785 (m_statEtaBins.size() + 1) * ratio_eta_stride + 0 * ratio_stride,
786 ratio_stride,
787 ratio_labels,
788 categories)
789 .columnWidth(10)
790 // only dump the footer for the last eta bin i.e. total
791 .minLabelWidth(max_label_width)
792 .dumpFooter(false);
793
794 // also dump a table for final tracks over used seeds (ratio_i==6 or 4) showing one row per eta bin
795 eta_labels.erase(eta_labels.end() - 1); // drop last line of table which shows again all eta bins summed.
796 std::size_t ratio_i = m_showResolvedStats ? 6 : 4;
797 table_out << makeTable(ratio,
798 ratio_i * ratio_stride,
799 ratio_eta_stride,
800 eta_labels,
801 categories)
802 .columnWidth(10)
803 .dumpHeader(false)
804 // only dump the footer for the last eta bin i.e. total
805 .dumpFooter(!m_dumpAllStatEtaBins.value() || ratio_i + 1 == ratio_labels.size())
806 .separateLastRow(false)
807 .minLabelWidth(max_label_width)
808 .labelPrefix(ratio_labels.at(ratio_i));
809 }
810
811 ATH_MSG_INFO("Ratios:\n"
812 << table_out.str());
813 }
814 }
#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 368 of file AthCommonDataStore.h.

369 {
370 h.renounce();
372 }
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 352 of file AthCommonDataStore.h.

352 {
354 }

◆ seedCollectionStride()

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

Definition at line 332 of file TrackFindingBaseAlg.h.

332 {
333 return m_statEtaBins.size() + 1;
334 }

◆ selectCounts()

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

Definition at line 518 of file TrackFindingBaseAlg.cxx.

518 {
519 bool enoughMeasurements = true, tooManyHoles = false, tooManyOutliers = false;
520 const auto &trackSelectorCfg = trackFinder().trackSelector.config();
521 std::size_t etaBin = (std::abs(eta) < trackSelectorCfg.absEtaEdges.front()) ? 0
522 : (std::abs(eta) >= trackSelectorCfg.absEtaEdges.back()) ? trackSelectorCfg.absEtaEdges.size() - 1
523 : trackSelectorCfg.binIndex(eta);
524 auto cutMin = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
525 return !cutSet.empty() && (val < (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
526 };
527 auto cutMax = [etaBin](std::size_t val, const std::vector<std::size_t> &cutSet) {
528 return !cutSet.empty() && (val > (etaBin < cutSet.size() ? cutSet[etaBin] : cutSet.back()));
529 };
530
531 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nPixelHits(track), m_minPixelHits);
532 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nStripHits(track), m_minStripHits);
533 enoughMeasurements = enoughMeasurements && !cutMin(s_branchState.nHgtdHits(track), m_minHgtdHits);
534 tooManyHoles = tooManyHoles || cutMax(s_branchState.nPixelHoles(track), m_maxPixelHoles);
535 tooManyHoles = tooManyHoles || cutMax(s_branchState.nStripHoles(track), m_maxStripHoles);
536 tooManyHoles = tooManyHoles || cutMax(s_branchState.nHgtdHoles(track), m_maxHgtdHoles);
537 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nPixelOutliers(track), m_maxPixelOutliers);
538 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nStripOutliers(track), m_maxStripOutliers);
539 tooManyOutliers = tooManyOutliers || cutMax(s_branchState.nHgtdOutliers(track), m_maxHgtdOutliers);
540
541 return {enoughMeasurements, tooManyHoles, tooManyOutliers};
542 }
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 837 of file TrackFindingBaseAlg.cxx.

837 {
838 if (not m_addCounts) return true;
839 double eta = -std::log(std::tan(0.5 * track.theta()));
840 auto [enoughMeasurementsPS, tooManyHolesPS, tooManyOutliersPS] = selectCounts(track, eta);
841 return enoughMeasurementsPS && !tooManyHolesPS && !tooManyOutliersPS;
842 }
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:

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

209 {
210
211 std::unique_ptr<ActsTrk::IMeasurementSelector> measurementSelector = ActsTrk::detail::getMeasurementSelector(
212 ctx,
213 m_pixelCalibTool.isEnabled() ? &(*m_pixelCalibTool) : nullptr,
214 m_stripCalibTool.isEnabled() ? &(*m_stripCalibTool) : nullptr,
215 m_hgtdCalibTool.isEnabled() ? &(*m_hgtdCalibTool) : nullptr,
216 measurements.measurementRanges(),
218 m_measurementSelectorConfig.m_chi2CutOffOutlier,
220 m_edgeHoleBorderWidth.value());
221
222 measurementSelector->connect(&options.extensions.createTrackStates);
223
224 return measurementSelector;
225 }
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::array< std::pair< float, float >, 4 > > &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 330 of file TrackFindingBaseAlg.cxx.

337 {
338 if (m_addCounts) {
339 updateCounts(track, trackState.typeFlags(),
340 measurementType(trackState));
341 if (m_checkCounts) {
342 checkCounts(track);
343 }
344 }
345
346 if (m_trackStatePrinter.isSet()) {
347 m_trackStatePrinter->printTrackState(tgContext, trackState,
348 measurementIndex, true);
349 }
350
351 if (!m_doBranchStopper) {
352 return BranchStopperResult::Continue;
353 }
354
355 const auto &parameters = trackState.hasFiltered() ? trackState.filtered()
356 : trackState.predicted();
357 double eta = -std::log(std::tan(0.5 * parameters[Acts::eBoundTheta]));
358 const auto &cutSet = getCuts(eta);
359
360 if (typeIndex < m_ptMinMeasurements.size() &&
361 !(track.nMeasurements() < m_ptMinMeasurements[typeIndex])) {
362 double pT = std::sin(parameters[Acts::eBoundTheta]) /
363 parameters[Acts::eBoundQOverP];
364 if (std::abs(pT) < cutSet.ptMin * m_branchStopperPtMinFactor) {
365 ++event_stat_category_i[kNStoppedTracksMinPt];
366 ATH_MSG_DEBUG("CkfBranchStopper: drop branch with q*pT="
367 << pT << " after " << track.nMeasurements()
368 << " measurements");
369 return BranchStopperResult::StopAndDrop;
370 }
371 }
372
373 if (typeIndex < m_absEtaMaxMeasurements.size() &&
374 !(track.nMeasurements() < m_absEtaMaxMeasurements[typeIndex]) &&
375 !(std::abs(eta) < trackSelectorCfg.absEtaEdges.back() +
377 ++event_stat_category_i[kNStoppedTracksMaxEta];
378 ATH_MSG_DEBUG("CkfBranchStopper: drop branch with eta="
379 << eta << " after " << track.nMeasurements()
380 << " measurements");
381 return BranchStopperResult::StopAndDrop;
382 }
383
384
385 // In the pixel endcap regions relax the requirement for minMeasurements before cutting the branch off
386 auto minMeasurementsBranchStop = std::abs(eta) > m_branchStopperAbsEtaMeasCut ? cutSet.minMeasurements - m_branchStopperMeasCutReduce : cutSet.minMeasurements;
387 bool enoughMeasurements = (track.nMeasurements() >= minMeasurementsBranchStop);
388 bool tooManyHoles = (track.nHoles() > cutSet.maxHoles);
389 bool tooManyOutliers = (track.nOutliers() > cutSet.maxOutliers);
390 bool tooManyHolesAndOutliers = (track.nHoles() + track.nOutliers() > cutSet.maxHolesAndOutliers);
391
392 if (m_addCounts) {
393 auto [enoughMeasurementsPS, tooManyHolesPS, tooManyOutliersPS] =
394 selectCounts(track, eta);
395 enoughMeasurements = enoughMeasurements && enoughMeasurementsPS;
396 tooManyHoles = tooManyHoles || tooManyHolesPS;
397 tooManyOutliers = tooManyOutliers || tooManyOutliersPS;
398 }
399
400 if (!(tooManyHoles || tooManyOutliers || tooManyHolesAndOutliers)) {
401 return BranchStopperResult::Continue;
402 }
403
404 if (!enoughMeasurements) {
405 ++event_stat_category_i[kNStoppedTracksMaxHoles];
406 }
407
408 if (m_addCounts) {
409 ATH_MSG_DEBUG("CkfBranchStopper: stop and "
410 << (enoughMeasurements ? "keep" : "drop")
411 << " branch with nHoles=" << track.nHoles() << " ("
412 << s_branchState.nPixelHoles(track) << " pixel+"
413 << s_branchState.nStripHoles(track) << " strip+"
414 << s_branchState.nHgtdHoles(track)
415 << " hgtd), nOutliers=" << track.nOutliers() << " ("
416 << s_branchState.nPixelOutliers(track) << "+"
417 << s_branchState.nStripOutliers(track) << "+"
418 << s_branchState.nHgtdOutliers(track)
419 << "), nMeasurements=" << track.nMeasurements() << " ("
420 << s_branchState.nPixelHits(track) << "+"
421 << s_branchState.nStripHits(track) << "+"
422 << s_branchState.nHgtdHits(track) << ")");
423 } else {
424 ATH_MSG_DEBUG("CkfBranchStopper: stop and "
425 << (enoughMeasurements ? "keep" : "drop")
426 << " branch with nHoles=" << track.nHoles()
427 << ", nOutliers=" << track.nOutliers()
428 << ", nMeasurements=" << track.nMeasurements());
429 }
430
431 return enoughMeasurements ? BranchStopperResult::StopAndKeep
432 : BranchStopperResult::StopAndDrop;
433 }
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, 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 25 of file TrackFindingBaseAlg.cxx.

25{ return *m_trackFinder; }
std::unique_ptr< CKF_pimpl > m_trackFinder

◆ trackFinder() [2/2]

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

Definition at line 26 of file TrackFindingBaseAlg.cxx.

26{ return *m_trackFinder; }

◆ updateCounts()

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

Definition at line 465 of file TrackFindingBaseAlg.cxx.

467 {
469 if (typeFlags.isHole()) {
470 s_branchState.nPixelHoles(track)++;
471 } else if (typeFlags.isOutlier()) {
472 s_branchState.nPixelOutliers(track)++;
473 } else if (typeFlags.isMeasurement()) {
474 s_branchState.nPixelHits(track)++;
475 }
476 } else if (detType == xAOD::UncalibMeasType::StripClusterType) {
477 if (typeFlags.isHole()) {
478 s_branchState.nStripHoles(track)++;
479 } else if (typeFlags.isOutlier()) {
480 s_branchState.nStripOutliers(track)++;
481 } else if (typeFlags.isMeasurement()) {
482 s_branchState.nStripHits(track)++;
483 }
484 } else if (detType == xAOD::UncalibMeasType::HGTDClusterType) {
485 if (typeFlags.isHole()) {
486 s_branchState.nHgtdHoles(track)++;
487 } else if (typeFlags.isOutlier()) {
488 s_branchState.nHgtdOutliers(track)++;
489 } else if (typeFlags.isMeasurement()) {
490 s_branchState.nHgtdHits(track)++;
491 }
492 }
493 }

◆ updateVHKA()

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

Definition at line 298 of file AthCommonDataStore.h.

298 {
299 for (auto &a : m_vhka) {
301 for (auto k : keys) {
302 k->setOwner(this);
303 }
304 }
305 }

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

341{};

◆ ATLAS_THREAD_SAFE [2/2]

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

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

129{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 107 of file TrackFindingBaseAlg.h.

107{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 128 of file TrackFindingBaseAlg.h.

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

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

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

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

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

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

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

◆ m_chi2CutOffSigma

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_chi2CutOffSigma
protected
Initial value:
{this, "chi2CutOffSigma", {},
"MeasurementSelector: maximum local chi2 contribution for each eta bin, in sigma of a 1D Gaussian. The p-value of a two-sided cut at this number of sigma gives the chi2 cut for each measurement dimension."}

Definition at line 100 of file TrackFindingBaseAlg.h.

100 {this, "chi2CutOffSigma", {},
101 "MeasurementSelector: maximum local chi2 contribution for each eta bin, in sigma of a 1D Gaussian. The p-value of a two-sided cut at this number of sigma gives the chi2 cut for each measurement dimension."};

◆ m_chi2OutlierCutOffSigma

Gaudi::Property<std::vector<double> > ActsTrk::TrackFindingBaseAlg::m_chi2OutlierCutOffSigma
protected
Initial value:
{this, "chi2OutlierCutOffSigma", {},
"MeasurementSelector: maximum local chi2 contribution for outliers for each eta bin, in sigma, see chi2CutOffSigma. If empty, there is no cut."}

Definition at line 102 of file TrackFindingBaseAlg.h.

102 {this, "chi2OutlierCutOffSigma", {},
103 "MeasurementSelector: maximum local chi2 contribution for outliers for each eta bin, in sigma, see chi2CutOffSigma. If empty, there is no cut."};

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

95{this};

◆ m_d0Max

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

Definition at line 133 of file TrackFindingBaseAlg.h.

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

◆ m_d0Min

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

Definition at line 132 of file TrackFindingBaseAlg.h.

132{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 381 of file AthCommonDataStore.h.

◆ m_doBranchStopper

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

Definition at line 108 of file TrackFindingBaseAlg.h.

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

◆ m_doSecondPassBranchStopper

Gaudi::Property<bool> ActsTrk::TrackFindingBaseAlg::m_doSecondPassBranchStopper {this, "doSecondPassBranchStopper", false, "run the branch stopper in the second pass, keeping instead of dropping"}
protected

Definition at line 111 of file TrackFindingBaseAlg.h.

111{this, "doSecondPassBranchStopper", false, "run the branch stopper in the second pass, keeping instead of dropping"};

◆ m_doTwoWay

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

Definition at line 112 of file TrackFindingBaseAlg.h.

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

◆ m_dumpAllStatEtaBins

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

Definition at line 165 of file TrackFindingBaseAlg.h.

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

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

160{this, "EndOfTheWorldVolumeIds", {}, ""};

◆ m_envelopeConstraint

Gaudi::Property<std::uint32_t> ActsTrk::TrackFindingBaseAlg::m_envelopeConstraint
protected
Initial value:
{this, "EnvelopeConstaint",
Acts::toUnderlying(ActsTrk::SystemEnvelope::ITkExit)}

Stop the propagation as soon as the envelope is left (Gen3 only).

Definition at line 157 of file TrackFindingBaseAlg.h.

157 {this, "EnvelopeConstaint",
158 Acts::toUnderlying(ActsTrk::SystemEnvelope::ITkExit)};

◆ m_etaBins

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

Definition at line 98 of file TrackFindingBaseAlg.h.

98{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 127 of file TrackFindingBaseAlg.h.

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

◆ m_etaMin

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

Definition at line 126 of file TrackFindingBaseAlg.h.

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

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

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

◆ m_maxChi2

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

Definition at line 142 of file TrackFindingBaseAlg.h.

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

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

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

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

◆ m_maxHolesAndOutliers

Gaudi::Property<std::vector<std::size_t> > ActsTrk::TrackFindingBaseAlg::m_maxHolesAndOutliers {this, "maxHolesAndOutliers", {}, "TrackSelector: maximum number of holes and outliers together"}
protected

Definition at line 140 of file TrackFindingBaseAlg.h.

140{this, "maxHolesAndOutliers", {}, "TrackSelector: maximum number of holes and outliers together"};

◆ m_maxOutliers

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

Definition at line 139 of file TrackFindingBaseAlg.h.

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

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

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

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

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

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

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

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

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

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

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

◆ m_monTool

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

Definition at line 84 of file TrackFindingBaseAlg.h.

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

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

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

◆ m_phiMin

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

Definition at line 124 of file TrackFindingBaseAlg.h.

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

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

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

◆ m_ptMin

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

Definition at line 130 of file TrackFindingBaseAlg.h.

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

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

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

164{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 117 of file TrackFindingBaseAlg.h.

117{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 339 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 163 of file TrackFindingBaseAlg.h.

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

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

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

◆ m_trackFinder

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

Definition at line 80 of file TrackFindingBaseAlg.h.

◆ m_trackingGeometrySvc

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

Definition at line 85 of file TrackFindingBaseAlg.h.

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

◆ m_tracksBackendHandlesHelper

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

Definition at line 92 of file TrackFindingBaseAlg.h.

92{this};

◆ m_trackStatePrinter

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

Definition at line 86 of file TrackFindingBaseAlg.h.

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

115{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 338 of file TrackFindingBaseAlg.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 387 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
privateinherited

Definition at line 386 of file AthCommonDataStore.h.

◆ m_z0Max

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

Definition at line 135 of file TrackFindingBaseAlg.h.

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

◆ m_z0Min

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

Definition at line 134 of file TrackFindingBaseAlg.h.

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

◆ s_branchState

BranchState ActsTrk::TrackFindingBaseAlg::s_branchState {}
staticconstexprprotected

Definition at line 303 of file TrackFindingBaseAlg.h.

303{};

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