recob::TrackTrajectory Class Reference

A trajectory in space reconstructed from hits. More...

#include "TrackTrajectory.h"

Inheritance diagram for recob::TrackTrajectory:

Public Types
using	Coord_t = tracking::Coord_t
	Type used for coordinates and values in general. More...
using	Point_t = tracking::Point_t
	Type for representation of position in physical 3D space. More...
using	Vector_t = tracking::Vector_t
	Type for representation of momenta in 3D space. More...
using	PointFlags_t = recob::TrajectoryPointFlags
	Type for flags of a point/hit. More...
using	flag = PointFlags_t::flag
	Flag traits (including the definition of flag mnemonics). More...
using	Positions_t = tracking::Positions_t
	Type of trajectory point list. More...
using	Momenta_t = tracking::Momenta_t
	Type of momentum list. More...
using	Flags_t = std::vector< PointFlags_t >
	Type of point flag list. More...
using	Ends_t = Trajectory_t::Ends_t
	Mnemonics for the access to begin and end of trajectory. More...
using	TrajectoryPoint_t = tracking::TrajectoryPoint_t
	A point in the trajectory, with position and momentum. More...
using	Rotation_t = tracking::Rotation_t
	Type for representation of space rotations. More...

Public Member Functions
	TrackTrajectory ()=default
	Default constructor; do not use it! it's needed by ROOT I/O. More...
	TrackTrajectory (Positions_t &&positions, Momenta_t &&momenta, Flags_t &&flags, bool hasMomenta)
	Constructor: specifies all the data for the trajectory. More...
	TrackTrajectory (const Trajectory &traj, Flags_t &&flags)
	Constructor: copies positions and momenta from an existing Trajecotry, adds the flags. More...
Trajectory_t const &	Trajectory () const
	Returns the plain trajectory of this object. More...
PointFlags_t const &	FlagsAtPoint (size_t i) const
	Returns the flags for the specified trajectory point. More...
Flags_t const &	Flags () const
	Returns all flags. More...
bool	HasValidPoint (size_t i) const
	Returns whether the specified point has NoPoint flag unset. More...
size_t	FirstValidPoint () const
	Returns the index of the first valid point in the trajectory. More...
size_t	NextValidPoint (size_t index) const
	Returns the index of the next valid point in the trajectory. More...
size_t	PreviousValidPoint (size_t index) const
	Returns the index of the previous valid point in the trajectory. More...
size_t	LastValidPoint () const
	Returns the index of the last valid point in the trajectory. More...
unsigned int	CountValidPoints () const
	Computes and returns the number of points with valid location. More...
Point_t const &	Vertex () const
	Returns the position of the first valid point of the trajectory [cm]. More...
Point_t const &	Start () const
	Returns the position of the first valid point of the trajectory [cm]. More...
Point_t const &	End () const
	Returns the position of the last valid point of the trajectory [cm]. More...
void	Extent (std::vector< double > &start, std::vector< double > &end) const
	Fills the first and last valid point in the trajectory. More...
std::pair< Point_t, Point_t >	Extent () const
	Returns a copy of the first and last valid point in the trajectory. More...
double	Length (size_t startAt=0) const
	Returns the approximate length of the trajectory. More...
size_t	NumberFitMomentum () const
Vector_t	VertexDirection () const
	Returns the direction of the trajectory at the first point. More...
Vector_t	StartDirection () const
	Returns the direction of the trajectory at the first point. More...
Vector_t	EndDirection () const
	Returns the direction of the trajectory at the last point. More...
double	Theta (size_t p) const
	Trajectory angle at point, with respect to positive z direction. More...
double	Theta () const
	Trajectory angle at start, with respect to positive z direction. More...
double	Phi (size_t p) const
	Azimuthal angle at a point on the trajectory, with respect to z. More...
double	Phi () const
	Azimuthal angle at a first valid point, with respect to z. More...
double	ZenithAngle (size_t p) const
	"Zenith" angle of trajectory, with respect to the vertical axis. More...
double	ZenithAngle () const
	"Zenith" angle of trajectory, with respect to the vertical axis. More...
double	AzimuthAngle (size_t p) const
	"Azimuth" angle of trajectory, with respect to the sky. More...
double	AzimuthAngle () const
	"Azimuth" angle of trajectory, with respect to the sky. More...
Vector_t const &	VertexMomentumVector () const
	Returns the momentum of the trajectory at the first valid point [GeV/c]. More...
Vector_t const &	StartMomentumVector () const
	Returns the momentum of the trajectory at the first valid point [GeV/c]. More...
Vector_t const &	EndMomentumVector () const
	Returns the momentum of the trajectory at the last valid point [GeV/c]. More...
double	VertexMomentum () const
double	StartMomentum () const
double	EndMomentum () const
void	Direction (double *start, double *end) const
	Fills the starting and ending direction of the trajectory. More...
std::pair< Vector_t, Vector_t >	Direction () const
	Returns the trajectory directions at first and last valid points. More...
template<typename Stream >
void	Dump (Stream &&out, unsigned int verbosity, std::string indent, std::string indentFirst) const
	Prints trajectory content into a stream. More...
template<typename Stream >
void	Dump (Stream &&out, unsigned int verbosity=1, std::string indent={}) const
	Prints trajectory content into a stream. More...
template<typename Stream >
void	LowLevelDump (Stream &&out, std::string indent, std::string indentFirst) const
	Prints low-level trajectory content into a stream. More...

Static Public Attributes
static constexpr unsigned int	MaxDumpVerbosity = 7
	Largest verbosity level supported by Dump(). More...
static constexpr size_t	InvalidIndex = std::numeric_limits<size_t>::max()
	Value returned on failed index queries. More...

Private Types
using	Trajectory_t = recob::Trajectory

Private Member Functions
template<int Dir>
size_t	ToValidPoint (size_t index) const
	Returns the index of the first valid point from index on. More...
bool	AtLeastValidTrajectoryPoints (unsigned int left) const
	Returns whether there are at least min valid points in the trajectory. More...
Access to trajectory information
size_t	NumberTrajectoryPoints () const
	Returns the number of stored trajectory points. More...
size_t	NPoints () const
	Returns the number of stored trajectory points. More...
