Fix of parsing bug related to the reading of the calib header. Added consistency...

[u/mrichter/AliRoot.git] / ITS / AliITStrackerSA.h

diff --git a/ITS/AliITStrackerSA.h b/ITS/AliITStrackerSA.h
index 41480f58da827417c6b3424bb435b17a6f11298b..80b5da1da22e56f5bd9d591110c043bee118f2ed 100644 (file)
--- a/ITS/AliITStrackerSA.h
+++ b/ITS/AliITStrackerSA.h
@@ -3,7 +3,8 @@
 
 
 
-#include "AliITStrackerV2.h"
+#include "AliITSgeomTGeo.h"
+#include "AliITStrackerMI.h"
 
 /* Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
@@ -12,50 +13,52 @@
 //  Stand alone tracker class                     //
 //  Origin:  Elisabetta Crescio                   //
 //  e-mail:  crescio@to.infn.it                   //
+//  adapted for cosmics by A.Dainese              //
 ////////////////////////////////////////////////////
 
 class AliITSclusterTable;
 class AliITStrackSA;
 class AliESDVertex;
+class AliESDEvent;
 class AliITSVertexer;
 class TTree;
+class TArrayD;
 
-class AliITStrackerSA : public AliITStrackerV2 {
+class AliITStrackerSA : public AliITStrackerMI {
 
 
  public:
 
   AliITStrackerSA();
-  AliITStrackerSA(AliITSgeom *geom);
-  AliITStrackerSA(AliITSgeom *geom,AliESDVertex *vert);
-  AliITStrackerSA(AliITSgeom *geom,AliITSVertexer *vertexer);
-  AliITStrackerSA(AliITStrackerSA& tracker);
+  AliITStrackerSA(const Char_t *geom);
+  AliITStrackerSA(const Char_t *geom,AliESDVertex *vert);
+  AliITStrackerSA(const Char_t *geom,AliITSVertexer *vertexer);
+  AliITStrackerSA(const AliITStrackerSA& tracker);
+  AliITStrackerSA& operator=(const AliITStrackerSA& source);
   virtual ~AliITStrackerSA();  
-  virtual Int_t Clusters2Tracks(AliESD *event){Int_t rc = AliITStrackerV2::Clusters2Tracks(event); if(!rc) rc=FindTracks(event); return rc;}
-  Int_t FindTracks(AliESD* event);
-  void  FindTracks(TTree *out,Int_t evnumber=0);
-  AliITStrackV2* FitTrack(AliITStrackSA* tr,Double_t* primaryVertex,
-                          Double_t *errorprimvert);
-
-  AliITStrackV2* FindTrackLowChiSquare(TObjArray* tracklist, Int_t dim) const;
-  Int_t LoadClusters(TTree *cf) {Int_t rc=AliITStrackerV2::LoadClusters(cf); SetClusterTree(cf);SetSixPoints(kTRUE); return rc;}
+  virtual Int_t Clusters2Tracks(AliESDEvent *event);
+  Int_t FindTracks(AliESDEvent* event);
+
+  AliITStrackV2* FitTrack(AliITStrackSA* tr,Double_t* primaryVertex);
+  AliITStrackV2* FitShortTrack(AliITStrackSA* tr,Double_t *primaryVertex,
+                              Int_t innLay);
+  AliITStrackV2* FitJumpyTrack(AliITStrackSA* tr,Double_t *primaryVertex);
+  Int_t FindTrackLowChiSquare(TObjArray* tracklist, Int_t dim) const;
+  Int_t LoadClusters(TTree *cf) {Int_t rc=AliITStrackerMI::LoadClusters(cf); SetClusterTree(cf); SetSixPoints(kTRUE); return rc;}
   void SetVertex(AliESDVertex *vtx){fVert = vtx;}
   void SetClusterTree(TTree * itscl){fITSclusters = itscl;}
-  void SetSixPoints(Bool_t sp = kTRUE){fSixPoints = sp;}
+  void SetSixPoints(Bool_t sp = kFALSE){fSixPoints = sp;}
   Bool_t GetSixPoints() const {return fSixPoints;}
+  void SetOuterStartLayer(Int_t osl = 0) {if(osl>(AliITSgeomTGeo::GetNLayers()-2)) osl=AliITSgeomTGeo::GetNLayers()-2; fOuterStartLayer = osl;}
+  Int_t GetOuterStartLayer() const {return fOuterStartLayer;}
+  void SetSAFlag(Bool_t fl){fITSStandAlone=fl;}  // StandAlone flag setter
+  Bool_t GetSAFlag() const {return fITSStandAlone;} // StandAlone flag getter
   void SetWindowSizes(Int_t n=46, Double_t *phi=0, Double_t *lam=0);
-  void UseFoundTracksV2(Int_t evnum,TTree* treev2);
-  void UseFoundTracksV2(AliESD *event);
 
   enum {kSAflag=0x8000}; //flag to mark clusters used in the SA tracker
 
  protected:
 
-  // copy constructor (NO copy allowed: the constructor is protected
-  // to avoid misuse)
-  AliITStrackerSA(const AliITStrackerSA& trkr);
-  // assignment operator (NO assignment allowed)
-  AliITStrackerSA& operator=(const AliITStrackerSA& /* trkr */);
   //Initialization
   void Init();
   void ResetForFinding();
