including a function to set CPV method

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

diff --git a/ITS/AliITSMeanVertexer.h b/ITS/AliITSMeanVertexer.h
index fe30545e7b0d5fb349628a64b2eb28fbcf90c0c6..99d35d4be3f69a36eacf514345c341b1777320ca 100644 (file)
--- a/ITS/AliITSMeanVertexer.h
+++ b/ITS/AliITSMeanVertexer.h
@@ -32,9 +32,10 @@ class AliITSMeanVertexer : public TObject {
  
     void   SetFilterOnContributors(Int_t nc=1){fFilterOnContributors = nc;}
     Bool_t Reconstruct(AliRawReader *rawReader);
-    void SetCutOnErrX(Double_t cut=3.){fErrXCut = cut;}
-    void SetCutOnR(Double_t cut=0.2){fRCut = cut;}
-    void SetCutOnCls(UInt_t vmin=50, UInt_t vmax=4000){fLowSPD0=vmin; fHighSPD0=vmax;}
+    void SetCutOnErrX(Double_t cut=9.5){fErrXCut = cut;}
+    void SetCutOnR(Double_t cut=0.5){fRCut = cut;}
+    void SetCutOnCls(UInt_t vmin=50, UInt_t vmax=7000){fLowSPD0=vmin; fHighSPD0=vmax;}
+    void SetZFiducialRegion(Double_t zcut=40.){fZCutDiamond=zcut;}
     void   WriteVertices(const char *filename);
 
     const TH2F*GetVertexXY() const { return fVertexXY; }
@@ -55,7 +56,9 @@ class AliITSMeanVertexer : public TObject {
     Bool_t Filter(AliESDVertex *vert,UInt_t mult);
     void   AddToMean(AliESDVertex *vert);
     Bool_t ComputeMean(Bool_t killOutliers);
-    void Reset();
+    void Reset(Bool_t redefine2D,Bool_t complete);
+    void ResetArray(){fAccEvents.ResetAllBits(kTRUE); fVertArray.Clear();
+      fIndex=0; for(Int_t i=0;i<fgkMaxNumOfEvents;i++)fClu0[i]=0; }
 
     static const Int_t fgkMaxNumOfEvents;   // max. number of events 
     AliITSDetTypeRec *fDetTypeRec; //! ITS reco class
@@ -84,8 +87,16 @@ class AliITSMeanVertexer : public TObject {
     Int_t fIndex;               //! current index on the arrays
     Double_t fErrXCut;          //! cut on error on X (error*1000<fErrXCut)
     Double_t fRCut;             //| cut on distance from first estimate (mm)
+    Double_t fZCutDiamond;      //! -/+ fZCutDiamond is the Z fiducial region
     UInt_t fLowSPD0;            //! low SPD0 cls value to accept event
     UInt_t fHighSPD0;           //! high SPD0 cls value to accept event
+    TH1F *fMultH;               //! debug hist: mult. on SPD0 before Filter
+    TH1F *fErrXH;               //! debug hist: error on X before Filter
+    TH1F *fMultHa;              //! debug hist: mult. on SPD0 after Filter
+    TH1F *fErrXHa;              //! debug hist: error on X after Filter
+    TH1F *fDistH;               //! debug hist: distance from peak
+    TH1F *fContrH;              //! debug hist: number of contributors
+    TH1F *fContrHa;             //! debug hist: number of contributors - after filter
 
     ClassDef(AliITSMeanVertexer,0)
 };