bj/lardataobj/RecoBase/Track.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////////
//
// \brief Definition of track object for LArSoft
//
// \author brebel@fnal.gov
//
////////////////////////////////////////////////////////////////////////////


#include "lardataobj/RecoBase/Track.h"
#include "lardataobj/RecoBase/TrackingPlane.h"

#include <iomanip>
#include <iostream>

namespace recob {

  //----------------------------------------------------------------------
  // ostream operator.
  //
  std::ostream& operator<< (std::ostream& stream, Track const& a)
  {
    auto const& start = a.VertexDirection();
    auto const& end   = a.EndDirection();
    stream << std::setiosflags(std::ios::fixed) << std::setprecision(3)
           << "\n Track ID "       << std::setw(4) << std::right << a.ID()
           << " Theta = "            << std::setw(6) << std::right << a.Theta()
           << " Phi = "              << std::setw(6) << std::right << a.Phi()
           << "\n  StartCosines : ( " << start.X() << " ; " << start.Y() << " ; " << start.Z()
           << ")  EndCosines : ( "   << end.X() << " ; " << end.Y() << " ; " << end.Z()
           << ")"
           << "\n  #Position and Direction = "
           << std::setw(5) << std::right << a.NumberTrajectoryPoints();
    stream << std::endl;

    return stream;
  }

  //----------------------------------------------------------------------------
  bool operator < (const Track & a, const Track & b)
  {
    if(a.ID() != b. ID())
      return a.ID() < b.ID();

    return false; //They are equal
  }

  //----------------------------------------------------------------------------
  Track::SVector6 Track::VertexParametersGlobal6D() const {
    Track::SVector6 result;
    result[0] = Vertex().X();
    result[1] = Vertex().Y();
    result[2] = Vertex().Z();
    result[3] = VertexDirection().X();
    result[4] = VertexDirection().Y();
    result[5] = VertexDirection().Z();
    return result;
  }

  Track::SVector6 Track::EndParametersGlobal6D() const {
    Track::SVector6 result;
    result[0] = End().X();
    result[1] = End().Y();
    result[2] = End().Z();
    result[3] = EndDirection().X();
    result[4] = EndDirection().Y();
    result[5] = EndDirection().Z();
    return result;
  }

  Track::SVector5 Track::VertexParametersLocal5D() const {
    //return tracking::Plane::Global6DToLocal5DParameters(VertexParametersGlobal6D(), fTraj.Vertex(), fTraj.VertexDirection());
    return Track::SVector5(0.,0.,0.,0.,(HasMomentum() ? 1./VertexMomentum() : 1. ));
  }

  Track::SVector5 Track::EndParametersLocal5D() const {
    //return tracking::Plane::Global6DToLocal5DParameters(EndParametersGlobal6D(), fTraj.End(), fTraj.EndDirection());
    return Track::SVector5(0.,0.,0.,0.,(HasMomentum() ? 1./EndMomentum() : 1. ));
  }

  Track::SMatrixSym66 Track::VertexCovarianceGlobal6D() const {
    return tracking::Plane::Local5DToGlobal6DCovariance(fCovVertex,HasMomentum(),fTraj.VertexMomentumVector(),fTraj.VertexDirection());
  }

  Track::SMatrixSym66 Track::EndCovarianceGlobal6D() const {
    return tracking::Plane::Local5DToGlobal6DCovariance(fCovEnd,HasMomentum(),fTraj.EndMomentumVector(),fTraj.EndDirection());
  }

}