// ITS Track Class
//
// Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
+// dEdx analysis by: Boris Batyunya, JINR, Boris.Batiounia@cern.ch
//-------------------------------------------------------------------------
#include "AliITSrecoV2.h"
-class AliTPCtrack;
+class AliESDtrack;
+class AliStrLine;
//_____________________________________________________________________________
-class AliITStrackV2 : public AliKalmanTrack {
+class AliITStrackV2 : public AliKalmanTrack {
public:
- AliITStrackV2():AliKalmanTrack(){}
- AliITStrackV2(const AliTPCtrack& t) throw (const Char_t *);
+ AliITStrackV2();
+ AliITStrackV2(AliESDtrack& t,Bool_t c=kFALSE) throw (const Char_t *);
AliITStrackV2(const AliITStrackV2& t);
- Int_t
- PropagateToVertex(Double_t x0=36.66,Double_t rho=1.2e-3,Double_t pm=0.139);
- void SetdEdx(Float_t dedx) {fdEdx=dedx;}
+ Int_t PropagateToVertex(Double_t d=0., Double_t x0=0.);
+ Int_t Propagate(Double_t alpha, Double_t xr);
+ virtual Int_t CorrectForMaterial(Double_t d, Double_t x0=21.82);
+ Int_t PropagateTo(Double_t xr, Double_t d, Double_t x0=21.82);
+ Double_t PropagateToDCA(AliKalmanTrack *p, Double_t d=0., Double_t x0=0.);
+ Int_t Update(const AliCluster* cl,Double_t chi2,UInt_t i);
+ Int_t Improve(Double_t x0,Double_t xyz[3],Double_t ers[3]);
+ void SetdEdx(Double_t dedx) {fdEdx=dedx;}
+ void SetSampledEdx(Float_t q, Int_t i);
+ void CookdEdx(Double_t low=0., Double_t up=0.51);
void SetDetectorIndex(Int_t i) {SetLabel(i);}
+ void ResetCovariance();
+ void ResetClusters() { SetChi2(0.); SetNumberOfClusters(0); }
+ void UpdateESDtrack(ULong_t flags) const;
- void *operator new(size_t s,void *p) { return p; }
- void *operator new(size_t s) { return ::operator new(s); }
-
+ AliESDtrack *GetESDtrack() const {return fESDtrack;}
+
Int_t GetDetectorIndex() const {return GetLabel();}
Double_t GetX() const {return fX;}
Double_t GetAlpha()const {return fAlpha;}
- Float_t GetdEdx() const {return fdEdx;}
+ Double_t GetdEdx() const {return fdEdx;}
+ Double_t GetPIDsignal() const {return GetdEdx();}
Double_t GetY() const {return fP0;}
Double_t GetZ() const {return fP1;}
Double_t GetSnp() const {return fP2;}
Double_t GetTgl() const {return fP3;}
- Double_t Get1Pt() const {return fP4*GetConvConst();}
- Double_t GetD() const;
+ Double_t GetC() const {return fP4;}
+ Double_t Get1Pt() const;
+ Double_t GetD(Double_t x=0, Double_t y=0) const;
+ Double_t GetZat(Double_t x=0) const;
+
Double_t GetSigmaY2() const {return fC00;}
Double_t GetSigmaZ2() const {return fC11;}
-
Int_t Compare(const TObject *o) const;
-
void GetExternalParameters(Double_t& xr, Double_t x[5]) const ;
void GetExternalCovariance(Double_t cov[15]) const ;
-
Int_t GetClusterIndex(Int_t i) const {return fIndex[i];}
Int_t GetGlobalXYZat(Double_t r,Double_t &x,Double_t &y,Double_t &z) const;
-
- Int_t Propagate(Double_t alpha,
- Double_t xr,Double_t x0,Double_t rho,Double_t pm=0.139);
-
+ void ApproximateHelixWithLine(Double_t xk, AliStrLine *line);
Double_t GetPredictedChi2(const AliCluster *cluster) const;
- Int_t Update(const AliCluster* cl,Double_t chi2,UInt_t i);
-
- Double_t GetPredictedChi2(const AliCluster *cluster, Double_t *m,
- Double_t x0, Double_t pm=0.139) const;
- Int_t Update(const Double_t *m, Double_t chi2, UInt_t i);
- Int_t Improve(Double_t x0,Double_t yv,Double_t zv);
- void ResetCovariance();
- void ResetClusters() { SetChi2(0.); SetNumberOfClusters(0); }
-
Int_t Invariant() const;
- //protected:
-Int_t
-PropagateTo(Double_t xr,Double_t x0=21.82,Double_t rho=2.33,Double_t pm=0.139);
+ void SetExtraCluster(Int_t i, Int_t idx) {fIndex[kMaxLayer+i]=idx;}
+ Int_t GetExtraCluster(Int_t i) const {return fIndex[kMaxLayer+i];}
-private:
+protected:
+ void GetXYZ(Float_t r[3]) const;
+
Double_t fX; // X-coordinate of this track (reference plane)
Double_t fAlpha; // rotation angle
Double_t fC30, fC31, fC32, fC33; // track
Double_t fC40, fC41, fC42, fC43, fC44; // parameters
- UInt_t fIndex[kMaxLayer]; // indices of associated clusters
+ Int_t fIndex[2*kMaxLayer]; // indices of associated clusters
+ Float_t fdEdxSample[4]; // array of dE/dx samples b.b.
- ClassDef(AliITStrackV2,1) //ITS reconstructed track
+ AliESDtrack *fESDtrack; //! pointer to the connected ESD track
+ ClassDef(AliITStrackV2,5) //ITS reconstructed track
};
inline
void AliITStrackV2::GetExternalParameters(Double_t& xr, Double_t x[5]) const {
//---------------------------------------------------------------------
- // This function return external TPC track representation
+ // This function return external ITS track representation
//---------------------------------------------------------------------
xr=fX;
- x[0]=GetY(); x[1]=GetZ(); x[2]=GetSnp(); x[3]=GetTgl(); x[4]=Get1Pt();
+ x[0]=GetY(); x[1]=GetZ(); x[2]=GetSnp(); x[3]=GetTgl();
+ x[4]=(TMath::Sign(1e-9,fP4) + fP4)*GetLocalConvConst();
+}
+
+inline
+void AliITStrackV2::SetSampledEdx(Float_t q, Int_t i) {
+ //----------------------------------------------------------------------
+ // This function stores dEdx sample corrected for the track segment length
+ // Origin: Boris Batyunya, JINR, Boris.Batiounia@cern.ch
+ //----------------------------------------------------------------------
+ if (i<0) return;
+ if (i>3) return;
+ Double_t s=GetSnp(), t=GetTgl();
+ q *= TMath::Sqrt((1-s*s)/(1+t*t));
+ fdEdxSample[i]=q;
+}
+
+inline void AliITStrackV2::GetXYZ(Float_t r[3]) const {
+ //---------------------------------------------------------------------
+ // Returns the position of the track in the global coord. system
+ //---------------------------------------------------------------------
+ Double_t cs=TMath::Cos(fAlpha), sn=TMath::Sin(fAlpha);
+ r[0]=fX*cs - fP0*sn; r[1]=fX*sn + fP0*cs; r[2]=fP1;
}
#endif