MarlinTrkUtils.h

#ifndef INCLUDE_MarlinTrkUtils
#define INCLUDE_MarlinTrkUtils 1

#include <vector>

#include <cfloat>

#include <LCIOSTLTypes.h>

namespace IMPL {
  class TrackImpl ;
  class TrackStateImpl ;
}

namespace EVENT {
  class TrackerHit;
  class TrackState;
}

namespace UTIL {
  class BitField64;
}


namespace MarlinTrk{
  class IMarlinTrack ;
}


namespace MarlinTrk{

  int createFinalisedLCIOTrack(
      IMarlinTrack* marlinTrk,
      std::vector<EVENT::TrackerHit*>& hit_list,
      IMPL::TrackImpl* track,
      bool fit_direction,
      EVENT::TrackState* pre_fit,
      float bfield_z,
      double maxChi2Increment=DBL_MAX);
  
  int createFinalisedLCIOTrack(
      IMarlinTrack* marlinTrk,
      std::vector<EVENT::TrackerHit*>& hit_list,
      IMPL::TrackImpl* track,
      bool fit_direction,
      const EVENT::FloatVec& initial_cov_for_prefit,
      float bfield_z,
      double maxChi2Increment=DBL_MAX);
  
  int createPrefit( std::vector<EVENT::TrackerHit*>& hit_list, IMPL::TrackStateImpl* pre_fit, float bfield_z, bool fit_direction );

  int createFit( std::vector<EVENT::TrackerHit*>& hit_list, IMarlinTrack* marlinTrk, EVENT::TrackState* pre_fit, float bfield_z, bool fit_direction, double maxChi2Increment=DBL_MAX );

  int finaliseLCIOTrack(
      IMarlinTrack* marlinTrk,
      IMPL::TrackImpl* track,
      std::vector<EVENT::TrackerHit*>& hit_list,
      bool fit_direction,
      IMPL::TrackStateImpl* atLastHit=0,
      IMPL::TrackStateImpl* atCaloFace=0);
  
  void addHitNumbersToTrack(IMPL::TrackImpl* track, std::vector<EVENT::TrackerHit*>& hit_list, bool hits_in_fit, UTIL::BitField64& cellID_encoder);

  void addHitNumbersToTrack(IMPL::TrackImpl* track, std::vector<std::pair<EVENT::TrackerHit* , double> >& hit_list, bool hits_in_fit, UTIL::BitField64& cellID_encoder);
  
}

#endif