MdtCalibInput.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MDTCALIBINTEFACES_MDTCALIBINPUT_H
#define MDTCALIBINTEFACES_MDTCALIBINPUT_H
 
#include "GeoPrimitives/GeoPrimitives.h"
 
#include <MuonReadoutGeometry/MdtReadoutElement.h>
#include <MuonReadoutGeometry/MuonDetectorManager.h>
#include <MuonReadoutGeometryR4/MdtReadoutElement.h>
#include <MuonReadoutGeometryR4/MuonDetectorManager.h>
#include <MuonDigitContainer/MdtDigit.h>
#include <MuonPrepRawData/MdtPrepData.h>
#include <TrkSurfaces/StraightLineSurface.h>
 
 
#include <GaudiKernel/PhysicalConstants.h>
#include <CxxUtils/CachedUniquePtr.h>
 
#include <iostream>
 
 
class MdtCalibInput {
public:
    MdtCalibInput(const MdtDigit& digit,
                  const MuonGM::MuonDetectorManager& detMgr);
    
    MdtCalibInput(const MdtDigit& digit,
                  const MuonGMR4::MuonDetectorManager& detMgr,
                  const ActsGeometryContext& gctx);
    
    MdtCalibInput(const Muon::MdtPrepData& prd);
    MdtCalibInput(const Identifier& id, 
                  const int adc, 
                  const int tdc, 
                  const Amg::Vector3D& globalPos);
 
    const Identifier& identify() const;
    int tdc() const;
    int adc() const;
    bool isMasked() const;
    const MuonGM::MdtReadoutElement* legacyDescriptor() const;
    const MuonGMR4::MdtReadoutElement* decriptor() const;
    
    const Amg::Vector3D& globalPos() const;
    const Amg::Vector3D& closestApproach() const;
    void setClosestApproach(const Amg::Vector3D& approach);
 
    const Amg::Vector3D& trackDirection() const;
    void setTrackDirection(const Amg::Vector3D& trackDir);
    
 
    double timeOfFlight() const;
    void setTimeOfFlight(const double toF);
 
    double triggerTime() const;
    void setTriggerTime(const double trigTime);
    double distanceToTrack() const;
    enum BFieldComp{
      alongWire = 0,
      alongTrack = 1,
    };
    Amg::Vector2D projectMagneticField(const Amg::Vector3D& fieldInGlob) const;
    double signalPropagationDistance() const;
    const Trk::SaggedLineSurface& idealSurface() const;
 
    const Trk::StraightLineSurface& saggedSurface() const;
    std::unique_ptr<Trk::StraightLineSurface> releaseSurface();
 
    const Amg::Vector3D& surfaceCenter() const;
    const Amg::Vector3D& saggedSurfCenter() const;
    double tubeLength() const;
    double readOutSide() const;
private:
  
  Identifier m_id{0};
  bool m_isMasked{false};
  int m_adc{0};
  int m_tdc{0};
  const MuonGM::MdtReadoutElement* m_legRE{nullptr};
 
  const ActsGeometryContext* m_gctx{nullptr};
  const MuonGMR4::MdtReadoutElement* m_RE{nullptr};
  IdentifierHash m_hash{};
  
  Amg::Transform3D m_globToLoc{Amg::Transform3D::Identity()};
  Amg::Vector3D m_globPos{Amg::Vector3D::Zero()};
  Amg::Vector3D m_approach{m_globPos};
  Amg::Vector3D m_trackDir{Amg::Vector3D::Zero()};
  
  static constexpr double s_inverseSpeed{1. / Gaudi::Units::c_light};
  double m_ToF{m_globPos.mag() * s_inverseSpeed};
  double m_trigTime{0.};
  double m_distToTrack{0.};
 
  CxxUtils::CachedUniquePtrT<Trk::StraightLineSurface> m_saggedSurf{};
 
};
 
std::ostream& operator<<(std::ostream& ostr, const MdtCalibInput& input);
 
 
#endif