//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Implementation of class TBSIRecoHit
//      see PndGemRecoHit.h for details
//
// Environment:
//
// Author List:
//      Sebastian Neubert     TUM            (original author)
//      Maxence Vandenbroucke TUM            (application to TBSICluster)
//
//-----------------------------------------------------------

// C/C++ Headers ----------------------
// root Headers ----------------------
#include "TMatrixT.h"
#include "TMath.h"
// Collaborating Class Headers --------
#include "RKTrackRep.h"
#include "GFDetPlane.h"
// This Class' Header ------------------
#include "TBSIRecoHit.h"
#include "TBSICluster.h"
#include "FairRootManager.h"
#include "TpcAlignmentManager.h"

// Class Member definitions -----------

ClassImp(TBSIRecoHit);


TBSIRecoHit::~TBSIRecoHit()
{}

TBSIRecoHit::TBSIRecoHit()
  : GFRecoHitIfc<GFPlanarHitPolicy>(fNparHitRep)
{}

TBSIRecoHit::TBSIRecoHit(TBSICluster* hit)
  : GFRecoHitIfc<GFPlanarHitPolicy>(fNparHitRep)
{ 

  Double_t hitX = hit->u();//hit->....

  fHitCoord[0][0] = hitX;
  //   fHitCoord[1][0] = hitY;
  
   fHitCov[0][0] = hit->xerr()*hit->xerr();
   //   fHitCov[1][1] = hit->yerr()*hit->yerr();
   
   TVector3 uu(0,0,0), vv(0,0,0), oo(0,0,0);
   hit->GetPlane(uu,vv,oo);
   fPolicy.setDetPlane(GFDetPlane(oo,uu,vv));
   //============================================================================
}

TMatrixT<double> TBSIRecoHit::getHMatrix(const GFAbsTrackRep* stateVector)
{
  
  if (dynamic_cast<const RKTrackRep*>(stateVector) != NULL) {
    // Uses TrackParP (q/p,v',w',v,w)
    // coordinates are defined by detplane!
    TMatrixT<double> HMatrix(1,5);
    
    HMatrix[0][0] = 0.;
    HMatrix[0][1] = 0.;
    HMatrix[0][2] = 0.;
    HMatrix[0][3] = 1.;
    HMatrix[0][4] = 0.;
    return HMatrix;
  }
  else {
    std::cerr << "TBSIRecoHit can only handle RK tracks"
	      << std::endl;
    throw;
  }
  
}

Double_t
TBSIRecoHit::residualScalar(GFAbsTrackRep* stateVector,
			     const TMatrixT<Double_t>& state)
{
  throw;
}