// -----                      PndSttHelixHit source file               -----
// -----                                                               -----
// -------------------------------------------------------------------------

#include "PndSttHelixHit.h"

#include <iostream>
#include "TMath.h"
#include "TVector3.h"
#include "PndSttHit.h"

using namespace std;

/** Default constructor **/
PndSttHelixHit::PndSttHelixHit() 
{
    Clear();
}

/** Standard constructor **/
// DONT USE THIS! delete it!
PndSttHelixHit::PndSttHelixHit(Int_t detID, 
			       TVector3& pos, TVector3& dpos,
			       Int_t mcindex, Int_t hitindex, Double_t edep)
  : FairHit(detID, pos, dpos, mcindex) 
{
  fHitIndex = hitindex;
  fELoss = edep;
}

// DONT USE THIS! add halflength to make it OK
PndSttHelixHit::PndSttHelixHit(Int_t detID, 
			       TVector3& pos, TVector3& dpos, 
			       Int_t mcindex, Int_t hitindex, Double_t edep,
			       TVector3& cpos, TVector3& dcpos, 
			       Double_t isochrone,
			       Double_t isochroneError, TVector3 wireDir)			        
 : FairHit(detID, pos, dpos, mcindex) 
{

  fHitIndex = hitindex;
  fELoss = edep;

  fXcen = cpos.X();  fYcen = cpos.Y();  fZcen = cpos.Z(); 
  fDxcen = dcpos.X();   fDycen = dcpos.Y();   fDzcen = dcpos.Z(); 
  fIsochrone = isochrone;
  fIsochroneError = isochroneError;
  fWireDirection = wireDir;
}


// this copies the PndSttHit to the PndSttHelixHit
void PndSttHelixHit::CopyHitToHelixHit(PndSttHit *aHit, Int_t hitindex)
{
//   FairHit::SetPosition(TVector3(0,0,0));
//   FairHit::SetPositionError(TVector3(0,0,0));
  FairHit::SetDetectorID(aHit->GetDetectorID());
  FairHit::SetRefIndex(aHit->GetRefIndex());

  fHitIndex = hitindex;
  fELoss = aHit->GetEnergyLoss();
  fDepCharge = aHit->GetDepCharge();
  //  fdEdx = aHit->GetdEdx();
  
  TVector3 centerPosition, centerPositionError;
  aHit->Position(centerPosition);
  aHit->PositionError(centerPositionError);

  fXcen  = centerPosition.X(); fYcen = centerPosition.Y(); fZcen = centerPosition.Z(); 
  fDxcen = centerPositionError.X(); fDycen = centerPositionError.Y(); fDzcen = centerPositionError.Z(); 

  fIsochrone      = aHit->GetIsochrone();
  fIsochroneError = aHit->GetIsochroneError();
  fWireDirection  = aHit->GetWireDirection();
  fHalfLength     = aHit->GetTubeHalfLength();
}



/** Public method Clear **/
void PndSttHelixHit::Clear() 
{
  fHitIndex = 0;
  fDepCharge = 0.; 
  fdEdx = 0.;      
  fELoss = 0.;   
  fXcen = 0.;
  fYcen = 0.;
  fZcen = 0.; 
  fDxcen = 0.; 
  fDycen = 0.; 
  fDzcen = 0.;
  fWireDirection = TVector3(0., 0., 0.);
  fIsochrone = 0.; 
  fIsochroneError = 0.;
  fHalfLength = 0.; 
 
}  

 
/** Destructor **/
PndSttHelixHit::~PndSttHelixHit() 
{
} 

void PndSttHelixHit::Print()
{
  // reco & center positions
  cout << "position " << GetX() << " " << GetY() << " " << GetZ() << endl;
  cout << "center position " << GetXcen()  << " " << GetYcen() << " " << GetZcen() << endl;

  // index
  cout << "mc, hit indexes " << GetRefIndex() << " " << GetHitIndex() << endl;

  // drift
  cout << "drift radius " << GetIsochrone() << endl;

  // wire
  cout << "wire dir " << GetWireDirection().X() << " " << GetWireDirection().Y() << " " << GetWireDirection().Z() << endl;
  cout << "wire half length " << GetTubeHalfLength() << endl;
}



ClassImp(PndSttHelixHit)