next up previous
Next: Vertex Fitting with Error Up: Vertex Fitting Previous: Vertex Fitting

Vertex Fitting without Error Matrix

In order to supply an initial vertex position, we first calculate an approximate vertex position by calculating the intersection in the xy plane of two helices arbitrarily chosen from the given set of tracks. There are two such intersections in general. We then compares the differences of z coordinates and take the intersection corresponding to the smaller distance. The mid-point in the z direction is our first guess of the common vertex ${\bf x}_v = (x_v, y_v, z_v)^T$.We approximate the helices to tangential lines at the points on the helices that are the closest in the xy plane to the trial vertex position:
\begin{displaymath}
\begin{array}
{lll}
 {\bf x} & = & {\bf x}_v^i + t ~{\bf e}_i\end{array}\end{displaymath} (26)
where
\begin{displaymath}
\begin{array}
{lll}
\begin{array}
{llll}
 {\bf x}_v^i & = & ...
 ...to the $i$-th helix at ${\bf x}_v^i$} \cr\end{array}\end{array}\end{displaymath} (27)
and $\phi_i$ is given by
\begin{displaymath}
\begin{array}
{lll}
 \phi_i = \arctan \left(
 \frac{ y_v - y...
 ...) - \phi_0^i
 + \frac{\pi}{2} ( 1 + Q/\vert Q\vert )\end{array}\end{displaymath} (28)
with (xci,yci) being the center of the i-th helix:
\begin{displaymath}
\begin{array}
{lll}
 x_c^i & = & x_0^i + 
 \left(d_\rho^i + ...
 ...ho^i + \frac{\alpha}{\kappa^i}\right) \sin\phi_0^i. \end{array}\end{displaymath} (29)
Then, what to minimize is the sum of the distance squared of each tangential line from a new trial vertex ${\bf x'}_v$:
\begin{displaymath}
\begin{array}
{lll}
 \chi^2 = \sum_i d_i^2\end{array}\end{displaymath} (30)
with
\begin{displaymath}
\begin{array}
{lll}
 d_i^2({\bf x'}_v) = ({\bf x'}_v - {\bf ...
 ...f e}_i \cdot ( {\bf x'}_v - {\bf x}_v^i ) \right)^2.\end{array}\end{displaymath} (31)
Notice that the $\chi^2$ is a quadratic funciton of ${\bf x'}_v$ so that the minimization condition
\begin{displaymath}
\begin{array}
{lll}
 \frac{\partial \chi^2}{\partial {\bf x'}_v} = 0\end{array}\end{displaymath} (32)
is a linear equation which can be solved by a single matrix inversion. The solution ${\bf x'}_v$ is our improved guess of the common vertex. We repeat this process until the vertex position converges.


next up previous
Next: Vertex Fitting with Error Up: Vertex Fitting Previous: Vertex Fitting
Keisuke Fujii
12/4/1998