Coordinates of SPD vertex
authormbroz <Michal.Broz@cern.ch>
Mon, 20 Oct 2014 12:10:40 +0000 (14:10 +0200)
committermbroz <Michal.Broz@cern.ch>
Mon, 20 Oct 2014 12:10:40 +0000 (14:10 +0200)
PWGUD/UPC/AliAnalysisTaskUpcPsi2s.cxx
PWGUD/UPC/AliAnalysisTaskUpcPsi2s.h

index 0c3a426..4a32c74 100644 (file)
@@ -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();
index 5589397..6341293 100644 (file)
@@ -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;