ATLAS Offline Software
Loading...
Searching...
No Matches
MuonR4::MsTrackFindingAlg Class Reference

#include <MsTrackFindingAlg.h>

Inheritance diagram for MuonR4::MsTrackFindingAlg:
Collaboration diagram for MuonR4::MsTrackFindingAlg:

Public Types

using OptBoundPars_t = Acts::Result<Acts::BoundTrackParameters>
using MeasVec_t = std::vector<const xAOD::UncalibratedMeasurement*>

Public Member Functions

virtual ~MsTrackFindingAlg ()
virtual StatusCode initialize () override final
 Standard algorithm hook to setup the extrapolator, retrieve the tools and declare algorithm's data dependencies.
virtual StatusCode execute (const EventContext &ctx) const override final
 Standard algorithm execution hook.
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
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
virtual void setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
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 Member Functions

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.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

std::unique_ptr< MsTrackSeedContainerfindTrackSeeds (const EventContext &ctx, const xAOD::MuonSegmentContainer &segments) const
 Iterates over the search tree and combines close-by segments to a track seed.
bool fitSeedCandidate (const Acts::GeometryContext &gCtx, const Acts::MagneticFieldContext &mCtx, const Acts::CalibrationContext &cCtx, const MsTrackSeed &seed, ActsTrk::MutableTrackContainer &outContainer) const
 Attempts to fit the track seed candidate to a full track and returns whether the fit succeeded.
std::pair< OptBoundPars_t, MeasVec_tprepareFit (const Acts::GeometryContext &tgContext, const Acts::MagneticFieldContext &mfContext, const Acts::CalibrationContext &calContext, const MsTrackSeed &seed) const
 Prepares the input by the fit by collecting the measurements on the segment &.
const MsTrackSeederbaselineSeeder () const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadHandleKey< xAOD::MuonSegmentContainerm_segmentKey {this, "SegmentContainer", "MuonSegmentsFromR4" }
 Declare the data dependency on the standard Mdt+Rpc+Tgc segment container & on the NSW segment container.
ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
 IdHelperSvc to decode the Identifiers.
const MuonGMR4::MuonDetectorManagerm_detMgr {nullptr}
 Pointer to the MuonDetectorManager.
SG::WriteHandleKey< MsTrackSeedContainerm_msTrkSeedKey {this, "MsTrkSeedKey", "MsTrackSeeds"}
 Temporary container write handle to push the seeds to store gate for later efficiency analysis.
ToolHandle< ISegmentSelectionToolm_segSelector {this, "SegmentSelectionTool" , "" }
 Segment selection tool to pick the good quality segments.
ToolHandle< ActsTrk::IFitterToolm_trackFitTool {this, "FittingTool", ""}
 Track fitting tool.
ToolHandle< ISpacePointCalibratorm_calibTool {this, "Calibrator", ""}
 Calibration tool to fill the track states.
PublicToolHandle< ActsTrk::ITrackingGeometryToolm_trackingGeometryTool {this, "TrackingGeometryTool", ""}
 Tracking geometry tool.
ToolHandle< ActsTrk::IExtrapolationToolm_extrapolationTool {this, "ExtrapolationTool" ,"" }
 Track extrapolation tool.
ToolHandle< MuonValR4::ITrackVisualizationToolm_visualizationTool {this, "VisualizationTool", ""}
 Visualization tool to debug the track finding.
Gaudi::Property< double > m_seedHalfLength {this, "SeedHalfLength", 50.*Gaudi::Units::cm}
 Maximum search window to search segments for.
SG::WriteHandleKey< ActsTrk::TrackContainerm_writeKey {this, "TrackWriteKey", "MsTracks"}
 Key to the output track container.
Gaudi::Property< bool > m_useMlSeeder {this, "UseMlSeeder", false, "Use segment-edge ML candidate ids to split seeding"}
 Use ML-guided segment grouping before baseline seeding.
Gaudi::Property< std::string > m_mlCandidateDecoration {this, "MlCandidateDecoration", "trackCandidateIds", "Segment vector<unsigned> decoration with ML track-candidate ids"}
 Segment decoration containing vector<unsigned> candidate IDs.
Gaudi::Property< unsigned > m_mlMinSegmentsPerCandidate {this, "MlMinSegmentsPerCandidate", 2}
Gaudi::Property< bool > m_mlFallbackToBaselineIfUndecorated
Gaudi::Property< bool > m_mlFallbackToBaselineIfNoCandidates
Gaudi::Property< bool > m_mlRunCandidatesInParallel
SG::ReadDecorHandleKey< xAOD::MuonSegmentContainerm_mlCandidateDecorKey
std::unique_ptr< MsTrackSeederm_seeder {}
 Pointer to the actual seeder implementation.
std::unique_ptr< MlMsTrackSeederm_mlSeeder {}
 Optional ML-wrapper seeder implementation.
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 39 of file MsTrackFindingAlg.h.

Member Typedef Documentation

◆ MeasVec_t

◆ OptBoundPars_t

using MuonR4::MsTrackFindingAlg::OptBoundPars_t = Acts::Result<Acts::BoundTrackParameters>

Definition at line 50 of file MsTrackFindingAlg.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ~MsTrackFindingAlg()

