Methods to get the decay length and pointing angle of V0 daughter (Annalisa)

diff --git a/PWGHF/vertexingHF/AliAODRecoCascadeHF.cxx b/PWGHF/vertexingHF/AliAODRecoCascadeHF.cxx
index 876b13525a7f1c943bcc0178d239aec6599ef5b8..443b494b92afac147cf3dea74e654248604725db 100644 (file)
--- a/PWGHF/vertexingHF/AliAODRecoCascadeHF.cxx
+++ b/PWGHF/vertexingHF/AliAODRecoCascadeHF.cxx
@@ -321,3 +321,20 @@ Bool_t AliAODRecoCascadeHF::TrigonometricalCut() const {
   return kTRUE;
 }
 
+//-----------------------------------------------------------------------------
+Double_t AliAODRecoCascadeHF::DecayLengthV0() const
+{
+  //
+  // Returns V0 decay length wrt primary vertex
+  //
+
+  AliAODv0 *v0 = (AliAODv0*)Getv0();
+
+  if (!v0) 
+    return -1.;
+  AliAODVertex *vtxPrimary = GetPrimaryVtx();
+  Double_t posVtx[3] = {0.,0.,0.};
+  vtxPrimary->GetXYZ(posVtx);
+  return v0->DecayLengthV0(posVtx);
+
+}

diff --git a/PWGHF/vertexingHF/AliAODRecoCascadeHF.h b/PWGHF/vertexingHF/AliAODRecoCascadeHF.h
index 59e0035b200b917b3752aa3fdbce7f3c954be40f..98492419b3d2b742a1903eae46aa402a1010082b 100644 (file)
--- a/PWGHF/vertexingHF/AliAODRecoCascadeHF.h
+++ b/PWGHF/vertexingHF/AliAODRecoCascadeHF.h
@@ -87,6 +87,11 @@ class AliAODRecoCascadeHF : public AliAODRecoDecayHF2Prong {
                   Int_t *pdgDg,Int_t *pdgDg2prong,
                   TClonesArray *mcArray, Bool_t isV0=kFALSE) const;
 
+  Double_t CosV0PointingAngle() const 
+  { if (Getv0()) return Getv0()->CosPointingAngle(GetPrimaryVtx()); else return -999.;}
+
+  Double_t DecayLengthV0() const;
+
  protected:
 
   ClassDef(AliAODRecoCascadeHF, 2); // heavy-flavour cascade class