recob::Seed Class Reference

#include "Seed.h"

Public Member Functions
	Seed ()
	Seed (double *Pt, double *Dir, double *PtErr, double *DirErr)
void	GetDirection (double *Dir, double *Err) const
void	GetPoint (double *Pt, double *Err) const
double	GetLength () const
void	Print () const
void	SetDirection (double *Dir, double *Err)
void	SetPoint (double *Pt, double *Err)
double	GetAngle (Seed const &AnotherSeed) const
double	GetProjDiscrepancy (Seed const &AnotherSeed) const
double	GetProjAngleDiscrepancy (Seed const &AnotherSeed) const
double	GetDistance (Seed const &AnotherSeed) const
Seed	Reverse ()
void	GetVectorBetween (Seed const &AnotherSeed, double *xyz) const
double	GetDistanceFrom (SpacePoint const &SomePoint) const
int	GetPointingSign (Seed const &AnotherSeed) const
	Seed (double *Pt, double *Dir)
void	SetDirection (double *Dir)
void	SetPoint (double *Pt)
bool	IsValid () const
void	SetValidity (bool Validity)

Private Attributes
double	fSeedPoint [3]
double	fSeedDirection [3]
double	fSeedPointError [3]
double	fSeedDirectionError [3]
bool	fIsValid

Friends
std::ostream &	operator<< (std::ostream &stream, Seed const &a)
bool	operator< (const Seed &a, const Seed &b)

Detailed Description

Definition at line 20 of file Seed.h.

Constructor & Destructor Documentation

recob::Seed::Seed

(

)

Definition at line 20 of file Seed.cxx.

References fIsValid.

Referenced by Reverse().

  {
    fIsValid = false;
  }

recob::Seed::Seed	(	double *	Pt,
		double *	Dir,
		double *	PtErr,
		double *	DirErr
	)

Definition at line 38 of file Seed.cxx.

References fIsValid, fSeedDirection, fSeedDirectionError, fSeedPoint, and fSeedPointError.

  {
    for (int i = 0; i != 3; i++) {
      fSeedPoint[i] = Pt[i];
      fSeedDirection[i] = Dir[i];
      fSeedPointError[i] = PtErr[i];
      fSeedDirectionError[i] = DirErr[i];
    }
    fIsValid = true;
  }

recob::Seed::Seed	(	double *	Pt,
		double *	Dir
	)

Definition at line 26 of file Seed.cxx.

References fIsValid, fSeedDirection, fSeedDirectionError, fSeedPoint, and fSeedPointError.

  {
    for (int i = 0; i != 3; i++) {
      fSeedPoint[i] = Pt[i];
      fSeedDirection[i] = Dir[i];
      fSeedPointError[i] = 0;
      fSeedDirectionError[i] = 0;
    }
    fIsValid = true;
  }

Member Function Documentation

double recob::Seed::GetAngle

(

Seed const &

AnotherSeed

)

const

Definition at line 165 of file Seed.cxx.

References fSeedDirection, GetDirection(), and GetLength().

  {
    double OtherDir[3];
    double OtherDirErr[3];
    AnotherSeed.GetDirection(OtherDir, OtherDirErr);

    double OtherMag = AnotherSeed.GetLength();
    double ThisMag = GetLength();

    /*
    std::cout<<"Seed angle calc: " <<
      OtherMag<< " " << ThisMag << " " <<
      ( (OtherDir[0]*fSeedDirection[0] +
   OtherDir[1]*fSeedDirection[1] +
   OtherDir[2]*fSeedDirection[2] ) / (ThisMag*OtherMag)) <<
      std::endl;
    */

    // Need a tiny offset, as acos(1.0) gives unpredictable results due to
    // floating point precision
    double eta = 0.00001;

    return std::acos((OtherDir[0] * fSeedDirection[0] + OtherDir[1] * fSeedDirection[1] +
                      OtherDir[2] * fSeedDirection[2] - eta) /
                     (ThisMag * OtherMag));
  }

void recob::Seed::GetDirection	(	double *	Dir,
		double *	Err
	)		const

Definition at line 80 of file Seed.cxx.

References fSeedDirection, and fSeedDirectionError.

Referenced by trkf::SeedAna::analyze(), trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), GetAngle(), trkf::SeedFinderAlgorithm::GetHitDistAndProj(), trkf::Track3DKalmanHitAlg::makeSurface(), and TrackProducerFromPFParticle::produce().

  {
    for (int i = 0; i != 3; i++) {
      rDir[i] = fSeedDirection[i];
      if (rDirErr) rDirErr[i] = fSeedDirectionError[i];
    }
  }