MuonR4::MsTrackFindingAlg::~MsTrackFindingAlg ( )
virtualdefault

Member Function Documentation

◆ baselineSeeder()

const MsTrackSeeder & MuonR4::MsTrackFindingAlg::baselineSeeder ( ) const
private

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

62{
63 return 0;
64}

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

◆ 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 MuonR4::MsTrackFindingAlg::execute ( const EventContext & ctx) const
finaloverridevirtual

Standard algorithm execution hook.

Attach the number of the parent seed to the output container

Definition at line 76 of file MsTrackFindingAlg.cxx.

76 {
77 ATH_MSG_VERBOSE("Run track finding in event "<<ctx.eventID().event_number());
78
79 const xAOD::MuonSegmentContainer* allEventSegs{nullptr};
80 ATH_CHECK(SG::get(allEventSegs, m_segmentKey, ctx));
81
82 auto seedContainer = findTrackSeeds(ctx, *allEventSegs);
83
84 const Acts::GeometryContext tgContext = m_trackingGeometryTool->getGeometryContext(ctx).context();
85 const Acts::MagneticFieldContext mfContext = m_extrapolationTool->getMagneticFieldContext(ctx);
86 const Acts::CalibrationContext calContext{ActsTrk::getCalibrationContext(ctx)};
87
88
89 Acts::VectorTrackContainer trackBackend{};
90 Acts::VectorMultiTrajectory trackStateBackend{};
91 ActsTrk::MutableTrackContainer cacheTrkContainer{std::move(trackBackend),
92 std::move(trackStateBackend)};
94 cacheTrkContainer.addColumn<std::size_t>("parentSeed");
95 unsigned seedIdx{0};
96 for (const MsTrackSeed& seed : *seedContainer) {
97 if (!fitSeedCandidate(tgContext, mfContext, calContext, seed,
98 cacheTrkContainer)) {
99 ++seedIdx;
100 continue;
101 }
102 auto lastTrack = cacheTrkContainer.getTrack(cacheTrkContainer.size() -1);
103 lastTrack.component<std::size_t, Acts::hashString("parentSeed")>() = seedIdx;
104 ++seedIdx;
105 }
106 SG::WriteHandle writeHandleSeed{m_msTrkSeedKey, ctx};
107 ATH_CHECK(writeHandleSeed.record(std::move(seedContainer)));
108
109 // Constant declination
110 Acts::ConstVectorTrackContainer ctrackBackend{std::move(cacheTrkContainer.container())};
111 Acts::ConstVectorMultiTrajectory ctrackStateBackend{std::move(cacheTrkContainer.trackStateContainer())};
112 auto ctc = std::make_unique<ActsTrk::TrackContainer>(std::move(ctrackBackend),
113 std::move(ctrackStateBackend));
114
115 SG::WriteHandle writeHandle{m_writeKey, ctx};
116 ATH_CHECK(writeHandle.record(std::move(ctc)));
117 return StatusCode::SUCCESS;
118 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_VERBOSE(x)
std::unique_ptr< MsTrackSeedContainer > findTrackSeeds(const EventContext &ctx, const xAOD::MuonSegmentContainer &segments) const
Iterates over the search tree and combines close-by segments to a track seed.
SG::WriteHandleKey< ActsTrk::TrackContainer > m_writeKey
Key to the output track container.
SG::WriteHandleKey< MsTrackSeedContainer > m_msTrkSeedKey
Temporary container write handle to push the seeds to store gate for later efficiency analysis.
PublicToolHandle< ActsTrk::ITrackingGeometryTool > m_trackingGeometryTool
Tracking geometry tool.
ToolHandle< ActsTrk::IExtrapolationTool > m_extrapolationTool
Track extrapolation tool.
SG::ReadHandleKey< xAOD::MuonSegmentContainer > m_segmentKey
Declare the data dependency on the standard Mdt+Rpc+Tgc segment container & on the NSW segment contai...
bool fitSeedCandidate(const Acts::GeometryContext &gCtx, const Acts::MagneticFieldContext &mCtx, const Acts::CalibrationContext &cCtx, const MsTrackSeed &seed, ActsTrk::MutableTrackContainer &outContainer) const
Attempts to fit the track seed candidate to a full track and returns whether the fit succeeded.
StatusCode record(std::unique_ptr< T > data)
Record a const object to the store.
Acts::TrackContainer< MutableTrackBackend, MutableTrackStateBackend, Acts::detail::ValueHolder > MutableTrackContainer
Acts::CalibrationContext getCalibrationContext(const EventContext &ctx)
The Acts::Calibration context is piped through the Acts fitters to (re)calibrate the Acts::SourceLink...
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
MuonSegmentContainer_v1 MuonSegmentContainer
Definition of the current "MuonSegment container version".

◆ 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 & AthCommonReentrantAlgorithm< 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 94 of file AthCommonReentrantAlgorithm.cxx.

88{
89 // If we didn't find any symlinks to add, just return the collection
90 // from the base class. Otherwise, return the extended collection.
91 if (!m_extendedExtraObjects.empty()) {
93 }
95}
An algorithm that can be simultaneously executed in multiple threads.

◆ filterPassed()

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

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ findTrackSeeds()

std::unique_ptr< MsTrackSeedContainer > MuonR4::MsTrackFindingAlg::findTrackSeeds ( const EventContext & ctx,
const xAOD::MuonSegmentContainer & segments ) const
private

Iterates over the search tree and combines close-by segments to a track seed.

Seeds with the same segments as other seeds are deduplicated

Parameters
ctxThe event's context to access StoreGate & Conditions
segmentsFull segment container

Definition at line 121 of file MsTrackFindingAlg.cxx.

122 {
123
124 auto seedContainer = m_seeder->findTrackSeeds(ctx, m_trackingGeometryTool->getGeometryContext(ctx), segments);
125
126 if (!m_visualizationTool.empty()) {
127 m_visualizationTool->displaySeeds(ctx, *m_seeder, segments, *seedContainer);
128 }
129 return seedContainer;
130 }
std::unique_ptr< MsTrackSeeder > m_seeder
Pointer to the actual seeder implementation.
ToolHandle< MuonValR4::ITrackVisualizationTool > m_visualizationTool
Visualization tool to debug the track finding.

◆ fitSeedCandidate()

bool MuonR4::MsTrackFindingAlg::fitSeedCandidate ( const Acts::GeometryContext & gCtx,
const Acts::MagneticFieldContext & mCtx,
const Acts::CalibrationContext & cCtx,
const MsTrackSeed & seed,
ActsTrk::MutableTrackContainer & outContainer ) const
private

Attempts to fit the track seed candidate to a full track and returns whether the fit succeeded.

Parameters
gCtxGeometry context to access the alignment of the surfaces
mCtxcMagnetic field context to access the field map during the fit
cCtxCalibration context to access the calibration constants from Store gate during the track state filling
seedThe seed of interest to fit
outContainerMutable track container to which the output track is written

Add the links to the segments making up this track as an extra column. Use the indices of the segment objects which can later be transformed into a full ElementLink as there is only one SegmentContainer from which the seeds are built

Definition at line 317 of file MsTrackFindingAlg.cxx.

321 {
322
323 ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Attempt to fit a new track seed \n"<<seed);
324 const EventContext& ctx{*calContext.get<const EventContext*>()};
325 const auto [initialPars, measurements] = prepareFit(tgContext, mfContext, calContext, seed);
326
327 if (!initialPars.ok()) {
328 ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" - Failed to construct valid parameters for seed \n"<<seed);
329 if (m_visualizationTool.isEnabled()) {
330 m_visualizationTool->displayTrackSeedObj(ctx, seed, initialPars, "FailedStartPars");
331 }
332 return false;
333 }
334 auto fitTraject = m_trackFitTool->fit(measurements, *initialPars,
335 tgContext, mfContext, calContext,
336 &(*initialPars).referenceSurface());
337 if (!fitTraject || fitTraject->size() == 0) {
338 ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Fit failed ");
339 if (m_visualizationTool.isEnabled()) {
340 m_visualizationTool->displayTrackSeedObj(ctx, seed, initialPars, "FailedFit");
341 }
342 return false;
343 }
344
345 auto track = fitTraject->getTrack(0);
346 const Amg::Vector3D trkP4 = ActsTrk::convertMomFromActs(track.fourMomentum()).first;
347 double pt = trkP4.perp() / 1000; //in GeV
348 if(pt < 2 ) {
349 double chi2PerDoF = track.chi2() / (std::max(track.nDoF(), 1u));
350 ATH_MSG_DEBUG(" ===cat dog: found low pt track candidate with pt "<<pt<<" GeV chi2/ndof "<< chi2PerDoF << " chi2 "<< track.chi2() << " nDOF "<< track.nDoF() <<"eta: "<<trkP4.eta());
351
352 track.container().trackStateContainer().visitBackwards(track.tipIndex(), [&](const auto& state) {
353 if(state.hasUncalibratedSourceLink()){
354 const auto* uncalib = dynamic_cast<const xAOD::MuonMeasurement*>(ActsTrk::detail::xAODUncalibMeasCalibrator::unpack(state.getUncalibratedSourceLink()));
355 if(uncalib){
356 ATH_MSG_DEBUG(" has meas: "<<m_idHelperSvc->toString(xAOD::identify(uncalib)) << " state " << state.typeFlags());
357 }
358 }
359 }
360 );
361
362 }
367 {
368 fitTraject->addColumn<std::vector<const xAOD::MuonSegment*>>("muonSegLinks");
369 fitTraject->getTrack(0).component<std::vector<const xAOD::MuonSegment*>>("muonSegLinks") = seed.segments();
370 }
371 outContainer.ensureDynamicColumns(*fitTraject);
372 auto destProxy = outContainer.getTrack(outContainer.addTrack());
373 destProxy.copyFrom(fitTraject->getTrack(0));
374 ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Good track fit...");
375 if (m_visualizationTool.isEnabled()) {
376 m_visualizationTool->displayTrackSeedObj(ctx, seed,
377 destProxy.createParametersAtReference(), "GoodFit");
378 }
379 return true;
380 }
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
ToolHandle< ActsTrk::IFitterTool > m_trackFitTool
Track fitting tool.
std::pair< OptBoundPars_t, MeasVec_t > prepareFit(const Acts::GeometryContext &tgContext, const Acts::MagneticFieldContext &mfContext, const Acts::CalibrationContext &calContext, const MsTrackSeed &seed) const
Prepares the input by the fit by collecting the measurements on the segment &.
std::pair< Amg::Vector3D, double > convertMomFromActs(const Acts::Vector4 &actsMom)
Converts an Acts four-momentum vector into an pair of an Athena three-momentum and the paritcle's ene...
Eigen::Matrix< double, 3, 1 > Vector3D

