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 <MuonRecToolInterfacesR4/IPatternVisualizationTool.h>
#include <MuonRecToolInterfacesR4/IFastRecoVisualizationTool.h>
#include <MuonRecToolInterfacesR4/ITrackSeedingTool.h>
 
#include "MuonFastRecoHelpers/GlobalPatternFinder.h"
#include "MuonFastRecoHelpers/FastMuonSABuilder.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::WriteHandleKey<xAOD::MuonContainer> m_outMuons{this, "OutMuons", "R4FastRecoSAMuons", "List of fast muon containers to write"};              
            ActsTrk::GeoContextReadKey_t m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};        
            ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
            ToolHandle<ISpacePointCalibrator> m_calibTool{this, "Calibrator", "" };
            ToolHandle<MuonValR4::IFastRecoVisualizationTool> m_patVisionTool{this, "VisualizationTool", ""};
            ToolHandle<MuonValR4::IPatternVisualizationTool> m_segVisionTool{this, "SegmentVisualizationTool", ""};
            ToolHandle<ITrackSeedingTool> m_seedingTool{this, "SeedingTool", ""};

            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.06, "Size of the search window in theta to link hits to a pattern"};
            DoubleProperty m_baseResidualSigma  {this, "BaseResidualSigma", 35, "Effective isotropic position uncertainty [mm], including detector resolution and unmodelled effects"};
            DoubleProperty m_phiTolerance {this, "PhiTolerance", 0.05, "Maximum allowed phi difference between two hits to be considered compatible"};
            UnsignedIntegerProperty m_minTriggerLayers {this, "MinBendingTriggerLayers", 2, "Minimum number of trigger layers in the bending direction required to accept a pattern"};
            UnsignedIntegerProperty m_minPrecisionLayers {this, "MinBendingPrecisionLayers", 8, "Minimum number of precision layers in the bending direction required to accept a pattern"};
            UnsignedIntegerProperty m_minPhiLayers {this, "MinPhiLayers", 1, "Minimum number of phi layers required to accept a pattern"};
            UnsignedIntegerProperty m_minStationLayers {this, "MinStationLayers", 5, "Minimum number of layers in a station to be considered a good station"};
            DoubleProperty m_meanNormRes2Cut {this, "meanNormRes2Cut", 0.2, "Quality cut on pattern'mean squared normalized residual"};
            UnsignedIntegerProperty m_maxSeedAttempts {this, "MaxSeedAttempts", 3, " Maximum number of attempts to build a pattern from hits already used in existing patterns"};
            UnsignedIntegerProperty m_maxMissLayersInStation {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_minHitDistance4Line {this, "MinHitDistance4Line",40, "Minimum distance (in mm) between two hits for being used to compute a reliable pattern line. Use the beamspot otherwise."};
            
            BooleanProperty m_recalibInFit{this, "recalibInFit" , false};
            BooleanProperty m_hessianResidual{this, "useHessianResidual", false};
            DoubleProperty m_seedHitChi2{this, "ResoSeedHitAssoc", 5. };
            BooleanProperty m_recalibSeed{this, "SeedRecalibrate", false};
            DoubleProperty m_goodSegmentCut{this, "GoodSegmentCut", 2.};
            DoubleProperty m_outlierRemovalCut{this, "OutlierRemoval", 3.};
            DoubleProperty m_recoveryPull{this, "RecoveryPull", 3.};
            UnsignedIntegerProperty m_precHitCut{this, "PrecHitCut" , 3};
            BooleanProperty m_useFastFitter{this, "useFastFitter", false};
            BooleanProperty m_fastPreFitter{this, "useFastPreFitter", false};
            BooleanProperty m_ignoreFailedPreFit{this, "ignoreFailedPreFit", false};
            UnsignedIntegerProperty m_maxIter{this, "maxIterations", 50};
 
 

            std::unique_ptr<FastReco::GlobalPatternFinder> m_globPatFinder{};
            std::unique_ptr<FastReco::FastMuonSABuilder> m_saBuilder{};
 
            
            
    };
}
 
 
#endif