double recob::Seed::GetDistance

(

Seed const &

AnotherSeed

)

const

Definition at line 209 of file Seed.cxx.

References fSeedPoint, and GetPoint().

  {
    double OtherPt[3];
    double OtherPtErr[3];

    AnotherSeed.GetPoint(OtherPt, OtherPtErr);

    return pow(pow(OtherPt[0] - fSeedPoint[0], 2) + pow(OtherPt[1] - fSeedPoint[1], 2) +
                 pow(OtherPt[2] - fSeedPoint[2], 2),
               0.5);
  }

double recob::Seed::GetDistanceFrom

(

recob::SpacePoint const &

SomePoint

)

const

Definition at line 222 of file Seed.cxx.

References recob::CrossProd(), fSeedDirection, fSeedPoint, and recob::SpacePoint::XYZ().

  {
    double SPxyz[3];
    SPxyz[0] = SomePoint.XYZ()[0];
    SPxyz[1] = SomePoint.XYZ()[1];
    SPxyz[2] = SomePoint.XYZ()[2];

    //  std::cout<<"Seed Dir Vec " << fSeedDirection[0]<<" " <<fSeedDirection[1]<< " " << fSeedDirection[2]<<std::endl;

    double ThisSeedLength =
      pow(pow(fSeedDirection[0], 2) + pow(fSeedDirection[1], 2) + pow(fSeedDirection[2], 2), 0.5);

    double SPProjOnSeed = (fSeedDirection[0] * (SPxyz[0] - fSeedPoint[0]) +
                           fSeedDirection[1] * (SPxyz[1] - fSeedPoint[1]) +
                           fSeedDirection[2] * (SPxyz[2] - fSeedPoint[2])) /
                          ThisSeedLength;

    // std::cout<<"proj : " <<SPProjOnSeed<<std::endl;
    // std::cout<<"Seed len :" << ThisSeedLength<<std::endl;

    if (SPProjOnSeed > (ThisSeedLength)) {
      //  std::cout<<"Seed over end"<<std::endl;
      return pow(pow(fSeedPoint[0] + fSeedDirection[0] - SPxyz[0], 2) +
                   pow(fSeedPoint[1] + fSeedDirection[1] - SPxyz[1], 2) +
                   pow(fSeedPoint[2] + fSeedDirection[2] - SPxyz[2], 2),
                 0.5);
    }
    else if (SPProjOnSeed < (0 - ThisSeedLength)) {
      //  std::cout<<"Seed under end"<<std::endl;
      return pow(pow(fSeedPoint[0] - fSeedDirection[0] - SPxyz[0], 2) +
                   pow(fSeedPoint[1] - fSeedDirection[1] - SPxyz[1], 2) +
                   pow(fSeedPoint[2] - fSeedDirection[2] - SPxyz[2], 2),
                 0.5);
    }
    else {
      //  std::cout<<"Seed valid region"<<std::endl;
      double crossprod[3];
      CrossProd(fSeedPoint[0] + fSeedDirection[0] - SPxyz[0],
                fSeedPoint[1] + fSeedDirection[1] - SPxyz[1],
                fSeedPoint[2] + fSeedDirection[2] - SPxyz[2],
                SPxyz[0] - fSeedPoint[0],
                SPxyz[1] - fSeedPoint[1],
                SPxyz[2] - fSeedPoint[2],
                crossprod[0],
                crossprod[1],
                crossprod[2]);

      return pow(pow(crossprod[0], 2) + pow(crossprod[1], 2) + pow(crossprod[2], 2), 0.5) /
             pow(pow(fSeedDirection[0], 2) + pow(fSeedDirection[1], 2) + pow(fSeedDirection[2], 2),
                 0.5);
    }
  }

double recob::Seed::GetLength

(

)

const

Definition at line 132 of file Seed.cxx.

References fSeedDirection.

Referenced by trkf::SeedFinderAlgorithm::FindSeeds(), and GetAngle().

  {
    return pow(pow(fSeedDirection[0], 2) + pow(fSeedDirection[1], 2) + pow(fSeedDirection[2], 2),
               0.5);
  }

void recob::Seed::GetPoint	(	double *	Pt,
		double *	Err
	)		const

