]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliExternalTrackParam.h
Add missing file to last entry.
[u/mrichter/AliRoot.git] / STEER / AliExternalTrackParam.h
index 0c7bf036f659c12c40e17f0c3379cff0c622006c..f9f100a6f2ccd0a1c67921019ba1765c02c0ec92 100644 (file)
  *****************************************************************************/
 #include "TObject.h"
 
-const Double_t kB2C=0.299792458e-3;
 const Double_t kAlmost1=0.999;
 const Double_t kAlmost0=1e-33;
 const Double_t kVeryBig=1./kAlmost0;
 
+const Double_t kB2C=0.299792458e-3;
+const Double_t kAlmost0Field=1.e-13;
+const Double_t kVeryBigConvConst=1/kB2C/kAlmost0Field;
+const Double_t kMostProbableMomentum=0.35;
+
 class AliKalmanTrack;
+class AliESDVertex;
 
 class AliExternalTrackParam: public TObject {
  public:
@@ -33,6 +38,7 @@ class AliExternalTrackParam: public TObject {
   AliExternalTrackParam(Double_t x, Double_t alpha, 
                        const Double_t param[5], const Double_t covar[15]);
   AliExternalTrackParam(const AliKalmanTrack& track);
+  virtual ~AliExternalTrackParam(){}
 
   void Reset();
   void Set(const AliKalmanTrack& track);
@@ -42,9 +48,14 @@ class AliExternalTrackParam: public TObject {
   Double_t GetSigmaY2() const {return fC[0];}
   Double_t GetSigmaZ2() const {return fC[2];}
   Double_t GetX() const {return fX;}
+  Double_t GetY()    const {return fP[0];}
+  Double_t GetZ()    const {return fP[1];}
   Double_t GetAlpha() const {return fAlpha;}
   Double_t GetSign() const {return (fP[4]>0) ? 1 : -1;}
   Double_t GetP() const;
+  Double_t Get1P() const;
+  Double_t GetC(Double_t b) const {return fP[4]*b*kB2C;}
+  void GetDZ(Double_t x,Double_t y,Double_t z,Double_t b,Float_t dz[2]) const; 
   Double_t GetD(Double_t xv, Double_t yv, Double_t b) const; 
   Double_t GetLinearD(Double_t xv, Double_t yv) const; 
   Bool_t CorrectForMaterial(Double_t d, Double_t x0, Double_t mass);
@@ -62,17 +73,18 @@ class AliExternalTrackParam: public TObject {
   Double_t GetDCA(const AliExternalTrackParam *p, Double_t b,
     Double_t &xthis,Double_t &xp) const;
   Double_t PropagateToDCA(AliExternalTrackParam *p, Double_t b);
+  Bool_t PropagateToDCA(const AliESDVertex *vtx, Double_t b, Double_t maxd);
 
   Bool_t GetPxPyPz(Double_t *p) const;
   Bool_t GetXYZ(Double_t *p) const;
   Bool_t GetCovarianceXYZPxPyPz(Double_t cv[21]) const;
   Bool_t GetPxPyPzAt(Double_t x, Double_t b, Double_t p[3]) const;
   Bool_t GetXYZAt(Double_t x, Double_t b, Double_t r[3]) const;
+  Bool_t GetYAt(Double_t x,  Double_t b,  Double_t &y) const;
   void Print(Option_t* option = "") const;
   // MI
-  virtual Bool_t   CorrectForMaterial(Double_t d, Double_t x0, Double_t rho, Double_t mass);
-  virtual Bool_t   PropagateTo(Double_t x, Double_t mass, Double_t maxStep, Bool_t rotateTo=kTRUE);
-
+  virtual Bool_t   PropagateTo(Double_t x, Double_t b, Double_t mass, Double_t maxStep, Bool_t rotateTo=kTRUE, Double_t maxSnp=0.8);
+  Double_t GetSnpAt(Double_t x,Double_t b) const;
 private:
   Double_t             fX;     // X coordinate for the point of parametrisation
   Double_t             fAlpha; // Local <-->global coor.system rotation angle