Store in AliVertex::fSigma the value of the DeltaPhi cut used in VertexerZ - can...
authorprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Dec 2009 13:53:35 +0000 (13:53 +0000)
committerprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Dec 2009 13:53:35 +0000 (13:53 +0000)
ITS/AliITSVertexer3D.cxx
ITS/AliITSVertexerZ.cxx
ITS/AliITSVertexerZ.h

index 00d5472..a7def28 100644 (file)
@@ -144,6 +144,7 @@ AliESDVertex* AliITSVertexer3D::FindVertexForCurrentEvent(TTree *itsClusterTree)
       Double_t chi2=99999.;
       Int_t    nContr=vtxz->GetNContributors();
       fCurrentVertex = new AliESDVertex(position,covmatrix,chi2,nContr);    
+      fCurrentVertex->SetDispersion(vtxz->GetDispersion());
       fCurrentVertex->SetTitle("vertexer: Z");
       fCurrentVertex->SetName("SPDVertexZ");
       delete vtxz;
index c476a72..81c1a5b 100644 (file)
@@ -106,7 +106,7 @@ AliITSVertexerZ::~AliITSVertexerZ() {
 void AliITSVertexerZ::ConfigIterations(Int_t noiter,Float_t *ptr){
   // configure the iterative procedure to gain efficiency for
   // pp events with very low multiplicity
-  Float_t defaults[5]={0.05,0.1,0.2,0.3,0.5};
+  Float_t defaults[5]={0.02,0.05,0.1,0.2,0.3};
   fMaxIter=noiter;
   if(noiter>5){
     Error("ConfigIterations","Maximum number of iterations is 5\n");
@@ -357,6 +357,7 @@ void AliITSVertexerZ::VertexZFinder(TTree *itsClusterTree){
   } while(niter<10 && TMath::Abs((zm-lim1)-(lim2-zm))>fTolerance);
   fCurrentVertex = new AliESDVertex(zm,ezm,ncontr);
   fCurrentVertex->SetTitle("vertexer: Z");
+  fCurrentVertex->SetDispersion(fDiffPhiMax);
   fNoVertices=1;
   points.Clear();
   if(ncontr>fMinTrackletsForPilup){ 
index b31e63a..cf492fe 100644 (file)
@@ -30,7 +30,7 @@ class AliITSVertexerZ : public AliITSVertexer {
   virtual AliESDVertex* FindVertexForCurrentEvent(TTree *itsClusterTree);
   virtual void PrintStatus() const;
   void SetDiffPhiMax(Float_t pm = 0.01){fDiffPhiMax = pm;}
-  void ConfigIterations(Int_t noiter=3,Float_t *ptr=0);
+  void ConfigIterations(Int_t noiter=4,Float_t *ptr=0);
   void SetFirstLayerModules(Int_t m1 = 0, Int_t m2 = 79){fFirstL1 = m1; fLastL1 = m2;}
   void SetSecondLayerModules(Int_t m1 = 80, Int_t m2 = 239){fFirstL2 = m1; fLastL2 = m2;}
   void SetLowLimit(Float_t lim=-40.){fLowLim = lim;}