◆ initialize()

StatusCode MuonR4::MsTrackFindingAlg::initialize ( )
finaloverridevirtual

Standard algorithm hook to setup the extrapolator, retrieve the tools and declare algorithm's data dependencies.

Definition at line 45 of file MsTrackFindingAlg.cxx.

45 {
46 ATH_CHECK(m_segmentKey.initialize());
47 ATH_CHECK(m_idHelperSvc.retrieve());
48 ATH_CHECK(detStore()->retrieve(m_detMgr));
49 ATH_CHECK(m_segSelector.retrieve());
50 ATH_CHECK(m_msTrkSeedKey.initialize());
51
52 ATH_CHECK(m_visualizationTool.retrieve(EnableTool{!m_visualizationTool.empty()}));
55 ATH_CHECK(m_trackFitTool.retrieve());
56 ATH_CHECK(m_calibTool.retrieve());
57 ATH_CHECK(m_writeKey.initialize());
58
59 if (m_trackingGeometryTool->trackingGeometry()->geometryVersion() !=
60 Acts::TrackingGeometry::GeometryVersion::Gen3){
61 ATH_MSG_ERROR("The MS track fit requires the Gen 3 geometry format");
62 return StatusCode::FAILURE;
63 }
64
65 MsTrackSeeder::Config seederCfg{};
66 seederCfg.seedHalfLength = m_seedHalfLength;
67 seederCfg.selector = m_segSelector.get();
68 seederCfg.detMgr = m_detMgr;
69
70 m_seeder = std::make_unique<MsTrackSeeder>(name(), std::move(seederCfg));
71 return StatusCode::SUCCESS;
72 }
#define ATH_MSG_ERROR(x)
const ServiceHandle< StoreGateSvc > & detStore() const
Gaudi::Property< double > m_seedHalfLength
Maximum search window to search segments for.
ToolHandle< ISpacePointCalibrator > m_calibTool
Calibration tool to fill the track states.
const MuonGMR4::MuonDetectorManager * m_detMgr
Pointer to the MuonDetectorManager.
ToolHandle< ISegmentSelectionTool > m_segSelector
Segment selection tool to pick the good quality segments.
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
IdHelperSvc to decode the Identifiers.

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

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

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

