MuonSegmentMatchingTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUON_MUONSEGMENTMATCHINGTOOL_H
#define MUON_MUONSEGMENTMATCHINGTOOL_H
 
#include <atomic>
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
#include "MuonRecHelperTools/IMuonEDMHelperSvc.h"
#include "MuonRecHelperTools/MuonEDMPrinterTool.h"
#include "MuonSegmentMakerToolInterfaces/IMuonSegmentInOverlapResolvingTool.h"
#include "MuonSegmentMakerToolInterfaces/IMuonSegmentMatchingTool.h"
#include "MuonSegmentMakerToolInterfaces/IMuonSegmentPairMatchingTool.h"
#include "TrkGeometry/MagneticFieldProperties.h"
 
namespace Muon {
 
class MuonSegment;

class MuonSegmentMatchingTool :  public extends<AthAlgTool, IMuonSegmentMatchingTool> {
  public:
    MuonSegmentMatchingTool(const std::string&, const std::string&, const IInterface*);

    virtual ~MuonSegmentMatchingTool() = default;

    StatusCode initialize();

    StatusCode finalize();

    bool match(const EventContext& ctx, const MuonSegment& seg1, const MuonSegment& seg2) const;
 
  private:
    bool straightLineMatch(const MuonSegment& seg1, const MuonSegment& seg2) const;

    bool curvedMatch(const MuonSegment& seg1, const MuonSegment& seg2) const;

    bool overlapMatch(const EventContext& ctx, const MuonSegment& seg1, const MuonSegment& seg2) const;

    bool initializeField() const;

    bool hasStereoAngle(const Identifier& id1, const Identifier& id2) const;

    bool isSLMatch(const Identifier& chid1, const Identifier& chid2) const;

    bool suppressNoise(const MuonSegment& seg1, const MuonSegment& seg2, const bool& useTightCuts) const;

    bool suppressNoisePhi(const MuonSegment& seg1, const MuonSegment& seg2, const bool& useTightCuts) const;

    void simpleEndcapExtrapolate(double x_segment, double y_segment, double z_segment, double theta_segment,
                                 double z_extrapolated, double& r_expected, double& theta_expected,
                                 double& rhoInv) const;
 

    bool endcapExtrapolationMatch(const MuonSegment& seg1, const MuonSegment& seg2, bool useTightCuts) const;
 
    ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{
        this,
        "MuonIdHelperSvc",
        "Muon::MuonIdHelperSvc/MuonIdHelperSvc",
    };
    ServiceHandle<Muon::IMuonEDMHelperSvc> m_edmHelperSvc{
        this,
        "edmHelper",
        "Muon::MuonEDMHelperSvc/MuonEDMHelperSvc",
        "Handle to the service providing the IMuonEDMHelperSvc interface",
    };  
 
    PublicToolHandle<Muon::MuonEDMPrinterTool> m_printer{
        this,
        "Printer",
        "Muon::MuonEDMPrinterTool/MuonEDMPrinterTool",
    };  
    ToolHandle<Muon::IMuonSegmentInOverlapResolvingTool> m_overlapResolvingTool{
        this,
        "MuonSegmentInOverlapResolvingTool",
        "Muon::MuonSegmentInOverlapResolvingTool/MuonSegmentInOverlapResolvingTool",
    };  
    ToolHandle<Muon::IMuonSegmentPairMatchingTool> m_pairMatchingTool{
        this,
        "MuonSegmentPairMatchingTool",
        "Muon::MuonSegmentPairMatchingTool/MuonSegmentPairMatchingTool",
    };  
 
    Gaudi::Property<bool> m_toroidOn{this, "ToroidOn", true, "Status of toroidal B-Field"};
 
    mutable std::atomic_uint m_straightLineMatches{0};
    mutable std::atomic_uint m_straightLineMatchesGood{0};
    mutable std::atomic_uint m_overlapMatches{0};
    mutable std::atomic_uint m_overlapMatchesGood{0};
    mutable std::atomic_uint m_curvedMatches{0};
    mutable std::atomic_uint m_curvedMatchesGood{0};
    mutable std::atomic_uint m_duplicateHitUses{0};
 
    bool m_isCosmics;
    bool m_doOverlapMatch;
    bool m_doStraightLineMatch;
    bool m_doCurvedMatch;
    bool m_thetaMatch;
    bool m_phiMatch;
    bool m_useLocalAngles{};
 
    // cuts for straight line match
    double m_straightLineMatchAngleCut;
    double m_straightLineMatchPositionCut;
 
    // cuts for sl overlap match
    double m_overlapMatchAngleDPhiCut;  
    double
           m_overlapMatchAngleDYZCut;  
    double m_overlapMatchPositionCut;          
    double m_overlapMatchPositionResidualCut;  
    double m_overlapMatchPhiHitPullCut;  
 
    // cuts for matching segments from different stations
    double m_angleABCut{};
    double m_maxDistSegments{};         
    double m_minDistSegmentsCosmics{};  
    double m_matchingbibm_lphisec{};
    double m_matchingbibo_lphisec{};
    double m_matchingbmbo_lphisec{};
    double m_matchingeiem_lphisec{};
    double m_matchingeieo_lphisec{};
    double m_matchingemeo_lphisec{};
    double m_matchingbibm_sphisec{};
    double m_matchingbibo_sphisec{};
    double m_matchingbmbo_sphisec{};
    double m_matchingeiem_sphisec{};
    double m_matchingeieo_sphisec{};
    double m_matchingemeo_sphisec{};
 
    double m_matchingbee_sphisec{};
 
    bool m_onlySameSectorIfTight{};  
    bool m_useTightCuts{};           
 
    bool   m_useEndcapExtrapolationMatching{};
    double m_drExtrapRMS{};
    double m_drExtrapAlignmentOffset{};
    double m_dthetaExtrapRMS{};
};
 
}  // namespace Muon
 
#endif