HLT/TPCLib/AliHLTTPCConfMapper.h

   1 // @(#) $Id$
   2 // Original: AliHLTConfMapper.h,v 1.11 2004/07/05 09:03:11 loizides
   3
   4 #ifndef ALIHLTTPCCONFMAPPER_H
   5 #define ALIHLTTPCCONFMAPPER_H
   6 //* This file is property of and copyright by the ALICE HLT Project        *
   7 //* ALICE Experiment at CERN, All rights reserved.                         *
   8 //* See cxx source for full Copyright notice                               *
   9
  10 /** @file   AliHLTTPCConfMapper.h
  11     @author Anders Vestbo, Matthias Richter
  12     @date   Conformal mapping base class.
  13     @brief
  14 */
  15
  16 #include <vector>
  17 #include "AliHLTTPCConfMapPoint.h"
  18
  19 class AliHLTTPCConfMapTrack;
  20 class AliHLTTPCVertex;
  21 class AliHLTTPCTrackArray;
  22 class AliHLTTPCSpacePointData;
  23
  24
  25 /**
  26  * @class AliHLTTPCConfMapper
  27  * Conformal mapping base class.
  28  *
  29  * @ingroup alihlt_tpc
  30  */
  31 class AliHLTTPCConfMapper {
  32
  33  public:
  34
  35   AliHLTTPCConfMapper();
  36   //  AliHLTTPCConfMapper(AliTPCParam *param,AliHLTTPCVertex *vertex,Bool_t bench=(Bool_t)false);
  37   /** destructor */
  38   virtual ~AliHLTTPCConfMapper();
  39
  40   /**
  41    * Init and calculate bounds for the internal arrays.
  42    */
  43   void InitVolumes();
  44
  45   /**
  46    * Reset internal arrays and free memory.
  47    */
  48   void Reset();
  49
  50   void InitSector(Int_t sector,Int_t *rowrange=0,Float_t *etarange=0);
  51   void SetVertex(AliHLTTPCVertex *vertex){fVertex = vertex;}
  52   void MainVertexTrackingA();
  53   void MainVertexTrackingB();
  54   void MainVertexTracking();
  55   void NonVertexTracking();
  56   void MainVertexSettings(Int_t trackletlength, Int_t tracklength,
  57                           Int_t rowscopetracklet, Int_t rowscopetrack,Double_t maxphi=0.1,Double_t maxeta=0.1);
  58   void NonVertexSettings(Int_t trackletlength, Int_t tracklength,
  59                          Int_t rowscopetracklet, Int_t rowscopetrack);
  60
  61   /**
  62    * Read an array of space point data.
  63    * @param count    array size
  64    * @param hits     array
  65    * @return kTRUE if success.
  66    */
  67   Bool_t ReadHits(UInt_t count, AliHLTTPCSpacePointData* hits );
  68   Bool_t ReadHitsChecked(UInt_t count, AliHLTTPCSpacePointData* hits,unsigned int sizeInByte );
  69   void ClusterLoop();
  70   void CreateTrack(AliHLTTPCConfMapPoint *hit);
  71   AliHLTTPCConfMapPoint *GetNextNeighbor(AliHLTTPCConfMapPoint *start_hit,AliHLTTPCConfMapTrack *track=NULL);
  72   Int_t EvaluateHit(AliHLTTPCConfMapPoint *start_hit,AliHLTTPCConfMapPoint *hit,AliHLTTPCConfMapTrack *track);
  73
  74   Double_t CalcDistance(const AliHLTTPCConfMapPoint *hit1,const AliHLTTPCConfMapPoint *hit2) const;
  75   Double_t TrackletAngle(AliHLTTPCConfMapTrack *track,Int_t n=3) const;
  76   Bool_t VerifyRange(const AliHLTTPCConfMapPoint *hit1,const AliHLTTPCConfMapPoint *hit2) const;
  77   Int_t FillTracks();
  78
  79   //getters
  80   Int_t GetNumberOfTracks()    const {return fNTracks;}
  81   AliHLTTPCTrackArray *GetTracks() const {return fTrack;}
  82   Double_t GetMaxDca()         const {return fMaxDca;}
  83   AliHLTTPCVertex* GetVertex()     const {return fVertex;}
  84   Double_t CpuTime();
  85
  86   //setters
  87   void SetTrackCuts(Double_t hitChi2Cut, Double_t goodHitChi2, Double_t trackChi2Cut, Int_t maxdist,Bool_t vertexconstraint);
  88   void SetTrackletCuts(Double_t maxangle,Double_t goodDist,Bool_t vc);    //Set cut of tracklet for the given vertex_constraint
  89   void SetNSegments(Int_t f,Int_t g) {fNumPhiSegment=f,fNumEtaSegment=g;} //Set number of subvolumes (#segments in (phi,eta)
  90   void SetMaxDca(Double_t f) {fMaxDca = f;}
  91
  92   //setter:
  93   void SetMinPoints(Int_t f,Bool_t vconstraint) {fMinPoints[(Int_t)vconstraint] = f; }
  94   void SetVertexConstraint(Bool_t f) {fVertexConstraint =f;}
  95
  96   void SetHitChi2Cut(Double_t f,Bool_t vert) {fHitChi2Cut[(Int_t)vert]=f;}
  97   void SetGoodHitChi2(Double_t f,Bool_t vert) {fGoodHitChi2[(Int_t)vert]=f;}
  98   void SetTrackChi2Cut(Double_t f,Bool_t vert) {fTrackChi2Cut[(Int_t)vert]=f;}
  99   void SetMaxDist(Int_t f,Bool_t vert) {fMaxDist[(Int_t)vert]=f;}
 100   void SetTrackletLength(Int_t f,Bool_t vert) {fTrackletLength[(Int_t)vert]=f;}
 101   void SetRowScopeTrack(Int_t f, Bool_t vc){fRowScopeTrack[(Int_t)vc] = f;}
 102   void SetRowScopeTracklet(Int_t f, Bool_t vc){fRowScopeTracklet[(Int_t)vc] = f;}
 103   void SetMaxAngleTracklet(Double_t f, Bool_t vc){fMaxAngleTracklet[(Int_t)vc] = f;}
 104
 105   void SetPointers();
 106   void SetParamDone(Bool_t vconstraint) {fParamSet[(Int_t)vconstraint] = kTRUE;}
 107
 108   void SetClusterCutZ(Int_t c){fClusterCutZ=c;}
 109
 110  private:
 111   /** copy constructor prohibited */
 112   AliHLTTPCConfMapper(const AliHLTTPCConfMapper&);
 113   /** assignment operator prohibited */
 114   AliHLTTPCConfMapper& operator=(const AliHLTTPCConfMapper&);
 115
 116   struct AliHLTTPCConfMapContainer
 117   {
 118     void *first; // first track
 119     void *last;  // last track
 120   };
 121
 122   Bool_t fBench; //run-time measurements
 123   Int_t fNTracks; //number of tracks build.
 124
 125   AliHLTTPCVertex *fVertex; //!
 126   Bool_t fParamSet[2];  //!
 127   Bool_t fVertexFinder; //Include vertexfinding or not
 128                         //(latter case vertex=(0,0,0))
 129
 130   /** the list of hits */
 131   vector<AliHLTTPCConfMapPoint> fHit;  //!
 132   AliHLTTPCTrackArray *fTrack;  //!
 133   Double_t fMaxDca;      //cut value for momentum fit
 134
 135   AliHLTTPCConfMapContainer *fVolume;  //!  Segment volume
 136   AliHLTTPCConfMapContainer *fRow;     //!  Row volume
 137
 138    //Number of cells (segments)
 139   Int_t  fNumRowSegment;          // Total number of padrows
 140   Int_t  fNumPhiSegment;          // number of phi segments
 141   Int_t  fNumEtaSegment;          // number of eta segments
 142   Int_t  fNumRowSegmentPlusOne;   // row+1
 143   Int_t  fNumPhiSegmentPlusOne;   // phi+1
 144   Int_t  fNumEtaSegmentPlusOne;   // eta+1
 145   Int_t  fNumPhiEtaSegmentPlusOne;// phieta+1
 146   Int_t  fBounds;                 // bounds
 147   Int_t  fPhiHitsOutOfRange;      // phi hits out of range
 148   Int_t  fEtaHitsOutOfRange;      // eta hits out of range
 149
 150   //tracking range:
 151   Float_t fPhiMin; //MinPhi angle to consider
 152   Float_t fPhiMax; //MaxPhi angle to consider
 153   Float_t fEtaMin; //MinEta to consider
 154   Float_t fEtaMax; //MaxEta to consider
 155   Int_t fRowMin;   //Minimum row to consider
 156   Int_t fRowMax;   //Maximum row to consider
 157
 158   Bool_t fVertexConstraint;       //vertex constraint (true or false)
 159   Int_t fTrackletLength[2];       //minimal length of tracks
 160   Int_t fRowScopeTracklet[2];     //number of row segments to look for the next point of a tracklet
 161   Int_t fRowScopeTrack[2];        //number of row segments to look for the next point of a track
 162   Int_t fMinPoints[2];            //minimum number of points on one track
 163
 164   // Cuts
 165   Double_t fMaxAngleTracklet[2];  //limit of angle between to pieces of a tracklet
 166   Int_t fMaxDist[2];              //maximum distance between two hits
 167   Double_t fHitChi2Cut[2];        //Maximum hit chi2
 168   Double_t fGoodHitChi2[2];       //Chi2 to stop looking for next hit
 169   Double_t fTrackChi2Cut[2];      //Maximum track chi2
 170   Double_t fGoodDist;             //In segment building, distance consider good enough
 171   Double_t fMaxPhi;               //Maximum phi
 172   Double_t fMaxEta;               //Maximum eta
 173
 174   // Tracking informtion
 175   Int_t fMainVertexTracks; //number of tracks coming from the main vertex
 176   Int_t fClustersUnused;   //number of unused clusters
 177
 178   Int_t fClusterCutZ;      //Cut for taking out clusters with large Z
 179
 180   ClassDef(AliHLTTPCConfMapper,0) //Base class for conformal mapping tracking
 181 };
 182
 183 #endif