Algorithm executing Phase-2 fast reconstruction. More...

#include <MuonFastReconstructionAlg.h>

Inheritance diagram for MuonR4::FastReconstructionAlg:

Collaboration diagram for MuonR4::FastReconstructionAlg:

Public Member Functions
virtual	~FastReconstructionAlg ()=default
virtual StatusCode	initialize () 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 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
using	SpacePointContainerVec = FastReco::GlobalPatternFinder::SpacePointContainerVec
	Abrivation for a vector of space-point containers.
using	PatternVec = FastReco::GlobalPatternFinder::PatternVec
	Abrivation for a vector of global patterns.
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
SG::ReadHandleKeyArray< SpacePointContainer >	m_inSpacePoints {this, "InSpacePoints", {"MuonSpacePoints"}, "List of space point containers to read"}
	Keys of SpacePoint containers to read.
SG::WriteHandleKeyArray< SpacePointContainer >	m_outSpacePoints {this, "OutSpacePoints", {}, "List of space point containers to write"}
	Write handle key for the output buckets.
StringProperty	m_outSpacePointSuffix {this, "OutSpacePointSuffix", "FastReco", "Suffix to add to input space point container names to create the output ones, when not provided"}
	Suffix to add to the input space point container names to create the output container names, when not provided.
SG::WriteHandleKey< GlobalPatternContainer >	m_outPatterns {this, "OutPatternContainer", "R4MuonGlobalPatterns"}
	Write handle key for the output global patterns.
SG::ReadHandleKey< ActsTrk::GeometryContext >	m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
	Geometry context key.
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
	Handle to the MuonIdHelper service.
ToolHandle< MuonValR4::IFastRecoVisualizationTool >	m_visionTool {this, "VisualizationTool", ""}
	Handle to the visualization tool.
BooleanProperty	m_useMdtHits {this, "UseMdtHits", true, "Activate the utilization of MDT hits to build patterns"}
	--------------— Configuration options for the global pattern finder --------------—
BooleanProperty	m_seedFromMdt {this, "SeedFromMdt", false, "Activate the seeding from MDT hits"}
	Toggle the seeding from MDT hits.
BooleanProperty	m_seedFromInner {this, "SeedFromInner", false, "Activate the seeding from Inner station"}
	Activate the seeding from Inner station.
DoubleProperty	m_thetaSearchWindow {this, "ThetaWindowSearch", 0.045, "Size of the search window in theta to link hits to a pattern"}
	Size of theta window [rad] to search for compatible hits with a seed, tailored to the target pt cutoff.
DoubleProperty	m_baseRWindow {this, "BaseRWindow", 140, "Çonstant term [mm] of the acceptance window for the residual between a test hit and the pattern line"}
	Çonstant term [mm] of the acceptance window for the residual between a test hit and the pattern line.
DoubleProperty	m_phiTolerance {this, "PhiTolerance", 0.1, "Maximum allowed phi difference between two hits to be considered compatible"}
	Maximum phi difference [rad] allowed between two hits belonging to the same pattern.
UnsignedIntegerProperty	m_minBendingTriggerHits {this, "MinBendingTriggerHits", 2, "Minimum number of trigger hits in the bending direction required to accept a pattern"}
	Requirement on trigger hits in the bending direction to accept a pattern.
UnsignedIntegerProperty	m_minBendingPrecisionHits {this, "MinBendingPrecisionHits", 8, "Minimum number of precision hits in the bending direction required to accept a pattern"}
	Requirement on precision hits in the bending direction to accept a pattern.
UnsignedIntegerProperty	m_minPhiHits {this, "MinPhiHits", 1, "Minimum number of phi measurements required to accept a pattern"}
	Minimum number of phi measurements required to accept a pattern.
DoubleProperty	m_meanNormRes2Cut {this, "meanNormRes2Cut", 0.2, "Quality cut on pattern'mean squared normalized residual"}
	Quality cut on pattern'mean squared normalized residual.
UnsignedIntegerProperty	m_maxSeedAttempts {this, "MaxSeedAttempts", 6, " Maximum number of attempts to build a pattern from hits already used in existing patterns"}
	Maximum number of attempts to build a pattern from hits already used in existing patterns.
UnsignedIntegerProperty	m_maxMissedLayerHits {this, "MaxMissedLayerHits", 3, "Maximum number of missed candidate hits in different measurement layers during pattern building allowed for a pattern branch before it is discarded"}
	Maximum number of missed candidate hits in different measurement layers during pattern building allowed for a pattern branch before it is discarded.