◆ prepareFit()

std::pair< MsTrackFindingAlg::OptBoundPars_t, MsTrackFindingAlg::MeasVec_t > MuonR4::MsTrackFindingAlg::prepareFit ( const Acts::GeometryContext & tgContext,
const Acts::MagneticFieldContext & mfContext,
const Acts::CalibrationContext & calContext,
const MsTrackSeed & seed ) const
private

Prepares the input by the fit by collecting the measurements on the segment &.

Parameters
gCtxGeometry context to access the alignment of the surfaces
mCtxcMagnetic field context to access the field map during the fit
cCtxCalibration context to access the calibration constants from Store gate during the track state filling
seedThe seed of interest to fit

Fetch the measurements from the segment & find the first segment which has phi hits

Ensure that the drift signs from the fit are stamped onto the Uncalibrated measurements

The middle or outer segment provide the phi information. Not so easy becasue we want to Take the y0 & precision direction from the inner segment but the phi & x0 from a straight line extrapolation onto the plane

Extrapolate the seed segment onto the inner plane. We want to take the precision intercept from the inner segment and the non-precision intercept from the extrapolated segment

the tube form the back extrapolated segment

The extrapolation needs to go backwards and stay within the surface boundaries

Definition at line 133 of file MsTrackFindingAlg.cxx.