Definition at line 89 of file Seed.cxx.

References fSeedPoint, and fSeedPointError.

Referenced by trkf::SeedAna::analyze(), trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), GetDistance(), trkf::SeedFinderAlgorithm::GetHitDistAndProj(), GetPointingSign(), GetProjAngleDiscrepancy(), GetProjDiscrepancy(), GetVectorBetween(), trkf::Track3DKalmanHitAlg::makeSurface(), and TrackProducerFromPFParticle::produce().

  {
    for (int i = 0; i != 3; i++) {
      rPt[i] = fSeedPoint[i];
      if (rPtErr) rPtErr[i] = fSeedPointError[i];
    }
  }

int recob::Seed::GetPointingSign

(

Seed const &

AnotherSeed

)

const

Definition at line 276 of file Seed.cxx.

References tca::DotProd(), fSeedDirection, fSeedPoint, and GetPoint().

  {
    double OtherPos[3], OtherErr[3];
    AnotherSeed.GetPoint(OtherPos, OtherErr);
    double DotProd = (OtherPos[0] - fSeedPoint[0]) * fSeedDirection[0] +
                     (OtherPos[1] - fSeedPoint[1]) * fSeedDirection[1] +
                     (OtherPos[2] - fSeedPoint[2]) * fSeedDirection[2];
    return ((DotProd > 0) - (DotProd < 0));
  }

double recob::Seed::GetProjAngleDiscrepancy

(

Seed const &

AnotherSeed

)

const

Definition at line 139 of file Seed.cxx.

References fSeedDirection, fSeedPoint, and GetPoint().

  {
    double OtherPt[3];
    double OtherPtErr[3];

    AnotherSeed.GetPoint(OtherPt, OtherPtErr);

    TVector3 OtherPtV(OtherPt[0], OtherPt[1], OtherPt[2]);
    TVector3 ThisDirV(fSeedDirection[0], fSeedDirection[1], fSeedDirection[2]);
    TVector3 ThisPtV(fSeedPoint[0], fSeedPoint[1], fSeedPoint[2]);

    return (OtherPtV - ThisPtV).Angle(ThisDirV.Unit());
  }

double recob::Seed::GetProjDiscrepancy

(

Seed const &

AnotherSeed

)

const

Definition at line 193 of file Seed.cxx.

References fSeedDirection, fSeedPoint, and GetPoint().

  {
    double OtherPt[3];
    double OtherPtErr[3];

    AnotherSeed.GetPoint(OtherPt, OtherPtErr);

    TVector3 OtherPtV(OtherPt[0], OtherPt[1], OtherPt[2]);
    TVector3 ThisDirV(fSeedDirection[0], fSeedDirection[1], fSeedDirection[2]);
    TVector3 ThisPtV(fSeedPoint[0], fSeedPoint[1], fSeedPoint[2]);

    return ((OtherPtV - ThisPtV) - ThisDirV.Unit() * (ThisDirV.Unit().Dot(OtherPtV - ThisPtV)))
      .Mag();
  }

void recob::Seed::GetVectorBetween	(	Seed const &	AnotherSeed,
		double *	xyz
	)		const

Definition at line 154 of file Seed.cxx.

References fSeedPoint, and GetPoint().

  {
    double xyzother[3], err[3];
    AnotherSeed.GetPoint(xyzother, err);

    xyz[0] = xyzother[0] - fSeedPoint[0];
    xyz[1] = xyzother[1] - fSeedPoint[1];
    xyz[2] = xyzother[2] - fSeedPoint[2];
  }

bool recob::Seed::IsValid

(

)

const

Definition at line 58 of file Seed.cxx.

References fIsValid.

Referenced by trkf::SeedFinderAlgorithm::FindSeeds().

  {
    return fIsValid;
  }

void recob::Seed::Print

(

)

const

Definition at line 50 of file Seed.cxx.

References fSeedDirection, and fSeedPoint.

  {
    std::cout << "Printing seed contents : " << fSeedPoint[0] << " " << fSeedPoint[1] << " "
              << fSeedPoint[2] << ", " << fSeedDirection[0] << " " << fSeedDirection[1] << " "
              << fSeedDirection[2] << std::endl;
  }

Seed recob::Seed::Reverse

(

)

Definition at line 64 of file Seed.cxx.

