#include <TRandom.h>
#include <Riostream.h>
+using std::endl;
+using std::cout;
/// \cond CLASSIMP
ClassImp(AliMUONClusterSplitterMLEM)
/// \endcond
delta /= pad->Charge();
chi2 += delta;
} // for (Int_t j=0;
- if (iflag == 0) qAver = qTot / npads;
+ if (iflag == 0 && npads) qAver = qTot / npads;
+ if (!npads)
+ {
+ AliError(Form("Got npads=0. Please check"));
+ }
f = chi2 / qAver;
}
param[j] = parOk[j];
}
param[fNpar] = 0.6;
- parmin[fNpar] = 0;
+ parmin[fNpar] = 1E-9;
parmax[fNpar++] = 1;
}
// Flag clusters for fit
nForFit = 0;
- while (minGroup[nForFit] >= 0 && nForFit < 3)
+ while (nForFit < 3 && minGroup[nForFit] >= 0)
{
if (fDebug) cout << clustNumb[minGroup[nForFit]] << " ";
clustFit[nForFit] = clustNumb[minGroup[nForFit]];