Added reference from vertex to candidate (Andrea)

diff --git a/PWG3/AliAnalysisTaskSESelectHF.cxx b/PWG3/AliAnalysisTaskSESelectHF.cxx
index c559e783f6aac0881872475620364af9369c7eaa..009be1e3f567d6cc2026d677f77d0ff8438382fb 100644 (file)
--- a/PWG3/AliAnalysisTaskSESelectHF.cxx
+++ b/PWG3/AliAnalysisTaskSESelectHF.cxx
@@ -159,6 +159,7 @@ void AliAnalysisTaskSESelectHF::UserExec(Option_t */*option*/)
       AliAODRecoDecayHF2Prong *dOut=new(aodD0toKpiRef[iOutD0toKpi++]) 
        AliAODRecoDecayHF2Prong(v,px,py,pz,d0,d0err,dIn->GetDCA());
       dOut->SetOwnPrimaryVtx((AliAODVertex*)((dIn->GetOwnPrimaryVtx())->Clone()));
+      v->SetParent(dOut);
     }
     if(unsetvtx) dIn->UnsetOwnPrimaryVtx();
   } // end loop on D0->Kpi

diff --git a/PWG3/AliAnalysisVertexingHF.cxx b/PWG3/AliAnalysisVertexingHF.cxx
index 4ecde50295f4308147e9a368d4b42692749f0dd1..d71e3cf5266768b87a0ff0c09830d566ccf0dbca 100644 (file)
--- a/PWG3/AliAnalysisVertexingHF.cxx
+++ b/PWG3/AliAnalysisVertexingHF.cxx
@@ -285,12 +285,14 @@ void AliAnalysisVertexingHF::FindCandidatesESDtoAOD(AliESDEvent *esd,
                AliAODRecoDecayHF2Prong(v,px,py,pz,d0,d0err,dcap1n1);
              if(fRecoPrimVtxSkippingTrks || fRmTrksFromPrimVtx) rd->SetOwnPrimaryVtx(io2Prong->GetOwnPrimaryVtx());
              rd->SetProngIDs(2,id);
+             v->SetParent(rd);
            }
            if(okJPSI) {
              AliAODRecoDecayHF2Prong *rd=new(aodJPSItoEleRef[iJPSItoEle++]) 
                AliAODRecoDecayHF2Prong(v,px,py,pz,d0,d0err,dcap1n1);
              if(fRecoPrimVtxSkippingTrks || fRmTrksFromPrimVtx) rd->SetOwnPrimaryVtx(io2Prong->GetOwnPrimaryVtx());
              rd->SetProngIDs(2,id);
+             if(!okD0) v->SetParent(rd); // do something better here...
            }
            //printf("DCA: %f\n",rd->GetDCA());
          } else {
@@ -349,6 +351,7 @@ void AliAnalysisVertexingHF::FindCandidatesESDtoAOD(AliESDEvent *esd,
                AliAODRecoDecayHF3Prong(v,px,py,pz,d0,d0err,dcas,io3Prong->GetSigmaVert(),io3Prong->GetDist12toPrim(),io3Prong->GetDist23toPrim(),io3Prong->GetCharge());
              if(fRecoPrimVtxSkippingTrks || fRmTrksFromPrimVtx) rd->SetOwnPrimaryVtx(io3Prong->GetOwnPrimaryVtx());
              rd->SetProngIDs(3,id);
+             v->SetParent(rd);
            } else {
              new(aodD0toKpiRef[i3Prong++]) AliAODRecoDecayHF3Prong(*io3Prong);
            }
@@ -389,6 +392,7 @@ void AliAnalysisVertexingHF::FindCandidatesESDtoAOD(AliESDEvent *esd,
                  AliAODRecoDecayHF4Prong(v,px,py,pz,d0,d0err,dcas,io4Prong->GetDist12toPrim(),io4Prong->GetDist23toPrim(),io4Prong->GetDist14toPrim(),io4Prong->GetDist34toPrim(),io4Prong->GetCharge());
                if(fRecoPrimVtxSkippingTrks || fRmTrksFromPrimVtx) rd->SetOwnPrimaryVtx(io4Prong->GetOwnPrimaryVtx());
                rd->SetProngIDs(4,id);
+               v->SetParent(rd);
              } else {
                new(aodD0toKpiRef[i4Prong++]) AliAODRecoDecayHF4Prong(*io4Prong);
              }
@@ -442,6 +446,7 @@ void AliAnalysisVertexingHF::FindCandidatesESDtoAOD(AliESDEvent *esd,
                AliAODRecoDecayHF3Prong(v,px,py,pz,d0,d0err,dcas,io3Prong->GetSigmaVert(),io3Prong->GetDist12toPrim(),io3Prong->GetDist23toPrim(),io3Prong->GetCharge());
              if(fRecoPrimVtxSkippingTrks || fRmTrksFromPrimVtx) rd->SetOwnPrimaryVtx(io3Prong->GetOwnPrimaryVtx());
              rd->SetProngIDs(3,id);
+             v->SetParent(rd);
            } else {
              new(aodD0toKpiRef[i3Prong++]) AliAODRecoDecayHF3Prong(*io3Prong);
            }