MdtCalibInput.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MDTCALIBINTEFACES_MDTCALIBINPUT_H
#define MDTCALIBINTEFACES_MDTCALIBINPUT_H
 
#include "GeoPrimitives/GeoPrimitives.h"
 
#include <GaudiKernel/PhysicalConstants.h>
#include <xAODMuonPrepData/MdtDriftCircleFwd.h>
#include <Identifier/Identifier.h>
#include <Identifier/IdentifierHash.h>
#include <variant>
namespace MuonGM{
  class MuonDetectorManager;
  class MdtReadoutElement;
}
namespace MuonGMR4{
  class MuonDetectorManager;
  class MdtReadoutElement;
}
namespace Muon{
  class MdtPrepData;
}
 
namespace Trk {
  class StraightLineSurface;
}
 
namespace ActsTrk {
  class GeometryContext;
}
class MdtDigit;
 
class MdtCalibInput {
    public:
      MdtCalibInput(const Identifier& id,
                    const int16_t adc,
                    const int16_t tdc,
                    const MuonGMR4::MdtReadoutElement* reEle,
                    const ActsTrk::GeometryContext& gctx);
    
      MdtCalibInput(const Identifier& id,
                    const int16_t adc,
                    const int16_t tdc,
                    const MuonGM::MdtReadoutElement* reEle);

      MdtCalibInput(const MdtDigit& digit,
                    const MuonGM::MuonDetectorManager& detMgr);
      
      MdtCalibInput(const MdtDigit& digit,
                    const MuonGMR4::MuonDetectorManager& detMgr,
                    const ActsTrk::GeometryContext& gctx);
      MdtCalibInput(const Muon::MdtPrepData& prd);
      MdtCalibInput(const xAOD::MdtDriftCircle& prd,
                    const ActsTrk::GeometryContext& gctx);
 
 
      MdtCalibInput(MdtCalibInput&& other) = default;
      MdtCalibInput& operator=(MdtCalibInput&& other) = default;
 
      ~MdtCalibInput();
      const Identifier& identify() const;
      int16_t tdc() const;
      int16_t adc() const;
      const MuonGM::MdtReadoutElement* legacyDescriptor() const;
      const MuonGMR4::MdtReadoutElement* decriptor() const;
    
      const Amg::Vector3D& closestApproach() const;
      void setClosestApproach(const Amg::Vector3D& approach);

      const Amg::Vector3D& trackDirection() const;
      void setTrackDirection(const Amg::Vector3D& trackDir,
                             bool hasPhi);
      bool trackDirHasPhi() const;

      double timeOfFlight() const;
      void setTimeOfFlight(const double toF);

      double triggerTime() const;
      void setTriggerTime(const double trigTime);
      double distanceToTrack() const;
      enum class BFieldComp{
        alongWire = 0,
        alongTrack = 1,
      };

      Amg::Vector2D projectMagneticField(const Amg::Vector3D& fieldInGlob) const;
      double signalPropagationDistance() const;
      const Trk::StraightLineSurface& legacySurface() const;

      const Amg::Vector3D& surfaceCenter() const;
      double tubeLength() const;
      double readOutSide() const;
      double innerTubeR() const;
  private:
    const Amg::Transform3D& localToGlobal() const;
    Amg::Vector3D center() const;
    Identifier m_id{};
    int16_t m_adc{0};
    int16_t m_tdc{0};

    const ActsTrk::GeometryContext* m_gctx{nullptr};
    using ReadoutEle_t = std::variant<const MuonGM::MdtReadoutElement*, const MuonGMR4::MdtReadoutElement*>;
    ReadoutEle_t m_RE{};
    IdentifierHash m_hash{};  
    Amg::Vector3D m_approach{center()};
    Amg::Vector3D m_trackDir{Amg::Vector3D::Zero()};
    bool m_trackHasPhi{false};
    static constexpr double s_inverseSpeed{1. / Gaudi::Units::c_light};
    double m_ToF{center().mag() * s_inverseSpeed};
    double m_trigTime{0.};
    double m_distToTrack{0.};
 
};
 
std::ostream& operator<<(std::ostream& ostr, const MdtCalibInput& input);
 
 
#endif