void esd_kink_fill_pointset(TEvePointSet* ps) { IlcESDEvent* esd = IlcEveEventManager::AssertESD(); for (Int_t n=0; n<esd->GetNumberOfTracks(); ++n) { IlcESDtrack* track = esd->GetTrack(n); if(track->GetKinkIndex(0)<0){ IlcESDkink *kink = esd->GetKink(TMath::Abs(track->GetKinkIndex(0))-1); const TVector3 Position(kink->GetPosition()); ps->SetNextPoint(Position.X(), Position.Y(), Position.Z()); ps->SetPointId(kink); } } }