136 {
137 const EventContext& ctx{*calContext.get<const EventContext*>()};
138 MeasVec_t measurements{};
139 measurements.reserve(100);
141 const xAOD::MuonSegment* refSeg{nullptr};
142 for (const xAOD::MuonSegment* segment : seed.segments()) {
144 m_calibTool->stampSignsOnMeasurements(*segment);
145 MeasVec_t segMeasurements = collectMeasurements(*segment, /*skipOutlier:*/ true);
146 if (msgLvl(MSG::VERBOSE)) {
147 std::stringstream sstr{};
148 for (const xAOD::UncalibratedMeasurement* m : segMeasurements) {
149 const Acts::Surface& surf{xAOD::muonSurface(m)};
150 sstr<<" *** "<<m_idHelperSvc->toString(xAOD::identify(m))
151 <<", "<<m->numDimensions()<<", "
152 <<", "<<surf.geometryId()<<" @ "<<Amg::toString(surf.localToGlobalTransform(tgContext))<<std::endl;
153 }
154 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Fetch measurements from segment: "<<Amg::toString(segment->position())
155 <<", direction: "<<Amg::toString(segment->direction()) << " eta " << segment->direction().eta() << " phi " << segment->direction().phi() <<"\n"<<sstr.str());
156 }
157 measurements.insert(measurements.end(),
158 std::make_move_iterator(segMeasurements.begin()),
159 std::make_move_iterator(segMeasurements.end()));
160
161 // Ususally we would like to take the first segment with a sufficient amount of phi hits to set the initial position and direction of the track fit. However in some cases the segment from the NSW has a missreconstructed phi direction which causes the track fit to loose all BW and OW hits in the first iteration. Therefore if the first segment is a NSW segment we first try use a non-NSW segments with enough phi hits. If we don't find any segment with enough phi hits we will use the NSW segment as reference as long as it passes the seeding quality criteria.
162 if (!refSeg && !isNswSegment(*segment) &&
163 m_segSelector->passSeedingQuality(ctx, *detailedSegment(*segment))) {
164 refSeg = segment;
165 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Set reference segment");
166 }
167 }
168 //if we did not find a reference segment let's try the NSW one before we give up on the track
169 if(!refSeg){
170 for (const xAOD::MuonSegment* segment : seed.segments()) {
171 if (isNswSegment(*segment) &&
172 m_segSelector->passSeedingQuality(ctx, *detailedSegment(*segment))) {
173 refSeg = segment;
174 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - NSW is the best what we have apparently....");
175 break;
176 }
177 }
178 }
179
180 if (!refSeg || measurements.empty()) {
181 ATH_MSG_WARNING(__func__<<"() "<<__LINE__
182 <<" - No reference segment passing seeding quality "<<
183 (refSeg != nullptr)<<" was found. #"<<measurements.size()<<" measurements. ");
184 return std::make_pair(OptBoundPars_t::failure(std::make_error_code(std::errc::invalid_argument)),
185 std::vector<const xAOD::UncalibratedMeasurement_v1*>{});
186 }
187 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - "<<measurements.size()<<" measurements");
188 Amg::Vector3D seedPos{atFirstSurface(tgContext, *refSeg)};
189 Amg::Vector3D seedDir{refSeg->direction()};
190 ATH_MSG_DEBUG(__func__<<"() "<<__LINE__<<" - Initial seed pos: "<<Amg::toString(seedPos)
191 <<", dir: "<<Amg::toString(seedDir) << " eta " << seedDir.eta() << " phi "
192 << (seedDir.phi() /Gaudi::Units::degree) );
193
194 const Acts::GeometryIdentifier volId = volumeId(xAOD::muonSurface(measurements.front()));
195 const Acts::TrackingVolume* volume = m_trackingGeometryTool->trackingGeometry()->findVolume(volId);
196 if (!volume) {
197 ATH_MSG_WARNING(__func__<<"() "<<__LINE__
198 <<" - Failed to find tracking volume for any seed measurement");
199 return std::make_pair(OptBoundPars_t::failure(std::make_error_code(std::errc::invalid_argument)),
200 std::vector<const xAOD::UncalibratedMeasurement_v1*>{});
201 }
202 if (volume->motherVolume() && volume->motherVolume()->isAlignable()) {
203 volume = volume->motherVolume();
204 }
205 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__
206 <<" - Bounding volume "<<volume->volumeName()
207 <<", trf: "<<Amg::toString(volume->localToGlobalTransform(tgContext))
208 <<", bounds: "<<volume->volumeBounds());
212 if (const xAOD::MuonSegment* frontSegment = seed.segments().front(); frontSegment != refSeg) {
213 const Amg::Vector3D frontSegPos = atFirstSurface(tgContext, *frontSegment);
214 const Acts::Surface& firstSurf = xAOD::muonSurface(firstMeasurement(*frontSegment));
215 const Amg::Transform3D toFirstTrf = firstSurf.localToGlobalTransform(tgContext).inverse();
216 const Amg::Vector3D locFrontSegPos = toFirstTrf * frontSegPos;
217 if (!volume->inside(tgContext, frontSegPos)) {
218 ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" - Segment "<<printID(*frontSegment)
219 <<" not inside mother volume: "<<volume->volumeName()<<", "
220 <<Amg::toString(volume->globalToLocalTransform(tgContext)*frontSegPos)
221 <<", bounds: "<<volume->volumeBounds()<<", "
222 <<SegmentFit::localSegmentPars(*frontSegment)
223 <<"\n"<<print(detailedSegment(*frontSegment)->measurements()));
224 }
225
229 const Acts::MultiIntersection firstIsect = firstSurf.intersect(tgContext, seedPos, seedDir,
230 Acts::BoundaryTolerance::Infinite());
231 const Amg::Vector3D locAtFirst = toFirstTrf * firstIsect.at(0).position();
232 if (firstSurf.type() == Acts::Surface::SurfaceType::Straw) {
233 const auto& bounds = static_cast<const Acts::LineBounds&>(firstSurf.bounds());
234 using enum Acts::LineBounds::BoundValues;
235 // we want the drift radius coordinate from the segment and the coordinate along
237 const Amg::Vector3D locStartPos{locFrontSegPos.x(), locFrontSegPos.y(),
238 std::clamp(locAtFirst.z(), -bounds.get(eHalfLengthZ), bounds.get(eHalfLengthZ))};
239 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - The first surface is a straw "
240 <<bounds<<", "<<Amg::toString(locAtFirst)<<" vs. "<<Amg::toString(locFrontSegPos));
241 seedPos = firstSurf.localToGlobalTransform(tgContext) * locStartPos;
242 } else if (firstSurf.type() == Acts::Surface::SurfaceType::Plane) {
243 if (isNswSegment(*frontSegment)) {
244 seedPos = frontSegPos;
245 }
246 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - The first surface is a straw "
247 <<firstSurf.bounds()<<", "<<Amg::toString(locAtFirst)<<" vs. "<<Amg::toString(locFrontSegPos));
248 }
249
250 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Updated seed position: "<<Amg::toString(seedPos));
251 }
252
253
254 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Using seed measurement "
255 << m_idHelperSvc->toString(xAOD::identify(measurements.front()))
256 << " with volume id " << volId);
257
258 auto boundSurf = MuonGMR4::bottomBoundary(*volume);
259 if (!boundSurf) {
260 ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" - Failed to find boundary surface for tracking volume");
261 return std::make_pair(OptBoundPars_t::failure(std::make_error_code(std::errc::invalid_argument)),
262 std::vector<const xAOD::UncalibratedMeasurement_v1*>{});
263 }
264 std::shared_ptr<const Acts::Surface> targetSurf{};
265
266 auto propagteToBoundary = [&](const Acts::Surface& volBoundary) -> Acts::Result<Amg::Vector3D> {
267
268 const Amg::Transform3D& trf{volBoundary.localToGlobalTransform(tgContext)};
269 auto pIsect = intersectPlane(seedPos, seedDir, trf.linear().col(2), trf.translation());
271 if (pIsect.pathLength() > Acts::s_epsilon || !pIsect.isValid()) {
272 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Intersection @"<<Amg::toString(pIsect.position())
273 <<" is forward "<<pIsect.pathLength()<<" or invalid "<<(!pIsect.isValid())
274 <<" within volume "<<volume->inside(tgContext, pIsect.position()));
275 return Acts::Result<Amg::Vector3D>::failure(std::make_error_code(std::errc::invalid_argument));
276 }
277 Acts::Result<Amg::Vector2D> locPos = volBoundary.globalToLocal(tgContext, pIsect.position(),
278 Amg::Vector3D::Zero());
279 if (!locPos.ok()){
280 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Intersection is not on surface "<<
281 Amg::toString(trf.inverse()*pIsect.position()));
282 return Acts::Result<Amg::Vector3D>::failure(std::make_error_code(std::errc::invalid_argument));
283 }
284 if (!volBoundary.insideBounds(*locPos)) {
285 ATH_MSG_VERBOSE(__func__<<"() "<<__LINE__<<" - Intersection is outside the boundaries: "<<
286 Amg::toString(*locPos)<<", bounds: "<<volBoundary.bounds());
287 return Acts::Result<Amg::Vector3D>::failure(std::make_error_code(std::errc::invalid_argument));
288 }
289 targetSurf = volBoundary.getSharedPtr();
290 return Acts::Result<Amg::Vector3D>::success(pIsect.position());
291 };
292
293
294 auto pIsect = propagteToBoundary(*boundSurf);
295 if (!pIsect.ok() && volume->isAlignable()) {
296 const Acts::VolumePlacementBase* placement = volume->volumePlacement();
297 for (std::size_t portal = 0; !pIsect.ok() && portal< placement->nPortalPlacements(); ++portal) {
298 pIsect = propagteToBoundary(placement->portalPlacement(portal)->surface());
299 }
300 }
301 if (!pIsect.ok()) {
302 ATH_MSG_WARNING(__func__<<"() "<<__LINE__<<" Cannot create valid start parameters from seed "<<seed);
303 // THROW_EXCEPTION("DIese kacke");
304 return std::make_pair(OptBoundPars_t::failure(std::make_error_code(std::errc::invalid_argument)),
305 std::vector<const xAOD::UncalibratedMeasurement_v1*>{});
306 }
307 auto fourPos = ActsTrk::convertPosToActs(*pIsect, (*pIsect).mag() / Gaudi::Units::c_light);
308 const double qOverP = 1./ ActsTrk::energyToActs(m_seeder->estimateQtimesP(*tgContext.get<const ActsTrk::GeometryContext*>(),
309 *mfContext.get<const AtlasFieldCacheCondObj*>(), seed));
310 auto initialPars = Acts::BoundTrackParameters::create(tgContext, targetSurf, fourPos,
311 seedDir, qOverP,
312 Acts::BoundMatrix::Identity(),
313 Acts::ParticleHypothesis::muon());
314 return std::make_pair(std::move(initialPars), std::move(measurements));
315
316 }
bool msgLvl(const MSG::Level lvl) const
std::vector< const xAOD::UncalibratedMeasurement * > MeasVec_t
const MeasVec & measurements() const
Returns the associated measurements.
Amg::Vector3D direction() const
Returns the direction as Amg::Vector.
constexpr double energyToActs(const double athenaE)
Converts an energy scalar from Athena to Acts units.
Acts::Vector4 convertPosToActs(const Amg::Vector3D &athenaPos, const double athenaTime=0.)
Converts a position vector & time from Athena units into Acts units.
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Eigen::Affine3d Transform3D
@ qOverP
perigee
const Acts::Surface * bottomBoundary(const Acts::TrackingVolume &volume)
Returns the boundary surface parallel to the x-y plane at negative local z.
Parameters localSegmentPars(const xAOD::MuonSegment &seg)
Returns the localSegPars decoration from a xAODMuon::Segment.
std::string printID(const xAOD::MuonSegment &seg)
Print the chamber ID of a segment, e.g.
std::vector< const xAOD::UncalibratedMeasurement * > collectMeasurements(const xAOD::MuonSegment &segment, bool skipOutlier=true)
Helper function to extract the measurements from the segment.
Acts::GeometryIdentifier volumeId(const Acts::Surface &surface)
Returns the identifier of the volume in which the surface is embedded.
const xAOD::UncalibratedMeasurement * firstMeasurement(const xAOD::MuonSegment &segment, const bool skipOutlier=true)
Retrieves the first measurement associated with the segment.
Amg::Vector3D atFirstSurface(const Acts::GeometryContext &gctx, const xAOD::MuonSegment &segment, const bool skipOutlier=true)
Expresses the segment position on the surface of the first measurement.
std::string print(const cont_t &container)
Print a space point container to string.
const Segment * detailedSegment(const xAOD::MuonSegment &seg)
Helper function to navigate from the xAOD::MuonSegment to the MuonR4::Segment.
UncalibratedMeasurement_v1 UncalibratedMeasurement
Define the version of the uncalibrated measurement class.
const Identifier & identify(const UncalibratedMeasurement *meas)
Returns the associated identifier from the muon measurement.
MuonSegment_v1 MuonSegment
Reference the current persistent version:
const Acts::Surface & muonSurface(const UncalibratedMeasurement *meas)
Returns the associated Acts surface to the measurement.

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

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

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< 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.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

