PmaSegment3D.h

Go to the documentation of this file.

#ifndef PmaSegment3D_h
#define PmaSegment3D_h

#include "larreco/RecoAlg/PMAlg/PmaElement3D.h"
#include "larreco/RecoAlg/PMAlg/PmaNode3D.h"
#include "larreco/RecoAlg/PMAlg/SortedObjects.h"

class TVector2;
#include "TVector3.h"
#include "larreco/RecoAlg/PMAlg/Utilities.h"

namespace pma {
  class Segment3D;
  class Track3D; // only declare here to keep "parent" of Segment3D
}

class pma::Segment3D : public pma::Element3D, public pma::SortedObjectBase {
public:
  Segment3D(void) : fParent(0) {}
  Segment3D(pma::Track3D* trk, pma::Node3D* vstart, pma::Node3D* vstop);

  Vector3D Start(void) const
  {
    auto const& p = static_cast<Node3D*>(Prev())->Point3D();
    return Vector3D(p.X(), p.Y(), p.Z());
  }
  Vector3D End(void) const
  {
    auto const& p = static_cast<Node3D*>(Next())->Point3D();
    return Vector3D(p.X(), p.Y(), p.Z());
  }

  double GetDistance2To(const TVector3& p3d) const override;

  double GetDistance2To(const TVector2& p2d, unsigned int view) const override;

  pma::Vector3D GetDirection3D(void) const override;

  TVector3 GetProjection(const TVector2& p, unsigned int view) const;

  TVector3 GetUnconstrainedProj3D(const TVector2& p2d, unsigned int view) const override;

  void SetProjection(pma::Hit3D& h) const override;

  double Length2(void) const override;

  pma::Track3D* Parent(void) const { return fParent; }

private:
  Segment3D(const pma::Segment3D& src);

  double SumDist2Hits(void) const override;

  pma::Track3D* fParent;

  static double GetDist2(const TVector3& psrc, const TVector3& p0, const TVector3& p1);
  static double GetDist2(const TVector2& psrc, const TVector2& p0, const TVector2& p1);
};

#endif