MuonFastReconstructionAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONR4_MUONFASTRECOALGS_MUONFASTRECONSTRUCTIONALG__H
#define MUONR4_MUONFASTRECOALGS_MUONFASTRECONSTRUCTIONALG__H
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "StoreGate/WriteHandleKey.h"
 
#include "MuonFastRecoHelpers/GlobalPatternFinder.h"
#include "MuonSpacePoint/SpacePointContainer.h"
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
 
namespace MuonR4{
    
 
    class FastReconstructionAlg: public AthReentrantAlgorithm{
        public:
            using AthReentrantAlgorithm::AthReentrantAlgorithm;
            virtual ~FastReconstructionAlg() = default;
 
            virtual StatusCode initialize() override;
            virtual StatusCode execute(const EventContext& ctx) const override;
 
        private:
            using SpacePointContainerVec = FastReco::GlobalPatternFinder::SpacePointContainerVec;
            using PatternVec = FastReco::GlobalPatternFinder::PatternVec;

            SG::ReadHandleKeyArray<SpacePointContainer> m_inSpacePoints{this, "InSpacePoints", {"MuonSpacePoints"}, "List of space point containers to read"};
            SG::WriteHandleKeyArray<SpacePointContainer> m_outSpacePoints{this, "OutSpacePoints", {}, "List of space point containers to write"};
            StringProperty m_outSpacePointSuffix{this, "OutSpacePointSuffix", "FastReco", "Suffix to add to input space point container names to create the output ones, when not provided"}; 
            SG::WriteHandleKey<GlobalPatternContainer> m_outPatterns{this, "OutPatternContainer", "R4MuonGlobalPatterns"};              
            SG::ReadHandleKey<ActsTrk::GeometryContext> m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};        
            ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
            ToolHandle<MuonValR4::IFastRecoVisualizationTool> m_visionTool{this, "VisualizationTool", ""};

            BooleanProperty m_useMdtHits {this, "UseMdtHits", true, "Activate the utilization of MDT hits to build patterns"};
            BooleanProperty m_seedFromMdt {this, "SeedFromMdt", false, "Activate the seeding from MDT hits"};
            BooleanProperty m_seedFromInner {this, "SeedFromInner", false, "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"};
            DoubleProperty m_baseRWindow {this, "BaseRWindow", 140, "Ç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"};
            UnsignedIntegerProperty m_minBendingTriggerHits {this, "MinBendingTriggerHits", 2, "Minimum number of trigger hits in the bending direction required to accept a pattern"};
            UnsignedIntegerProperty m_minBendingPrecisionHits {this, "MinBendingPrecisionHits", 8, "Minimum number of precision hits in the bending direction required to accept a pattern"};
            UnsignedIntegerProperty m_minPhiHits {this, "MinPhiHits", 1, "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"};
            UnsignedIntegerProperty m_maxSeedAttempts {this, "MaxSeedAttempts", 6, " 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"};
            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."};
            DoubleProperty m_precisionWeight {this, "PrecisionWeight", 0.75, "Weight of precision hits in the pattern score of a pattern"};
            DoubleProperty m_hitScoreSaturation {this, "HitScoreSaturation", 10.0, "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"};
            DoubleProperty m_phiBonusSaturation {this, "PhiBonusSaturation", 4.0, "Saturation for phi bonus — beyond this many phi hits the bonus is maxed"};
            
            std::unique_ptr<FastReco::GlobalPatternFinder> m_globPatFinder{};

            
    };
}
 
 
#endif