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 25 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 45 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 32 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 196 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 102 of file Seed.cxx.

References fSeedDirection, and fSeedDirectionError.

Referenced by trkf::SeedAna::analyze(), trkf::FeatureTracker::CheckSeedLineInt(), trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::FeatureTracker::ExtendSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), GetAngle(), trkf::BezierCurveHelper::GetBezierPointCubic(), trkf::BezierCurveHelper::GetBezierPointQuartic(), trkf::BezierCurveHelper::GetBezierPointsCubic(), trkf::BezierCurveHelper::GetBezierPointsQuartic(), trkf::SeedFinderAlgorithm::GetHitDistAndProj(), trkf::BezierTrackerAlgorithm::GetSeedDirProjected(), 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 243 of file Seed.cxx.

References fSeedPoint, and GetPoint().

Referenced by trkf::BezierTrackerAlgorithm::EvaluateOccupancy().

  {
    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 257 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 159 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 112 of file Seed.cxx.

References fSeedPoint, and fSeedPointError.

Referenced by trkf::SeedAna::analyze(), trkf::FeatureTracker::CheckSeedLineInt(), trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::FeatureTracker::ExtendSeed(), trkf::SeedFinderAlgorithm::FindSeeds(), trkf::BezierCurveHelper::GetBezierPointCubic(), trkf::BezierCurveHelper::GetBezierPointQuartic(), trkf::BezierCurveHelper::GetBezierPointsCubic(), trkf::BezierCurveHelper::GetBezierPointsQuartic(), trkf::BezierTrack::GetClosestApproach(), 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 322 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 168 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 226 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 183 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 74 of file Seed.cxx.

References fIsValid.

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

  {
    return fIsValid;
  }

void recob::Seed::Print

(

)

const

Definition at line 59 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 81 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 137 of file Seed.cxx.

References fIsValid, fSeedDirection, and fSeedDirectionError.

Referenced by trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::FeatureTracker::ExtendSeed(), 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 123 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 148 of file Seed.cxx.

References fIsValid, fSeedPoint, and fSeedPointError.

Referenced by trkf::SeedFinderAlgorithm::ConsolidateSeed(), trkf::FeatureTracker::ExtendSeed(), 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 130 of file Seed.cxx.

References SetPoint().

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

void recob::Seed::SetValidity

(

bool

Validity

)

Definition at line 94 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 361 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 345 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 31 of file Seed.h.

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

double recob::Seed::fSeedDirection[3]

private

Definition at line 28 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 30 of file Seed.h.

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

double recob::Seed::fSeedPoint[3]

private

Definition at line 27 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 29 of file Seed.h.

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

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/lardataobj/v07_06_00/source/lardataobj/RecoBase/Seed.h
• /cvmfs/larsoft.opensciencegrid.org/products/lardataobj/v07_06_00/source/lardataobj/RecoBase/Seed.cxx