DoubleProperty	m_minLayerSeparation {this, "MinLayerSeparation",40, "Minimum separation (in mm) between the measurement layers of two hits for being used to compute a reliable pattern line. Use the beamspot otherwise."}
	Minimum separation [mm] between the measurement layers of two hits for being used to compute a reliable pattern line.
DoubleProperty	m_precisionWeight {this, "PrecisionWeight", 0.75, "Weight of precision hits in the pattern score of a pattern"}
	Weight of precision hits in the score of a pattern, w.r.t trigger hits.
DoubleProperty	m_hitScoreSaturation {this, "HitScoreSaturation", 10.0, "Number of eta hits per station at which the pattern score saturates"}
	number of eta hits per station at which the pattern score saturates
DoubleProperty	m_residualPenalty {this, "ResidualPenalty", 2.0, "How strongly to penalize residual in the pattern score — higher = stricter quality requirement"}
	How strongly to penalize residual — higher = stricter quality requirement.
DoubleProperty	m_phiBonusSaturation {this, "PhiBonusSaturation", 4.0, "Saturation for phi bonus — beyond this many phi hits the bonus is maxed"}
	Saturation for phi bonus — beyond this many phi hits the bonus is maxed.
std::unique_ptr< FastReco::GlobalPatternFinder >	m_globPatFinder {}
	Pointer to the actual global pattern finder.
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

Algorithm executing Phase-2 fast reconstruction.

This algorithm performs the Phase-2 fast reconstruction stage before the precision tracking in the Phase-2 Event Filter. It receives as input the collection of space-point containers and produces """"to be defined""" muon objects, which will be used to perform trigger decisions and be combined with InnerDetector tracks in the Combined Fast Reconstruction step. It will optionally write the GlobalPattern objects into the event store. The algorithm builds global patterns through the standalone module GlobalPatternFinder, fits segments in each station and use them to estimate the muon momentum.

Definition at line 29 of file MuonFastReconstructionAlg.h.

Member Typedef Documentation

◆ PatternVec

using MuonR4::FastReconstructionAlg::PatternVec = FastReco::GlobalPatternFinder::PatternVec

private

Abrivation for a vector of global patterns.

Definition at line 41 of file MuonFastReconstructionAlg.h.

◆ SpacePointContainerVec

using MuonR4::FastReconstructionAlg::SpacePointContainerVec = FastReco::GlobalPatternFinder::SpacePointContainerVec

private

Abrivation for a vector of space-point containers.

Definition at line 39 of file MuonFastReconstructionAlg.h.

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ~FastReconstructionAlg()

virtual MuonR4::FastReconstructionAlg::~FastReconstructionAlg ( )

virtualdefault

Member Function Documentation

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

{
  return 0;
}

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

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

◆ declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

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

95{ return m_detStore; }

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

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

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

85{ return m_evtStore; }

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

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode MuonR4::FastReconstructionAlg::execute ( const EventContext & ctx ) const

overridevirtual

Write out the filtered spacepoint containers, i.e. those having a pattern

Write out the global patterns. Make this optional in the future

Definition at line 72 of file MuonFastReconstructionAlg.cxx.

                                                                       {
    SpacePointContainerVec inSpacePoints{};
    inSpacePoints.reserve(m_inSpacePoints.size());
    for (const auto& key: m_inSpacePoints) {
        auto& spc = inSpacePoints.emplace_back(nullptr);
        ATH_CHECK(SG::get(spc, key, ctx));
        ATH_MSG_DEBUG("Reading " << spc->size() << " SP buckets from collection: " << key);
    }
 
    const ActsTrk::GeometryContext* gctx{nullptr};
    ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
 
    using BucketPerContainer = GlobalPatternFinder::BucketPerContainer;
    BucketPerContainer outBuckets{};
    for (const SpacePointContainer* inSpc : inSpacePoints) {
        outBuckets.emplace(inSpc, BucketPerContainer::mapped_type{});
    }
    PatternVec patterns {m_globPatFinder->findPatterns(*gctx, inSpacePoints, outBuckets)};
    for (size_t idx = 0; idx < inSpacePoints.size(); ++idx) {
        const SpacePointContainer* inSpc = inSpacePoints.at(idx);
        const auto& outBucketVec = outBuckets.at(inSpc);
 
        SG::WriteHandle<SpacePointContainer> writeHandle{m_outSpacePoints.at(idx), ctx};
        ATH_CHECK(writeHandle.record(std::make_unique<SpacePointContainer>()));
 
        for (const SpacePointBucket* bucket : outBucketVec) {
            writeHandle->push_back(std::make_unique<SpacePointBucket>(*bucket));
        }
 
        ATH_MSG_DEBUG("Written " << writeHandle->size()
                        << " SP Buckets into StoreGate with key "
                        << m_outSpacePoints.at(idx));
    }
    SG::WriteHandle<GlobalPatternContainer> writeHandle{m_outPatterns, ctx};
    ATH_CHECK(writeHandle.record(std::make_unique<GlobalPatternContainer>()));
    for (const GlobalPattern& pat : patterns) {
        writeHandle->push_back(std::make_unique<GlobalPattern>(pat));
    }
    ATH_MSG_DEBUG("Written "<<writeHandle->size()<<" GlobalPatterns into StoreGate.");
 
    /* Consume the patterns to build muon candidates. WORK IN PROGRESS */
    return StatusCode::SUCCESS;
}

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

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

