MarlinTrackFit.h

#ifndef MARLINTRACKFIT_H
#define MARLINTRACKFIT_H 1
#include "HelixClass.h"

extern "C" {
  void tfithl_(int & NPT, double * XF,double * YF, float * RF, float *PF, double * WF,float * ZF,
               float * WZF, int & IOPT, float * VV0,
               float * EE0, float & CH2PH, float & CH2Z);
}

extern "C" {
  void trackfit_(int & nhits, int * idet, int * itype, 
                 float * x, float * y, float * z, float * phireso, float * zreso,
                 float * ref, int & ierr, float * rfit, float * rfite, float & chi2, int & ndf,
                 int & noutl, int * idoutl, int & fitcode);
  
}


class MarlinTrackFit {
 public:

  MarlinTrackFit();

  ~MarlinTrackFit();

  int DoFitting(int useExtraPoint, int fitOpt, // inputs
                int & nhits, float bField, int * idet, int * itype, // inputs 
                float chi2PrefitCut, // inputs 
                float * x, float * y, float * z, float * RPReso, float * ZReso, // inputs 
                float * param, float * eparam, float * RefPoint, float & chi2, int & ndf, // outputs
                float & chi2rphi, float & chi2z,int * lhits); // outputs

  void CrudeErrorEstimates(float BField, float * reso, float * param, float * eparam);
  void setParametersForIPErrors(float * par);

 private:    

  float _parIpReso[3];
    
  void CalculateChi2(HelixClass & helix, int & nhits, 
                     float * xhit, float * yhit, float * zhit,
                     int * iType, float * rphireso, float * zreso,
                     float * chi2hit, float & chi2rphi, float & chi2z);

  void ConvertLCtoTANAGRA(float * param, float & bField, float * ref);
  void ConvertTANAGRAtoLC(float * rfit, float & bField, int & fitCode, float * param);
  void ConvertCovMatrix(float * rfit, float * rfite, float * param, float & bField, int fitCode, float * eparam);

};

#endif