39 <<
" Flip called after segment initialization. " 40 <<
"Will lead to incorrect results!" <<
std::endl;
87 if (!
fFlipped) { x =
xx; y = yy; w = 1.0/(sigy*sigy); }
88 else { x = yy; y =
xx; w = 1.0/(sigx*sigx); }
125 if (d==0.0) d = 1.0E-9;
140 else { x = yy; y =
xx; }
144 double m = this->
BestM(x,y);
159 else { x = yy; y =
xx; }
162 double d1 = this->
GetD1(x,y);
163 if (d1>0.0) chi1 = 2.0*d1/
fResEX1;
167 double d2 = this->
GetD2(x,y);
168 if (d2>0.0) chi2 = 2.0*d2/
fResEX2;
171 if (chi1<chi2)
return chi1;
181 return (dx10*
fdX)+(dy10*
fdY);
190 return -(dx20*
fdX)-(dy20*
fdY);
double BestM(double x, double y)
void SetEndPoint1(double x, double y, double resin, double resex)
double GetD2(double x, double y)
void Flip(bool flip=true)
void AddPoint(double x, double sigx, double y, double sigy)
double Chi2(double x, double y)
double Chi2EndPoint(double x, double y)
double GetD1(double x, double y)
VertexDT finding algorithm.
Segment(bool flipped=false)
double Chi2Line(double x, double y)
void SetEndPoint2(double x, double y, double resin, double resex)