◆ filterPassed()

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

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

◆ initialize()

StatusCode MuonR4::FastReconstructionAlg::initialize ( )

overridevirtual

Definition at line 11 of file MuonFastReconstructionAlg.cxx.

                                             {
    if (m_outSpacePoints.size() != m_inSpacePoints.size()) {
        m_outSpacePoints.clear();
        for (const auto& key: m_inSpacePoints) {
            m_outSpacePoints.emplace_back(key.key() + m_outSpacePointSuffix);
        }
    }
    ATH_CHECK(m_inSpacePoints.initialize());
    ATH_CHECK(m_outSpacePoints.initialize());
    ATH_CHECK(m_outPatterns.initialize());  
 
    ATH_CHECK(m_geoCtxKey.initialize());
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_visionTool.retrieve(EnableTool{!m_visionTool.empty()}));
 
    GlobalPatternFinder::Config patCfg{};
    patCfg.useMdtHits = m_useMdtHits;
    patCfg.seedFromMdt = m_seedFromMdt;
    patCfg.thetaSearchWindow = m_thetaSearchWindow;
    patCfg.baseRWindow = m_baseRWindow;
    patCfg.phiTolerance = m_phiTolerance;
    patCfg.minBendingTriggerHits = m_minBendingTriggerHits;
    patCfg.minBendingPrecisionHits = m_minBendingPrecisionHits;
    patCfg.minPhiHits = m_minPhiHits;
    patCfg.meanNormRes2Cut = m_meanNormRes2Cut;
    patCfg.maxSeedAttempts = m_maxSeedAttempts;
    patCfg.maxMissedLayerHits = m_maxMissedLayerHits;
    patCfg.minLayerSeparation = m_minLayerSeparation;
    patCfg.precisionWeight = m_precisionWeight;
    patCfg.hitScoreSaturation = m_hitScoreSaturation;
    patCfg.residualPenalty = m_residualPenalty;
    patCfg.phiBonusSaturation = m_phiBonusSaturation;
 
    if (m_seedFromInner) {
        patCfg.layerSeedings.push_back(LayerIndex::Inner);
    }
    patCfg.visionTool = m_visionTool.get();
    patCfg.idHelperSvc = m_idHelperSvc.get();
    m_globPatFinder = std::make_unique<GlobalPatternFinder>(name(), std::move(patCfg));
 
    //Print Configuration
    ATH_MSG_DEBUG(" Configuration:\n"
            << " Theta search window [rad]: " << m_thetaSearchWindow << "\n"
            << " Base R window [mm]: " << m_baseRWindow << "\n"
            << " Max missed layer hits: " << m_maxMissedLayerHits << "\n"
            << " Min layer separation [mm]: " << m_minLayerSeparation << "\n"
            << " Phi tolerance [rad]: " << m_phiTolerance << "\n"
            << " Min bending trigger hits: " << m_minBendingTriggerHits << "\n"
            << " Min bending precision hits: " << m_minBendingPrecisionHits << "\n"
            << " Min phi hits: " << m_minPhiHits << "\n"
            << " Seed from inner: " << m_seedFromInner << "\n"
            << " Use MDT hits: " << m_useMdtHits << "\n"
            << " Seed from MDT: " << m_seedFromMdt << "\n"
            << " Max seed attempts: " << m_maxSeedAttempts << "\n"
            << " Hit score saturation: " << m_hitScoreSaturation << "\n"
            << " Residual penalty: " << m_residualPenalty << "\n"
            << " Phi bonus saturation: " << m_phiBonusSaturation);
 
    return StatusCode::SUCCESS;
}

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

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, 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, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

