Additional initialization

[u/mrichter/AliRoot.git] / STEER / AliTrackReference.h

diff --git a/STEER/AliTrackReference.h b/STEER/AliTrackReference.h
index 8b779c8e76ded325fe02671cb775d46edea20d56..4ab594004f779aa5a8dd39e21935b9c852355bc7 100644 (file)
--- a/STEER/AliTrackReference.h
+++ b/STEER/AliTrackReference.h
@@ -5,15 +5,21 @@
 
 /* $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 AliTrackReference : public TObject {
 
 public:
 
   AliTrackReference();
-  AliTrackReference(Int_t label, TVirtualMC *vMC);
+  AliTrackReference(Int_t label);
   virtual ~AliTrackReference() {}
 
   virtual Int_t GetTrack() const {return fTrack;}
@@ -27,16 +33,26 @@ public:
   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 (fZ==0)?TMath::Pi()/2:TMath::ACos(fZ/TMath::Sqrt(fX*fX+fY*fY+fZ*fZ)); }
+  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