- for (Int_t i = 0; i < fK; i++) {
- res = fA[i];
- AliKMeansClustering::SoftKMeans2(i+1, ic, phi, eta, res->GetMx(),res->GetMy(), res->GetSigma2(), res->GetRk());
- res->Sort(ic, phi, eta);
- Int_t j = (res->GetInd())[0];
- rk0[i] = (res->GetTarget())[j];
- }
-
- Double_t* mPhi = 0;
- Double_t* mEta = 0;
- Double_t* sigma2 = 0;
- Float_t rmax = -1.;
- Int_t imax = 0;
-
- for (Int_t i = 0; i < fK; i++) {
- if (rk0[i] > rmax) {
+ for (Int_t k = 0; k < 20; k++) {
+ Float_t rmax = -1.;
+ Int_t imax = 0;
+ for (Int_t i = 0; i < fK; i++) {
+ res = fA[i];
+ AliKMeansClustering::SoftKMeans2(i+1, ic, phi, eta, res->GetMx(),res->GetMy(), res->GetSigma2(), res->GetRk());
+ res->Sort(ic, phi, eta);
+ Int_t j = (res->GetInd())[0];
+ rk0[i] = (res->GetTarget())[j];
+ if (rk0[i] > rmax) {