3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 //-------------------------------------------------------
10 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
12 // The track parameterization is fixed in the following way:
13 // param0: local Y-coordinate of a track (cm)
14 // param1: local Z-coordinate of a track (cm)
15 // param2: local sine of the track momentum azimuth angle
16 // param3: tangent of the track momentum dip angle
17 // param4: 1/pt (1/(GeV/c))
19 //-------------------------------------------------------
21 #include <AliKalmanTrack.h>
24 #include "AliTPCreco.h"
25 #include "AliExternalTrackParam.h"
29 //_____________________________________________________________________________
30 class AliTPCtrack : public AliKalmanTrack {
33 AliTPCtrack(Double_t x, Double_t alpha, const Double_t p[5],
34 const Double_t cov[15], Int_t index);
35 AliTPCtrack(const AliESDtrack& t);
36 AliTPCtrack(const AliTPCtrack& t);
37 AliTPCtrack& operator=(const AliTPCtrack& o);
38 virtual ~AliTPCtrack() {}
40 Int_t Compare(const TObject *o) const;
42 void SetdEdx(Double_t dedx) {fdEdx=dedx;}
43 Double_t GetdEdx() const {return fdEdx;}
44 Double_t GetPIDsignal() const {return GetdEdx();}
46 Int_t GetClusterIndex(Int_t i) const {return fIndex[i];}
47 void SetClusterIndex(Int_t i, Int_t idx) {fIndex[i]=idx;}
49 Double_t GetC() const {return AliExternalTrackParam::GetC(GetBz());}
51 Double_t GetPredictedChi2(const AliCluster *cluster) const;
52 Bool_t PropagateTo(Double_t xr, Double_t rho=0.9e-3, Double_t x0=28.94);
53 Bool_t Update(const AliCluster *c, Double_t chi2, Int_t i);
54 Bool_t Rotate(Double_t alpha) {
55 return AliExternalTrackParam::Rotate(GetAlpha()+alpha);
58 Bool_t PropagateToVertex(const AliESDVertex *v,
59 Double_t rho=1.2e-3, Double_t x0=36.66);
60 void ResetClusters() {SetNumberOfClusters(0); SetChi2(0.);}
61 void UpdatePoints();//update points
62 Float_t* GetPoints() {return fPoints;}
64 Float_t Density(Int_t row0, Int_t row1); //calculate cluster density
65 Float_t Density2(Int_t row0, Int_t row1); //calculate cluster density
66 Double_t GetD(Double_t x=0, Double_t y=0) const {
67 return AliExternalTrackParam::GetD(x,y,GetBz());
69 AliExternalTrackParam & GetReference(){ return fReference;}
70 void UpdateReference(){ new (&fReference) AliExternalTrackParam(*this);}
71 Int_t GetKinkIndex(Int_t i) const{ return fKinkIndexes[i];}
72 Int_t* GetKinkIndexes() { return fKinkIndexes;}
73 Int_t GetV0Index(Int_t i) const{ return fV0Indexes[i];}
74 Int_t* GetV0Indexes() { return fV0Indexes;}
76 void SetFirstPoint(Int_t f) {fFirstPoint=f;}
77 void SetLastPoint(Int_t f) {fLastPoint=f;}
78 void SetRemoval(Int_t f) {fRemoval=f;}
79 void SetLab2(Int_t f) {fLab2=f;}
80 void SetKinkIndex(Int_t i, Int_t idx) {fKinkIndexes[i]=idx;}
81 void SetBConstrain(Bool_t b) {fBConstrain=b;}
82 void SetNFoundable(Int_t n) {fNFoundable=n;}
83 void SetNShared(Int_t s) {fNShared=s;}
85 Int_t GetFirstPoint() const {return fFirstPoint;}
86 Int_t GetLastPoint() const {return fLastPoint;}
87 Int_t GetRemoval() const {return fRemoval;}
88 Int_t GetLab2() const {return fLab2;}
89 Bool_t GetBConstrain() const {return fBConstrain;}
90 Int_t GetNShared() const {return fNShared;}
91 Int_t GetNFoundable() const {return fNFoundable;}
94 Double_t GetBz() const;
95 Double_t fdEdx; // dE/dx
97 Int_t fIndex[kMaxRow]; // indices of associated clusters
98 Float_t fPoints[4]; //first, max dens row end points of the track and max density
100 Float_t fSdEdx; // sigma of dedx
102 Int_t fNFoundable; //number of foundable clusters - dead zone taken to the account
103 Bool_t fBConstrain; // indicate seeding with vertex constrain
104 Int_t fLastPoint; // last cluster position
105 Int_t fFirstPoint; // first cluster position
106 Int_t fRemoval; // removal factor
107 Int_t fTrackType; // track type - 0 - normal - 1 - kink - 2 -V0 3- double found
108 Int_t fLab2; // index of corresponding track (kink, V0, double)
109 Int_t fNShared; // number of shared points
110 AliExternalTrackParam fReference; // track parameters at the middle of the chamber
111 Float_t fKinkPoint[12]; //radius, of kink, dfi and dtheta
112 Int_t fKinkIndexes[3]; // kink indexes - minus = mother + daughter
113 Int_t fV0Indexes[3]; // kink indexes - minus = mother + daughter
115 ClassDef(AliTPCtrack,4) // Time Projection Chamber reconstructed tracks