Automatic switch to ITSpureSA vertexing in absence of TPC

If the TPC is absent then the ITS tracker runs only in ITSpureSA mode, therefore the AliReconstruction will switch to the VertexerTracks to ITSpureSA mode if no TPC clusters were found (done per event basis)
sevdokim · Jan 29, 2016 · 1aea71a · 1aea71a
1 parent bec60a7
commit 1aea71a
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 4 deletions.
diff --git a/STEER/ESD/AliVertexerTracks.cxx b/STEER/ESD/AliVertexerTracks.cxx
@@ -196,8 +196,8 @@ AliESDVertex* AliVertexerTracks::FindPrimaryVertex(const AliVEvent *vEvent)
   //
   int bcRound = fBCSpacing/25;   // profit from larger than 25ns spacing and set correct BC
   TDirectory * olddir = gDirectory;
-  TFile *f = 0;
-  if(nTrks>500) f = new TFile("VertexerTracks.root","recreate");
+  //  TFile *f = 0;
+  //  if(nTrks>500) f = new TFile("VertexerTracks.root","recreate");
   TObjArray trkArrayOrig(nTrks);
   UShort_t *idOrig = new UShort_t[nTrks];
   Double_t *zTr = new Double_t[nTrks];
@@ -276,12 +276,13 @@ AliESDVertex* AliVertexerTracks::FindPrimaryVertex(const AliVEvent *vEvent)
   delete [] zTr; zTr=NULL;
   delete [] err2zTr; err2zTr=NULL;
 
+  /*
   if(f) {
     f->Close(); delete f; f = NULL;
     gSystem->Unlink("VertexerTracks.root");
     olddir->cd();
   }
-
+  */
   // set vertex ID for tracks used in the fit
   // (only for ESD)
   if(!inputAOD && fCurrentVertex) {

diff --git a/STEER/ESD/AliVertexerTracks.h b/STEER/ESD/AliVertexerTracks.h
@@ -95,7 +95,8 @@ 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; }
+  void  SetITSpureSA(Bool_t v=kTRUE) { fITSpureSA=v; return; }
+  Bool_t GetITSpureSA() { return fITSpureSA; }
   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; }

diff --git a/STEER/STEER/AliReconstruction.cxx b/STEER/STEER/AliReconstruction.cxx
@@ -2390,7 +2390,13 @@ Bool_t AliReconstruction::ProcessEvent(Int_t iEvent)
 	   }
 	 } 
        }
+       Bool_t usePureSA = ftVertexer->GetITSpureSA();
+       if (!fesd->GetNumberOfTPCClusters()) {
+	 AliInfo("No TPC clusters: VertexerTracks will run with ITSpureSA tracks");
+	 ftVertexer->SetITSpureSA(kTRUE);
+       }
        AliESDVertex *pvtx=ftVertexer->FindPrimaryVertex(fesd);
+       ftVertexer->SetITSpureSA(usePureSA); // restore use settings
        if (pvtx) {
 	 if(constrSPD){
 	   TString title=pvtx->GetTitle();