CDCTangent.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <tracking/trackingUtilities/eventdata/hits/CDCRLWireHitPair.h>
 
#include <tracking/trackingUtilities/geometry/ParameterLine2D.h>
#include <framework/geometry/VectorUtil.h>
 
#include <Math/Vector2D.h>
 
#include <iosfwd>
 
namespace Belle2 {
  namespace TrackingUtilities {
    class CDCRLWireHit;
    class CDCRecoHit2D;


    class CDCTangent : public CDCRLWireHitPair {
    public:
      CDCTangent() = default;

      explicit CDCTangent(const CDCRLWireHitPair& rlWireHitPair);

      CDCTangent(const CDCRLWireHit& fromRLWireHit,
                 const CDCRLWireHit& toRLWireHit);

      CDCTangent(const CDCRLWireHitPair& rlWireHitPair,
                 const ParameterLine2D& line);

      CDCTangent(const CDCRLWireHit& fromRLWireHit,
                 const CDCRLWireHit& toRLWireHit,
                 const ParameterLine2D& line);

      const ROOT::Math::XYVector& getFromRecoPos2D() const
      { return getLine().support(); }

      ROOT::Math::XYVector getFromRecoDisp2D() const;

      ROOT::Math::XYVector getToRecoPos2D() const
      { return getLine().at(1); }

      ROOT::Math::XYVector getToRecoDisp2D() const;

      const ROOT::Math::XYVector& getFlightVec2D() const
      { return getLine().tangential(); }

      double getCosFlightDifference(const CDCTangent& tangent) const
      { return VectorUtil::CosPhi(getFlightVec2D(), tangent.getFlightVec2D()); }

      CDCRecoHit2D getFromRecoHit2D() const;

      CDCRecoHit2D getToRecoHit2D() const;

      void adjustLine();

      void adjustRLInfo();

      void reverse() override;

      CDCTangent reversed() const;

      const ParameterLine2D& getLine() const
      { return m_line; }
 
    private:
      ParameterLine2D m_line;
    };

    std::ostream& operator<<(std::ostream& output, const CDCTangent& tangent);
  }
}