73 if( NDaughters<2 )
return;
74 float constrP[3]={fP[0], fP[1], fP[2]};
75 float constrC[6]={fC[0], fC[1], fC[2], fC[3], fC[4], fC[5]};
79 SetVtxGuess( fVtxGuess[0], fVtxGuess[1], fVtxGuess[2] );
81 for(
int i=0;
i<NDaughters;
i++ ) vtxFlag[
i] = 1;
83 Int_t nRest = NDaughters;
87 Int_t worstDaughter = 0;
88 for( Int_t it=0; it<NDaughters; it++ ){
89 if( !vtxFlag[it] )
continue;
90 const AliKFParticle &
p = *(vDaughters[it]);
92 float chi = p.GetDeviationFromVertex( tmp );
98 if( worstChi < ChiCut )
break;
100 vtxFlag[worstDaughter] = 0;
101 *
this -= *(vDaughters[worstDaughter]);
106 SetVtxGuess( fP[0], fP[1], fP[2] );
111 for(
int i=0;
i<6;
i++ ) fC[
i] = constrC[
i];
114 const AliKFParticle **vDaughtersNew=
new const AliKFParticle *[NDaughters];
115 for(
int i=0;
i<NDaughters;
i++ ){
116 if( vtxFlag[
i] ) vDaughtersNew[nDaughtersNew++] = vDaughters[
i];
119 delete[] vDaughtersNew;
122 if( nRest<=2 && GetChi2() > ChiCut*ChiCut*GetNDF() ) {
123 for(
int i=0;
i<NDaughters;
i++ ) vtxFlag[
i] = 0;