size_t	FirstPoint () const
	Returns the index of the first point in the trajectory (yep, it's 0). More...
size_t	LastPoint () const
	Returns the index of the last point in the trajectory. More...
bool	HasPoint (size_t i) const
	Returns whether the specified trajectory point is available. More...
const Positions_t &	Positions () const
	Returns reference to stored vector of positions. More...
const Momenta_t &	Momenta () const
	Returns reference to stored vector of momenta. More...
bool	TrajectoryAtPoint (size_t i, TVector3 &pos, TVector3 &dir) const
	Fills position and direction at the specified trajectory point. More...
TrajectoryPoint_t	TrajectoryPoint (size_t i) const
	Returns position and momentum at the specified trajectory point. More...
Point_t const &	LocationAtPoint (size_t i) const
	Returns the position at the specified trajectory point. More...
Access to direction and momentum information
Vector_t	DirectionAtPoint (size_t i) const
	Computes and returns the direction of the trajectory at a point. More...
bool	HasMomentum () const
	Returns whether information about the momentum is available. More...
double	MomentumAtPoint (size_t i) const
	Computes and returns the modulus of the momentum at a point. More...
Vector_t const &	MomentumVectorAtPoint (size_t i) const
	Returns the momentum vector at a point. More...
Rotation_t	GlobalToLocalRotationAtPoint (size_t p) const
	Returns a rotation matrix that brings trajectory direction along z. More...
void	GlobalToLocalRotationAtPoint (unsigned int p, TMatrixD &rot) const
	Fills a rotation matrix that brings trajectory direction along z. More...
Rotation_t	LocalToGlobalRotationAtPoint (size_t p) const
	Returns a rotation matrix bringing relative directions to global. More...
void	LocalToGlobalRotationAtPoint (unsigned int p, TMatrixD &rot) const
	Fills a rotation matrix bringing relative directions to global. More...

Private Attributes
Flags_t	fFlags
	Flags of each of the points in trajectory. More...

Detailed Description

A trajectory in space reconstructed from hits.

---

See also

recob::Trajectory, recob::Track, recob::trackutil::makeTrajectory()

The track trajectory class contains a trajectory in 6D space representing the path walked by a particle. A trajectory point is made of a 3D position component (measured in centimeters) and a momentum component (measured in GeV/c). The associated hits are integral part of the track trajectory. To store additional point-by-point information, the track trajectory augments recob::Trajectory, of which it presents most of the interface, with point metadata called for convenience "flags".

Each point is formally associated to a reconstructed hit, as for recob::Trajectory requirements. One flag set is provided for each point, whose flags describe the point and/or the hit.

The meaning of the flags is documented also in the flags class recob::TrajectoryPointFlagTraits, which can be accessed as typedef flags in this class (e.g. recob::TrackTrajectory::flag::NoPoint).

Invariants

The same as recob::Trajectory, plus:

• there must be one flag set per trajectory point
• there must be at least two points with the flag NoPoint not set

Definition at line 62 of file TrackTrajectory.h.

Member Typedef Documentation

using recob::TrackTrajectory::Coord_t = tracking::Coord_t

Type used for coordinates and values in general.

Definition at line 67 of file TrackTrajectory.h.

using recob::TrackTrajectory::Ends_t = Trajectory_t::Ends_t

Mnemonics for the access to begin and end of trajectory.

Definition at line 91 of file TrackTrajectory.h.

using recob::TrackTrajectory::flag = PointFlags_t::flag

Flag traits (including the definition of flag mnemonics).

Definition at line 79 of file TrackTrajectory.h.

using recob::TrackTrajectory::Flags_t = std::vector<PointFlags_t>

Type of point flag list.

Definition at line 88 of file TrackTrajectory.h.

using recob::TrackTrajectory::Momenta_t = tracking::Momenta_t

Type of momentum list.

Definition at line 85 of file TrackTrajectory.h.

using recob::TrackTrajectory::Point_t = tracking::Point_t

Type for representation of position in physical 3D space.

Definition at line 70 of file TrackTrajectory.h.

using recob::TrackTrajectory::PointFlags_t = recob::TrajectoryPointFlags

Type for flags of a point/hit.

Definition at line 76 of file TrackTrajectory.h.

using recob::TrackTrajectory::Positions_t = tracking::Positions_t

Type of trajectory point list.

Definition at line 82 of file TrackTrajectory.h.

using recob::TrackTrajectory::Rotation_t = tracking::Rotation_t

Type for representation of space rotations.

Definition at line 97 of file TrackTrajectory.h.

using recob::TrackTrajectory::Trajectory_t = recob::Trajectory

private

Definition at line 63 of file TrackTrajectory.h.

using recob::TrackTrajectory::TrajectoryPoint_t = tracking::TrajectoryPoint_t

A point in the trajectory, with position and momentum.

Definition at line 94 of file TrackTrajectory.h.

using recob::TrackTrajectory::Vector_t = tracking::Vector_t

Type for representation of momenta in 3D space.

Definition at line 73 of file TrackTrajectory.h.

Constructor & Destructor Documentation

recob::TrackTrajectory::TrackTrajectory ( )

default

Default constructor; do not use it! it's needed by ROOT I/O.

recob::TrackTrajectory::TrackTrajectory	(	Positions_t &&	positions,
		Momenta_t &&	momenta,
		Flags_t &&	flags,
		bool	hasMomenta
	)

Constructor: specifies all the data for the trajectory.

Parameters

positions	(moved) trajectory as a sorted list of points
momenta	(moved) momentum along the trajectory, one per point
flags	(moved) flag sets, one flag set per point
hasMomenta	whether the information on momentum modulus is provided

Exceptions

std::runtime_error

if the invariants are violated

See also

recob::trackutil::makeTrackTrajectory()

The most convenient way to create a recob::Trajectory is to use recob::trackutil::makeTrackTrajectory().

Requirements

• one momentum is required for each trajectory point
• one flag is required for each trajectory point
• at least two points must be provided

Definition at line 21 of file TrackTrajectory.cxx.

References AtLeastValidTrajectoryPoints(), CountValidPoints(), fFlags, recob::Trajectory::NPoints(), and util::flags::to_string().

  : Trajectory_t(std::move(positions), std::move(momenta), hasMomenta)
  , fFlags(std::move(flags))
{
  // additional invariant check
  if (fFlags.size() != NPoints()) {
    throw std::runtime_error("recob::TrackTrajectory constructed with "
      + std::to_string(NPoints()) + " points "
      + std::to_string(fFlags.size())
      + " point flags! it requires the same number for both."
      );
  }
  if (!AtLeastValidTrajectoryPoints(2U)) {
    throw std::runtime_error("recob::TrackTrajectory constructed with only "
      + std::to_string(CountValidPoints())
      + " valid positions! at least 2 are required."
      );
  }
} // recob::TrackTrajectory::TrackTrajectory()

recob::TrackTrajectory::TrackTrajectory ( const Trajectory &  traj, Flags_t &&  flags  )

inline

Constructor: copies positions and momenta from an existing Trajecotry, adds the flags.

Parameters

traj	existing Trajectory
flags	(moved) flag sets, one flag set per point

Exceptions

std::runtime_error

if the invariants are violated

Definition at line 139 of file TrackTrajectory.h.

      : TrackTrajectory(Positions_t(traj.Positions()),Momenta_t(traj.Momenta()),std::move(flags),traj.HasMomentum()) {}

Member Function Documentation

bool recob::TrackTrajectory::AtLeastValidTrajectoryPoints ( unsigned int  left ) const

private

Returns whether there are at least min valid points in the trajectory.

Definition at line 107 of file TrackTrajectory.cxx.

References HasValidPoint(), art::left(), min, and recob::Trajectory::NPoints().

Referenced by Direction(), and TrackTrajectory().

{
  if (min == 0) return true;
  unsigned int left = min;
  for (size_t i = 0; i < NPoints(); ++i) {
    if (!HasValidPoint(i)) continue;
    if (--left == 0) return true;
  } // for
  return false;
  
} // moreThanTwoValidTrajectoryPoints()

double recob::TrackTrajectory::AzimuthAngle ( size_t  p ) const

inline

"Azimuth" angle of trajectory, with respect to the sky.

Parameters

p	the index point to extract the angle from (no default!)

Returns

the azimuth angle, in [radians]

See also

AzimuthAngle(), ZenithAngle(size_t), Phi(size_t)

The azimuth is returned, as defined in recob::Trajectory::Azimuth(), for point with the specified index p.

If the point is invalid, the behaviour is undefined.

Note

This function has no default value for p; if p is not specified at all, the method Azimuth() is called instead.

Definition at line 448 of file TrackTrajectory.h.

References recob::Trajectory::AzimuthAngle(), and Trajectory().

Referenced by recob::Track::AzimuthAngle().

      { return Trajectory().AzimuthAngle(p); }

double recob::TrackTrajectory::AzimuthAngle ( ) const

inline

"Azimuth" angle of trajectory, with respect to the sky.

Returns

the azimuth angle, in [radians]

See also

AzimuthAngle(size_t), ZenithAngle(), Phi()

Zenith(size_t), Theta()

The azimuth angle is returned, as defined in recob::Trajectory::Azimuth(), for the first valid point in the trajectory.

Note

This is not equivalent to Azimuth(0), but instead to Azimuth(FirstValidPoint()).

Definition at line 464 of file TrackTrajectory.h.

References FirstValidPoint().

      { return AzimuthAngle(FirstValidPoint()); }

unsigned int recob::TrackTrajectory::CountValidPoints

(

)

const

Computes and returns the number of points with valid location.

Returns

number of points in the trajectory with valid location

This method is slow, taking O(NPoints()) time.

Definition at line 45 of file TrackTrajectory.cxx.

References Extent(), HasValidPoint(), and recob::Trajectory::NPoints().

Referenced by recob::Track::CountValidPoints(), LastValidPoint(), and TrackTrajectory().

                                                          {
  
  unsigned int count = 0;
  for (size_t index = 0; index < NPoints(); ++index) {
    if (HasValidPoint(index)) ++count;
  } // for
  return count;
  
} // recob::TrackTrajectory::CountValidPoints()

void recob::TrackTrajectory::Direction	(	double *	start,
		double *	end
	)		const

Fills the starting and ending direction of the trajectory.

Parameters

start	(output) direction at the beginning of the trajectory
end	(output) direction at the end of the trajectory

The two arguments are expected to point each one to an area with room for at least three double numbers. The two filled vectors have norm 1.

The labelling of start and end is consistent within the trajectory but is not guaranteed to be physically correct.

Deprecated:

Use physics vectors instead (see Direction())

Definition at line 100 of file TrackTrajectory.cxx.

References AtLeastValidTrajectoryPoints(), EndDirection(), recob::details::legacy::FillTwoVectors(), and StartDirection().

Referenced by recob::Track::Direction().

                                                                     {
  details::legacy::FillTwoVectors(StartDirection(), EndDirection(), start, end); 
} // recob::TrackTrajectory::Direction()

std::pair<Vector_t, Vector_t> recob::TrackTrajectory::Direction ( ) const

inline

Returns the trajectory directions at first and last valid points.

Returns

a pair with the first and last direction

The two returned vectors have norm 1. The labelling of start and end is consistent within the trajectory but is not guaranteed to be physically correct.

Example:

recob::Trajectory::Vector_t startDir, endDir;

std::tie(startDir, endDir) = traj.Direction(); // assign start and end

Definition at line 539 of file TrackTrajectory.h.

References Dump(), EndDirection(), recob::Trajectory::GlobalToLocalRotationAtPoint(), art::detail::indent(), recob::Trajectory::LocalToGlobalRotationAtPoint(), and StartDirection().

Referenced by EndMomentum().

      { return { StartDirection(), EndDirection() }; }

recob::Trajectory::Vector_t recob::Trajectory::DirectionAtPoint ( size_t  i ) const

inherited

Computes and returns the direction of the trajectory at a point.

Parameters

i	index of the point in the trajectory

Returns

the direction at that point

The direction is computed as unit vector parallel to the momentum at that trajectory point. If the index is not contained in the trajectory, the result is undefined.

Definition at line 157 of file Trajectory.cxx.

References recob::Trajectory::HasMomentum(), and recob::Trajectory::MomentumVectorAtPoint().

Referenced by recob::Trajectory::AzimuthAngle(), recob::Track::DirectionAtPoint(), EndDirection(), recob::Trajectory::EndDirection(), recob::Trajectory::EndMomentum(), EndMomentum(), trkf::BezierTrack::FillSeedVector(), recob::Trajectory::GlobalToLocalRotationAtPoint(), trkf::TrajectoryMCSFitter::linearRegression(), recob::Trajectory::LocalToGlobalRotationAtPoint(), trkf::KalmanFilterTrajectoryFitter::restoreInputPoints(), trkf::KalmanFilterFinalTrackFitter::restoreInputPoints(), StartDirection(), recob::Trajectory::StartDirection(), and recob::Trajectory::ZenithAngle().

{
  
  auto const& mom = MomentumVectorAtPoint(i);
  return HasMomentum()? (mom / mom.R()): mom;
  
} // recob::Trajectory::DirectionAtPoint()

template<typename Stream >

void recob::TrackTrajectory::Dump	(	Stream &&	out,
		unsigned int	verbosity,
		std::string	indent,
		std::string	indentFirst
	)		const

Prints trajectory content into a stream.

Template Parameters

Stream

type of the output stream

Parameters

out	stream to output the information into
verbosity	verbosity level (default: 1)
indent	indentation string (default: none)
indentFirst	indentation for first output line (default: as indent)

The amount of information dumped to screen is regulated by the Indentation string is prepended to each line, and the first line has its own special indentation string (indentFirst).

The output can be multi-line, it ends with no end-of-line and it does not inserts an end-of-line at its beginning (unless that is explicitly inside indentFirst). The lowest verbosity is guaranteed to be on a single line.

Information printed out (verbosity argument)

• level 0: start position, direction, momentum modulus and number of points
• level 1: also end position, direction and momentum modulus
• level 2: also trajectory length
• level 3: also angles at start
• level 4: also 9 intermediate valid trajectory points
• level 5: also 10 more intermediate valid trajectory points (19 total)
• level 6: all valid trajectory points
• level 7: all trajectory points

Referenced by Direction(), and Dump().

template<typename Stream >

void recob::TrackTrajectory::Dump ( Stream &&  out, unsigned int  verbosity = 1, std::string  indent = {}  ) const

inline

Prints trajectory content into a stream.

Template Parameters

Stream

type of the output stream

Parameters

out	stream to output the information into
verbosity	verbosity level (default: 1)
indent	indentation string (default: none)

See also

Dump(Stream&&, unsigned int, std::string, std::string)

Implementation detail for Dump(Stream&&, unsigned int, std::string).

Definition at line 601 of file TrackTrajectory.h.

References Dump(), and LowLevelDump().

                                                                    {})
      const
      { Dump(std::forward<Stream>(out), verbosity, indent, indent); }