75{
76 return BaseAlg::sysExecute (ctx);
77}

◆ sysInitialize()

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

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

105 {
107
108 if (sc.isFailure()) {
109 return sc;
110 }
111
112 ServiceHandle<ICondSvc> cs("CondSvc",name());
113 for (auto h : outputHandles()) {
114 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
115 // do this inside the loop so we don't create the CondSvc until needed
116 if ( cs.retrieve().isFailure() ) {
117 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
118 return StatusCode::SUCCESS;
119 }
120 if (cs->regHandle(this,*h).isFailure()) {
122 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
123 << " with CondSvc");
124 }
125 }
126 }
127 return sc;
128}
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.

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_calibTool

ToolHandle<ISpacePointCalibrator> MuonR4::MsTrackFindingAlg::m_calibTool {this, "Calibrator", ""}
private

Calibration tool to fill the track states.

Definition at line 101 of file MsTrackFindingAlg.h.

101{this, "Calibrator", ""};

◆ m_detMgr

const MuonGMR4::MuonDetectorManager* MuonR4::MsTrackFindingAlg::m_detMgr {nullptr}
private

Pointer to the MuonDetectorManager.

Definition at line 93 of file MsTrackFindingAlg.h.

93{nullptr};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_extrapolationTool

ToolHandle<ActsTrk::IExtrapolationTool> MuonR4::MsTrackFindingAlg::m_extrapolationTool {this, "ExtrapolationTool" ,"" }
private

