Coordinates of SPD vertex

diff --git a/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.cxx b/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.cxx
index 0c3a42694bd2e3969fefd53927481fe49ebba73f..4a32c7409c526346f5da07d9eae5eb639b3edcce 100644 (file)
--- a/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.cxx
+++ b/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.cxx
@@ -142,6 +142,7 @@ void AliAnalysisTaskUpcPsi2s::Init()
        fVtxPos[i] = -666; 
        fVtxErr[i] = -666;
        fKfVtxPos[i] = -666;
+       fSpdVtxPos[i] = -666;
        }
 
 }//Init
@@ -225,6 +226,7 @@ void AliAnalysisTaskUpcPsi2s::UserCreateOutputObjects()
   fJPsiTree ->Branch("fVtxNDF", &fVtxNDF, "fVtxNDF/D");
   
   fJPsiTree ->Branch("fKfVtxPos", &fKfVtxPos[0], "fKfVtxPos[3]/D");
+  fJPsiTree ->Branch("fSpdVtxPos", &fSpdVtxPos[0], "fSpdVtxPos[3]/D");
   
   fJPsiTree ->Branch("fZDCAenergy", &fZDCAenergy, "fZDCAenergy/D");
   fJPsiTree ->Branch("fZDCCenergy", &fZDCCenergy, "fZDCCenergy/D");
@@ -281,6 +283,7 @@ void AliAnalysisTaskUpcPsi2s::UserCreateOutputObjects()
   fPsi2sTree ->Branch("fVtxNDF", &fVtxNDF, "fVtxNDF/D");
   
   fPsi2sTree ->Branch("fKfVtxPos", &fKfVtxPos[0], "fKfVtxPos[3]/D");
+  fPsi2sTree ->Branch("fSpdVtxPos", &fSpdVtxPos[0], "fSpdVtxPos[3]/D");
   
   fPsi2sTree ->Branch("fZDCAenergy", &fZDCAenergy, "fZDCAenergy/D");
   fPsi2sTree ->Branch("fZDCCenergy", &fZDCCenergy, "fZDCCenergy/D");
@@ -832,6 +835,14 @@ void AliAnalysisTaskUpcPsi2s::RunAODtree()
   fVtxErr[2] = CovMatx[2];
   fVtxChi2 = fAODVertex->GetChi2();
   fVtxNDF = fAODVertex->GetNDF();
+  
+  //SPD primary vertex
+  AliAODVertex *fSPDVertex = aod->GetPrimaryVertexSPD();
+  if(fSPDVertex){
+       fSpdVtxPos[0] = fSPDVertex->GetX();
+       fSpdVtxPos[1] = fSPDVertex->GetY();
+       fSpdVtxPos[2] = fSPDVertex->GetZ();
+       }
 
   //Tracklets
   fNtracklets = aod->GetTracklets()->GetNumberOfTracklets();
@@ -1404,6 +1415,14 @@ void AliAnalysisTaskUpcPsi2s::RunESDtree()
   fVtxErr[2] = CovMatx[2];
   fVtxChi2 = fESDVertex->GetChi2();
   fVtxNDF = fESDVertex->GetNDF();
+    
+  //SPD primary vertex
+  AliESDVertex *fSPDVertex = (AliESDVertex*) esd->GetPrimaryVertexSPD();
+  if(fSPDVertex){
+       fSpdVtxPos[0] = fSPDVertex->GetX();
+       fSpdVtxPos[1] = fSPDVertex->GetY();
+       fSpdVtxPos[2] = fSPDVertex->GetZ();
+       }
 
   //Tracklets
   fNtracklets = esd->GetMultiplicity()->GetNumberOfTracklets();

diff --git a/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.h b/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.h
index 55893976ab541f6d7b03292330081377f15056dd..634129365916f4a1dc1ea360948e65e357f05acf 100644 (file)
--- a/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.h
+++ b/PWGUD/UPC/AliAnalysisTaskUpcPsi2s.h
@@ -82,6 +82,7 @@ class AliAnalysisTaskUpcPsi2s : public AliAnalysisTaskSE {
   Double_t fVtxErr[3];
   Double_t fVtxChi2,fVtxNDF;
   Double_t fKfVtxPos[3];
+  Double_t fSpdVtxPos[3];
   UShort_t fBCrossNum, fNtracklets, fNLooseTracks;
   //vzero, zdc
   Double_t fZDCAenergy, fZDCCenergy;