MuonSegmentMatchingTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 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 : virtual public IMuonSegmentMatchingTool, public AthAlgTool {
  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",
    };  
 
    ToolHandle<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_dumpAngles;             
 
    bool   m_useEndcapExtrapolationMatching;
    double m_drExtrapRMS;
    double m_drExtrapAlignmentOffset;
    double m_dthetaExtrapRMS;
};
 
}  // namespace Muon
 
#endif