Moving the efficiency calculations before the cleanup (Guiseppe)

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

diff --git a/STEER/AliTrackReference.h b/STEER/AliTrackReference.h
index 8b779c8e76ded325fe02671cb775d46edea20d56..198f6c4f0633ae6cbe8b3cdfaf50aca960d985f7 100644 (file)
--- a/STEER/AliTrackReference.h
+++ b/STEER/AliTrackReference.h
@@ -5,17 +5,33 @@
 
 /* $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:
 
+    enum constants {kDisappeared = -1,
+                   kITS   = 0,
+                   kTPC   = 1,
+                   kFRAME = 2,
+                   kTRD   = 3,
+                   kTOF   = 4,
+                   kMUON  = 5,
+                   kHMPID = 6
+    };
   AliTrackReference();
-  AliTrackReference(Int_t label, TVirtualMC *vMC);
+  AliTrackReference(Int_t label, Int_t id = -999);
+  AliTrackReference(const AliTrackReference &tr);
   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;}
@@ -27,28 +43,52 @@ 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 Int_t   UserId() const {return fUserId;}
+  virtual Int_t   DetectorId() const {return fDetectorId;}
+  virtual void SetDetectorId(Int_t id){fDetectorId = id;}
   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;}
+  virtual void SetUserId(Int_t userId){fUserId=userId;}
+ 
+  // 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());}
 
+  Bool_t IsSortable() const {return kTRUE;}
+  Int_t Compare(const TObject *obj) const {
+    Int_t ll = ((AliTrackReference*)obj)->GetTrack();
+    if (ll < fTrack) return 1;
+    if (ll > fTrack) return -1;
+    return 0;
+  }
 
+  virtual void Print(Option_t* opt="") const;
+    
 protected:
-  Int_t     fTrack;  // Track number
-  Float_t   fX;      // X reference position of the track
-  Float_t   fY;      // Y reference position of the track
-  Float_t   fZ;      // Z reference position of the track
-  Float_t   fPx;     // momentum
-  Float_t   fPy;     // momentum
-  Float_t   fPz;     // momentum
-  Float_t   fLength; // track lenght from its origin in cm
-  Float_t   fTime;   // time of flight in cm  
-
-  ClassDef(AliTrackReference,3)  //Base class for all Alice track references
+  Int_t     fTrack;      // Track number
+  Float_t   fX;          // X reference position of the track
+  Float_t   fY;          // Y reference position of the track
+  Float_t   fZ;          // Z reference position of the track
+  Float_t   fPx;         // momentum
+  Float_t   fPy;         // momentum
+  Float_t   fPz;         // momentum
+  Float_t   fLength;     // track lenght from its origin in cm
+  Float_t   fTime;       // time of flight in cm  
+  Int_t     fUserId;     // optional Id defined by user
+  Int_t     fDetectorId; // Detector Id
+  ClassDef(AliTrackReference,6)  //Base class for all Alice track references
 };
 #endif