DV MollerTrackRecon

From New IAC Wiki
Jump to navigation Jump to search

Moller events WITH Solenoid

LUND file with Moller events (with origin of coordinates occurring at each event)

2       1       1       1       1       0       0.000563654     3.53715 0       6.2002
1 -1 1 11 0 0 0.69 -2.4999 10993.7998 10993.80 0.000511 0 0 0
2 -1 1 11 0 0 -0.69 2.4999 6.5852 7.08 0.000511 0 0 0


From a GEMC run WITH the Solenoid ced is used to obtain the information from the eg12_rec.ev file.

      Event29.png


We take the phi angle from the Simulated Event momentum as the initial phi angle. The obtain the final phi angle, we can look at the final position of the electron with in the drift chambers.

     Detector position.png

Examining the position from Timer Based Tracking, we can see that after rotations about first the y-axis, then the z-axis transforms from the detector frame of reference to the lab frame of reference.

Euler Angles

We can use the Euler angles to perform the rotations.

For the rotation about the y axis.

Euler1.png

And the rotation about the z axis.

Euler2.png

Transformation Matrix

The Euler angles can be applied using a transformation matrix

[math]\left( \begin{array}{ccc} \cos (\theta ) & 0 & -\sin (\theta ) \\ 0 & 1 & 0 \\ \sin (\theta ) & 0 & \cos (\theta ) \\ \end{array} \right).\left( \begin{array}{c} x \\ y \\ z \\ \end{array} \right)[/math]


[math]=\left( \begin{array}{c} x \cos (\theta )-z \sin (\theta ) \\ y \\ z \cos (\theta )+x \sin (\theta ) \\ \end{array} \right)[/math]



For event #29, in sector 3, the location of the first interaction is given by

Conversions.png


Converting -25 degrees to radians, [math]\theta =-0.436332[/math] which is the rotation the detectors are rotated from the y axis.

[math]\left( \begin{array}{ccc} \cos (\theta ) & 0 & -\sin (\theta ) \\ 0 & 1 & 0 \\ \sin (\theta ) & 0 & \cos (\theta ) \\ \end{array} \right).\left( \begin{array}{c} -15.76 \\ 0 \\ 237.43 \\ \end{array} \right)[/math]

[math]=\left( \begin{array}{c} 86.0588 \\ 0. \\ 221.845 \\ \end{array} \right)[/math]

Finding [math]\phi =\frac{120\ 2 \pi }{360};[/math] since "sector -1" =3-1=2*60=120 degrees

[math]\left( \begin{array}{ccc} \cos (\phi ) & -\sin (\phi ) & 0 \\ \sin (\phi ) & \cos (\phi ) & 0 \\ 0 & 0 & 1 \\ \end{array} \right).\left( \begin{array}{c} 86.0588 \\ 0. \\ 221.845 \\ \end{array} \right)[/math]

[math]\left( \begin{array}{c} -43.0294 \\ 74.5291 \\ 221.845 \\ \end{array} \right)[/math]

This shows how the coordinates are transformed and explains the validity of using the TBTracking information to obtain a phi angle in the lab frame.

Phi shifts

This shows that the final position within the drift chambers is equal to the eTBTracking position. Normalizing the position with respect to the vertex location for the even, we can use the angle phi that this position makes with respect to origin to compare to the original phi from the momentum term to see if phi has shifted.

Phi shift line.png

Collecting the Phi shifts as discussed above, for both with and without the Solenoid, and writing to separate data files, a root macro is written to plot the phi shifts.

Phi shifts.png

Large Phi shifts

The large shifts in phi occur for two reasons: Small Energy, and large angles in theta with respect to the vertex position as shown in this sample of phi angles greater than 20°

Info.png

DV_RunGroupC_Moller#Moller_Track_Reconstruction