#ifndef _COMBINEDRESIDFUNC_H_
#define _COMBINEDRESIDFUNC_H_

#define PIOVER2 1.5707963
#include "FDC/DFDCPseudo.h"
#include "CDC/DCDCTrackHit.h"

class combinedResidFunc:public residFunc {
 public:
  combinedResidFunc(vector<const DFDCPseudo*> *pseudopoints,
		    vector<const DCDCTrackHit*> *trackhits,
		    MyTrajectory *trajectory);
  void resid(const HepVector *x, void *data, HepVector *f);
  void deriv(const HepVector *x, void *data, HepMatrix *J);
  void residAndDeriv(const HepVector *x, void *data, HepVector *f,
		     HepMatrix *J);
  inline unsigned int getn() {return n + n_cdc;};
  inline unsigned int getp() {return trajPtr->getNumberOfParams();};
 private:
  unsigned int n, n_cdc;
  vector<const DFDCPseudo*> *ppPtr;
  vector<const DCDCTrackHit*> *trkhitPtr;
  MyTrajectory *trajPtr;
  HepVector pseudo2HepVector(const DFDCPseudo &pseudopoint);
  DLine trackhit2line(const DCDCTrackHit &trackhit);
  vector<double> delta;
  int debug_level;
};

#endif // _COMBINEDRESIDFUNC_H_

// end of C++ source