3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //-------------------------------------------------------------------------
10 // that is the base for AliTPCtracker, AliITStrackerV2 and AliTRDtracker
11 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
12 //-------------------------------------------------------------------------
14 #include <TGeoGlobalMagField.h>
18 #include "AliRecoParam.h"
19 #include "AliPlaneEff.h"
25 class AliExternalTrackParam;
29 class AliTracker : public TObject {
32 virtual ~AliTracker(){}
33 virtual Int_t Clusters2Tracks(AliESDEvent *event)=0;
34 virtual Int_t PropagateBack(AliESDEvent *event)=0;
35 virtual Int_t RefitInward(AliESDEvent *event)=0;
36 virtual Int_t PostProcess(AliESDEvent */*event*/) {return 0;}
37 void SetVertex(const Double_t *xyz, const Double_t *ers=0) {
38 fX=xyz[0]; fY=xyz[1]; fZ=xyz[2];
39 if (ers) { fSigmaX=ers[0]; fSigmaY=ers[1]; fSigmaZ=ers[2]; }
43 virtual Int_t LoadClusters(TTree *)=0;
44 virtual void UnloadClusters()=0;
45 virtual void FillClusterArray(TObjArray* array) const;
46 virtual AliCluster *GetCluster(Int_t index) const=0;
47 virtual AliPlaneEff *GetPlaneEff() {return NULL;}
48 virtual Bool_t GetTrackPoint(Int_t /* index */ , AliTrackPoint& /* p */) const { return kFALSE;}
49 virtual Bool_t GetTrackPointTrackingError(Int_t /* index */,
50 AliTrackPoint& /* p */, const AliESDtrack* /* t */) { return kFALSE;}
51 virtual void UseClusters(const AliKalmanTrack *t, Int_t from=0) const;
52 virtual void CookLabel(AliKalmanTrack *t,Float_t wrong) const;
53 Double_t GetX() const {return fX;}
54 Double_t GetY() const {return fY;}
55 Double_t GetZ() const {return fZ;}
56 Double_t GetSigmaX() const {return fSigmaX;}
57 Double_t GetSigmaY() const {return fSigmaY;}
58 Double_t GetSigmaZ() const {return fSigmaZ;}
60 static Double_t GetTrackPredictedChi2(AliExternalTrackParam *track,
61 Double_t mass, Double_t step,
62 const AliExternalTrackParam *backup);
64 Double_t MeanMaterialBudget(const Double_t *start, const Double_t *end, Double_t *mparam);
66 Bool_t PropagateTrackTo(AliExternalTrackParam *track, Double_t x, Double_t m,
67 Double_t maxStep, Bool_t rotateTo=kTRUE, Double_t maxSnp=0.8, Double_t sign=1.);
68 static Bool_t PropagateTrackToBxByBz(AliExternalTrackParam *track, Double_t x,
70 Double_t maxStep, Bool_t rotateTo=kTRUE, Double_t maxSnp=0.8,Double_t sign=1.);
72 static Double_t GetBz(const Double_t *r);
73 static void GetBxByBz(const Double_t r[3], Double_t b[3]);
74 static Double_t GetBz();
75 static Bool_t UniformField();
77 static void FillResiduals(const AliExternalTrackParam *t,
78 Double_t *p, Double_t *cov,
79 UShort_t id, Bool_t updated=kTRUE);
80 static void FillResiduals(const AliExternalTrackParam *t,
81 const AliCluster *c, Bool_t updated=kTRUE);
82 static void SetFillResiduals(AliRecoParam::EventSpecie_t es, Bool_t flag=kTRUE) { fFillResiduals=flag; fEventSpecie = es ;}
83 static void SetResidualsArray(TObjArray **arr) { fResiduals=arr; }
84 static TObjArray ** GetResidualsArray() { return fResiduals; }
86 void SetEventInfo(AliEventInfo *evInfo) {fEventInfo = evInfo;}
87 const AliEventInfo* GetEventInfo() const {return fEventInfo;}
90 AliTracker(const AliTracker &atr);
92 AliTracker & operator=(const AliTracker & atr);
93 static Bool_t fFillResiduals; // Fill residuals flag
94 static TObjArray **fResiduals; //! Array of histograms with residuals
96 Double_t fX; //X-coordinate of the primary vertex
97 Double_t fY; //Y-coordinate of the primary vertex
98 Double_t fZ; //Z-coordinate of the primary vertex
100 Double_t fSigmaX; // error of the primary vertex position in X
101 Double_t fSigmaY; // error of the primary vertex position in Y
102 Double_t fSigmaZ; // error of the primary vertex position in Z
104 static AliRecoParam::EventSpecie_t fEventSpecie ; //! event specie, see AliRecoParam
105 AliEventInfo* fEventInfo; //! pointer to the event info object
107 ClassDef(AliTracker,5) //abstract tracker
110 //__________________________________________________________________________
111 inline Bool_t AliTracker::UniformField()
113 AliMagF* fld = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
114 return fld ? fld->IsUniform():kTRUE;