Track.hh

#ifndef Track_hh
#define Track_hh

#include <vector>

#include "TLinearFitter.h"

#include "StHit.h"

class Track : public vector<StHit*> {
private:
  struct LessHit {
    bool operator()(const StHit* hit1, const StHit* hit2) const
    {
      return hit1->position().z() < hit2->position().z();
    }
  };

public:
  StHit* firstHit() const;

  StHit* lastHit() const;

  void merge(Track* track);

  bool fit();

  double chi2zx();

  double chi2zy();

  int ndf() const;

  bool ok() const;

  bool accept(StHit* hit) const;

  double length() const;

  double x0() const;

  double y0() const;

  double dxdz() const;

  double dydz() const;

  double x0error() const;

  double y0error() const;

  double dxdzError() const;

  double dydzError() const;

private:
  static TLinearFitter mXfitter;
  static TLinearFitter mYfitter;

  double mLength;
  double mChiSquareX;
  double mChiSquareY;
  int    mNumberFreeParameters;
  double mX0;
  double mY0;
  double mdxdz;
  double mdydz;
  double mX0error;
  double mY0error;
  double mdxdzError;
  double mdydzError;
};

inline StHit* Track::firstHit() const { return front(); }
inline StHit* Track::lastHit() const { return back(); }
inline double Track::length() const { return mLength; }
inline double Track::chi2zx() { return mChiSquareX; }
inline double Track::chi2zy() { return mChiSquareY; }
inline int    Track::ndf() const { return mNumberFreeParameters; }
inline double Track::x0() const { return mX0; }
inline double Track::y0() const { return mY0; }
inline double Track::dxdz() const { return mdxdz; }
inline double Track::dydz() const { return mdydz; }
inline double Track::x0error() const { return mX0error; }
inline double Track::y0error() const { return mY0error; }
inline double Track::dxdzError() const { return mdxdzError; }
inline double Track::dydzError() const { return mdydzError; }

ostream& operator<<(ostream&, const Track&);

#endif