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"

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