Point_t const& recob::TrackTrajectory::End ( ) const

inline

Returns the position of the last valid point of the trajectory [cm].

Definition at line 260 of file TrackTrajectory.h.

References evd::details::end(), Extent(), LastValidPoint(), and recob::Trajectory::LocationAtPoint().

Referenced by recob::Track::End(), Extent(), trkf::TrackKalmanFitter::fitTrack(), and trkf::KalmanFilterFinalTrackFitter::setDirFlip().

      { return LocationAtPoint(LastValidPoint()); }

Vector_t recob::TrackTrajectory::EndDirection ( ) const

inline

Returns the direction of the trajectory at the last point.

Definition at line 332 of file TrackTrajectory.h.

References recob::Trajectory::DirectionAtPoint(), and LastValidPoint().

Referenced by Direction(), recob::Track::EndDirection(), and trkf::TrackKalmanFitter::fitTrack().

      { return DirectionAtPoint(LastValidPoint()); }

double recob::TrackTrajectory::EndMomentum ( ) const

inline

Computes and returns the modulus of momentum at the last point [GeV/c].

See also

EndMomentumVector()

Definition at line 493 of file TrackTrajectory.h.

References Direction(), recob::Trajectory::DirectionAtPoint(), EndMomentumVector(), recob::Trajectory::HasMomentum(), recob::Trajectory::MomentumAtPoint(), and recob::Trajectory::MomentumVectorAtPoint().

