Difference between revisions of "Wire angle correspondance"
(152 intermediate revisions by the same user not shown) | |||
Line 517: | Line 517: | ||
<center>[[File:Layer1bins.png]]</center> | <center>[[File:Layer1bins.png]]</center> | ||
+ | <center>[[File:Layer1bins_Isotropic.png]]</center> | ||
<pre> | <pre> | ||
Line 528: | Line 529: | ||
− | <center>[[File:Layer1binWeighted.png]]</center> | + | <center>[[File:Layer1binWeighted.png]][[File:MolThetaLabWeighted.png]]</center> |
+ | |||
+ | |||
+ | <center>[[File:Layer1binWeighted_Isotropic.png]][[File:MolThetaLabWeighted2.png]]</center> | ||
Using the expression for n in terms of Theta: | Using the expression for n in terms of Theta: | ||
Line 538: | Line 542: | ||
− | <center>[[File:MolThetaWireWeighted.png]]</center> | + | <center>[[File:MolThetaWireWeighted.png]][[File:MolThetaWireWeightedIsotropic.png]]</center> |
− | <center>[[File:MolThetaWireWeightedAdjusted.png]]</center> | + | <center>[[File:MolThetaWireWeightedAdjusted.png]][[File:MolThetaWireWeightedAdjustedIsotropic.png]]</center> |
− | <center>[[File:TheoryDCbinsWire.png]]</center> | + | <center>[[File:TheoryDCbinsWire.png]][[File:TheoryDCbinsWireIsotropic.png]]</center> |
=Moller Occupancy Rate= | =Moller Occupancy Rate= | ||
+ | |||
+ | |||
+ | ==Comparing Rates== | ||
+ | Using the values from Whitney: | ||
+ | |||
+ | <center>[[File:Whitney_values.png]]</center> | ||
+ | |||
+ | For a 5cm target of LH2: | ||
+ | |||
+ | <center><math>\rho_{target}\times l_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6.022\times10^{23} molecules\ LH_2}{1 mole} \frac{2 atoms}{1\ molecule\ LH_2} \times \frac{1m^3}{(100 cm)^3} \times \frac{5 cm}{ }=2.11\times 10^{23} \times \frac{10^{-24} cm^{2}}{barn}=2.11\times 10^{-1}\ barns^{-1}</math></center> | ||
+ | |||
+ | |||
+ | Using the beam current of 100nA, | ||
+ | |||
+ | <center><math>\Phi_{beam}=\frac{100\times 10^{-9}\ A}{1} \times \frac{1\ C}{1\ A} \times \frac{1\ e^{-}}{1\ s} \times \frac{1}{1.602\times 10^{-19}C}\Rightarrow 6.2422\times 10^{11}\ \frac{1}{s}=\Phi_{e^{-}}</math></center> | ||
+ | |||
+ | |||
+ | Given the beam Luminosity of: | ||
+ | <center><math>\mathcal {L}=\frac{1.32\times 10^{35}}{cm^2 \cdot s}\frac{10^{-24} cm^{2}}{barn}=\frac{1.32\times 10^{11}}{barn\cdot s}=\Phi_{beam}\ \rho\ l_{target}</math></center> | ||
+ | |||
+ | We can check to make sure the density makes sense | ||
+ | |||
+ | <center><math>\Rightarrow \rho_{target}=\frac{1.32\times 10^{35}}{ \Phi_{beam}\ l_{target}cm^2 \cdot s}=\frac{1.32\times 10^{35}}{ 6.24\times 10^{11}\ \cdot 5cm^3}=\frac{4.23\times 10^{22}}{cm^3}</math></center> | ||
+ | |||
+ | |||
+ | <center><math>\Rightarrow \rho_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6.022\times10^{23} molecules\ LH_2}{1 mole} \times \frac{1m^3}{(100 cm)^3}\frac{2 atoms}{1\ molecule\ LH_2}=\frac{4.23\times 10^{22}}{cm^3}</math></center> | ||
+ | |||
+ | |||
+ | Using the values from Whitney : | ||
+ | |||
+ | |||
+ | https://wiki.iac.isu.edu/index.php/CLAS12_RateEst_byWA | ||
+ | |||
+ | {| class="wikitable" align="center" border=1 | ||
+ | |+ '''Rates''' | ||
+ | |- | ||
+ | | Energy || 6 GeV || 11 GeV | ||
+ | |- | ||
+ | |Process || (nb) || (nb) | ||
+ | |- | ||
+ | |Moller || 22773001||75008636 | ||
+ | |- | ||
+ | | DIS + radiative tail || 128 || 83 | ||
+ | |- | ||
+ | | Elastic e-p || 5511220 ||3670740 | ||
+ | |- | ||
+ | |Elastic radiative tail || 24705 || 12944 | ||
+ | |- | ||
+ | | π0 electro-production || 14802 ||17908 | ||
+ | |- | ||
+ | | π0 photo-production || 569 || 852 | ||
+ | |- | ||
+ | | π+ electro-production || 4032 || 5536 | ||
+ | |- | ||
+ | | π+ photo-production || 282 || 487 | ||
+ | |- | ||
+ | | π− electro-production || 2806 || 3843 | ||
+ | |- | ||
+ | | π− photo-production || 199 || 342 | ||
+ | |- | ||
+ | | Total || 2.83317E7 || 7.87214E7 | ||
+ | |} | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <center><math>R_{events}= \sigma_{events} \mathcal{L}=7.87\times 10^{-2}\ barn\ \cdot \frac{1.33\times 10^{11}}{barn\cdot s}=\frac{1.05\times 10^{10}events}{s}</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>R_{Moller}= \sigma_{Moller} \mathcal{L}=7.50\times 10^{-2}\ barn\ \cdot \frac{1.33\times 10^{11}}{barn\cdot s}=\frac{9.97\times 10^{9}Moller}{s}</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | [[File:CLAS12ExpWS_WA.pdf]] | ||
+ | |||
+ | <center>[[File:Whitney_2.png]][[File:Whitney_3.png]]</center> | ||
+ | |||
+ | |||
+ | <center>[[File:Whitney_4.png]]</center> | ||
+ | |||
+ | |||
+ | <center><math>\sigma=\frac{R}{\mathcal{L}}=\frac{dN}{dt}\cdot \frac{1}{\mathcal{L}}</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>\Rightarrow \int\limits_{0}^{t}\, dt=\int\limits_{0}^{N}\frac{1}{\sigma \mathcal{L}}\, dN</math></center> | ||
+ | |||
+ | |||
+ | <center><math>t_{simulated}=\frac{N_{events}}{\sigma_{events} \Phi \rho \ell}=\frac{1000000\ barn \cdot s}{7.87\times 10^{-2} \cdot 1.33\times 10^{11}\ barn}=9.54\times 10^{-5} s</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <center><math>N_0=\Delta t \cdot R_{events}=\Delta t \cdot \frac{N_{events}}{t_{simulated}}=500\times 10^{-9}\ s \cdot \frac{1\times 10^6}{0.0954\times 10^{-3}\ s}=5240</math></center> | ||
+ | |||
+ | |||
+ | <center><math>Occupancy=\frac{N_{hits}}{N_0}=\frac{N_{hits}}{\Delta t \cdot R_{events}}=\frac{t_{simulated}\cdot N_{hits}}{N_{events}\cdot \Delta t}=</math></center> | ||
+ | |||
+ | Similarly, | ||
+ | |||
+ | |||
+ | <center><math>N_{events} =R_{events}\cdot t_{simulated}=\frac{1.05\times 10^{10}events}{s} \cdot 9.54\times 10^{-5} s=1001700\ events</math></center> | ||
+ | |||
+ | |||
+ | <center><math>N_{Moller} =R_{Moller}\cdot t_{simulated}=\frac{9.97\times 10^{9}Moller}{s} \cdot 9.54\times 10^{-5} s=951138\ Moller</math></center> | ||
+ | |||
+ | |||
+ | <center><math>\sigma=\frac{N_{events}}{N_{incident}}\frac{1}{\rho \ell}\Rightarrow N_{incident}=\frac{N_{events}}{\sigma \rho \ell}=\frac{951138}{.075\cdot 2.11\times 10^{-1}}=6\times 10^7\ incident </math></center> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>\left ( \frac{Number\ of\ hits}{Moller\ electron}\right ) \left (\frac{Moller\ electrons}{incidents\ electron} \right) \left (\frac{incident\ electrons}{sec} \right )</math></center> | ||
+ | |||
+ | ==GEANT Simulations== | ||
https://wiki.iac.isu.edu/index.php/Converting_to_barns | https://wiki.iac.isu.edu/index.php/Converting_to_barns | ||
Line 553: | Line 673: | ||
Converting the number of electrons scattered per angle theta to barns, we can use the relation | Converting the number of electrons scattered per angle theta to barns, we can use the relation | ||
− | <center><math>\mathcal{L}=\frac{ | + | <center><math>\mathcal{L}=\frac{R_{scattered}}{\sigma} = \Phi_{beam}\ \rho\ \ell </math></center> |
+ | |||
+ | |||
+ | |||
+ | <center><math>\sigma=\frac{R_{scattered}}{\mathcal{L}} =\frac{R_{scattered}}{\Phi_{beam}\ \rho\ \ell}= \frac{R_{scattered}}{R_{incident}\ \rho\ \ell}=\frac{N_{scattered}}{\Delta t}\cdot \frac{\Delta t}{N_{incident}}\cdot \frac{1}{\rho \ell}</math></center> | ||
+ | |||
+ | If the time is taken to be the same for the amount scattered as for the amount incident (the time simulated), this can be viewed as the probability of one incident electron producing a Moller event. | ||
− | <center><math>\sigma=\frac{N_{scattered}} | + | <center><math>\sigma=\frac{N_{scattered}}{N_{incident}}\frac{1}{\rho \ell}</math></center> |
− | + | While this expression has no explicit dependancies on energy, the ratio is a function of the energy, as well as the physical makeup of the target. | |
This gives, for LH2: | This gives, for LH2: | ||
− | <center><math>\rho_{target}\times l_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6\times10^{23} | + | <center><math>\rho_{target}\times l_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6\times10^{23} molecule}{1 mole} \times \frac{2\ atoms}{molecule}\times \frac{1m^3}{(100 cm)^3} \times \frac{1 cm}{ } \times \frac{10^{-24} cm^{2}}{barn} =4.2\times 10^{-2} barns^{-1}</math></center> |
+ | |||
Line 576: | Line 703: | ||
− | + | <center><math>\sigma = \frac{N_{events}}{N_{incident}\ \rho\ \ell}=\frac{975593}{40000000\ \cdot 4.2\times 10^{-2} barns^{-1}}=\frac{0.024}{4.2\times 10^{-2} barns^{-1}}=0.58 barns</math></center> | |
+ | |||
+ | |||
+ | <center><math>\sigma=\frac{R_{events}}{\mathcal{L}} \Rightarrow \mathcal{L}=\frac{R_{events}}{\sigma}</math></center> | ||
+ | |||
+ | |||
+ | <center><math>\mathcal{L}=\frac{dN_{events}}{dt}\frac{1}{ \sigma}\Rightarrow \int_{0}^{t_{simulated}}\mathcal {L}\, dt= \int_{0}^{N_{events}}\frac{1}{\sigma}\, dN</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>\mathcal{L} \cdot t_{simulated}=\frac{N_{events}}{\sigma}</math></center> | ||
+ | |||
+ | |||
+ | For a Luminosity of <math>\mathcal{L}=\frac{1.3\times 10^{11}}{barn\cdot s}</math> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>\frac{1.3\times 10^{11}}{barn\cdot s} \cdot t_{simulated}=\frac{975593}{.58\ barn}</math></center> | ||
+ | |||
+ | |||
+ | |||
+ | <center><math>t=1.3\times 10^{-5}\ s</math></center> | ||
+ | |||
+ | |||
+ | ===Comparison of GEANT Simulation to Whitney Data=== | ||
+ | |||
+ | |||
+ | {| class="wikitable" align="center" border=1 | ||
+ | |+ '''Rates''' | ||
+ | |- | ||
+ | | Simulation || GEANT4 || Whitney | ||
+ | |- | ||
+ | | Event Cross-Section (nb) || .61 || 0.079 | ||
+ | |- | ||
+ | |Moller Cross-Section (nb) || 0.58||0.075 | ||
+ | |- | ||
+ | | Length of Target (cm) || 1 || 5 | ||
+ | |- | ||
+ | | t_{simulated} (s) || 1.3E-5 ||9.54E-5 | ||
+ | |- | ||
+ | | N_{events} || 1026940 || 1001700 | ||
+ | |- | ||
+ | | N_{Moller} || 975593 ||951138 | ||
+ | |- | ||
+ | | N_{incident} || 4E7 || 6E7 | ||
+ | |} | ||
+ | |||
+ | |||
+ | Moller events occur for about 2.5% of the incident electrons on a LH2 target. We can assume the number of Moller events that occur within the DC range to be around 30% of the total Moller events occuring for the number of incident electrons for LH2 as well. Since the differential cross-section over the angel theta is proportional to the differtial cross-section over wire number we can dividing the Moller differential cross-section by the product of the density and length of the target material | ||
+ | |||
+ | ==Isotropic Weighted Distribution== | ||
+ | ===Isotropic in CM=== | ||
+ | Assuming a 1cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron, | ||
<pre> | <pre> | ||
TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); | TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); | ||
− | Combo->Add(Layer1binsWeighted, | + | Combo->Add(Layer1binsWeighted,23.8); |
Combo->Draw(); | Combo->Draw(); | ||
</pre> | </pre> | ||
Line 590: | Line 769: | ||
− | <center><math>\frac{ | + | <center><math>\frac{2003912\ Moller\ Events\ in\ Range}{.30\ Moller\ Events\ in\ Range\ per\ Moller\ Electron}=6.68\times 10^6\ Moller\ electrons</math></center> |
− | <center><math>\frac{ | + | <center><math>\frac{6.68\times 10^6\ Moller\ electrons}{.02\ Moller\ Events\ \ per\ Incident\ Electron}=3.34\times 10^8\ Incident\ electrons</math></center> |
− | |||
− | |||
− | + | <center><math>\sigma=\frac{N_{Moller}}{N_{incident}\cdot \rho \ell}=\frac{.02\ barn}{4.2\times 10^{-2}}=.48\ barn</math></center> | |
− | |||
− | + | <center><math>\mathcal{L}=N_{incident}\ \rho \ell=3.34\times 10^8\ Incident\ electrons \cdot 4.2\times 10^{-2}=\frac{1.40\times 10^{7}\ electrons}{barn \cdot s}</math></center> | |
− | |||
+ | Assuming a 5cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron, | ||
− | < | + | <pre> |
+ | TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); | ||
+ | Combo->Add(Layer1binsWeighted,4.76); | ||
+ | Combo->Draw(); | ||
+ | </pre> | ||
+ | <center>[[File:Layer1binsMollerRate_2.png]]</center> | ||
− | |||
− | + | <center><math>\frac{400782\ Moller\ Events\ in\ Range}{.30\ Moller\ Events\ in\ Range\ per\ Moller\ Electron}=1.34\times 10^6\ Moller\ electrons</math></center> | |
− | |||
+ | <center><math>\frac{1.34\times 10^6\ Moller\ electrons}{.02\ Moller\ Events\ \ per\ Incident\ Electron}=6.67\times 10^7\ Incident\ electrons</math></center> | ||
− | |||
− | <center><math> | + | <center><math>\sigma=\frac{N_{Moller}}{N_{incident}\cdot \rho \ell}=\frac{.02\ barn}{2.1\times 10^{-1}}=.095\ barn</math></center> |
− | <center><math>\ | + | <center><math>\mathcal{L}=N_{incident}\ \rho \ell=6.67\times 10^7\ Incident\ electrons \cdot 2.1\times 10^{-1}=\frac{1.40\times 10^{7}\ electrons}{barn \cdot s}</math></center> |
+ | ===Isotropic in Lab=== | ||
+ | Assuming a 1cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron, | ||
+ | <pre> | ||
+ | TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); | ||
+ | Combo->Add(Layer1binsWeighted,23.8); | ||
+ | Combo->Draw(); | ||
+ | </pre> | ||
− | <center> | + | <center>[[File:Layer1binsMollerRateIsotropic.png]]</center> |
+ | <center><math>\frac{833000\ Moller\ Events\ in\ Range}{.30\ Moller\ Events\ in\ Range\ per\ Moller\ Electron}=2.78\times 10^6\ Moller\ electrons</math></center> | ||
+ | |||
+ | |||
+ | <center><math>\frac{2.78\times 10^6\ Moller\ electrons}{.02\ Moller\ Events\ \ per\ Incident\ Electron}=1.39\times 10^8\ Incident\ electrons</math></center> | ||
− | |||
− | |||
+ | <center><math>\sigma=\frac{N_{Moller}}{N_{incident}\cdot \rho \ell}=\frac{.02\ barn}{4.2\times 10^{-2}}=.48\ barn</math></center> | ||
− | |||
− | |||
+ | <center><math>\mathcal{L}=N_{incident}\ \rho \ell=3.34\times 10^8\ Incident\ electrons \cdot 4.2\times 10^{-2}=\frac{5.83\times 10^{6}\ electrons}{barn \cdot s}</math></center> | ||
− | + | Assuming a 5cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron, | |
+ | <pre> | ||
+ | TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); | ||
+ | Combo->Add(Layer1binsWeighted,4.76); | ||
+ | Combo->Draw(); | ||
+ | </pre> | ||
− | + | <center>[[File:Layer1binsMollerRate_2Isotropic.png]]</center> | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
+ | <center><math>\frac{166600\ Moller\ Events\ in\ Range}{.30\ Moller\ Events\ in\ Range\ per\ Moller\ Electron}=5.55\times 10^5\ Moller\ electrons</math></center> | ||
+ | <center><math>\frac{5.55\times 10^5\ Moller\ electrons}{.02\ Moller\ Events\ \ per\ Incident\ Electron}=2.78\times 10^7\ Incident\ electrons</math></center> | ||
− | |||
+ | <center><math>\sigma=\frac{N_{Moller}}{N_{incident}\cdot \rho \ell}=\frac{.02\ barn}{2.1\times 10^{-1}}=.095\ barn</math></center> | ||
− | |||
− | <center><math>\ | + | <center><math>\mathcal{L}=N_{incident}\ \rho \ell=6.67\times 10^7\ Incident\ electrons \cdot 2.1\times 10^{-1}=\frac{5.83\times 10^{6}\ electrons}{barn \cdot s}</math></center> |
=[[DV_RunGroupC_Moller#Wire_angle_correspondance|back]]= | =[[DV_RunGroupC_Moller#Wire_angle_correspondance|back]]= |
Latest revision as of 07:27, 30 December 2016
Straight Track Moller Rates
Goal: Confirm that the moller Rate (Hz) seen in on a wire agrees with the Rate calculated from the cross section. The torus is turned off so the wire is more easily correlated with the moller scattering angle. No need for track reconstruction.
Reproduce the plot below using Torus off
https://wiki.iac.isu.edu/index.php/File:TheoryDCHits.png
Then calculate Rate and compare to the reference.
Determining wire-theta correspondence
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)
DC Geometry(geom)
![Wires parallel all types.png](/./images/3/39/Wires_parallel_all_types.png)
Examining the geometry file, we can see that each plane, not just the plane of the sense wires, is separated by a distance of D=.3861 cm. We can use the geometry of the wire placements to find
Finding the separation distance between two adjacent sense wires
![Hex.png](/./images/2/2c/Hex.png)
Since the separation between adjacent sense wires is uniform and at a set angle of 25 degrees with respect to the beam line, we can use this fact to determine the angle theta each wire makes when measured from the vertex.
![DC geom.png](/./images/c/ce/DC_geom.png)
The new x and z coordinates for wire 2 can be found using the change in the components
This can be extended to any point along the same wire plane, starting from the coordinates for wire 1
The angle theta that the wire makes with the vertex is given by
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.
![ParallelWireZoom.png](/./images/9/90/ParallelWireZoom.png)
Coordinates(cm) | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | Layer 6 |
---|---|---|---|---|---|---|
x | 20.42 | 21.51 | 21.39 | 22.49 | 22.37 | 23.47 |
y | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
z | 243.41 | 244.18 | 245.51 | 246.28 | 247.62 | 248.38 |
Coordinates(cm) | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | Layer 6 |
---|---|---|---|---|---|---|
x | 21.64 | 22.73 | 22.61 | 23.71 | 23.58 | 24.69 |
y | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
z | 242.86 | 243.62 | 244.94 | 245.71 | 247.05 | 247.81 |
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 |
62 | 24.42 | |||||
63 | 24.76 | |||||
64 | 25.10 | |||||
77 | 29.46 | 29.60 | 29.41 | 29.55 | 29.37 | 29.51 |
78 | 29.79 | 29.93 | 29.74 | 29.88 | 29.69 | 29.83 |
79 | 30.12 | 30.26 | 30.07 | 30.21 | 30.02 | 30.16 |
110 | 40.19 | 40.28 | 40.05 | 40.13 | 39.91 | 39.99 |
111 | 40.51 | 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 |
77 | 29.46 | 29.60 | 29.41 | 29.56 | 29.37 | 29.51 |
78 | 29.79 | 29.93 | 29.74 | 29.88 | 29.70 | 29.84 |
79 | 30.13 | 30.27 | 30.07 | 30.21 | 30.02 | 30.16 |
110x | 40.19 | 40.28 | 40.05 | 40.13 | 39.91 | 39.99 |
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
Coordinates(cm) | Layer 1 | Layer 2 | Layer 3 | Layer 4 | Layer 5 | Layer 6 |
---|---|---|---|---|---|---|
1.22 | 1.22 | 1.22 | 1.22 | 1.21 | 1.22 | |
0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | |
.55 | .56 | .57 | .57 | .57 | .57 |
Using the geometric construction for determining angle theta to wire 2:
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.
where
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
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
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,
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 .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,
Similarly, finding the difference between wires 111 and 112,
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
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
As discussed earlier, taking the 1st derivative of this function will give us the spacing of the bins as a function of wire number.
The derivative of this function then will tell us where the bin spacing is at a minimum and a maximum.
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
![Layer1bins.png](/./images/f/f1/Layer1bins.png)
![Layer1bins Isotropic.png](/./images/a/ae/Layer1bins_Isotropic.png)
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)
![MolThetaLabWeighted.png](/./images/8/84/MolThetaLabWeighted.png)
![Layer1binWeighted Isotropic.png](/./images/c/ca/Layer1binWeighted_Isotropic.png)
![MolThetaLabWeighted2.png](/./images/0/06/MolThetaLabWeighted2.png)
Using the expression for n in terms of Theta:
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)
![MolThetaWireWeightedIsotropic.png](/./images/9/94/MolThetaWireWeightedIsotropic.png)
![MolThetaWireWeightedAdjusted.png](/./images/7/72/MolThetaWireWeightedAdjusted.png)
![MolThetaWireWeightedAdjustedIsotropic.png](/./images/9/9c/MolThetaWireWeightedAdjustedIsotropic.png)
![TheoryDCbinsWire.png](/./images/a/aa/TheoryDCbinsWire.png)
![TheoryDCbinsWireIsotropic.png](/./images/0/08/TheoryDCbinsWireIsotropic.png)
Moller Occupancy Rate
Comparing Rates
Using the values from Whitney:
![Whitney values.png](/./images/8/80/Whitney_values.png)
For a 5cm target of LH2:
Using the beam current of 100nA,
Given the beam Luminosity of:
We can check to make sure the density makes sense
Using the values from Whitney :
https://wiki.iac.isu.edu/index.php/CLAS12_RateEst_byWA
Energy | 6 GeV | 11 GeV |
Process | (nb) | (nb) |
Moller | 22773001 | 75008636 |
DIS + radiative tail | 128 | 83 |
Elastic e-p | 5511220 | 3670740 |
Elastic radiative tail | 24705 | 12944 |
π0 electro-production | 14802 | 17908 |
π0 photo-production | 569 | 852 |
π+ electro-production | 4032 | 5536 |
π+ photo-production | 282 | 487 |
π− electro-production | 2806 | 3843 |
π− photo-production | 199 | 342 |
Total | 2.83317E7 | 7.87214E7 |
![Whitney 2.png](/./images/3/37/Whitney_2.png)
![Whitney 3.png](/./images/3/30/Whitney_3.png)
![Whitney 4.png](/./images/2/22/Whitney_4.png)
Similarly,
GEANT Simulations
https://wiki.iac.isu.edu/index.php/Converting_to_barns
https://wiki.iac.isu.edu/index.php/Check_Differential_Cross-Section
Converting the number of electrons scattered per angle theta to barns, we can use the relation
If the time is taken to be the same for the amount scattered as for the amount incident (the time simulated), this can be viewed as the probability of one incident electron producing a Moller event.
While this expression has no explicit dependancies on energy, the ratio is a function of the energy, as well as the physical makeup of the target.
This gives, for LH2:
From earlier simulations for random angle Phi, we know that the full range of Theta is limited depending on the target material.
![MollerThetaLab 4e7 LH2 11GeV.png](/./images/a/aa/MollerThetaLab_4e7_LH2_11GeV.png)
![MollerThetaLab 4e7 NH3 11GeV.png](/./images/1/14/MollerThetaLab_4e7_NH3_11GeV.png)
![MollerThetaLab 4e7 LH2 11GeV Detector.png](/./images/d/d8/MollerThetaLab_4e7_LH2_11GeV_Detector.png)
![MollerThetaLab 4e7 NH3 11GeV Detector.png](/./images/e/ee/MollerThetaLab_4e7_NH3_11GeV_Detector.png)
For a Luminosity of
Comparison of GEANT Simulation to Whitney Data
Simulation | GEANT4 | Whitney |
Event Cross-Section (nb) | .61 | 0.079 |
Moller Cross-Section (nb) | 0.58 | 0.075 |
Length of Target (cm) | 1 | 5 |
t_{simulated} (s) | 1.3E-5 | 9.54E-5 |
N_{events} | 1026940 | 1001700 |
N_{Moller} | 975593 | 951138 |
N_{incident} | 4E7 | 6E7 |
Moller events occur for about 2.5% of the incident electrons on a LH2 target. We can assume the number of Moller events that occur within the DC range to be around 30% of the total Moller events occuring for the number of incident electrons for LH2 as well. Since the differential cross-section over the angel theta is proportional to the differtial cross-section over wire number we can dividing the Moller differential cross-section by the product of the density and length of the target material
Isotropic Weighted Distribution
Isotropic in CM
Assuming a 1cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron,
TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); Combo->Add(Layer1binsWeighted,23.8); Combo->Draw();
![Layer1binsMollerRate.png](/./images/6/67/Layer1binsMollerRate.png)
Assuming a 5cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron,
TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); Combo->Add(Layer1binsWeighted,4.76); Combo->Draw();
![Layer1binsMollerRate 2.png](/./images/7/70/Layer1binsMollerRate_2.png)
Isotropic in Lab
Assuming a 1cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron,
TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); Combo->Add(Layer1binsWeighted,23.8); Combo->Draw();
![Layer1binsMollerRateIsotropic.png](/./images/0/0a/Layer1binsMollerRateIsotropic.png)
Assuming a 5cm long target of LH2 where around 2% of the incident electrons will result in a Moller electron,
TH1F *Combo=new TH1F("Layer1binsMollerRate","Layer1binsMollerRate",113,0,113); Combo->Add(Layer1binsWeighted,4.76); Combo->Draw();
![Layer1binsMollerRate 2Isotropic.png](/./images/6/6a/Layer1binsMollerRate_2Isotropic.png)