EtaHoughTransformAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONR4_MUONPATTERNRECOGNITIONALGS_ETAHOUGHTRANSFORMALG__H
#define MUONR4_MUONPATTERNRECOGNITIONALGS_ETAHOUGHTRANSFORMALG__H
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "MuonRecToolInterfacesR4/IPatternVisualizationTool.h"
 
#include "StoreGate/WriteHandleKey.h"
 
#include <MuonPatternEvent/MuonPatternContainer.h>
#include <MuonPatternEvent/HoughEventData.h>
#include <MuonSpacePoint/SpacePointContainer.h>
 
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
 
// muon includes
 
 
namespace MuonR4{
    
    class EtaHoughTransformAlg: public AthReentrantAlgorithm{
        public:
            using AthReentrantAlgorithm::AthReentrantAlgorithm;    
 
            virtual ~EtaHoughTransformAlg() = default;
            virtual StatusCode initialize() override;
            virtual StatusCode execute(const EventContext& ctx) const override;
 
        private:
            
            using HoughSetupForBucket = HoughEventData::HoughSetupForBucket;
            void preProcess(const EventContext& ctx,
                            const ActsTrk::GeometryContext& gctx,
                            const SpacePointContainer & spacePoints,
                            HoughEventData & data) const; 

            void prepareHoughPlane(HoughEventData & data) const; 
            
            void processBucket(const EventContext& ctx,
                               HoughEventData & data, 
                               HoughSetupForBucket& currentBucket) const; 

            void fillFromSpacePoint(HoughEventData & data, 
                                    const MuonR4::HoughHitType & SP) const; 

            bool passSeedQuality (const HoughSetupForBucket& currentBucket, const MuonR4::ActsPeakFinderForMuon::Maximum & maximum) const; 
            
            void extendWithPhiHits(std::vector<HoughHitType> & hitList, 
                                   HoughSetupForBucket& bucket,
                                   const double tanBeta,
                                   const double interceptY) const ;  
            static bool isPrecisionHit(const HoughHitType& hit);

            void visualizeBucketAccumulator(const EventContext& ctx, 
                                            HoughEventData& data, 
                                            HoughSetupForBucket& bucket, 
                                            const std::vector<ActsPeakFinderForMuon::Maximum>& maxima) const;

            void visualizeBucketSeed(const EventContext& ctx,
                                     HoughSetupForBucket& bucket,
                                     const SegmentSeed& seed,
                                     bool seedPassedSelection) const;
 
 
            // target resolution in the angle
            DoubleProperty m_targetResoTanTheta{this, "ResolutionTargetTanTheta", 0.05};
            // target resolution in the y intercept
            DoubleProperty m_targetResoIntercept{this, "ResolutionTargetIntercept", 10.};
            // minimum search window half width, tan(theta) 
            // - in multiples of the target resolution
            DoubleProperty m_minSigmasSearchTanTheta{this, "minSigmasSearchTanTheta", 2.0};
            // minimum search window half width, intercept 
            // - in multiples of the target resolution
            DoubleProperty m_minSigmasSearchIntercept{this, "minSigmasSearchIntercept", 2.0};
            // Cut on the number of weighted hits on the maximum
            DoubleProperty m_peakThreshold{this, "peakThreshold", 2.5};
            // Minimum distance in tanBeta between two maxima
            DoubleProperty m_minMaxDistTheta{this, "MaximumSeparationTheta", 0.};
            // Minimum distance in the intercept between two maxima
            DoubleProperty m_minMaxDistIntercept{this, "MaximumSeparationIntercept", 15.};
            // Fraction of weighted counts around the peak to be associated to the maximum
            DoubleProperty m_peakFractionCutOff{this, "PeakFractionCutOff", 0.6};
            // How many valid precision hits have to be on the pattern
            UnsignedIntegerProperty m_nPrecHitCut{this, "nMinPrecHits", 3};
            // Phi strip safety margin. Phi hits are only appended if the maximum kind
            // of crosses the strip in theory
            DoubleProperty m_phiStripSafety{this, "PhiSafetyMargin", 15.*Gaudi::Units::cm};
            
            // number of accumulator bins for the angle 
            IntegerProperty m_nBinsTanTheta{this, "nBinsTanTheta", 7};
            // number of accumulator bins for the intercept 
            IntegerProperty m_nBinsIntercept{this, "nBinsIntercept", 15};
            // Flag to steer whether space points shall be downweighted according to their instance
            // multiplicity of the phi measurement such that it effectively contributes with weight 1
            BooleanProperty m_downWeightMultiplePrd{this, "downWeightPrdMultiplicity", false};
            ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
 
            // input space points from SG
            SG::ReadHandleKey<SpacePointContainer> m_spacePointKey{this, "SpacePointContainer", "MuonSpacePoints"};
            // output maxima for downstram processing
            SG::WriteHandleKey<EtaHoughMaxContainer> m_maxima{this, "EtaHoughMaxContainer", "MuonHoughStationMaxima"};
            // ACTS geometry context
            ActsTrk::GeoContextReadKey_t m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};
            ToolHandle<MuonValR4::IPatternVisualizationTool> m_visionTool{this, "VisualizationTool", ""};
 
    };
}
 
 
#endif