for (Int_t i=0;i<6;i++) fParams[i] = 0;
for (Int_t i=0;i<9;i++) fSumXY[i] = 0;
for (Int_t i=0;i<9;i++) fSumXZ[i] = 0;
- for (Int_t i=0;i<6;i++) {
+ for (Int_t i=0;i<5;i++) {
fSumPolY[i]=0;
fSumPolZ[i]=0;
}
for (Int_t i=0;i<6;i++) fParams[i] = 0;
for (Int_t i=0;i<9;i++) fSumXY[i] = 0;
for (Int_t i=0;i<9;i++) fSumXZ[i] = 0;
- for (Int_t i=0;i<6;i++) {
+ for (Int_t i=0;i<5;i++) {
fSumPolY[i]=0;
fSumPolZ[i]=0;
}
for (Int_t i=0;i<6;i++) fParams[i] = rieman.fParams[i];
for (Int_t i=0;i<9;i++) fSumXY[i] = rieman.fSumXY[i];
for (Int_t i=0;i<9;i++) fSumXZ[i] = rieman.fSumXZ[i];
- for (Int_t i=0;i<6;i++) {
+ for (Int_t i=0;i<5;i++) {
fSumPolY[i]=rieman.fSumPolY[i];
fSumPolZ[i]=rieman.fSumPolZ[i];
}
for (Int_t i=0;i<6;i++) fParams[i] = 0;
for (Int_t i=0;i<9;i++) fSumXY[i] = 0;
for (Int_t i=0;i<9;i++) fSumXZ[i] = 0;
- for (Int_t i=0;i<6;i++) {
+ for (Int_t i=0;i<5;i++) {
fSumPolY[i]=0;
fSumPolZ[i]=0;
}
Double_t x0 = -fParams[1]/fParams[0];
if (-fParams[2]*fParams[0]+fParams[1]*fParams[1]+1<0) return 0;
Double_t rm1 = fParams[0]/TMath::Sqrt(-fParams[2]*fParams[0]+fParams[1]*fParams[1]+1);
- if ( 1./(rm1*rm1)-(x-x0)*(x-x0)<=0) return 0;
- Double_t res = (x-x0)/TMath::Sqrt(1./(rm1*rm1)-(x-x0)*(x-x0));
+ Double_t arg = (1./rm1-(x-x0))*(1./rm1+(x-x0));
+ if ( arg <= 0) return 0;
+ Double_t res = (x-x0)/TMath::Sqrt(arg);
if (fParams[0]<0) res*=-1.;
return res;
}