Possibility to reconstruct a primary vertex using ITS-standalone tracks only (A....
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 10 Mar 2010 17:32:05 +0000 (17:32 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 10 Mar 2010 17:32:05 +0000 (17:32 +0000)
STEER/AliVertexerTracks.cxx
STEER/AliVertexerTracks.h

index 5030484..d75af21 100644 (file)
@@ -67,6 +67,7 @@ fMaxd0z0(0.5),
 fMinDetFitter(100.),
 fMaxTgl(1000.),
 fITSrefit(kTRUE),
+fITSpureSA(kFALSE),
 fFiducialR(3.),
 fFiducialZ(30.),
 fnSigmaForUi00(1.5),
@@ -101,6 +102,7 @@ fMaxd0z0(0.5),
 fMinDetFitter(100.),
 fMaxTgl(1000.),
 fITSrefit(kTRUE),
+fITSpureSA(kFALSE),
 fFiducialR(3.),
 fFiducialZ(30.),
 fnSigmaForUi00(1.5),
@@ -176,8 +178,10 @@ AliESDVertex* AliVertexerTracks::FindPrimaryVertex(AliVEvent *vEvent)
     }
     if(skipThis) continue;
 
-    // skip pure ITS SA tracks
-    if(track->GetStatus()&AliESDtrack::kITSpureSA) continue;
+    // skip pure ITS SA tracks (default)
+    if(!fITSpureSA && (track->GetStatus()&AliESDtrack::kITSpureSA)) continue;
+    // or use only pure ITS SA tracks
+    if(fITSpureSA && !(track->GetStatus()&AliESDtrack::kITSpureSA)) continue;
 
     // kITSrefit
     if(fMode==0 && fITSrefit && !(track->GetStatus()&AliESDtrack::kITSrefit)) continue;
index a5e5f4f..96f008f 100644 (file)
@@ -86,6 +86,7 @@ class AliVertexerTracks : public TObject {
   void  SetDCAcutIter0(Double_t maxdca) { fDCAcutIter0=maxdca; return; }
   void  SetFinderAlgorithm(Int_t opt=1) { fAlgo=opt; return; }
   void  SetITSrefitRequired() { fITSrefit=kTRUE; return; }
+  void  SetITSpureSA() { fITSpureSA=kTRUE; return; }
   Bool_t GetITSrefitRequired() const { return fITSrefit; }
   void  SetITSrefitNotRequired() { fITSrefit=kFALSE; return; }
   void  SetFiducialRZ(Double_t r=3,Double_t z=30) { fFiducialR=r; fFiducialZ=z; return; }
@@ -159,6 +160,8 @@ class AliVertexerTracks : public TObject {
   Double_t  fMaxTgl;          // maximum tgl of tracks
   Bool_t    fITSrefit;        // if kTRUE (default), use only kITSrefit tracks
                               // if kFALSE, use all tracks (also TPC only)
+  Bool_t    fITSpureSA;       // if kFALSE (default) skip ITSpureSA tracks
+                              // if kTRUE use only those
   Double_t  fFiducialR;       // radius of fiducial cylinder for tracks 
   Double_t  fFiducialZ;       // length of fiducial cylinder for tracks
   Double_t  fnSigmaForUi00;   // n. sigmas from finder in TrackToPoint
@@ -182,7 +185,7 @@ class AliVertexerTracks : public TObject {
   AliVertexerTracks(const AliVertexerTracks & source);
   AliVertexerTracks & operator=(const AliVertexerTracks & source);
 
-  ClassDef(AliVertexerTracks,12) // 3D Vertexing with tracks 
+  ClassDef(AliVertexerTracks,14) // 3D Vertexing with tracks 
 };
 
 #endif