x1=x2=-b/(2*a);
return 1;
}
+ if(a==0) {
+ x1 = -c/b;
+ return 1;
+ }
// delta>0
x1 = (-b+TMath::Sqrt(delta))/(2*a);
x2 = (-b-TMath::Sqrt(delta))/(2*a);
Double_t chi2 = 0;
for(Int_t i=0;i<nPhiBins;i++) {
Double_t theo = (Double_t)npeff/(Double_t)nPhiBins;
+ if(theo==0) continue;
chi2+= (iPhiBin[i] - theo)*(iPhiBin[i] - theo)/theo;
}