|
|
Line 36: |
Line 36: |
| ===[[Determining wire-theta correspondence]]=== | | ===[[Determining wire-theta correspondence]]=== |
| | | |
− | ===CED Verification=== | + | ===[[CED Verification of DC Angle Theta and Wire Correspondance]]=== |
− | 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.
| |
− |
| |
− | <center>[[File:ParallelWireZoom.png]]</center>
| |
− |
| |
− |
| |
− | {| border=1 align=center
| |
− | |+ Table 1: Superlayer 1 Wire 1 Coordinates
| |
− | |- style="font-weight:bold; text-align:center;"
| |
− | ! style=" border:1px solid gray;"|Coordinates(cm)
| |
− | ! style=" border:1px solid gray;"|Layer 1
| |
− | ! style=" border:1px solid gray;"|Layer 2
| |
− | ! style=" border:1px solid gray;"|Layer 3
| |
− | ! style=" border:1px solid gray;"|Layer 4
| |
− | ! style=" border:1px solid gray;"|Layer 5
| |
− | ! style=" border:1px solid gray;"|Layer 6
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''x'''
| |
− | | style="border:1px solid gray;"|20.42
| |
− | | style="border:1px solid gray;"|21.51
| |
− | | style="border:1px solid gray;"|21.39
| |
− | | style="border:1px solid gray;"|22.49
| |
− | | style="border:1px solid gray;"|22.37
| |
− | | style="border:1px solid gray;"|23.47
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''y'''
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''z'''
| |
− | | style="border:1px solid gray;"|243.41
| |
− | | style="border:1px solid gray;"|244.18
| |
− | | style="border:1px solid gray;"|245.51
| |
− | | style="border:1px solid gray;"|246.28
| |
− | | style="border:1px solid gray;"|247.62
| |
− | | style="border:1px solid gray;"|248.38
| |
− | |}
| |
− |
| |
− | {| border=1 align=center
| |
− | |+ Table 2: Superlayer 1 Wire 2 Coordinates
| |
− | |- style="font-weight:bold; text-align:center;"
| |
− | ! style=" border:1px solid gray;"|Coordinates(cm)
| |
− | ! style=" border:1px solid gray;"|Layer 1
| |
− | ! style=" border:1px solid gray;"|Layer 2
| |
− | ! style=" border:1px solid gray;"|Layer 3
| |
− | ! style=" border:1px solid gray;"|Layer 4
| |
− | ! style=" border:1px solid gray;"|Layer 5
| |
− | ! style=" border:1px solid gray;"|Layer 6
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''x'''
| |
− | | style="border:1px solid gray;"|21.64
| |
− | | style="border:1px solid gray;"|22.73
| |
− | | style="border:1px solid gray;"|22.61
| |
− | | style="border:1px solid gray;"|23.71
| |
− | | style="border:1px solid gray;"|23.58
| |
− | | style="border:1px solid gray;"|24.69
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''y'''
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | | style="border:1px solid gray;"|0.00
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''z'''
| |
− | | style="border:1px solid gray;"|242.86
| |
− | | style="border:1px solid gray;"|243.62
| |
− | | style="border:1px solid gray;"|244.94
| |
− | | style="border:1px solid gray;"|245.71
| |
− | | style="border:1px solid gray;"|247.05
| |
− | | style="border:1px solid gray;"|247.81
| |
− | |}
| |
− |
| |
− | Corresponding theta angles can be found for other wires, in Region 1, Superlayers 1 and 2.
| |
− |
| |
− | {| border=1 align=center
| |
− | |+ Table 3: Superlayer 1 Wire-Angle Theta Correspondence in Degrees
| |
− | |- style="font-weight:bold; text-align:center;"
| |
− | ! style=" border:1px solid gray;"|Wire Number
| |
− | ! style=" border:1px solid gray;"|Layer 1
| |
− | ! style=" border:1px solid gray;"|Layer 2
| |
− | ! style=" border:1px solid gray;"|Layer 3
| |
− | ! style=" border:1px solid gray;"|Layer 4
| |
− | ! style=" border:1px solid gray;"|Layer 5
| |
− | ! style=" border:1px solid gray;"|Layer 6
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''1'''
| |
− | | style="border:1px solid gray;"|4.79
| |
− | | style="border:1px solid gray;"|5.03
| |
− | | style="border:1px solid gray;"|4.98
| |
− | | style="border:1px solid gray;"|5.22
| |
− | | style="border:1px solid gray;"|5.16
| |
− | | style="border:1px solid gray;"|5.40
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 2'''
| |
− | | style="border:1px solid gray;"|5.09
| |
− | | style="border:1px solid gray;"|5.33
| |
− | | style="border:1px solid gray;"|5.27
| |
− | | style="border:1px solid gray;"|5.51
| |
− | | style="border:1px solid gray;"|5.45
| |
− | | style="border:1px solid gray;"|5.69
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 62'''
| |
− | | style="border:1px solid gray;"|24.42
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 63'''
| |
− | | style="border:1px solid gray;"|24.76
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 64'''
| |
− | | style="border:1px solid gray;"|25.10
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | | style="border:1px solid gray;"|
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''77'''
| |
− | | style="border:1px solid gray;"|29.46
| |
− | | style="border:1px solid gray;"|29.60
| |
− | | style="border:1px solid gray;"|29.41
| |
− | | style="border:1px solid gray;"|29.55
| |
− | | style="border:1px solid gray;"|29.37
| |
− | | style="border:1px solid gray;"|29.51
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''78'''
| |
− | | style="border:1px solid gray;"|29.79
| |
− | | style="border:1px solid gray;"|29.93
| |
− | | style="border:1px solid gray;"|29.74
| |
− | | style="border:1px solid gray;"|29.88
| |
− | | style="border:1px solid gray;"|29.69
| |
− | | style="border:1px solid gray;"|29.83
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''79'''
| |
− | | style="border:1px solid gray;"|30.12
| |
− | | style="border:1px solid gray;"|30.26
| |
− | | style="border:1px solid gray;"|30.07
| |
− | | style="border:1px solid gray;"|30.21
| |
− | | style="border:1px solid gray;"|30.02
| |
− | | style="border:1px solid gray;"|30.16
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''110'''
| |
− | | style="border:1px solid gray;"|40.19
| |
− | | style="border:1px solid gray;"|40.28
| |
− | | style="border:1px solid gray;"|40.05
| |
− | | style="border:1px solid gray;"|40.13
| |
− | | style="border:1px solid gray;"|39.91
| |
− | | style="border:1px solid gray;"|39.99
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''111'''
| |
− | | style="border:1px solid gray;"|40.51
| |
− | | style="border:1px solid gray;"|40.59
| |
− | | style="border:1px solid gray;"|40.36
| |
− | | style="border:1px solid gray;"|40.44
| |
− | | style="border:1px solid gray;"|40.21
| |
− | | style="border:1px solid gray;"|40.29
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 112'''
| |
− | | style="border:1px solid gray;"|40.82
| |
− | | style="border:1px solid gray;"|40.90
| |
− | | style="border:1px solid gray;"|40.67
| |
− | | style="border:1px solid gray;"|40.75
| |
− | | style="border:1px solid gray;"|40.52
| |
− | | style="border:1px solid gray;"|40.60
| |
− | |}
| |
− |
| |
− |
| |
− |
| |
− | {| border=1 align=center
| |
− | |+ Table 4: Superlayer 2 Wire-Angle Theta Correspondence in Degrees
| |
− | |- style="font-weight:bold; text-align:center;"
| |
− | ! style=" border:1px solid gray;"|Wire Number
| |
− | ! style=" border:1px solid gray;"|Layer 1
| |
− | ! style=" border:1px solid gray;"|Layer 2
| |
− | ! style=" border:1px solid gray;"|Layer 3
| |
− | ! style=" border:1px solid gray;"|Layer 4
| |
− | ! style=" border:1px solid gray;"|Layer 5
| |
− | ! style=" border:1px solid gray;"|Layer 6
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''1'''
| |
− | | style="border:1px solid gray;"|4.79
| |
− | | style="border:1px solid gray;"|5.03
| |
− | | style="border:1px solid gray;"|4.98
| |
− | | style="border:1px solid gray;"|5.22
| |
− | | style="border:1px solid gray;"|5.16
| |
− | | style="border:1px solid gray;"|5.40
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 2'''
| |
− | | style="border:1px solid gray;"|5.09
| |
− | | style="border:1px solid gray;"|5.33
| |
− | | style="border:1px solid gray;"|5.27
| |
− | | style="border:1px solid gray;"|5.51
| |
− | | style="border:1px solid gray;"|5.45
| |
− | | style="border:1px solid gray;"|5.69
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''77'''
| |
− | | style="border:1px solid gray;"|29.46
| |
− | | style="border:1px solid gray;"|29.60
| |
− | | style="border:1px solid gray;"|29.41
| |
− | | style="border:1px solid gray;"|29.56
| |
− | | style="border:1px solid gray;"|29.37
| |
− | | style="border:1px solid gray;"|29.51
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''78'''
| |
− | | style="border:1px solid gray;"|29.79
| |
− | | style="border:1px solid gray;"|29.93
| |
− | | style="border:1px solid gray;"|29.74
| |
− | | style="border:1px solid gray;"|29.88
| |
− | | style="border:1px solid gray;"|29.70
| |
− | | style="border:1px solid gray;"|29.84
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''79'''
| |
− | | style="border:1px solid gray;"|30.13
| |
− | | style="border:1px solid gray;"|30.27
| |
− | | style="border:1px solid gray;"|30.07
| |
− | | style="border:1px solid gray;"|30.21
| |
− | | style="border:1px solid gray;"|30.02
| |
− | | style="border:1px solid gray;"|30.16
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''110x'''
| |
− | | style="border:1px solid gray;"|40.19
| |
− | | style="border:1px solid gray;"|40.28
| |
− | | style="border:1px solid gray;"|40.05
| |
− | | style="border:1px solid gray;"|40.13
| |
− | | style="border:1px solid gray;"|39.91
| |
− | | style="border:1px solid gray;"|39.99
| |
− | |-
| |
− | | style="border:1px solid gray;"|'''111'''
| |
− | | style="border:1px solid gray;"|40.51
| |
− | | style="border:1px solid gray;"|40.59
| |
− | | style="border:1px solid gray;"|40.36
| |
− | | style="border:1px solid gray;"|40.44
| |
− | | style="border:1px solid gray;"|40.22
| |
− | | style="border:1px solid gray;"|40.30
| |
− | |-
| |
− | | style="border:1px solid gray;"|''' 112'''
| |
− | | style="border:1px solid gray;"|40.82
| |
− | | style="border:1px solid gray;"|40.90
| |
− | | style="border:1px solid gray;"|40.67
| |
− | | style="border:1px solid gray;"|40.75
| |
− | | style="border:1px solid gray;"|40.52
| |
− | | style="border:1px solid gray;"|40.60
| |
− | |}
| |
| | | |
| ====Super Layer 1:Layer 1==== | | ====Super Layer 1:Layer 1==== |
Introduction
Moller Scattering
Moller Scattering Definition
GEANT4 Simulation of Moller Scattering
LH2 Target
NH2 Target
LH2 Vs. NH3
Effects Due to Target Material
Differential Cross-Section Offset
Weighted Isotropic Distribution in Lab Frame
GEMC Simulation
Drift Chamber
Super Layer 1:Layer 1
Table 1: Superlayer 1 Change of Coordinates From Wire 1 to Wire 2(cm)
Coordinates(cm)
|
Layer 1
|
Layer 2
|
Layer 3
|
Layer 4
|
Layer 5
|
Layer 6
|
[math]\Delta x[/math]
|
1.22
|
1.22
|
1.22
|
1.22
|
1.21
|
1.22
|
[math]\Delta y[/math]
|
0.00
|
0.00
|
0.00
|
0.00
|
0.00
|
0.00
|
[math]\Delta z[/math]
|
.55
|
.56
|
.57
|
.57
|
.57
|
.57
|
Using the geometric construction for determining angle theta to wire 2:
[math]\tan{\theta}=\frac{x+\Delta x}{z+\Delta z}\Rightarrow \theta=\arctan{\frac{x+\Delta x}{z+\Delta z}}[/math]
[math]\arctan{\frac{20.42+1.212}{243.39-0.5652}}=\arctan{\frac{21.632}{242.8248}}=\arctan{0.0891}=5.09^{\circ}[/math]
This equation can be solved for a hypothetical wire 0, which will allow the wire number to be the multiplicative factor for the change from the starting position.
[math]\tan{\theta}=\frac{x-\Delta x+n\Delta x}{z-\Delta z+n\Delta z}\Rightarrow \theta=\arctan{\frac{x'+n\Delta x}{z'+n\Delta z}}[/math]
where
[math]x'=20.42-1.212\approx 19.21cm[/math]
[math]z'=243.39+0.5652\approx 243.96cm[/math]
[math]\theta=\arctan{\frac{19.21+n1.212}{243.39-n0.5652}}[/math]
[math]\tan{\theta}=\frac{19.21+n1.212}{243.39-n0.5652}\Rightarrow \tan{ \theta}=-2.14437 + \frac{957.412}{430.626 - n}\Rightarrow 430.626 - n = \frac{957.412}{\tan{ \theta}+2.14437}[/math]
[math]\Rightarrow n = \frac{-957.412}{\tan{ \theta}+2.14437}+430.626[/math]
Using Mathematica, a series expansion about n=0 can be found:
In[1]:= Series[ArcTan[(19.21+1.212 x)/(243.39-0.5652 n)],{n,0,4}]
Out[1]= 0.0787635+0.00513098 n+9.83726*10^-6 n^2-2.61675*10^-8 n^3-2.22826*10^-10 n^4+O[n]^5
[math]\theta\approx 0.0787635+0.00513098 n+9.83726\times 10^{-6} n^2-2.61675\times 10^{-8} n^3[/math]
This expression will find the angle theta in radians given the wire number. To convert from radians to degrees, we can multiply by 180 and divide by Pi.
In[2]:= 180*(0.07876354504106763`+0.005130982891104289` n+9.837257652182922`*^-6 n^2-2.616751708660921`*^-8 n^3-2.2282595631625944`*^-10 n^4+O[n]^5)/3.14159625359
Out[2]= 4.51281+0.293983 n+0.000563633 n^2-1.49929*10^-6 n^3-1.2767*10^-8 n^4+O[n]^5
[math]\theta\approx 4.51281+0.293983 n+0.000563633 n^2-1.49929\times10^{-6} n^3[/math]
This tells us that the expression for theta will follow a function that comes from a series expansion. To account for instrument and measurement errors, angle measurements from ced can be used to find a better fit.
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},{62,24.42},{63,24.76},{64,25.10},{78,29.79},{111,40.50},{112,40.82}}
Out[1]= {{1,4.79},{2,5.09},{62,24.42},{63,24.76},{64,25.1},{78,29.79},{111,40.5},{112,40.82}}
Testing for a linear fit:
In[2]:= line1=Fit[data1,{1,n},n]
Out[2]= 4.39975 + 0.32469 n
Examing the range limits for the angle theta for layer 1:
Taking the difference of the upper and lower limits in theta,
[math]\Delta degrees=40.82^{\circ}-4.79^{\circ}=36.03^{\circ}[/math]
Dividing by the change in wire numbers (112-1=111), we find
[math]\frac{\Delta degrees}{\Delta wire\ number}=\frac{36.03^{\circ}}{111\ wires}\approx\ \frac{0.325^{\circ}}{wire\ number}[/math]
This would imply that if the wires were evenly placed, their change in angle theta would increase by the factor of .325 degrees for each increase in wire number, starting obviously with wire 1 at 4.79 degrees. In addition, this implies that the bin spacing for each wire would be around .325 degrees in width.
Chercking this, we can find the difference between wires 1 and 2,
[math]\frac{\Delta degree}{\Delta wire\ number}=\frac{5.09^{\circ}-4.79^{\circ} }{wire\ number}=\frac{.3^{\circ}}{wire\ number} [/math]
Similarly, finding the difference between wires 111 and 112,
[math]\frac{\Delta degree}{\Delta wire\ number}=\frac{40.82^{\circ}-40.51^{\circ}}{wire\ number}=\frac{.31^{\circ}}{wire\ number}[/math]
These differing values show that the bin width is not uniform in length, therefore a first order, linear fit, will not suffice.
Testing for a quadratic fit:
In[3]:= quad1=Fit[data1,{1,n,n^2},n]
Out[3]= 4.45564 +0.32015 n+0.0000417787 n^2
The quadratic fit does not work since it's first derivative
[math]\theta'\approx \Delta \theta=.325072-6.83308\ n[/math]
does not give the same spacing between low and high values of n as seen in the ced data. The coefficient near the first order will have to smaller than .31 to find a correlation that would agree.
Testing for a polynomial of degree 3 fit:
In[4]:= polynomial1=Fit[data1,{1,n,n^2,n^3},n]
Out[4]= 4.49876 +0.293001 n+0.000679074 n^2-3.57132*10^-6 n^3
This fit best matches the data found in ced
[math]\theta\equiv 4.49876 +0.293001 n+0.000679074 n^2-3.57132\times 10^{-6} n^3[/math]
As discussed earlier, taking the 1st derivative of this function will give us the spacing of the bins as a function of wire number.
[math]\theta '\equiv 0.293001 + 0.00135815 n - 0.000010714 n^2[/math]
The derivative of this function then will tell us where the bin spacing is at a minimum and a maximum.
[math]\theta ''\equiv .001021976-1.61514\times 10^{-5} n[/math]
[math]0.00135815 - 0.0000214279 n=0\Rightarrow Bin\ Spacing\ Maximum\ at\ n=63.3822[/math]
Comparing this maximum, we can see from Table 3 that near the maximum of 63 are seperated by larger distances than at n=1 or n=112. This is the midpoint of the plane as seen in the geometry file.
Binning based on wire numbers
The bin size based on wire number will need to be a uniform width of 1, as in an increment of 1 between the integer values of the wires. This uniformity in bin size based on wire numbers is not uniform when viewed by the angle theta due to the Drift Chamber geometry discussed earlier.
Modifying evioreader
gStyle->SetStripDecimals(kTRUE);
TF1 *fit_function=new TF1("fit_function","[0]+[1]*x+[2]*x*x+[3]*x*x*x",4.49876,41.12592);
fit_function->SetParameters(4.49876,0.293001,0.000679074,-0.00000357132);
TGaxis *A1 = new TGaxis(0,5000,113,5000,"fit_function",510,"-");
A1->SetTitle("Angle Theta(degrees)");
A1->Draw();
![Layer1binWeighted.png](/./images/8/8b/Layer1binWeighted.png)
![Layer1binWeighted Isotropic.png](/./images/c/ca/Layer1binWeighted_Isotropic.png)
Using the expression for n in terms of Theta:
[math]n = \frac{-957.412}{\tan{ \theta}+2.14437}+430.626[/math]
This relationship can be used to multiply each Moller Scattering angle theta in the lab frame, with it's differential cross-section weight, to find the Moller differential cross-section as a function of wire number in the lab frame.
![MolThetaWireWeighted.png](/./images/c/cd/MolThetaWireWeighted.png)
![MolThetaWireWeightedAdjusted.png](/./images/7/72/MolThetaWireWeightedAdjusted.png)
![TheoryDCbinsWire.png](/./images/a/aa/TheoryDCbinsWire.png)
Detector Occupancy