References fSeedDirection, fSeedDirectionError, fSeedPoint, fSeedPointError, n, and Seed().

  {
    double NewSeedDir[3];
    for (size_t n = 0; n != 3; ++n) {
      NewSeedDir[n] = -fSeedDirection[n];
    }
    return Seed(fSeedPoint, NewSeedDir, fSeedPointError, fSeedDirectionError);
  }

void recob::Seed::SetDirection	(	double *	Dir,
		double *	Err
	)

Definition at line 112 of file Seed.cxx.

References fIsValid, fSeedDirection, and fSeedDirectionError.

Referenced by trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), and SetDirection().

  {
    for (int i = 0; i != 3; i++) {
      fSeedDirection[i] = Dir[i];
      fSeedDirectionError[i] = Err[i];
    }
    fIsValid = true;
  }

void recob::Seed::SetDirection

(

double *

Dir

)

Definition at line 98 of file Seed.cxx.

References SetDirection().

  {
    double Empty[3] = {0, 0, 0};
    SetDirection(Dir, Empty);
  }

void recob::Seed::SetPoint	(	double *	Pt,
		double *	Err
	)

Definition at line 122 of file Seed.cxx.

References fIsValid, fSeedPoint, and fSeedPointError.

Referenced by trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), and SetPoint().

  {
    for (int i = 0; i != 3; i++) {
      fSeedPoint[i] = Pt[i];
      fSeedPointError[i] = Err[i];
    }
    fIsValid = true;
  }

void recob::Seed::SetPoint

(

double *

Pt

)

Definition at line 105 of file Seed.cxx.

References SetPoint().

  {
    double Empty[3] = {0, 0, 0};
    SetPoint(Pt, Empty);
  }

void recob::Seed::SetValidity

(

bool

Validity

)

Definition at line 74 of file Seed.cxx.

References fIsValid.

Referenced by trkf::SeedFinderAlgorithm::ConsolidateSeed(), and trkf::SeedFinderAlgorithm::FindSeeds().

  {
    fIsValid = Validity;
  }

Friends And Related Function Documentation

bool operator< ( const Seed &  a, const Seed &  b  )

friend

Definition at line 315 of file Seed.cxx.

  {
    if (a.fSeedPoint[2] != b.fSeedPoint[2]) { return a.fSeedPoint[2] < b.fSeedPoint[2]; }
    else if (a.fSeedPoint[1] != b.fSeedPoint[1]) {
      return a.fSeedPoint[1] < b.fSeedPoint[1];
    }

    return a.fSeedPoint[0] < b.fSeedPoint[0];
  }

std::ostream& operator<< ( std::ostream &  stream, Seed const &  a  )

friend

Definition at line 303 of file Seed.cxx.

  {
    o << "Printing seed contents : " << a.fSeedPoint[0] << " " << a.fSeedPoint[1] << " "
      << a.fSeedPoint[2] << ", " << a.fSeedDirection[0] << " " << a.fSeedDirection[1] << " "
      << a.fSeedDirection[2];

    return o;
  }

Member Data Documentation

bool recob::Seed::fIsValid

private

Definition at line 30 of file Seed.h.

Referenced by IsValid(), Seed(), SetDirection(), SetPoint(), and SetValidity().

double recob::Seed::fSeedDirection[3]

private

Definition at line 27 of file Seed.h.

Referenced by GetAngle(), GetDirection(), GetDistanceFrom(), GetLength(), GetPointingSign(), GetProjAngleDiscrepancy(), GetProjDiscrepancy(), recob::operator<<(), Print(), Reverse(), Seed(), and SetDirection().

double recob::Seed::fSeedDirectionError[3]

private

Definition at line 29 of file Seed.h.

Referenced by GetDirection(), Reverse(), Seed(), and SetDirection().

double recob::Seed::fSeedPoint[3]

private

Definition at line 26 of file Seed.h.

Referenced by GetDistance(), GetDistanceFrom(), GetPoint(), GetPointingSign(), GetProjAngleDiscrepancy(), GetProjDiscrepancy(), GetVectorBetween(), recob::operator<(), recob::operator<<(), Print(), Reverse(), Seed(), and SetPoint().

double recob::Seed::fSeedPointError[3]

private

Definition at line 28 of file Seed.h.

Referenced by GetPoint(), Reverse(), Seed(), and SetPoint().

---

The documentation for this class was generated from the following files:

• lardataobj/v10_01_00/source/lardataobj/RecoBase/Seed.h
• lardataobj/v10_01_00/source/lardataobj/RecoBase/Seed.cxx