MuonSegmentInOverlapResolvingTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUON_MUONSEGMENTSOVERLAPRESOLVINGTOOL_H
#define MUON_MUONSEGMENTSOVERLAPRESOLVINGTOOL_H
 
#include <string>
#include <vector>
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
#include "MuonRecHelperTools/IMuonEDMHelperSvc.h"
#include "MuonRecHelperTools/MuonEDMPrinterTool.h"
#include "MuonSegmentMakerToolInterfaces/IMuonSegmentInOverlapResolvingTool.h"
#include "TrkExInterfaces/IPropagator.h"
#include "TrkGeometry/MagneticFieldProperties.h"
#include "TrkToolInterfaces/IResidualPullCalculator.h"
 
class MsgStream;
 
namespace Trk {
class MagneticFieldProperties;
class MeasurementBase;
}  // namespace Trk
 
namespace MuonGM {
class MdtReadoutElement;
}
 
namespace Muon {
 
class MuonSegment;
class MuonSegmentInOverlapResolvingTool : virtual public IMuonSegmentInOverlapResolvingTool, public AthAlgTool {
  public:
    MuonSegmentInOverlapResolvingTool(const std::string&, const std::string&, const IInterface*);

    virtual ~MuonSegmentInOverlapResolvingTool(){};

    StatusCode initialize();

    SegmentMatchResult matchResult(const EventContext& ctx, const MuonSegment& seg1, const MuonSegment& seg2) const;

    SegmentPositionMatchResult bestPositionAlongTubeMatch(const MuonSegment& seg1, const MuonSegment& seg2,
                                                          const Amg::Vector3D& segDir1Min) const;

    SegmentPhiMatchResult bestPhiMatch(const MuonSegment& seg1, const MuonSegment& seg2) const;

    SegmentPhiMatchResult bestPhiMatchAnalytic(const MuonSegment& seg1, const MuonSegment& seg2) const;

    Amg::Vector3D updateSegmentDirection(const MuonSegment& seg, double phi) const;

    SegmentGeometrySummary segmentGeometrySummary(const MuonSegment& seg) const;
 
  private:
    double checkPhiHitConsistency(const EventContext& ctx, const Muon::MuonSegment& segment, SegmentPhiMatchResult& phiMatchResult,
                                  SegmentPositionMatchResult& matchResult) const;
 
    Amg::Vector3D estimateSegmentDirection(const MuonSegment& seg1, const MuonSegment& seg2, double& phi,
                                           double& stereoangle) const;
 
    ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{
        this,
        "MuonIdHelperSvc",
        "Muon::MuonIdHelperSvc/MuonIdHelperSvc",
    };
    ServiceHandle<IMuonEDMHelperSvc> m_edmHelperSvc{
        this,
        "edmHelper",
        "Muon::MuonEDMHelperSvc/MuonEDMHelperSvc",
        "Handle to the service providing the IMuonEDMHelperSvc interface",
    };  
 
    PublicToolHandle<MuonEDMPrinterTool> m_printer{
        this,
        "Printer",
        "Muon::MuonEDMPrinterTool/MuonEDMPrinterTool",
    };  
    ToolHandle<Trk::IPropagator> m_propagator{
        this,
        "AtlasRungeKuttaPropagator",
        "Trk::RungeKuttaPropagator/AtlasRungeKuttaPropagator",
    };
    ToolHandle<Trk::IResidualPullCalculator> m_pullCalculator{
        this,
        "ResidualPullCalculator",
        "Trk::ResidualPullCalculator/ResidualPullCalculator",
    };
 
    Trk::MagneticFieldProperties m_magFieldProperties;  
};
 
}  // namespace Muon
 
#endif