◆ msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

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

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

  {
    h.renounce();
    PBASE::renounce (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.

                                                          {
    handlesArray.renounce();
  }

◆ setFilterPassed()

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

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

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

{
  return BaseAlg::sysExecute (ctx);
}

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

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

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

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

◆ m_baseRWindow

DoubleProperty MuonR4::FastReconstructionAlg::m_baseRWindow {this, "BaseRWindow", 140, "Çonstant term [mm] of the acceptance window for the residual between a test hit and the pattern line"}

private

Çonstant term [mm] of the acceptance window for the residual between a test hit and the pattern line.

Represents the minimum window size regardless of the geometry and absorbing unmodelled effects such as muon bending in the magnetic field and scattering

Definition at line 69 of file MuonFastReconstructionAlg.h.

69{this, "BaseRWindow", 140, "Çonstant term [mm] of the acceptance window for the residual between a test hit and the pattern line"};

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

SG::ReadHandleKey<ActsTrk::GeometryContext> MuonR4::FastReconstructionAlg::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}

private

Geometry context key.

Definition at line 52 of file MuonFastReconstructionAlg.h.

52{this, "AlignmentKey", "ActsAlignment", "cond handle key"};

◆ m_globPatFinder

std::unique_ptr<FastReco::GlobalPatternFinder> MuonR4::FastReconstructionAlg::m_globPatFinder {}

private

Pointer to the actual global pattern finder.

Definition at line 96 of file MuonFastReconstructionAlg.h.

96{};

◆ m_hitScoreSaturation

DoubleProperty MuonR4::FastReconstructionAlg::m_hitScoreSaturation {this, "HitScoreSaturation", 10.0, "Number of eta hits per station at which the pattern score saturates"}

private

number of eta hits per station at which the pattern score saturates

Definition at line 89 of file MuonFastReconstructionAlg.h.

89{this, "HitScoreSaturation", 10.0, "Number of eta hits per station at which the pattern score saturates"};

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonR4::FastReconstructionAlg::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Handle to the MuonIdHelper service.

Definition at line 54 of file MuonFastReconstructionAlg.h.

54{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ m_inSpacePoints

SG::ReadHandleKeyArray<SpacePointContainer> MuonR4::FastReconstructionAlg::m_inSpacePoints {this, "InSpacePoints", {"MuonSpacePoints"}, "List of space point containers to read"}

private

Keys of SpacePoint containers to read.

Definition at line 44 of file MuonFastReconstructionAlg.h.

44{this, "InSpacePoints", {"MuonSpacePoints"}, "List of space point containers to read"};

◆ m_maxMissedLayerHits

UnsignedIntegerProperty MuonR4::FastReconstructionAlg::m_maxMissedLayerHits {this, "MaxMissedLayerHits", 3, "Maximum number of missed candidate hits in different measurement layers during pattern building allowed for a pattern branch before it is discarded"}

private

Maximum number of missed candidate hits in different measurement layers during pattern building allowed for a pattern branch before it is discarded.

Definition at line 83 of file MuonFastReconstructionAlg.h.

83{this, "MaxMissedLayerHits", 3, "Maximum number of missed candidate hits in different measurement layers during pattern building allowed for a pattern branch before it is discarded"};

◆ m_maxSeedAttempts

UnsignedIntegerProperty MuonR4::FastReconstructionAlg::m_maxSeedAttempts {this, "MaxSeedAttempts", 6, " Maximum number of attempts to build a pattern from hits already used in existing patterns"}

private

Maximum number of attempts to build a pattern from hits already used in existing patterns.

Definition at line 81 of file MuonFastReconstructionAlg.h.

81{this, "MaxSeedAttempts", 6, " Maximum number of attempts to build a pattern from hits already used in existing patterns"};

◆ m_meanNormRes2Cut

DoubleProperty MuonR4::FastReconstructionAlg::m_meanNormRes2Cut {this, "meanNormRes2Cut", 0.2, "Quality cut on pattern'mean squared normalized residual"}

private

Quality cut on pattern'mean squared normalized residual.

Set to a large value to disable the cut, e.g. 10.

Definition at line 79 of file MuonFastReconstructionAlg.h.

79{this, "meanNormRes2Cut", 0.2, "Quality cut on pattern'mean squared normalized residual"};

◆ m_minBendingPrecisionHits

UnsignedIntegerProperty MuonR4::FastReconstructionAlg::m_minBendingPrecisionHits {this, "MinBendingPrecisionHits", 8, "Minimum number of precision hits in the bending direction required to accept a pattern"}

private

Requirement on precision hits in the bending direction to accept a pattern.

Definition at line 75 of file MuonFastReconstructionAlg.h.

75{this, "MinBendingPrecisionHits", 8, "Minimum number of precision hits in the bending direction required to accept a pattern"};

◆ m_minBendingTriggerHits

UnsignedIntegerProperty MuonR4::FastReconstructionAlg::m_minBendingTriggerHits {this, "MinBendingTriggerHits", 2, "Minimum number of trigger hits in the bending direction required to accept a pattern"}

private

Requirement on trigger hits in the bending direction to accept a pattern.

Definition at line 73 of file MuonFastReconstructionAlg.h.

73{this, "MinBendingTriggerHits", 2, "Minimum number of trigger hits in the bending direction required to accept a pattern"};

◆ m_minLayerSeparation

DoubleProperty MuonR4::FastReconstructionAlg::m_minLayerSeparation {this, "MinLayerSeparation",40, "Minimum separation (in mm) between the measurement layers of two hits for being used to compute a reliable pattern line. Use the beamspot otherwise."}

private

Minimum separation [mm] between the measurement layers of two hits for being used to compute a reliable pattern line.

Use the beamspot otherwise.

Definition at line 85 of file MuonFastReconstructionAlg.h.

85{this, "MinLayerSeparation",40, "Minimum separation (in mm) between the measurement layers of two hits for being used to compute a reliable pattern line. Use the beamspot otherwise."};

◆ m_minPhiHits

UnsignedIntegerProperty MuonR4::FastReconstructionAlg::m_minPhiHits {this, "MinPhiHits", 1, "Minimum number of phi measurements required to accept a pattern"}

private

Minimum number of phi measurements required to accept a pattern.

Definition at line 77 of file MuonFastReconstructionAlg.h.

77{this, "MinPhiHits", 1, "Minimum number of phi measurements required to accept a pattern"};

◆ m_outPatterns

SG::WriteHandleKey<GlobalPatternContainer> MuonR4::FastReconstructionAlg::m_outPatterns {this, "OutPatternContainer", "R4MuonGlobalPatterns"}

private

Write handle key for the output global patterns.

Definition at line 50 of file MuonFastReconstructionAlg.h.

50{this, "OutPatternContainer", "R4MuonGlobalPatterns"};

◆ m_outSpacePoints

SG::WriteHandleKeyArray<SpacePointContainer> MuonR4::FastReconstructionAlg::m_outSpacePoints {this, "OutSpacePoints", {}, "List of space point containers to write"}

private

Write handle key for the output buckets.

Definition at line 46 of file MuonFastReconstructionAlg.h.

46{this, "OutSpacePoints", {}, "List of space point containers to write"};

◆ m_outSpacePointSuffix

StringProperty MuonR4::FastReconstructionAlg::m_outSpacePointSuffix {this, "OutSpacePointSuffix", "FastReco", "Suffix to add to input space point container names to create the output ones, when not provided"}

private

Suffix to add to the input space point container names to create the output container names, when not provided.

Definition at line 48 of file MuonFastReconstructionAlg.h.

48{this, "OutSpacePointSuffix", "FastReco", "Suffix to add to input space point container names to create the output ones, when not provided"};

◆ m_phiBonusSaturation

DoubleProperty MuonR4::FastReconstructionAlg::m_phiBonusSaturation {this, "PhiBonusSaturation", 4.0, "Saturation for phi bonus — beyond this many phi hits the bonus is maxed"}

private

Saturation for phi bonus — beyond this many phi hits the bonus is maxed.

Definition at line 93 of file MuonFastReconstructionAlg.h.

93{this, "PhiBonusSaturation", 4.0, "Saturation for phi bonus — beyond this many phi hits the bonus is maxed"};

◆ m_phiTolerance

DoubleProperty MuonR4::FastReconstructionAlg::m_phiTolerance {this, "PhiTolerance", 0.1, "Maximum allowed phi difference between two hits to be considered compatible"}

private

Maximum phi difference [rad] allowed between two hits belonging to the same pattern.

Definition at line 71 of file MuonFastReconstructionAlg.h.

71{this, "PhiTolerance", 0.1, "Maximum allowed phi difference between two hits to be considered compatible"};

◆ m_precisionWeight

DoubleProperty MuonR4::FastReconstructionAlg::m_precisionWeight {this, "PrecisionWeight", 0.75, "Weight of precision hits in the pattern score of a pattern"}

private

Weight of precision hits in the score of a pattern, w.r.t trigger hits.

Definition at line 87 of file MuonFastReconstructionAlg.h.

87{this, "PrecisionWeight", 0.75, "Weight of precision hits in the pattern score of a pattern"};

◆ m_residualPenalty

DoubleProperty MuonR4::FastReconstructionAlg::m_residualPenalty {this, "ResidualPenalty", 2.0, "How strongly to penalize residual in the pattern score — higher = stricter quality requirement"}

private

How strongly to penalize residual — higher = stricter quality requirement.

Definition at line 91 of file MuonFastReconstructionAlg.h.

91{this, "ResidualPenalty", 2.0, "How strongly to penalize residual in the pattern score — higher = stricter quality requirement"};

◆ m_seedFromInner

BooleanProperty MuonR4::FastReconstructionAlg::m_seedFromInner {this, "SeedFromInner", false, "Activate the seeding from Inner station"}

private

Activate the seeding from Inner station.

Definition at line 64 of file MuonFastReconstructionAlg.h.

64{this, "SeedFromInner", false, "Activate the seeding from Inner station"};

◆ m_seedFromMdt

BooleanProperty MuonR4::FastReconstructionAlg::m_seedFromMdt {this, "SeedFromMdt", false, "Activate the seeding from MDT hits"}

private

Toggle the seeding from MDT hits.

Definition at line 62 of file MuonFastReconstructionAlg.h.

62{this, "SeedFromMdt", false, "Activate the seeding from MDT hits"};

◆ m_thetaSearchWindow

DoubleProperty MuonR4::FastReconstructionAlg::m_thetaSearchWindow {this, "ThetaWindowSearch", 0.045, "Size of the search window in theta to link hits to a pattern"}

private

Size of theta window [rad] to search for compatible hits with a seed, tailored to the target pt cutoff.

Definition at line 66 of file MuonFastReconstructionAlg.h.

66{this, "ThetaWindowSearch", 0.045, "Size of the search window in theta to link hits to a pattern"};

◆ m_useMdtHits

BooleanProperty MuonR4::FastReconstructionAlg::m_useMdtHits {this, "UseMdtHits", true, "Activate the utilization of MDT hits to build patterns"}

private

--------------— Configuration options for the global pattern finder --------------—

Toggle the utilization of MDT hits to build patterns

Definition at line 60 of file MuonFastReconstructionAlg.h.

60{this, "UseMdtHits", true, "Activate the utilization of MDT hits to build patterns"};

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

ToolHandle<MuonValR4::IFastRecoVisualizationTool> MuonR4::FastReconstructionAlg::m_visionTool {this, "VisualizationTool", ""}

private

Handle to the visualization tool.

Definition at line 56 of file MuonFastReconstructionAlg.h.

56{this, "VisualizationTool", ""};

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ PatternVec

◆ SpacePointContainerVec

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ ~FastReconstructionAlg()

Member Function Documentation

◆ cardinality()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ filterPassed()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_baseRWindow

◆ m_detStore

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_geoCtxKey

◆ m_globPatFinder

◆ m_hitScoreSaturation

◆ m_idHelperSvc

◆ m_inSpacePoints

◆ m_maxMissedLayerHits

◆ m_maxSeedAttempts

◆ m_meanNormRes2Cut

◆ m_minBendingPrecisionHits

◆ m_minBendingTriggerHits

◆ m_minLayerSeparation

◆ m_minPhiHits

◆ m_outPatterns

◆ m_outSpacePoints

◆ m_outSpacePointSuffix

◆ m_phiBonusSaturation

◆ m_phiTolerance

◆ m_precisionWeight

◆ m_residualPenalty

◆ m_seedFromInner

◆ m_seedFromMdt

◆ m_thetaSearchWindow

◆ m_useMdtHits

◆ m_varHandleArraysDeclared

◆ m_vhka

◆ m_visionTool