Track extrapolation tool.

Definition at line 105 of file MsTrackFindingAlg.h.

105{this, "ExtrapolationTool" ,"" };

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonR4::MsTrackFindingAlg::m_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

IdHelperSvc to decode the Identifiers.

Definition at line 91 of file MsTrackFindingAlg.h.

91{this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_mlCandidateDecoration

Gaudi::Property<std::string> MuonR4::MsTrackFindingAlg::m_mlCandidateDecoration {this, "MlCandidateDecoration", "trackCandidateIds", "Segment vector<unsigned> decoration with ML track-candidate ids"}
private

Segment decoration containing vector<unsigned> candidate IDs.

Definition at line 115 of file MsTrackFindingAlg.h.

115{this, "MlCandidateDecoration", "trackCandidateIds", "Segment vector<unsigned> decoration with ML track-candidate ids"};

◆ m_mlCandidateDecorKey

SG::ReadDecorHandleKey<xAOD::MuonSegmentContainer> MuonR4::MsTrackFindingAlg::m_mlCandidateDecorKey
private
Initial value:
{
this, "MlCandidateDecorationKey", "", "Scheduler dependency on ML candidate decoration"}

Definition at line 126 of file MsTrackFindingAlg.h.

126 {
127 this, "MlCandidateDecorationKey", "", "Scheduler dependency on ML candidate decoration"};

◆ m_mlFallbackToBaselineIfNoCandidates

Gaudi::Property<bool> MuonR4::MsTrackFindingAlg::m_mlFallbackToBaselineIfNoCandidates
private
Initial value:
{
this, "MlFallbackToBaselineIfNoCandidates", false,
"Run baseline seeder if ML grouping produced no seed candidates"}

Definition at line 120 of file MsTrackFindingAlg.h.

120 {
121 this, "MlFallbackToBaselineIfNoCandidates", false,
122 "Run baseline seeder if ML grouping produced no seed candidates"};

◆ m_mlFallbackToBaselineIfUndecorated

Gaudi::Property<bool> MuonR4::MsTrackFindingAlg::m_mlFallbackToBaselineIfUndecorated
private
Initial value:
{
this, "MlFallbackToBaselineIfUndecorated", true,
"Run baseline seeder if the input segment container has no ML decoration"}

Definition at line 117 of file MsTrackFindingAlg.h.

117 {
118 this, "MlFallbackToBaselineIfUndecorated", true,
119 "Run baseline seeder if the input segment container has no ML decoration"};

◆ m_mlMinSegmentsPerCandidate

Gaudi::Property<unsigned> MuonR4::MsTrackFindingAlg::m_mlMinSegmentsPerCandidate {this, "MlMinSegmentsPerCandidate", 2}
private

Definition at line 116 of file MsTrackFindingAlg.h.

116{this, "MlMinSegmentsPerCandidate", 2};

◆ m_mlRunCandidatesInParallel

Gaudi::Property<bool> MuonR4::MsTrackFindingAlg::m_mlRunCandidatesInParallel
private
Initial value:
{
this, "MlRunCandidatesInParallel", true,
"Run baseline seeding independently for ML candidate groups in parallel"}

Definition at line 123 of file MsTrackFindingAlg.h.

123 {
124 this, "MlRunCandidatesInParallel", true,
125 "Run baseline seeding independently for ML candidate groups in parallel"};

◆ m_mlSeeder

std::unique_ptr<MlMsTrackSeeder> MuonR4::MsTrackFindingAlg::m_mlSeeder {}
private

Optional ML-wrapper seeder implementation.

Definition at line 131 of file MsTrackFindingAlg.h.

131{};

◆ m_msTrkSeedKey

SG::WriteHandleKey<MsTrackSeedContainer> MuonR4::MsTrackFindingAlg::m_msTrkSeedKey {this, "MsTrkSeedKey", "MsTrackSeeds"}
private

Temporary container write handle to push the seeds to store gate for later efficiency analysis.

Definition at line 95 of file MsTrackFindingAlg.h.

95{this, "MsTrkSeedKey", "MsTrackSeeds"};

◆ m_seeder

std::unique_ptr<MsTrackSeeder> MuonR4::MsTrackFindingAlg::m_seeder {}
private

Pointer to the actual seeder implementation.

Definition at line 129 of file MsTrackFindingAlg.h.

129{};

◆ m_seedHalfLength

Gaudi::Property<double> MuonR4::MsTrackFindingAlg::m_seedHalfLength {this, "SeedHalfLength", 50.*Gaudi::Units::cm}
private

Maximum search window to search segments for.

Definition at line 109 of file MsTrackFindingAlg.h.

109{this, "SeedHalfLength", 50.*Gaudi::Units::cm};

◆ m_segmentKey

SG::ReadHandleKey<xAOD::MuonSegmentContainer> MuonR4::MsTrackFindingAlg::m_segmentKey {this, "SegmentContainer", "MuonSegmentsFromR4" }
private

Declare the data dependency on the standard Mdt+Rpc+Tgc segment container & on the NSW segment container.

Definition at line 89 of file MsTrackFindingAlg.h.

89{this, "SegmentContainer", "MuonSegmentsFromR4" };

◆ m_segSelector

ToolHandle<ISegmentSelectionTool> MuonR4::MsTrackFindingAlg::m_segSelector {this, "SegmentSelectionTool" , "" }
private

Segment selection tool to pick the good quality segments.

Definition at line 97 of file MsTrackFindingAlg.h.

97{this, "SegmentSelectionTool" , "" };

◆ m_trackFitTool

ToolHandle<ActsTrk::IFitterTool> MuonR4::MsTrackFindingAlg::m_trackFitTool {this, "FittingTool", ""}
private

Track fitting tool.

Definition at line 99 of file MsTrackFindingAlg.h.

99{this, "FittingTool", ""};

◆ m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> MuonR4::MsTrackFindingAlg::m_trackingGeometryTool {this, "TrackingGeometryTool", ""}
private

Tracking geometry tool.

Definition at line 103 of file MsTrackFindingAlg.h.

103{this, "TrackingGeometryTool", ""};

◆ m_useMlSeeder

Gaudi::Property<bool> MuonR4::MsTrackFindingAlg::m_useMlSeeder {this, "UseMlSeeder", false, "Use segment-edge ML candidate ids to split seeding"}
private

Use ML-guided segment grouping before baseline seeding.

Definition at line 113 of file MsTrackFindingAlg.h.

113{this, "UseMlSeeder", false, "Use segment-edge ML candidate ids to split seeding"};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_visualizationTool

ToolHandle<MuonValR4::ITrackVisualizationTool> MuonR4::MsTrackFindingAlg::m_visualizationTool {this, "VisualizationTool", ""}
private

Visualization tool to debug the track finding.

Definition at line 107 of file MsTrackFindingAlg.h.

107{this, "VisualizationTool", ""};

◆ m_writeKey

SG::WriteHandleKey<ActsTrk::TrackContainer> MuonR4::MsTrackFindingAlg::m_writeKey {this, "TrackWriteKey", "MsTracks"}
private

Key to the output track container.

Definition at line 111 of file MsTrackFindingAlg.h.

111{this, "TrackWriteKey", "MsTracks"};

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