+ if(im2==0) rd=0;
+ if(im2>0) rd=fRandom->Integer(im2);
+ index=triggers2[rd];
+ AliVParticle *tr1 = (AliVParticle*)list->At(index);
+
+
+
+ for(Int_t kk=0;kk<im;kk++){
+ //if(kk==rd) continue;
+ if(index==triggers[kk]) continue;
+ Int_t lab=triggers[kk];
+ AliVParticle *tr2 = (AliVParticle*)list->At(lab);
+
+ Double_t detat=tr1->Eta()-tr2->Eta();
+ Double_t dphit=RelativePhi(tr1->Phi(),tr2->Phi());
+ Double_t deltaRt=TMath::Sqrt(detat*detat+dphit*dphit);
+
+ if(deltaRt>0.4) number=number-1;}