Referenced by recob::Track::EndMomentum().

      { return EndMomentumVector().R(); }

Vector_t const& recob::TrackTrajectory::EndMomentumVector ( ) const

inline

Returns the momentum of the trajectory at the last valid point [GeV/c].

Definition at line 477 of file TrackTrajectory.h.

References LastValidPoint(), and recob::Trajectory::MomentumVectorAtPoint().

Referenced by EndMomentum(), and recob::Track::EndMomentumVector().

      { return MomentumVectorAtPoint(LastValidPoint()); }

void recob::TrackTrajectory::Extent	(	std::vector< double > &	start,
		std::vector< double > &	end
	)		const

Fills the first and last valid point in the trajectory.

Parameters

start	(output) position of the beginning of the trajectory
end	(output) position of the end of the trajectory

The labelling of start and end is consistent within the trajectory but is not guaranteed to be physically correct.

Deprecated:

Use physics vectors instead (see Extent())

Definition at line 58 of file TrackTrajectory.cxx.

References End(), recob::details::legacy::FillTwoVectors(), and Start().

Referenced by recob::Track::MomentumVectorAtPoint().

  { details::legacy::FillTwoVectors(Start(), End(), start, end); }

std::pair<Point_t, Point_t> recob::TrackTrajectory::Extent ( ) const

inline

Returns a copy of the first and last valid point in the trajectory.

Returns

a pair: the first and last point in the trajectory

The labelling of start and end is consistent within the trajectory but is not guaranteed to be physically correct.

Example:

recob::TrackTrajectory::Point_t start, end;

std::tie(start, end) = traj.Extent(); // assign both start and end

Definition at line 295 of file TrackTrajectory.h.

References End(), Length(), and Start().

Referenced by CountValidPoints(), and End().

      { return { Start(), End() }; }

size_t recob::Trajectory::FirstPoint ( ) const

inlineinherited

