+ // Look at the Mass Spectrum:
+ if(nmass==0){
+ AliWarning("*** No peak was found in Mass difference spectrum!");
+ softProblem=kTRUE;
+ }
+ else if(nmass>0){
+ Int_t massind[3];
+ Float_t *massPos = mass->GetPositionX();
+ //check the found peaks (sorted in ascending mass order)
+ TMath::Sort(nmass,massPos,massind,kFALSE);
+ for(Int_t imass=0;imass<TMath::Min(nmass,3);imass++){
+ AliInfo(Form("*** the Mass peak for %s found at = %f GeV/c^2",part[imass],massPos[massind[imass]]));
+ if(TMath::Abs( massPos[massind[imass]]-masses[imass])> fmassPeakMax){
+ AliWarning(Form("*** the Mass peak position for %s is not in the right place, found at = %f GeV/c^2",part[imass],massPos[massind[imass]]));
+ softProblem=kTRUE;
+ }
+ }
+ // harder check on the pion peak (more statistically significant)
+ if(TMath::Abs( massPos[massind[0]]-masses[0])> 3*fmassPeakMax){
+ hardProblem=kTRUE;
+ }
+ }