PndTrkTrack.h

Go to the documentation of this file.

#ifndef PNDTRKTRACK_H
#define PNDTRKTRACK_H 1

#include "TVector3.h"
#include "TVector2.h"
#include "PndTrkCluster.h"
#include "PndTrackCand.h"
#include "PndTrack.h"
#include "PndTrkHit.h"
#include "PndTrkTools.h"

class PndTrkTrack : public TObject
{
  
 public:    

  PndTrkTrack();
  PndTrkTrack(PndTrkCluster *cluster);
  PndTrkTrack(PndTrkCluster *cluster, double x, double y, double radius);
  PndTrkTrack(PndTrkHit *hit, PndTrkCluster *cluster, double x, double y, double radius);
  PndTrkTrack(double x, double y, double radius);
  PndTrkTrack(PndTrack *trk);
  PndTrkTrack(const PndTrkTrack& track);
 
 ~PndTrkTrack();    

  PndTrkTrack& operator=(const PndTrkTrack& track); 

  Bool_t operator==( PndTrkTrack track) ; // CHECK this needs to be changed

  void SetRadius(double radius) { fRadius = radius; }
  void SetCenter(double x, double y) { fCenterX = x; fCenterY = y; }
  void SetCenter(TVector3 center)  { fCenterX = center.X(); fCenterY = center.Y(); }
  void SetTanL(double tanl) { fTanL = tanl; }
  void SetZ0(double z0) { fZ0 = z0;}
  void SetCluster(PndTrkCluster *cluster) { fCluster = *cluster; }
  void SetRefHit(PndTrkHit *hit) { fRefHit = hit; }
  Double_t GetRadius() { return fRadius; }
  TVector2 GetCenter() { return TVector2(fCenterX, fCenterY); }
  Double_t GetTanL() { return fTanL; } 
  Double_t GetZ0() { return fZ0; }

  PndTrkCluster GetCluster() { return fCluster; }

  Double_t ComputePhi(TVector3 hit);
  Double_t ComputePhiFrom(TVector3 hit, TVector3 from); 

  TVector3 ComputeMomentumAtPosition(TVector3 position, TVector3 &newposition);
  void ComputeCharge();
  PndTrackCand ConvertToPndTrackCand(); 
  PndTrack ConvertToPndTrack();
  Int_t GetCharge() { return fCharge; }


  void Draw(Color_t color = kBlack);
  void LightUp();

  void Clear(Option_t* opt = "");

 protected:
  PndTrkCluster fCluster; 
  PndTrkHit *fRefHit; 
  double fCenterX, fCenterY, fRadius, fTanL, fZ0; 
  int fCharge; 
  double fPhiMin; 
  double fPhiMax; 
  // parameters:
  // pT = p cos(lam)
  // pL = p sin(lam)
  // ==> pL = pT * tan(lam)
  // z = z0 + s sin(lam) = z0 + tan(lam) [s cos(lam)]
  // with s cos(lam) = proj track length = -q phi_rad* R
  // ==> z = z0 - q tanl phi_rad R
  // If I have z = z0 + k phi_deg ==> k 180/pi = -q tanl R ==> tanl = -q k (180/pi) /R  


  ClassDef(PndTrkTrack,1);
};


#endif