Returns the index of the first point in the trajectory (yep, it's 0).

Definition at line 176 of file Trajectory.h.

Referenced by recob::Track::FirstPoint(), recob::Trajectory::Start(), recob::Trajectory::StartDirection(), recob::Trajectory::StartMomentumVector(), and Trajectory().

      { return 0U; }

size_t recob::TrackTrajectory::FirstValidPoint ( ) const

inline

Returns the index of the first valid point in the trajectory.

Returns

index of the first point in the trajectory, or InvalidIndex

Returns the index of the first point with the flag NoPoint unset. It never returns InvalidIndex unless the track trajectory is invalid.

Definition at line 199 of file TrackTrajectory.h.

References NextValidPoint().

Referenced by AzimuthAngle(), trkf::TrajectoryMCSFitter::breakTrajInSegments(), recob::Track::FirstValidPoint(), Phi(), Start(), StartDirection(), StartMomentumVector(), Theta(), and ZenithAngle().

      { return NextValidPoint(0U); }

Flags_t const& recob::TrackTrajectory::Flags ( ) const

inline

Returns all flags.

Definition at line 173 of file TrackTrajectory.h.

References fFlags.

Referenced by trkf::TrackKalmanFitter::fitTrack().

      { return fFlags; }

PointFlags_t const& recob::TrackTrajectory::FlagsAtPoint ( size_t  i ) const

inline

Returns the flags for the specified trajectory point.

Parameters

i	index of the point in the trajectory

Returns

flags for the specified trajectory point [cm]

If the point index is invalid, the result is undefined.

Definition at line 167 of file TrackTrajectory.h.

References fFlags.

Referenced by recob::Track::FlagsAtPoint(), and HasValidPoint().

      { return fFlags[i]; }

recob::Trajectory::Rotation_t recob::Trajectory::GlobalToLocalRotationAtPoint ( size_t  p ) const

inherited

Returns a rotation matrix that brings trajectory direction along z.

Parameters

p	index of the trajectory point where to apply the rotation

Returns

a rotation matrix suitable to point the trajectory along z

The returned rotation matrix, applied to the direction vector of the trajectory at point p, will make that direction point toward the positive z axis direction, in a sort of "first person view" of the trajectory at that point. If p does not denote a valid trajectory point, the result is undefined.

The return value can be used on a Vector_t to rotate it. For example:

auto rot = traj.GlobalToLocalRotationAtPoint(0);
auto local = rot * traj.DirectionAtPoint(0);

local will be a Vector_t object (like the direction at the right side of the product) that points to (0, 0, 1).

While the transformation that yields the rotation matrix is continuous, the direction of the new "local" x and y axes is not defined.

Definition at line 173 of file Trajectory.cxx.

References recob::Trajectory::DirectionAtPoint(), and recob::tracking::Plane::Global3DToLocal3DRotation().

Referenced by recob::Trajectory::Direction(), Direction(), recob::Track::GlobalToLocalRotationAtPoint(), and recob::Trajectory::GlobalToLocalRotationAtPoint().

{
  return recob::tracking::Plane::Global3DToLocal3DRotation(DirectionAtPoint(p));
} // recob::Trajectory::GlobalToLocalRotationAtPoint()

void recob::Trajectory::GlobalToLocalRotationAtPoint ( unsigned int  p, TMatrixD &  rot  ) const

inherited

Fills a rotation matrix that brings trajectory direction along z.

Parameters

p	index of the trajectory point where to apply the rotation
rot	(output) the rotation matrix to be filled

Deprecated:

Use GlobalToLocalRotationAtPoint(size_t) instead.

Definition at line 181 of file Trajectory.cxx.

References recob::Trajectory::GlobalToLocalRotationAtPoint(), and recob::Trajectory::LocalToGlobalRotationAtPoint().

{
  auto trackRot = GlobalToLocalRotationAtPoint(p);
  fillTMatrixD(trackRot, rot);
} // recob::Trajectory::GlobalToLocalRotationAtPoint(TMatrixD)

bool recob::Trajectory::HasMomentum ( ) const

inlineinherited

Returns whether information about the momentum is available.

Returns

whether information about the momentum is available

The trajectory may or may not store valid momentum information. If not, the methods returning momentum information will stick to a modulus 1 GeV/c and the momentum numerically matches the direction.

Definition at line 460 of file Trajectory.h.

References recob::Trajectory::fHasMomentum.

Referenced by recob::Trajectory::DirectionAtPoint(), EndMomentum(), recob::Track::FlagsAtPoint(), NumberFitMomentum(), recob::Trajectory::NumberFitMomentum(), and trkmkr::KalmanFilterFitTrackMaker::restoreInputPoints().

      { return fHasMomentum; }

bool recob::Trajectory::HasPoint ( size_t  i ) const

inlineinherited

Returns whether the specified trajectory point is available.

Parameters

i	index of the trajectory point

Returns

whether the specified trajectory point is available.

Definition at line 188 of file Trajectory.h.

References recob::Trajectory::NPoints().

Referenced by recob::Track::HasPoint(), HasValidPoint(), and Trajectory().

      { return i < NPoints(); }

bool recob::TrackTrajectory::HasValidPoint ( size_t  i ) const

inline

Returns whether the specified point has NoPoint flag unset.

Returns

whether the specified point has NoPoint flag unset

A point with flag NoPoint set is actually an invalid point, that the algorithm could not at all set, but it has still a hit associated with it.

If the point index is invalid, false is returned.

Definition at line 186 of file TrackTrajectory.h.

References FlagsAtPoint(), recob::Trajectory::HasPoint(), recob::TrajectoryPointFlags::isSet(), recob::TrajectoryPointFlagTraits::NoPoint, and Trajectory().

Referenced by AtLeastValidTrajectoryPoints(), CountValidPoints(), and recob::Track::HasValidPoint().

      {
        return Trajectory().HasPoint(i)
                    && !FlagsAtPoint(i).isSet(flag::NoPoint);
      }

size_t recob::Trajectory::LastPoint ( ) const

inlineinherited

Returns the index of the last point in the trajectory.

Definition at line 180 of file Trajectory.h.

References recob::Trajectory::NPoints().

Referenced by recob::Trajectory::End(), recob::Trajectory::EndDirection(), recob::Trajectory::EndMomentumVector(), recob::Track::LastPoint(), LastValidPoint(), recob::Trajectory::Length(), Length(), and Trajectory().

      { return NPoints() - 1; }

size_t recob::TrackTrajectory::LastValidPoint ( ) const

inline

Returns the index of the last valid point in the trajectory.

Returns

index of the last point in the trajectory, or InvalidIndex

Returns the index of the last point with the flag NoPoint unset. It never returns InvalidIndex unless the track trajectory is invalid.

Definition at line 237 of file TrackTrajectory.h.

References CountValidPoints(), recob::Trajectory::LastPoint(), PreviousValidPoint(), and recob::Trajectory::TrajectoryPoint().

Referenced by trkf::TrajectoryMCSFitter::breakTrajInSegments(), End(), EndDirection(), EndMomentumVector(), recob::Track::LastValidPoint(), and Length().

      { return PreviousValidPoint(LastPoint()); }

double recob::TrackTrajectory::Length

(

size_t

startAt = 0

)

const

Returns the approximate length of the trajectory.

Parameters

startAt

(default: 0, from beginning) point to start from

Returns

the approximate length of the trajectory [cm]

The residual length from the trajectory point startAt to the end of the trajectory is computed and returned. By default, the whole trajectory length is returned. All invalid points are skipped. If after skipping, less than two valid points are left, 0 is returned.

The length approximation is just the sum of Euclidean distances between each valid trajectory point and the next (starting from the first valid one with index startAt or larger).

This operation is slow, and the result should be stored in a variable.

Parameters

startAt

(default: 0, from beginning) point to start from

Returns

the approximate length of the trajectory [cm]

The residual length from the trajectory point startAt to the end of the trajectory is computed and returned. By default, the whole trajectory length is returned. If a non-existing point is specified, 0 is returned.

The length approximation is just the sum of Euclidean distances between all consecutive trajectory points (starting from the one with index startAt).

This operation is slow, and the result should be stored in a variable.

Definition at line 79 of file TrackTrajectory.cxx.

References recob::Trajectory::LastPoint(), LastValidPoint(), recob::Trajectory::LocationAtPoint(), and R.

Referenced by trkf::TrajectoryMCSFitter::breakTrajInSegments(), recob::Track::Covariances(), Extent(), trkmkr::KalmanFilterFitTrackMaker::setMomentum(), trkf::KalmanFilterTrajectoryFitter::setMomValue(), and trkmkr::KalmanFilterFitTrackMaker::setParticleID().

                                                          {
  
  // sanity check
  if (startAt >= LastPoint()) return 0.;
  
  // just sum the distance between all locations in the trajectory
  size_t iCurr = ToValidPoint<+1>(startAt);
  size_t iNext = iCurr;
  size_t iLast = LastValidPoint();
  Point_t const* curr = &(LocationAtPoint(iCurr));
  Coord_t length = 0.0;
  while ((iNext = ToValidPoint<+1>(++iNext)) <= iLast) {
    Point_t const* next = &LocationAtPoint(iNext);
    length += (*next - *curr).R();
    curr = next;
  } // while
  return length;
} // recob::TrackTrajectory::Length()

recob::Trajectory::Rotation_t recob::Trajectory::LocalToGlobalRotationAtPoint ( size_t  p ) const

inherited

Returns a rotation matrix bringing relative directions to global.

Parameters

p	index of the trajectory point where to apply the rotation

Returns

a rotation matrix suitable to convert directions to global

The returned rotation matrix, applied to the unit vector (0, 0, 1) (local z axis direction), will turn it into the trajectory direction at point p. If p does not denote a valid trajectory point, the result is undefined.

While the transformation that yields the rotation matrix is continuous, the conversion of the directions orthogonal to the local z is not defined.

Definition at line 190 of file Trajectory.cxx.

References recob::Trajectory::DirectionAtPoint(), and recob::tracking::Plane::Local3DToGlobal3DRotation().

Referenced by recob::Trajectory::Direction(), Direction(), recob::Trajectory::GlobalToLocalRotationAtPoint(), recob::Track::LocalToGlobalRotationAtPoint(), and recob::Trajectory::LocalToGlobalRotationAtPoint().

{
  return recob::tracking::Plane::Local3DToGlobal3DRotation(DirectionAtPoint(p));
} // recob::Trajectory::GlobalToLocalRotationAtPoint()

void recob::Trajectory::LocalToGlobalRotationAtPoint ( unsigned int  p, TMatrixD &  rot  ) const

inherited

Fills a rotation matrix bringing relative directions to global.

Parameters

p	index of the trajectory point where to apply the rotation
rot	(output) the rotation matrix to be filled

Deprecated:

Use LocalToGlobalRotationAtPoint(size_t) instead.

Definition at line 198 of file Trajectory.cxx.

References recob::Trajectory::LocalToGlobalRotationAtPoint().

{
  auto trackRot = LocalToGlobalRotationAtPoint(p);
  fillTMatrixD(trackRot, rot);
} // recob::Trajectory::LocalToGlobalRotationAtPoint(TMatrixD)

Point_t const& recob::Trajectory::LocationAtPoint ( size_t  i ) const

inlineinherited

Returns the position at the specified trajectory point.

Parameters

i	index of the point in the trajectory

Returns

position at the specified trajectory point [cm]

If the point index is invalid, the result is undefined.

Definition at line 255 of file Trajectory.h.

References evd::details::end(), recob::Trajectory::Extent(), and recob::Trajectory::fPositions.

Referenced by trkf::TrajectoryMCSFitter::breakTrajInSegments(), recob::Trajectory::End(), End(), trkf::BezierTrack::FillSeedVector(), recob::Trajectory::Length(), Length(), trkf::TrajectoryMCSFitter::linearRegression(), recob::Track::LocationAtPoint(), TrackProducerFromTrack::produce(), TrackProducerFromTrackTrajectory::produce(), TrackProducerFromPFParticle::produce(), trkf::KalmanFilterTrajectoryFitter::restoreInputPoints(), trkf::KalmanFilterFinalTrackFitter::restoreInputPoints(), trkmkr::KalmanFilterFitTrackMaker::restoreInputPoints(), recob::Trajectory::Start(), Start(), and recob::Trajectory::TrajectoryPoint().

      { return fPositions[i]; }

template<typename Stream >

void recob::TrackTrajectory::LowLevelDump	(	Stream &&	out,
		std::string	indent,
		std::string	indentFirst
	)		const

Prints low-level trajectory content into a stream.

Template Parameters

Stream

type of the output stream

Parameters

out	stream to output the information into
indent	indentation string (default: none)
indentFirst	indentation for first output line (default: as indent)

Referenced by Dump().

const Momenta_t& recob::Trajectory::Momenta ( ) const

inlineinherited

Returns reference to stored vector of momenta.

Returns

reference to stored vector of momenta

Definition at line 201 of file Trajectory.h.

References dir, recob::Trajectory::fMomenta, and recob::Trajectory::TrajectoryAtPoint().

{ return fMomenta; }

double recob::Trajectory::MomentumAtPoint ( size_t  i ) const

inlineinherited

Computes and returns the modulus of the momentum at a point.

Parameters

i	index of the point in the trajectory

Returns

modulus of the momentum at that point [GeV/c]

See also

HasMomentum()

The modulus of the momentum at the specified trajectory point is computed and returned. If the trajectory does not have momentum information, the value 1 GeV/c is always returned. This can be tested trajectory by trajectory by HasMomentum(). If the index is not valid in the trajectory, the result is undefined.

Definition at line 477 of file Trajectory.h.

References recob::Trajectory::MomentumVectorAtPoint().

Referenced by EndMomentum(), and recob::Track::MomentumAtPoint().

      { return MomentumVectorAtPoint(i).R(); }

Vector_t const& recob::Trajectory::MomentumVectorAtPoint ( size_t  i ) const

inlineinherited

Returns the momentum vector at a point.

Parameters

i	index of the point in the trajectory

Returns

the momentum at that point [GeV/c]

See also

HasMomentum(), MomentumAtPoint(), DirectionAtPoint()

The momentum at the specified trajectory point is returned. If the trajectory does not have momentum information, the returned value will represent the direction, that is a momentum with modulus 1 GeV/c. This can be tested trajectory by trajectory by HasMomentum(). If the index is not valid in the trajectory, the result is undefined.

Definition at line 492 of file Trajectory.h.

References recob::Trajectory::Direction(), and recob::Trajectory::fMomenta.

Referenced by recob::Trajectory::DirectionAtPoint(), EndMomentum(), recob::Trajectory::EndMomentumVector(), EndMomentumVector(), recob::Trajectory::MomentumAtPoint(), recob::Track::MomentumVectorAtPoint(), recob::Trajectory::Phi(), trkmkr::KalmanFilterFitTrackMaker::restoreInputPoints(), recob::Trajectory::StartMomentumVector(), StartMomentumVector(), recob::Trajectory::Theta(), and recob::Trajectory::TrajectoryPoint().

      { return fMomenta[i]; }

size_t recob::TrackTrajectory::NextValidPoint ( size_t  index ) const

inline

Returns the index of the next valid point in the trajectory.

Parameters

index

starting index

Returns

index of next valid point in the trajectory, or InvalidIndex

Returns the index of the first point with the flag NoPoint unset, starting with the point with the specified index (included), and moving forward toward the end of the trajectory. It returns InvalidIndex if point at index is invalid and there are no valid points left after it.

Definition at line 213 of file TrackTrajectory.h.

Referenced by trkf::TrajectoryMCSFitter::breakTrajInSegments(), FirstValidPoint(), trkf::TrajectoryMCSFitter::linearRegression(), and recob::Track::NextValidPoint().

      { return ToValidPoint<+1>(index); }

size_t recob::Trajectory::NPoints ( ) const

inlineinherited

Returns the number of stored trajectory points.

Returns

the number of stored trajectory points

See also

TrajectoryAtPoint(), PositionAtPoint(), MomentumVectorAtPoint(), DirectionAtPoint()

For each point, both position and momentum are available.

Definition at line 172 of file Trajectory.h.

References recob::Trajectory::fPositions.

Referenced by AtLeastValidTrajectoryPoints(), CountValidPoints(), trkf::TrajectoryMCSFitter::fitMcs(), recob::Trajectory::HasPoint(), trkmkr::TrackMaker::initEvent(), recob::Trajectory::LastPoint(), trkmkr::TrackMaker::makeTrack(), recob::Track::NPoints(), NumberFitMomentum(), recob::Trajectory::NumberFitMomentum(), recob::Trajectory::NumberTrajectoryPoints(), TrackTrajectory(), and Trajectory().

      { return fPositions.size(); }

size_t recob::TrackTrajectory::NumberFitMomentum ( ) const

inline

Definition at line 320 of file TrackTrajectory.h.

References recob::Trajectory::HasMomentum(), and recob::Trajectory::NPoints().

      { return HasMomentum()? NPoints(): 0U; }

size_t recob::Trajectory::NumberTrajectoryPoints ( ) const

inlineinherited

Returns the number of stored trajectory points.

Returns

the number of stored trajectory points

See also

TrajectoryAtPoint(), PositionAtPoint(), MomentumVectorAtPoint(), DirectionAtPoint()

For each point, both position and momentum are available.

Definition at line 161 of file Trajectory.h.

References recob::Trajectory::NPoints().

Referenced by trkf::BezierTrack::NSegments(), recob::Track::Trajectory(), and Trajectory().

      { return NPoints(); }

double recob::TrackTrajectory::Phi ( size_t  p ) const

inline

Azimuthal angle at a point on the trajectory, with respect to z.

Parameters

p	the index point to extract the angle from (no default!)

Returns

the azimuthal angle, in [radians]

See also

Phi(), Theta(size_t), ZenithAngle(size_t)

The angle is returned, as defined in recob::Trajectory::Phi(), for point with the specified index p.

If the point is invalid, the behaviour is undefined.

Note

This function has no default value for p; if p is not specified at all, the method Phi() is called instead.

Definition at line 382 of file TrackTrajectory.h.

References recob::Trajectory::Phi(), and Trajectory().

Referenced by recob::Track::Phi().

      { return Trajectory().Phi(p); }

double recob::TrackTrajectory::Phi ( ) const

inline

Azimuthal angle at a first valid point, with respect to z.

Returns

angle with respect to positive z, in [radians]

See also

Phi(size_t), Theta(), ZenithAngle()

The angle is returned, as defined in recob::Trajectory::Phi(), for the first valid point in the trajectory.

Note

This is not equivalent to Phi(0), but instead to Phi(FirstValidPoint()).

Definition at line 396 of file TrackTrajectory.h.

References FirstValidPoint().

      { return Phi(FirstValidPoint()); }

const Positions_t& recob::Trajectory::Positions ( ) const

inlineinherited

Returns reference to stored vector of positions.

Returns

reference to stored vector of positions

Definition at line 195 of file Trajectory.h.

References recob::Trajectory::fPositions.

{ return fPositions; }

size_t recob::TrackTrajectory::PreviousValidPoint ( size_t  index ) const

inline

Returns the index of the previous valid point in the trajectory.

Parameters

index

starting index

Returns

index of previous valid point in trajectory, or InvalidIndex

Returns the index of the first point with the flag NoPoint unset, starting with the point with the specified index (included), and moving backward toward the start of the trajectory. It returns InvalidIndex if point at index is invalid and there are no valid points before it.

Definition at line 227 of file TrackTrajectory.h.

References ToValidPoint().

Referenced by LastValidPoint(), and recob::Track::PreviousValidPoint().

      { return ToValidPoint<-1>(index); }

Point_t const& recob::TrackTrajectory::Start ( ) const

inline

Returns the position of the first valid point of the trajectory [cm].

Definition at line 256 of file TrackTrajectory.h.

References FirstValidPoint(), and recob::Trajectory::LocationAtPoint().

Referenced by trkf::Geometric3DVertexFitter::closestPointAlongTrack(), Extent(), trkf::Geometric3DVertexFitter::fitTwoTracks(), recob::Track::FlagsAtPoint(), trkf::Geometric3DVertexFitter::getParsCovsOnPlane(), trkf::Geometric3DVertexFitter::pDist(), and Vertex().

      { return LocationAtPoint(FirstValidPoint()); }

Vector_t recob::TrackTrajectory::StartDirection ( ) const

inline

Returns the direction of the trajectory at the first point.

Definition at line 328 of file TrackTrajectory.h.

References recob::Trajectory::DirectionAtPoint(), and FirstValidPoint().

Referenced by trkf::Geometric3DVertexFitter::closestPointAlongTrack(), Direction(), trkf::Geometric3DVertexFitter::fitTwoTracks(), recob::Track::FlagsAtPoint(), trkf::Geometric3DVertexFitter::getParsCovsOnPlane(), trkf::Geometric3DVertexFitter::pDist(), and VertexDirection().

      { return DirectionAtPoint(FirstValidPoint()); }

double recob::TrackTrajectory::StartMomentum ( ) const

inline

Computes and returns the modulus of momentum at the first point [GeV/c].

See also

StartMomentumVector()

Definition at line 488 of file TrackTrajectory.h.

References StartMomentumVector().

Referenced by recob::Track::StartMomentum(), and VertexMomentum().

      { return StartMomentumVector().R(); }

Vector_t const& recob::TrackTrajectory::StartMomentumVector ( ) const

inline

Returns the momentum of the trajectory at the first valid point [GeV/c].

Definition at line 473 of file TrackTrajectory.h.

References FirstValidPoint(), and recob::Trajectory::MomentumVectorAtPoint().

Referenced by StartMomentum(), recob::Track::StartMomentumVector(), and VertexMomentumVector().

      { return MomentumVectorAtPoint(FirstValidPoint()); }

double recob::TrackTrajectory::Theta ( size_t  p ) const

inline

Trajectory angle at point, with respect to positive z direction.

Parameters

p	the index point to extract the angle from (no default!)

Returns

angle with respect to positive z, in [radians]

The angle is returned, as defined in recob::Trajectory::Theta(), for point with the specified index p.

If the point is invalid, the behaviour is undefined.

Note

This function has no default value for p; if p is not specified at all, the method Theta() is called instead.

Definition at line 350 of file TrackTrajectory.h.

References recob::Trajectory::Theta(), and Trajectory().

Referenced by recob::Track::ParticleId(), and recob::Track::Theta().

      { return Trajectory().Theta(p); }

double recob::TrackTrajectory::Theta ( ) const

inline

Trajectory angle at start, with respect to positive z direction.

Returns

angle with respect to positive z, in [radians]

The angle is returned, as defined in recob::Trajectory::Theta(), for the first valid point in the trajectory.

Note

This is not equivalent to Theta(0), but instead to Theta(FirstValidPoint()).

Definition at line 364 of file TrackTrajectory.h.

References FirstValidPoint().

      { return Theta(FirstValidPoint()); }

template<int Dir>

size_t recob::TrackTrajectory::ToValidPoint ( size_t  index ) const

private

Returns the index of the first valid point from index on.

Template Parameters

Dir	the direction to move when an index has an invalid point

Parameters

index

the starting index

Returns

index of the first valid point from index on (or InvalidIndex)

The valid direction Dir values are only +1 and -1. The first point considered is always the one at index. If no valid point is found, InvalidIndex is returned. The invariant guarantees that all these calls return a valid index:

if (!(
     (ToValidPoint<+1>(0)                 != InvalidIndex) // FirstValidPoint()
  && (ToValidPoint<-1>(LastPoint())       != InvalidIndex) // LastValidPoint()
  && (ToValidPoint<+1>(FirstValidPoint()) != InvalidIndex) // LastValidPoint()
  && (ToValidPoint<-1>(LastValidPoint())  != InvalidIndex) // FirstValidPoint()
  )) throw std::logic_error("Invalid TrackTrajectory!");

Referenced by PreviousValidPoint().

Trajectory_t const& recob::TrackTrajectory::Trajectory ( ) const

inline

Returns the plain trajectory of this object.

Definition at line 144 of file TrackTrajectory.h.

References recob::Trajectory::FirstPoint(), recob::Trajectory::HasPoint(), recob::Trajectory::LastPoint(), recob::Trajectory::NPoints(), recob::Trajectory::NumberTrajectoryPoints(), and recob::Trajectory::TrajectoryAtPoint().

Referenced by AzimuthAngle(), HasValidPoint(), Phi(), trkf::KalmanFilterFinalTrackFitter::produce(), Theta(), and ZenithAngle().

      { return static_cast<Trajectory_t const&>(*this); }

bool recob::Trajectory::TrajectoryAtPoint ( size_t  i, TVector3 &  pos, TVector3 &  dir  ) const

inherited

Fills position and direction at the specified trajectory point.

Parameters

i	index of the trajectory point
pos	(output) filled with the position at the given point [cm]
dir	(output) filled with the direction at the given point

Returns

false if the specified point was invalid

Deprecated:

Use TrajectoryPoint() instead (and note that it returns momentum, not direction!)

Definition at line 66 of file Trajectory.cxx.

References recob::Trajectory::Extent().

Referenced by recob::Trajectory::Momenta(), recob::Trajectory::Trajectory(), Trajectory(), and recob::Track::TrajectoryAtPoint().

{
  if (!HasPoint(i)) return false;
  
  auto const& origPos = LocationAtPoint(i);
  decltype(auto) origDir = DirectionAtPoint(i);
  
  pos.SetXYZ(origPos.X(), origPos.Y(), origPos.Z());
  dir.SetXYZ(origDir.X(), origDir.Y(), origDir.Z());
  return true;
} // recob::Trajectory::TrajectoryAtPoint(TVector&)

TrajectoryPoint_t recob::Trajectory::TrajectoryPoint ( size_t  i ) const

inlineinherited

Returns position and momentum at the specified trajectory point.

Parameters

i	index of the trajectory point

Returns

the information at the specified trajectory point

Note that this method returns the momentum, not the direction.

auto trajPoint = traj.TrajectoryPoint(0);
std::cout << "Start of trajectory at " << trajPoint.position
  << " cm, with momentum " << trajPoint.momentum << " GeV/c"
  << std::endl;

If the specified index is not valid, result is undefined.

Definition at line 231 of file Trajectory.h.

References recob::Trajectory::LocationAtPoint(), and recob::Trajectory::MomentumVectorAtPoint().

Referenced by recob::Track::CountValidPoints(), and LastValidPoint().

      { return { LocationAtPoint(i), MomentumVectorAtPoint(i) }; }

Point_t const& recob::TrackTrajectory::Vertex ( ) const

inline

Returns the position of the first valid point of the trajectory [cm].

Definition at line 252 of file TrackTrajectory.h.

References Start().

Referenced by trkf::TrackKalmanFitter::fitTrack(), trkf::KalmanFilterFinalTrackFitter::setDirFlip(), and recob::Track::Vertex().

      { return Start(); }

Vector_t recob::TrackTrajectory::VertexDirection ( ) const

inline

Returns the direction of the trajectory at the first point.

Definition at line 324 of file TrackTrajectory.h.

References StartDirection().

Referenced by trkf::TrackKalmanFitter::fitTrack(), trkmkr::KalmanFilterFitTrackMaker::isFlipDirection(), trkf::KalmanFilterTrajectoryFitter::setDirFlip(), and recob::Track::VertexDirection().

      { return StartDirection(); }

double recob::TrackTrajectory::VertexMomentum ( ) const

inline

Computes and returns the modulus of momentum at the first point [GeV/c].

See also

StartMomentum()

Definition at line 483 of file TrackTrajectory.h.

References StartMomentum().

Referenced by recob::Track::VertexMomentum().

      { return StartMomentum(); }

Vector_t const& recob::TrackTrajectory::VertexMomentumVector ( ) const

inline

Returns the momentum of the trajectory at the first valid point [GeV/c].

Definition at line 469 of file TrackTrajectory.h.

References StartMomentumVector().

Referenced by recob::Track::VertexMomentumVector().

      { return StartMomentumVector(); }

double recob::TrackTrajectory::ZenithAngle ( size_t  p ) const

inline

"Zenith" angle of trajectory, with respect to the vertical axis.

Parameters

p	the index point to extract the angle from (no default!)

Returns

opposite of the actual zenith angle, in [radians]

See also

AzimuthAngle(size_t)

The zenith is returned, as defined in recob::Trajectory::Zenith(), for point with the specified index p.

If the point is invalid, the behaviour is undefined.

Note

This function has no default value for p; if p is not specified at all, the method Zenith() is called instead.

Definition at line 415 of file TrackTrajectory.h.

References Trajectory(), and recob::Trajectory::ZenithAngle().

Referenced by recob::Track::ZenithAngle().

      { return Trajectory().ZenithAngle(p); }

double recob::TrackTrajectory::ZenithAngle ( ) const

inline

"Zenith" angle of trajectory, with respect to the vertical axis.

Returns

opposite of the actual zenith angle, in [radians]

See also

Zenith(size_t), Theta()

The zenith angle is returned, as defined in recob::Trajectory::Zenith(), for the first valid point in the trajectory.

Note

This is not equivalent to Zenith(0), but instead to Zenith(FirstValidPoint()).

Definition at line 430 of file TrackTrajectory.h.

References FirstValidPoint().

      { return ZenithAngle(FirstValidPoint()); }

Member Data Documentation

Flags_t recob::TrackTrajectory::fFlags

private

Flags of each of the points in trajectory.

Definition at line 625 of file TrackTrajectory.h.

Referenced by Flags(), FlagsAtPoint(), and TrackTrajectory().

constexpr size_t recob::TrackTrajectory::InvalidIndex = std::numeric_limits<size_t>::max()

static

Value returned on failed index queries.

Definition at line 621 of file TrackTrajectory.h.

Referenced by trkf::TrajectoryMCSFitter::breakTrajInSegments().

constexpr unsigned int recob::TrackTrajectory::MaxDumpVerbosity = 7

static

Largest verbosity level supported by Dump().

Definition at line 618 of file TrackTrajectory.h.

---

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

• lardataobj/v1_36_00/source/lardataobj/RecoBase/TrackTrajectory.h
• lardataobj/v1_36_00/source/lardataobj/RecoBase/TrackTrajectory.cxx