1 #ifndef ALITRACKPOINTARRAY_H
2 #define ALITRACKPOINTARRAY_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 //////////////////////////////////////////////////////////////////////////////
7 // Class AliTrackPoint //
8 // This class represent a single track space-point. //
9 // It is used to access the points array defined in AliTrackPointArray. //
10 // Note that the space point coordinates are given in the global frame. //
12 // cvetan.cheshkov@cern.ch 3/11/2005 //
13 //////////////////////////////////////////////////////////////////////////////
16 #include <TMatrixDSym.h>
20 class AliTrackPoint : public TObject {
24 AliTrackPoint(Float_t x, Float_t y, Float_t z, const Float_t *cov, UShort_t volid, Float_t charge = 0, Float_t drifttime = 0,Float_t chargeratio = 0, Int_t clutype = 0);
25 AliTrackPoint(const Float_t *xyz, const Float_t *cov, UShort_t volid, Float_t charge = 0, Float_t drifttime = 0,Float_t chargeratio = 0, Int_t clutype=0);
26 AliTrackPoint(const AliTrackPoint &p);
27 AliTrackPoint& operator= (const AliTrackPoint& p);
28 virtual ~AliTrackPoint() {}
30 // Multiplication with TGeoMatrix and distance between points (chi2) to be implemented
32 void SetXYZ(Float_t x, Float_t y, Float_t z, const Float_t *cov = 0);
33 void SetXYZ(const Float_t *xyz, const Float_t *cov = 0);
34 void SetVolumeID(UShort_t volid) { fVolumeID = volid; }
35 void SetCharge(Float_t charge) { fCharge = charge; }
36 void SetDriftTime(Float_t time) { fDriftTime = time; }
37 void SetChargeRatio(Float_t ratio) { fChargeRatio= ratio; }
38 void SetClusterType(Int_t clutype) { fClusterType= clutype; }
39 void SetExtra(Bool_t flag=kTRUE) { fIsExtra = flag; }
41 Float_t GetX() const { return fX; }
42 Float_t GetY() const { return fY; }
43 Float_t GetZ() const { return fZ; }
44 void GetXYZ(Float_t *xyz, Float_t *cov = 0) const;
45 const Float_t *GetCov() const { return &fCov[0]; }
46 UShort_t GetVolumeID() const { return fVolumeID; }
47 Float_t GetCharge() const { return fCharge; }
48 Float_t GetDriftTime() const { return fDriftTime;}
49 Float_t GetChargeRatio() const { return fChargeRatio;}
50 Int_t GetClusterType() const { return fClusterType;}
51 Bool_t IsExtra() const { return fIsExtra;}
53 Float_t GetResidual(const AliTrackPoint &p, Bool_t weighted = kFALSE) const;
54 Bool_t GetPCA(const AliTrackPoint &p, AliTrackPoint &out) const;
56 Float_t GetAngle() const;
57 Bool_t GetRotMatrix(TGeoRotation& rot) const;
58 void SetAlignCovMatrix(const TMatrixDSym alignparmtrx);
60 AliTrackPoint& Rotate(Float_t alpha) const;
61 AliTrackPoint& MasterToLocal() const;
63 void Print(Option_t *) const;
67 Float_t fX; // X coordinate
68 Float_t fY; // Y coordinate
69 Float_t fZ; // Z coordinate
70 Float_t fCharge; // Cluster charge in arbitrary units
71 Float_t fDriftTime;// Drift time in SDD (in ns)
72 Float_t fChargeRatio; // Charge ratio in SSD
73 Int_t fClusterType; // Cluster Type (encoded info on size and shape)
74 Float_t fCov[6]; // Cov matrix
75 Bool_t fIsExtra; // attached by tracker but not used in fit
76 UShort_t fVolumeID; // Volume ID
78 ClassDef(AliTrackPoint,7)
81 //////////////////////////////////////////////////////////////////////////////
82 // Class AliTrackPointArray //
83 // This class contains the ESD track space-points which are used during //
84 // the alignment procedures. Each space-point consist of 3 coordinates //
85 // (and their errors) and the index of the sub-detector which contains //
86 // the space-point. //
87 // cvetan.cheshkov@cern.ch 3/11/2005 //
88 //////////////////////////////////////////////////////////////////////////////
90 class AliTrackPointArray : public TObject {
95 AliTrackPointArray(Int_t npoints);
96 AliTrackPointArray(const AliTrackPointArray &array);
97 AliTrackPointArray& operator= (const AliTrackPointArray& array);
98 virtual ~AliTrackPointArray();
100 // Bool_t AddPoint(Int_t i, AliCluster *cl, UShort_t volid);
101 Bool_t AddPoint(Int_t i, const AliTrackPoint *p);
103 Int_t GetNPoints() const { return fNPoints; }
104 Int_t GetCovSize() const { return fSize; }
105 Bool_t GetPoint(AliTrackPoint &p, Int_t i) const;
106 // Getters for fast access to the coordinate arrays
107 const Float_t* GetX() const { return &fX[0]; }
108 const Float_t* GetY() const { return &fY[0]; }
109 const Float_t* GetZ() const { return &fZ[0]; }
110 const Float_t* GetCharge() const { return &fCharge[0]; }
111 const Float_t* GetDriftTime() const { return &fDriftTime[0]; }
112 const Float_t* GetChargeRatio() const { return &fChargeRatio[0]; }
113 const Int_t* GetClusterType() const { return &fClusterType[0]; }
114 const Bool_t* GetExtra() const { return &fIsExtra[0]; }
115 const Float_t* GetCov() const { return &fCov[0]; }
116 const UShort_t* GetVolumeID() const { return &fVolumeID[0]; }
118 Bool_t HasVolumeID(UShort_t volid) const;
120 void Sort(Bool_t down=kTRUE);
123 Bool_t fSorted; // Sorted flag
125 Int_t fNPoints; // Number of space points
126 Float_t *fX; //[fNPoints] Array with space point X coordinates
127 Float_t *fY; //[fNPoints] Array with space point Y coordinates
128 Float_t *fZ; //[fNPoints] Array with space point Z coordinates
129 Float_t *fCharge; //[fNPoints] Array with clusters charge
130 Float_t *fDriftTime; //[fNPoints] Array with drift times
131 Float_t *fChargeRatio; //[fNPoints] Array with charge ratio
132 Int_t *fClusterType; //[fNPoints] Array with cluster type
133 Bool_t *fIsExtra; //[fNPoints] Array with extra flags
134 Int_t fSize; // Size of array with cov matrices = 6*N of points
135 Float_t *fCov; //[fSize] Array with space point coordinates cov matrix
136 UShort_t *fVolumeID; //[fNPoints] Array of space point volume IDs
138 ClassDef(AliTrackPointArray,7)