- /*
- if (iseed > 0) {
- // Find distance to the nearest neighbour
- dist[0] = dist[1] = TMath::Sqrt ((param0[min][0]-param0[min][2])*
- (param0[min][0]-param0[min][2])
- +(param0[min][1]-param0[min][3])*
- (param0[min][1]-param0[min][3]));
- if (iseed > 1) {
- dist[2] = TMath::Sqrt ((param0[min][0]-param0[min][5])*
- (param0[min][0]-param0[min][5])
- +(param0[min][1]-param0[min][6])*
- (param0[min][1]-param0[min][6]));
- rad = TMath::Sqrt ((param0[min][2]-param0[min][5])*
- (param0[min][2]-param0[min][5])
- +(param0[min][3]-param0[min][6])*
- (param0[min][3]-param0[min][6]));
- if (dist[2] < dist[0]) dist[0] = dist[2];
- if (rad < dist[1]) dist[1] = rad;
- if (rad < dist[2]) dist[2] = rad;
- }
- cout << dist[0] << " " << dist[1] << " " << dist[2] << endl;
- if (dist[TMath::LocMin(iseed+1,dist)] < 1.) { fNpar -= 3; break; }
- }
- */
-
- for (Int_t i=0; i<fNpar; i++) {
- parOk[i] = param0[min][i];
- //errOk[i] = fmin;
- errOk[i] = chi2n;
- // Bounded params
- parOk[i] = TMath::Max (parOk[i], parmin[i]);
- parOk[i] = TMath::Min (parOk[i], parmax[i]);
- }
+ for (Int_t i = 0; i < fNpar; ++i) {
+ parOk[i] = param0[min][i];
+ //errOk[i] = fmin;
+ errOk[i] = chi2n;
+ // Bounded params
+ parOk[i] = TMath::Max (parOk[i], parmin[i]);
+ parOk[i] = TMath::Min (parOk[i], parmax[i]);
+ }