vertexer->SetSkipTracks(nTrksToSkip,skipped);
vertexESD = (AliESDVertex*)vertexer->FindPrimaryVertex(event);
- } else if(fRmTrksFromPrimVtx) {
+ } else if(fRmTrksFromPrimVtx && nTrks>0) {
// removing the prongs tracks
TObjArray rmArray(nTrks);
d->SetOwnPrimaryVtx(vtx1);
}
Double_t field=aod->GetMagneticField();
- ReconstructKF(d,pdgs,field);
+ Bool_t outKF=ReconstructKF(d,pdgs,field);
+ if(!outKF) return 0;
}
//2track cuts
if(d->GetDist12toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]|| d->GetDist23toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]) return 0;
//------------------
Bool_t AliRDHFCutsLctopKpi::ReconstructKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field) const{
- const Int_t nprongs=d->GetNProngs();
+ const Int_t nprongs=d->GetNProngs();
+ if(nprongs==0) return kFALSE;
+
Int_t iprongs[nprongs];
for(Int_t i=0;i<nprongs;i++) iprongs[i]=i;