ITS/AliITSVertexerZ.h

   1 #ifndef ALIITSVERTEXERZ_H
   2 #define ALIITSVERTEXERZ_H
   3
   4 #include<AliITSVertexer.h>
   5
   6 ///////////////////////////////////////////////////////////////////
   7 //                                                               //
   8 // Class for primary vertex finding                              //
   9 //                                                               //
  10 ///////////////////////////////////////////////////////////////////
  11
  12 class TFile;
  13 class TString;
  14 class AliESDVertex;
  15 class TH1F;
  16
  17 class AliITSVertexerZ : public AliITSVertexer {
  18
  19  public:
  20
  21   AliITSVertexerZ();
  22   AliITSVertexerZ(TString filename,Float_t x0=0., Float_t y0=0.);
  23   virtual ~AliITSVertexerZ();
  24   // The number of contributors set in the AliESDVertex object is the
  25   // number of tracklets used to determine the vertex position
  26   // If this number is <1, the procedure could not find a vertex position
  27   // and by default the Z coordinate is set to 0
  28   // Number of contributors = -1  --> No tracklets
  29   // Number of contributors = -2  --> No SPD recpoints
  30   virtual AliESDVertex* FindVertexForCurrentEvent(Int_t evnumb);
  31   virtual void FindVertices();
  32   virtual void PrintStatus() const;
  33   void SetDiffPhiMax(Float_t pm = 0.01){fDiffPhiMax = pm;}
  34   void ConfigIterations(Int_t noiter=3,Float_t *ptr=0);
  35   void SetFirstLayerModules(Int_t m1 = 0, Int_t m2 = 79){fFirstL1 = m1; fLastL1 = m2;}
  36   void SetSecondLayerModules(Int_t m1 = 80, Int_t m2 = 239){fFirstL2 = m1; fLastL2 = m2;}
  37   void SetLowLimit(Float_t lim=-20.){fLowLim = lim;}
  38   void SetHighLimit(Float_t lim=20.){fHighLim = lim;}
  39   Float_t GetLowLimit() const {return fLowLim;}
  40   Float_t GetHighLimit() const {return fHighLim;}
  41   void SetBinWidthCoarse(Float_t bw=0.01){fStepCoarse = bw;}
  42   void SetBinWidthFine(Float_t bw=0.0005){fStepFine = bw;}
  43   void SetPPsetting(Float_t cl2=250., Float_t coarsebin=0.02){fPPsetting[0]=cl2; fPPsetting[1]=coarsebin;}
  44   Float_t GetBinWidthCoarse() const {return fStepCoarse;}
  45   Float_t GetBinWidthFine() const {return fStepFine;}
  46   void SetTolerance(Float_t tol = 20./10000.){fTolerance = tol;}
  47   Float_t GetTolerance() const {return fTolerance;}
  48   virtual void MakeTracklet(Double_t * /* pA */, Double_t * /*pB */, Int_t & /* nolines */) {} // implemented in a derived class
  49
  50  protected:
  51   AliITSVertexerZ(const AliITSVertexerZ& vtxr);
  52   AliITSVertexerZ& operator=(const AliITSVertexerZ& /* vtxr */);
  53   void ResetHistograms();
  54   void VertexZFinder(Int_t evnumber);
  55   Float_t GetPhiMaxIter(Int_t i) const {return fPhiDiffIter[i];}
  56
  57
  58   Int_t fFirstL1;          // first module of the first pixel layer used
  59   Int_t fLastL1;           // last module of the first pixel layer used
  60   Int_t fFirstL2;          // first module of the second pixel layer used
  61   Int_t fLastL2;           // last module of the second pixel layer used
  62   Float_t fDiffPhiMax;     // Maximum delta phi allowed among corr. pixels
  63   Float_t fX0;             // Nominal x coordinate of the vertex
  64   Float_t fY0;             // Nominal y coordinate of the vertex
  65   Float_t fZFound;         //! found value for the current event
  66   Float_t fZsig;           //! RMS of Z
  67   TH1F *fZCombc;           //! histogram with coarse z distribution
  68   TH1F *fZCombv;           //! histogram with very coarse z (step=3coarse)distribution
  69   TH1F *fZCombf;           //! histogram with fine z distribution
  70   Float_t fLowLim;         // low limit for fZComb histograms
  71   Float_t fHighLim;        // high limit for fZComb histograms
  72   Float_t fStepCoarse;     // bin width for fZCombc
  73   Float_t fStepFine;       // bin width for fZCombf
  74   Float_t fTolerance;      // tolerance on the symmetry of the Z interval
  75   Float_t fPPsetting[2];   // [0] is the max. number of clusters on L2 to use [1] as fStepCoarse
  76   Int_t fMaxIter;            // Maximum number of iterations (<=5)
  77   Float_t fPhiDiffIter[5];   // Delta phi used in iterations
  78
  79   ClassDef(AliITSVertexerZ,4);
  80 };
  81
  82 #endif