Difference between revisions of "Wire angle correspondance"
Line 156: | Line 156: | ||
<center><math>\frac{\Delta degrees}{\Delta wire\ number}=\frac{36.03^{\circ}}{111\ wires}\approx\ \frac{0.3246^{\circ}}{wire\ number}</math></center> | <center><math>\frac{\Delta degrees}{\Delta wire\ number}=\frac{36.03^{\circ}}{111\ wires}\approx\ \frac{0.3246^{\circ}}{wire\ number}</math></center> | ||
− | This would imply that if the wires were evenly placed, their change in angle theta would increase by the factor of .3246 degrees for each increase in wire number, starting obviously with | + | This would imply that if the wires were evenly placed, their change in angle theta would increase by the factor of .3246 degrees for each increase in wire number, starting obviously with wire 1 at 4.79 degrees. |
− | + | Chercking this, we can find the difference between wires 1 and 2, | |
<center><math>\frac{\Delta degree}{\Delta wire\ number}=\frac{5.09^{\circ}-4.79^{\circ} }{wire\ number}=\frac{.3^{\circ}}{wire\ number} </math></center> | <center><math>\frac{\Delta degree}{\Delta wire\ number}=\frac{5.09^{\circ}-4.79^{\circ} }{wire\ number}=\frac{.3^{\circ}}{wire\ number} </math></center> | ||
− | + | Similarly, finding the difference between wires 111 and 112, | |
<center><math>\frac{\Delta degree}{\Delta wire\ number}\frac{40.82^{\circ}-40.50^{\circ}}{wire\ number}=\frac{.32^{\circ}}{wire\ number}</math></center> | <center><math>\frac{\Delta degree}{\Delta wire\ number}\frac{40.82^{\circ}-40.50^{\circ}}{wire\ number}=\frac{.32^{\circ}}{wire\ number}</math></center> |
Revision as of 20:48, 28 November 2016
Determining wire-theta correspondance
To associate the hits with the Moller scattering angle theta, the occupancy plots of the drift chamber hits by means of wire numbers and layer must be translated using the physical constraints of the detector. Using the data released for the DC:
DC: Drift Chambers(specs)
This gives the detector with a working range of 5 to 40 degrees in Theta for the lab frame, with a resolution of 1m radian.
This sets the lower limit:
This sets the upper limit:
Taking the difference,
Dividing by 112, we find
CED Verification
Using CED to verify the angle and wire correlation,
Zooming in on the view paralell to the direction of the wires in ced, we can examine the wire corresponding theta angle in the drift chamber.
Corresponding theta angles can be found for other wires, in Region 1, Superlayers 1 and 2.
Wire Number | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | Layer 6 |
---|---|---|---|---|---|---|
1 | 4.79 | 5.03 | 4.98 | 5.22 | 5.16 | 5.40 |
2 | 5.09 | 5.33 | 5.27 | 5.51 | 5.45 | 5.69 |
78 | 29.79 | 29.93 | 29.74 | 29.88 | 29.69 | 29.83 |
111 | 40.50 | 40.59 | 40.36 | 40.44 | 40.21 | 40.29 |
112 | 40.82 | 40.90 | 40.67 | 40.75 | 40.52 | 40.60 |
Wire Number | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | Layer 6 |
---|---|---|---|---|---|---|
1 | 4.79 | 5.03 | 4.98 | 5.22 | 5.16 | 5.40 |
2 | 5.09 | 5.33 | 5.27 | 5.51 | 5.45 | 5.69 |
78 | 29.79 | 29.93 | 29.74 | 29.88 | 29.70 | 29.84 |
111 | 40.51 | 40.59 | 40.36 | 40.44 | 40.22 | 40.30 |
112 | 40.82 | 40.90 | 40.67 | 40.75 | 40.52 | 40.60 |
Super Layer 1:Layer 1
Examing the range limits for the angle theta:
Taking the difference of the upper and lower limits in theta,
Dividing by the change in wire numbers (112-1=111), we find
This would imply that if the wires were evenly placed, their change in angle theta would increase by the factor of .3246 degrees for each increase in wire number, starting obviously with wire 1 at 4.79 degrees.
Chercking this, we can find the difference between wires 1 and 2,
Similarly, finding the difference between wires 111 and 112,
Changing cell size
Using Mathematica, a line can be fitted to the data collected on the wire number to angle theta correspondence.
Declaring the data set:
In[1]:= data1={{1,4.79},{2,5.09},{78,29.79},{111,40.50},{112,40.82}} Out[1]= {{1,4.79},{2,5.09},{78,29.79},{111,40.5},{112,40.82}}
Testing for a linear fit:
In[2]:= line1=Fit[data1,{1,x},x] Out[2]= 4.45395 +0.324738 x
Testing for a quadratic fit:
In[3]:= quad1=Fit[data1,{1,x,x^2},x] Out[3]= 4.45253 +0.325072 x-3.05048*10^-6 x^2
Testing for a polynomial of degree 3 fit:
In[118]:= polynomial1=Fit[data1,{1,x,x^2,x^3},x] Out[118]= 4.48739 +0.300911 x+0.000510988 x^2-2.6919*10^-6 x^3