MuonPhiHoughTransformAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONR4__PHIHOUGHTRANSFORMALG__H
#define MUONR4__PHIHOUGHTRANSFORMALG__H
 
#include "MuonPatternEvent/StationHoughMaxContainer.h"
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "GaudiKernel/ToolHandle.h"
#include "StoreGate/ReadHandleKey.h"
#include "StoreGate/WriteHandleKey.h"
#include "StoreGate/ReadCondHandleKey.h"
#include "xAODMuonPrepData/MdtDriftCircleContainer.h"
#include <xAODMuonPrepData/RpcStripContainer.h>
#include <xAODMuonPrepData/TgcStripContainer.h>
#include <MuonSpacePoint/MuonSpacePointContainer.h>
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
#include "MuonPatternEvent/MuonHoughDefs.h"
#include "Gaudi/Property.h"
 
// muon includes
 
 
namespace MuonR4{
    
    class MuonPhiHoughTransformAlg: public AthReentrantAlgorithm{
        public:
                MuonPhiHoughTransformAlg(const std::string& name, ISvcLocator* pSvcLocator);
                virtual ~MuonPhiHoughTransformAlg() = default;
                virtual StatusCode initialize() override;
                virtual StatusCode execute(const EventContext& ctx) const override;
 
        private:
            
            template <class ContainerType> StatusCode retrieveContainer(const EventContext& ctx,
                                                                        const SG::ReadHandleKey<ContainerType>& key,
                                                                        const ContainerType* & contToPush) const;
 
            StatusCode prepareHoughPlane(MuonHoughEventData & data) const; 
 
            StatusCode preProcessMaximum(MuonHoughEventData & data, const MuonR4::HoughMaximum & maximum) const; 
 
            HoughSegmentSeed recoverSinglePhiMax(MuonHoughEventData & data, const MuonR4::HoughMaximum & maximum) const; 
            
            std::vector<MuonR4::ActsPeakFinderForMuon::Maximum> findRankedSegmentSeeds (MuonHoughEventData & data, const MuonR4::HoughMaximum & maximum) const; 
 
            int countIncompatibleEtaHits(const MuonR4::ActsPeakFinderForMuon::Maximum & phiMaximum, const MuonR4::HoughMaximum &  etaMaximum) const; 
            
            MuonR4::HoughSegmentSeed buildSegmentSeed(const HoughMaximum & etaMax,  const MuonR4::ActsPeakFinderForMuon::Maximum & phiMax) const; 
 
            // read handle key for the input maxima (from a previous eta-transform)
            SG::ReadHandleKey<StationHoughMaxContainer> m_maxima{this, "StationHoughMaxContainer", "MuonHoughStationMaxima"};
            // write handle key for the output segment seeds 
            SG::WriteHandleKey<StationHoughSegmentSeedContainer> m_segmentSeeds{this, "StationHoughSegmentSeedContainer", "MuonHoughStationSegmentSeeds"};
 
            // access to the ACTS geometry context 
            SG::ReadHandleKey<ActsGeometryContext> m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};
 
            // steers the target resolution in tan(phi) 
            DoubleProperty m_targetResoTanPhi{this, "ResolutionTargetTanAngle", 0.01};
            // steers the target resolution in the x-axis intercept
            DoubleProperty m_targetResoIntercept{this, "ResolutionTargetIntercept", 10.};
            // minimum search window half width, tan(phi) 
            // - in multiples of the target resolution
            DoubleProperty m_minSigmasSearchTanPhi{this, "minSigmasSearchTanPhi", 3.};
            // minimum search window half width, intercept 
            // - in multiples of the target resolution
            DoubleProperty m_minSigmasSearchIntercept{this, "minSigmasSearchIntercept", 3.};
            
            // number of accumulator bins in tan(phi)
            // target resolution in the angle
            IntegerProperty m_nBinsTanPhi{this, "nBinsTanAngle", 15};
            // number of accumulator bins in the x-axis intercept
            IntegerProperty m_nBinsIntercept{this, "nBinsIntercept", 100};
            // maximum number of eta measurements allowed to be discarded by 
            // a valid phi-extension
            IntegerProperty m_maxEtaHolesOnMax{this, "maxEtaHoles", 1};
            // flag to steer whether to recover maxima with a single phi measurement
            // using a beam spot projection. Should not be used in splashes or cosmics. 
            BooleanProperty m_recoverSinglePhiWithBS{this, "recoverSinglePhiHitsWithBS", true}; 
 
    };
}
 
 
#endif