Possibility to set the vertex smearing in AliITSVertexerFast via AliITSRecoParam...

diff --git a/ITS/AliITSRecoParam.cxx b/ITS/AliITSRecoParam.cxx
index 481caf8df636783c1315cbcfe990ab8f74257e0c..fcb813757f83397475e5486e37330e9d69102c1b 100644 (file)
--- a/ITS/AliITSRecoParam.cxx
+++ b/ITS/AliITSRecoParam.cxx
@@ -163,7 +163,10 @@ fTrackleterPhiOverlapCut(0.005),
 fTrackleterZetaOverlapCut(0.05),
 fUseCosmicRunShiftsSSD(kFALSE),
 fSPDRemoveNoisyFlag(kTRUE),
-fSPDRemoveDeadFlag(kTRUE)
+fSPDRemoveDeadFlag(kTRUE),
+fVertexerFastSmearX(0.005),
+fVertexerFastSmearY(0.005),
+fVertexerFastSmearZ(0.01)
 {
   //
   // constructor

diff --git a/ITS/AliITSRecoParam.h b/ITS/AliITSRecoParam.h
index 1d120d3c7f2234a49ad0d8332872bd085e178de1..397e2d584c591956d7307850c5c64b1fda58b06c 100644 (file)
--- a/ITS/AliITSRecoParam.h
+++ b/ITS/AliITSRecoParam.h
@@ -77,10 +77,16 @@ class AliITSRecoParam : public AliDetectorRecoParam
   void     SetVertexerZ() { SetVertexer(1); }
   void     SetVertexerCosmics() { SetVertexer(2); }
   void     SetVertexerIons() { SetVertexer(3); }
-  void     SetVertexerSmearMC() { SetVertexer(4); }
+  void     SetVertexerSmearMC(Float_t smearx=0.005, Float_t smeary=0.005, Float_t smearz=0.01) { 
+    fVertexerFastSmearX=smearx;  fVertexerFastSmearY=smeary; fVertexerFastSmearZ=smearz; SetVertexer(4); 
+  }
   void     SetVertexerFixedOnTDI() {SetVertexer(5);} // for injection tests
   void     SetVertexerFixedOnTED() {SetVertexer(6);} // for injection tests
   Int_t    GetVertexer() const { return fVertexer; }
+  Float_t  GetVertexerFastSmearX() const {return fVertexerFastSmearX;}
+  Float_t  GetVertexerFastSmearY() const {return fVertexerFastSmearY;}
+  Float_t  GetVertexerFastSmearZ() const {return fVertexerFastSmearZ;}
+
   void     SetClusterFinder(Int_t cf=0) { fClusterFinder=cf; }
   void     SetClusterFinderV2() { SetClusterFinder(0); }
   void     SetClusterFinderOrig() { SetClusterFinder(1); }
@@ -519,9 +525,14 @@ class AliITSRecoParam : public AliDetectorRecoParam
   // should be removed at the local reconstruction step (default and safe way is true for both)
   Bool_t  fSPDRemoveNoisyFlag;  // Flag saying whether noisy pixels should be removed
   Bool_t  fSPDRemoveDeadFlag;   // Flag saying whether dead pixels should be removed
+  
+  // VertexerFast configuration
+  Float_t fVertexerFastSmearX;  // gaussian sigma for x MC vertex smearing 
+  Float_t fVertexerFastSmearY;  // gaussian sigma for y MC vertex smearing
+  Float_t fVertexerFastSmearZ;  // gaussian sigma for z MC vertex smearing
 
   
-  ClassDef(AliITSRecoParam,19) // ITS reco parameters
+  ClassDef(AliITSRecoParam,20) // ITS reco parameters
 };
 
 #endif

diff --git a/ITS/AliITSReconstructor.cxx b/ITS/AliITSReconstructor.cxx
index a2ef77d02e72708a6fd421a461b31bd02d29065f..7b5c3cb4a08e0ce178bd73ff44393d6796ea5850 100644 (file)
--- a/ITS/AliITSReconstructor.cxx
+++ b/ITS/AliITSReconstructor.cxx
@@ -202,7 +202,9 @@ AliVertexer* AliITSReconstructor::CreateVertexer() const
     vptr =  new AliITSVertexerIons();
   }
   else if(vtxOpt==4){
-    Double_t smear[3]={0.005,0.005,0.01};
+    Double_t smear[3]={GetRecoParam()->GetVertexerFastSmearX(),
+                      GetRecoParam()->GetVertexerFastSmearY(),
+                      GetRecoParam()->GetVertexerFastSmearZ()};
     Info("CreateVertexer","a AliITSVertexerFast object has been selected\n"); 
     vptr = new AliITSVertexerFast(smear);
   }

diff --git a/ITS/AliITSVertexerFast.cxx b/ITS/AliITSVertexerFast.cxx
index 2e20c39dbde820b5a66ae4a849e4ad50bdc52c1f..293074f0f65907c35b240e7cb55c902d71a7afb6 100644 (file)
--- a/ITS/AliITSVertexerFast.cxx
+++ b/ITS/AliITSVertexerFast.cxx
@@ -53,7 +53,7 @@ fSmear(0)
   // Standard constructor
   fSmear = new Double_t[3];
   for(Int_t i=0;i<3;i++)fSmear[i]=smear[i];
-  Info("AliITSVertexerFast","Gaussian smaring of the generated vertex. Parameters %f12.5 , %f12.5 , %f12.5 \n",fSmear[0],fSmear[1],fSmear[2]);
+  AliInfo(Form("Gaussian smaring of the generated vertex. Sigmas (x,y,z) = %12.5f , %12.5f , %12.5f cm",fSmear[0],fSmear[1],fSmear[2]));
   AliRunLoader *rl =AliRunLoader::Instance();
   TTree *trK=(TTree*)rl->TreeK();
   if(!trK)AliFatal("This class should be used only with simulated events!!");