1 #ifndef ALIITSVERTEXERZ_H
2 #define ALIITSVERTEXERZ_H
4 #include<AliITSVertexer.h>
6 ///////////////////////////////////////////////////////////////////
8 // Class for primary vertex finding //
10 ///////////////////////////////////////////////////////////////////
17 class AliITSVertexerZ : public AliITSVertexer {
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 Int_t GetPeakRegion(TH1F*h, Int_t &binmin, Int_t &binmax) const;
45 Float_t GetBinWidthCoarse() const {return fStepCoarse;}
46 Float_t GetBinWidthFine() const {return fStepFine;}
47 void SetTolerance(Float_t tol = 20./10000.){fTolerance = tol;}
48 Float_t GetTolerance() const {return fTolerance;}
49 // virtual void MakeTracklet(Double_t * /* pA */, Double_t * /*pB */, Int_t & /* nolines */) {} // implemented in a derived class
52 AliITSVertexerZ(const AliITSVertexerZ& vtxr);
53 AliITSVertexerZ& operator=(const AliITSVertexerZ& vtxr );
54 void ResetHistograms();
55 void VertexZFinder(Int_t evnumber);
56 Float_t GetPhiMaxIter(Int_t i) const {return fPhiDiffIter[i];}
59 Int_t fFirstL1; // first module of the first pixel layer used
60 Int_t fLastL1; // last module of the first pixel layer used
61 Int_t fFirstL2; // first module of the second pixel layer used
62 Int_t fLastL2; // last module of the second pixel layer used
63 Float_t fDiffPhiMax; // Maximum delta phi allowed among corr. pixels
64 Float_t fX0; // Nominal x coordinate of the vertex
65 Float_t fY0; // Nominal y coordinate of the vertex
66 Float_t fZFound; //! found value for the current event
67 Float_t fZsig; //! RMS of Z
68 TH1F *fZCombc; //! histogram with coarse z distribution
69 TH1F *fZCombv; //! histogram with very coarse z (step=3coarse)distribution
70 TH1F *fZCombf; //! histogram with fine z distribution
71 Float_t fLowLim; // low limit for fZComb histograms
72 Float_t fHighLim; // high limit for fZComb histograms
73 Float_t fStepCoarse; // bin width for fZCombc
74 Float_t fStepFine; // bin width for fZCombf
75 Float_t fTolerance; // tolerance on the symmetry of the Z interval
76 Float_t fPPsetting[2]; // [0] is the max. number of clusters on L2 to use [1] as fStepCoarse
77 Int_t fMaxIter; // Maximum number of iterations (<=5)
78 Float_t fPhiDiffIter[5]; // Delta phi used in iterations
80 ClassDef(AliITSVertexerZ,4);