@@ -80,8 +83,16 @@ class AliITStrackerSA : public AliITStrackerV2 {
  
   Int_t SearchClusters(Int_t layer,Double_t phiwindow,Double_t lambdawindow, 
                        AliITStrackSA* trs,Double_t zvertex,Int_t flagp); 
- 
+
+  void GetCoorAngles(AliITSRecPoint* cl,Double_t &phi,Double_t &lambda,Float_t &x,Float_t &y,Float_t &z,Double_t* vertex);
+  void GetCoorErrors(AliITSRecPoint* cl,Float_t &sx,Float_t &sy, Float_t &sz);
+
+  AliITSclusterTable* GetClusterCoord(Int_t layer,Int_t n) const {return (AliITSclusterTable*)fCluCoord[layer]->UncheckedAt(n);}
+  void RemoveClusterCoord(Int_t layer, Int_t n) {fCluCoord[layer]->RemoveAt(n);fCluCoord[layer]->Compress();}
+
+
   Double_t fPhiEstimate; //Estimation of phi angle on next layer
+  Bool_t fITSStandAlone; //Tracking is performed in the ITS alone if kTRUE
   Float_t fPoint1[2];   //! coord. of 1-st point to evaluate the curvature
   Float_t fPoint2[2];   //! coord. of 2-nd point to evaluate the curvature
   Float_t fPoint3[2];   //! coord. of 3-rd point to evaluate the curvature
@@ -96,12 +107,15 @@ class AliITStrackerSA : public AliITStrackerV2 {
   Double_t *fLambdaWin; // lambda window sizes
   AliESDVertex *fVert;        //! primary vertex
   AliITSVertexer *fVertexer;  //! vertexer 
-  AliITSgeom *fGeom;          //! ITS geometry
   TObjArray *fListOfTracks;   //! container for found tracks 
   TTree *fITSclusters;        //! pointer to ITS tree of clusters
   Bool_t fSixPoints;          // If true 6/6 points are required (default). 5/6 otherwise
-  AliITSclusterTable* fTable; //  table with clusters
-  ClassDef(AliITStrackerSA,1)
+  Int_t fOuterStartLayer;     // Search for tracks with <6 points: outer layer to start from
+
+  TClonesArray** fCluLayer; //! array with clusters 
+  TClonesArray** fCluCoord; //! array with cluster info
+
+  ClassDef(AliITStrackerSA,6)
 };
 
 #endif