ca282c36801c52f8f860cd5727cb6f205e1c38b3
[u/mrichter/AliRoot.git] / ITS / AliITStrackerSA.h
1 #ifndef ALIITSTRACKERSA_H
2 #define ALIITSTRACKERSA_H 
3
4
5
6 #include "AliITStrackerV2.h"
7
8 /* Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. *
9  * See cxx source for full Copyright notice                               */
10
11 ////////////////////////////////////////////////////
12 //  Stand alone tracker class                     //
13 //  Origin:  Elisabetta Crescio                   //
14 //  e-mail:  crescio@to.infn.it                   //
15 ////////////////////////////////////////////////////
16
17 class AliITSclusterTable;
18 class AliITStrackSA;
19 class AliESDVertex;
20 class AliITSVertexer;
21 class TTree;
22
23 class AliITStrackerSA : public AliITStrackerV2 {
24
25
26  public:
27
28   AliITStrackerSA();
29   AliITStrackerSA(AliITSgeom *geom);
30   AliITStrackerSA(AliITSgeom *geom,AliESDVertex *vert);
31   AliITStrackerSA(AliITSgeom *geom,AliITSVertexer *vertexer);
32   AliITStrackerSA(AliITStrackerSA& tracker);
33   virtual ~AliITStrackerSA();  
34   virtual Int_t Clusters2Tracks(AliESD *event){Int_t rc = AliITStrackerV2::Clusters2Tracks(event); if(!rc)rc=FindTracks(event); return rc;}
35   Int_t FindTracks(AliESD* event);
36   void  FindTracks(TTree *out,Int_t evnumber=0);
37   AliITStrackV2* FitTrack(AliITStrackSA* tr,Double_t* primaryVertex,
38                           Double_t *errorprimvert);
39
40   AliITStrackV2* FindTrackLowChiSquare(TObjArray* tracklist, Int_t dim) const;
41   Int_t LoadClusters(TTree *cf) {Int_t rc=AliITStrackerV2::LoadClusters(cf); SetClusterTree(cf); return rc;}
42   void SetVertex(AliESDVertex *vtx){fVert = vtx;}
43   void SetClusterTree(TTree * itscl){fITSclusters = itscl;}
44   void SetSixPoints(Bool_t sp = kTRUE){fSixPoints = sp;}
45   Bool_t GetSixPoints() const {return fSixPoints;}
46   void SetWindowSizes(Int_t n=46, Double_t *phi=0, Double_t *lam=0);
47   void UseFoundTracksV2(Int_t evnum,TTree* treev2);
48   void UseFoundTracksV2(AliESD *event);
49
50  protected:
51
52   // copy constructor (NO copy allowed: the constructor is protected
53   // to avoid misuse)
54   AliITStrackerSA(const AliITStrackerSA& trkr);
55   // assignment operator (NO assignment allowed)
56   AliITStrackerSA& operator=(const AliITStrackerSA& /* trkr */);
57   //Initialization
58   void Init();
59   void ResetForFinding();
60   void UpdatePoints();
61
62   static Double_t Curvature(Double_t x1,Double_t y1,Double_t x2,Double_t y2,
63                      Double_t x3,Double_t y3);
64
65   Double_t ChoosePoint(Double_t p1, Double_t p2, Double_t pp); 
66
67   static Int_t   FindIntersection(Float_t a1, Float_t b1, Float_t c1, Float_t c2, 
68                            Float_t& x1,Float_t& y1, Float_t& x2, Float_t& y2);
69   static Int_t   FindEquation(Float_t x1, Float_t y1, Float_t x2, Float_t y2, 
70                        Float_t x3, Float_t y3,Float_t& a, Float_t& b, 
71                        Float_t& c);
72  
73   static Int_t FindLabel(Int_t l1, Int_t l2, Int_t l3, Int_t l4, Int_t l5, Int_t l6);
74   static Int_t Label(Int_t gl1, Int_t gl2, Int_t gl3, Int_t gl4, Int_t gl5, 
75               Int_t gl6,Int_t gl7, Int_t gl8, Int_t gl9, Int_t gl10,Int_t gl11,
76               Int_t gl12, Int_t gl13, Int_t gl14,Int_t gl15, Int_t gl16, 
77               Int_t gl17, Int_t gl18, Int_t numberofpoints=6);
78  
79   Int_t SearchClusters(Int_t layer,Double_t phiwindow,Double_t lambdawindow, 
80                        AliITStrackSA* trs,Double_t zvertex,Int_t flagp, AliITSclusterTable* table); 
81  
82
83   Double_t fPhiEstimate; //Estimation of phi angle on next layer
84   Float_t fPoint1[2];   //! coord. of 1-st point to evaluate the curvature
85   Float_t fPoint2[2];   //! coord. of 2-nd point to evaluate the curvature
86   Float_t fPoint3[2];   //! coord. of 3-rd point to evaluate the curvature
87   Float_t fPointc[2];   //! current point coord (for curvature eval.)
88   Double_t fLambdac;    //! current value of the Lambda angle in the window
89   Double_t fPhic;       //! current value of the Phi angle in the window
90   Float_t fCoef1;       //! param. of the equation of the circ. approx a layer
91   Float_t fCoef2;       //! param. of the equation of the circ. approx a layer
92   Float_t fCoef3;       //! param. of the equation of the circ. approx a layer
93   Int_t fNloop;         //  Number of iterqations on phi and lambda windows
94   Double_t *fPhiWin;    // phi window sizes
95   Double_t *fLambdaWin; // lambda window sizes
96   AliESDVertex *fVert;        //! primary vertex
97   AliITSVertexer *fVertexer;  //! vertexer 
98   AliITSgeom *fGeom;          //! ITS geometry
99   TObjArray *fListOfTracks;   //! container for found tracks 
100   TTree *fITSclusters;        //! pointer to ITS tree of clusters
101   Bool_t fSixPoints;          // If true 6/6 points are required (default). 5/6 otherwise
102   AliITSclusterTable* fTable; //  table with clusters
103   ClassDef(AliITStrackerSA,1)
104 };
105
106 #endif
107