/* $Id$ */
+//
+// Track Reference object is created every time particle is
+// crossing detector bounds.
+// The object is created by Step Manager
+//
+
#include "TObject.h"
-#include "TVirtualMC.h"
+#include "TMath.h"
+
+class AliExternalTrackParam;
class AliTrackReference : public TObject {
public:
AliTrackReference();
- AliTrackReference(Int_t label, TVirtualMC *vMC);
+ AliTrackReference(Int_t label);
virtual ~AliTrackReference() {}
+ static AliExternalTrackParam * MakeTrack(const AliTrackReference *ref, Double_t mass);
virtual Int_t GetTrack() const {return fTrack;}
virtual void SetTrack(Int_t track) {fTrack=track;}
virtual void SetLength(Float_t length){fLength=length;}
virtual void SetTime(Float_t time) {fTime = time;}
- virtual Float_t GetLength(){return fLength;}
- virtual Float_t GetTime(){return fTime;}
+ virtual Float_t GetLength() const {return fLength;}
+ virtual Float_t GetTime() const {return fTime;}
+ virtual Int_t Label() const {return fTrack;}
+ virtual void SetLabel(Int_t track) {fTrack=track;}
+ virtual Float_t R() const {return TMath::Sqrt(fX*fX+fY*fY);}
+ virtual Float_t Pt() const {return TMath::Sqrt(fPx*fPx+fPy*fPy);}
+ virtual Float_t Phi() const {return TMath::Pi()+TMath::ATan2(-fPy,-fPx);}
+ virtual Float_t Theta() const {return (fPz==0)?TMath::Pi()/2:TMath::ACos(fPz/P());}
virtual Float_t X() const {return fX;}
virtual Float_t Y() const {return fY;}
virtual Float_t Z() const {return fZ;}
virtual Float_t Px() const {return fPx;}
virtual Float_t Py() const {return fPy;}
virtual Float_t Pz() const {return fPz;}
+ virtual Float_t P() const {return TMath::Sqrt(fPx*fPx+fPy*fPy+fPz*fPz);}
virtual void SetPosition(Float_t x, Float_t y, Float_t z){fX=x; fY=y; fZ=z;}
virtual void SetMomentum(Float_t px, Float_t py, Float_t pz){fPx=px; fPy=py; fPz=pz;}
+ // Methods to get position of the track reference in
+ // in the TPC/TRD/TOF Tracking coordinate system
+
+ virtual Float_t PhiPos() const {return TMath::Pi()+TMath::ATan2(-fY, -fX);}
+ virtual Float_t Alpha() const
+ {return TMath::Pi()*(20*((((Int_t)(PhiPos()*180/TMath::Pi()))/20))+10)/180.;}
+ virtual Float_t LocalX() const {return fX*TMath::Cos(-Alpha()) - fY*TMath::Sin(-Alpha());}
+ virtual Float_t LocalY() const {return fX*TMath::Sin(-Alpha()) + fY*TMath::Cos(-Alpha());}
+
protected:
Int_t fTrack; // Track number