|
|
| Line 1: |
Line 1: |
| − | =[[File:LUND_Spread.C|Creating uniform LUND files]]=
| + | [[File:LUND_Spread.C|Creating uniform LUND files]] |
| | | | |
| | The LUND file is created by creating an isotropic distribution of particles within the Moller Center of Mass frame of reference. These particles are uniformly distributed through the angle theta with respect to the beam line in the range 90-180 in the Center of Mass frame. This is done at a set angle phi (10 degrees) with respect to the perpendicular components with respect to the beam line. | | The LUND file is created by creating an isotropic distribution of particles within the Moller Center of Mass frame of reference. These particles are uniformly distributed through the angle theta with respect to the beam line in the range 90-180 in the Center of Mass frame. This is done at a set angle phi (10 degrees) with respect to the perpendicular components with respect to the beam line. |
| | | | |
| | <center>[[File:MolThetaCM_spread.png]][[File:MolMomCM_spread.png]]</center> | | <center>[[File:MolThetaCM_spread.png]][[File:MolMomCM_spread.png]]</center> |
| − |
| |
| − |
| |
| − | ==Center of Mass for Stationary Target==
| |
| − | For an [[DV_Calculations_of_4-momentum_components#Center_of_Mass_Frame|incoming electron of 11GeV striking a stationary electron]] we would expect:
| |
| − |
| |
| − |
| |
| − | <center><math>\Longrightarrow \left( \begin{matrix} E^*_{2} \\ p^*_{2(x)} \\ p^*_{2(y)} \\ p^*_{2(z)}\end{matrix} \right)=\left(\begin{matrix}\gamma^* & 0 & 0 & -\beta^* \gamma^*\\0 & 1 & 0 & 0 \\ 0 & 0 & 1 &0 \\ -\beta^* \gamma^* & 0 & 0 & \gamma^* \end{matrix} \right) . \left( \begin{matrix}m\\ 0 \\ 0 \\ 0\end{matrix} \right)</math></center>
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | <center>[[File:Init_e_Mom_CM.png|600 px]][[File:Init_Mol_Mom_CM.png|600 px]]</center>
| |
| − |
| |
| − |
| |
| − |
| |
| − | <center>[[File:Init_e_Theta_CM.png|600 px]][[File:Init_Mol_Theta_CM.png|600 px]]</center>
| |
| − |
| |
| − | ==Phase space Limiting Particles==
| |
| − | Since the angle phi has been constrained to remain constant, the x and y components of the momentum will increase in the positive first quadrant. This implies that the z component of the momentum must decrease by the relation:
| |
| − |
| |
| − | <center><math>p^2=p_x^2+p_y^2+p_z^2</math></center>
| |
| − |
| |
| − | In the Center of Mass frame, this becomes:
| |
| − |
| |
| − | <center><math>p_x^{*2}+p_y^{*2} = p^{*2}-p_z^{*2}</math></center>
| |
| − |
| |
| − |
| |
| − | Since the momentum in the CM frame is a constant, this implies that pz must decrease. We can use the variable rapidity:
| |
| − |
| |
| − | <center><math>y \equiv \frac {1}{2} \ln \left(\frac{E+p_z}{E-p_z}\right)</math></center>
| |
| − |
| |
| − | where
| |
| − |
| |
| − | <center><math> P^+ \equiv E+p_z</math></center>
| |
| − |
| |
| − | <center><math> P^- \equiv E-p_z</math></center>
| |
| − |
| |
| − | this implies that as
| |
| − | <center><math>p_z \rightarrow 0 \Rrightarrow \frac{E+p_z}{E-p_z} \rightarrow 1 \Rrightarrow \ln 1 \rightarrow 0 \Rrightarrow y=0</math></center>
| |
| − |
| |
| − |
| |
| − | For forward travel in the light cone:
| |
| − | <center><math>p_z \rightarrow E \Rrightarrow \frac{E+p_z}{E-p_z} \rightarrow \infin \Rrightarrow \ln \infin \rightarrow \infin \Rrightarrow y \rightarrow \infin </math></center>
| |
| − |
| |
| − | This corresponds to the scattered electron proven [[DV_Calculations_of_4-momentum_components|earlier]].
| |
| − |
| |
| − | For backward travel in the light cone:
| |
| − | <center><math>p_z \rightarrow -E \Rrightarrow \frac{E+p_z}{E-p_z} \rightarrow 0 \Rrightarrow \ln 0 \rightarrow -\infin \Rrightarrow y \rightarrow -\infin </math></center>
| |
| − |
| |
| − |
| |
| − | Similarly, this corresponds to the Moller electron.
| |
| − |
| |
| − |
| |
| − | For a particle that transforms from the Lab frame to the CM frame where the particle is not within the light cone:
| |
| − |
| |
| − | <center><math>p_x^2+p_y^2=52.589054^2+9.272868^2=53.400MeV > 53.015 MeV (E) \therefore p_z \rightarrow imaginary</math></center>
| |
| − |
| |
| − | These particles are outside the light cone and are more timelike, thus not visible in normal space. This will reduce the number of particles that will be detected.
| |
| − |
| |
| − | ==Determining Momentum Components After Collision in CM Frame==
| |
| − | The energy and total momentum of the Moller and scattered electron remain the same under a rotation in any frame of reference. After the collision, these quantities remain the same, but the x, y, z components change.
| |
| − |
| |
| − | In this frame, we can cycle through values of theta from 90 to 180 degrees which physically correspond to a stationary electron being impinged by an electron.
| |
| − |
| |
| − |
| |
| − | ===Theta Dependent Components===
| |
| − | <center>[[File:xz_lab.png | 400 px]]</center>
| |
| − | <center>'''Figure 3: Definition of Moller electron variables in the CM Frame in the x-z plane.'''</center>
| |
| − |
| |
| − | <center>Using <math>\theta '_2=\arccos \left(\frac{p^'_{2(z)}}{p^'_{2}}\right)</math></center>
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>\Longrightarrow {p^'_{2(z)}=p^'_{2}\cos(\theta '_2)}</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | Checking on the sign resulting from the cosine function, we are limited to:
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: red" | <math>90^\circ \le \theta '_2 \le 180^\circ \equiv \frac{\pi}{2} \le \theta '_2 \le \pi Radians</math>
| |
| − | |}
| |
| − |
| |
| − | Since,
| |
| − | <center><math>\frac{p^'_{2(z)}}{p^'_{2}}=cos(\theta '_2)</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\Longrightarrow p^'_{2(z)}\ should\ always\ be\ negative</math></center>
| |
| − |
| |
| − |
| |
| − | ===Phi Dependent Components===
| |
| − | Since only the z direction is considered to be the relativistic direction of motion, this implies that the x and y components are not effected by a Lorentz transformation and remain the same in the CM and Lab frame. Holding the angle Phi constant at an initial value of 10 degrees, allows us to find the x and y components.
| |
| − |
| |
| − | <center>[[File:xy_lab.png | 400 px]]</center>
| |
| − | <center>'''Figure 4: Definition of Moller electron variables in the Lab Frame in the x-y plane.'''</center>
| |
| − |
| |
| − | <center>Similarly, <math>\phi '_2=\arccos \left( \frac{p^'_{2(x) Lab}}{p^'_{2(xy)}} \right)</math></center>
| |
| − |
| |
| − |
| |
| − | <center>where <math>p_{2(xy)}^'=\sqrt{(p_{2(x)}^')^2+(p^'_{2(y)})^2}</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>(p^'_{2(xy)})^2=(p^'_{2(x)})^2+(p^'_{2(y)})^2</math></center>
| |
| − |
| |
| − |
| |
| − | <center>and using <math>p^2=p_{(x)}^2+p_{(y)}^2+p_{(z)}^2</math></center>
| |
| − |
| |
| − |
| |
| − | <center>this gives <math>(p^'_{2})^2=(p^'_{2(xy)})^2+(p^'_{2(z)})^2</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\Longrightarrow (p'_{2})^2-(p'_{2(z)})^2 = (p'_{2(xy)})^2</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\Longrightarrow p_{2(xy)}^'=\sqrt{(p^'_{2})^2-(p^'_{2(z)})^2}</math></center>
| |
| − |
| |
| − |
| |
| − | <center>which gives<math>\phi '_2 = \arccos \left( \frac{p_{2(x)}'}{\sqrt{p_{2}^{'\ 2}-p_{2(z)}^{'\ 2}}}\right)</math></center>
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>\Longrightarrow p_{2(x)}'=\sqrt{p_{2}^{'\ 2}-p_{2(z)}^{'\ 2}} \cos(\phi)</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − | <center>Similarly, using <math>p_{2}^2=p_{2(x)}^2+p_{2(y)}^2+p_{2(z)}^2</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\Longrightarrow p_{2}^{'\ 2}-p_{2(x)}^{'\ 2}-p_{2(z)}^{'\ 2}=p_{2(y)}^{'\ 2}</math></center>
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>p_{2(y)}'=\sqrt{p_{2}^{'\ 2}-p_{2(x)}^{'\ 2}-p_{2(z)}^{'\ 2}}</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − | Checking on the sign from the cosine results for <math>\phi '_2</math>
| |
| − |
| |
| − |
| |
| − | We have the limiting range that <math>\phi</math> must fall within:
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: grey" | <math>-\pi \le \phi '_2 \le \pi\ Radians</math>
| |
| − | |}
| |
| − |
| |
| − | <center>[[File:xy_plane.png | 400px]]</center>
| |
| − |
| |
| − | Examining the signs of the components which make up the angle <math>\phi</math> in the 4 quadrants which make up the xy plane:
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: red" | <math>For\ 0 \ge \phi '_2 \ge \frac{-\pi}{2}\ Radians</math>
| |
| − | |-
| |
| − | | <center>p<sub>x</sub>=POSITIVE</center>
| |
| − | |-
| |
| − | | <center>p<sub>y</sub>=NEGATIVE</center>
| |
| − | |}
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: red" | <math>For\ 0 \le \phi '_2 \le \frac{\pi}{2}\ Radians</math>
| |
| − | |-
| |
| − | | <center>p<sub>x</sub>=POSITIVE</center>
| |
| − | |-
| |
| − | | <center>p<sub>y</sub>=POSITIVE</center>
| |
| − | |}
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: red" | <math>For\ \frac{-\pi}{2} \ge \phi '_2 \ge -\pi\ Radians</math>
| |
| − | |-
| |
| − | | <center>p<sub>x</sub>=NEGATIVE</center>
| |
| − | |-
| |
| − | | <center>p<sub>y</sub>=NEGATIVE</center>
| |
| − | |}
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: red" | <math>For\ \frac{\pi}{2} \le \phi '_2 \le \pi\ Radians</math>
| |
| − | |-
| |
| − | |<center>p<sub>x</sub>=NEGATIVE</center>
| |
| − | |-
| |
| − | | <center>p<sub>y</sub>=POSITIVE</center>
| |
| − | |}
| |
| − |
| |
| − | ===Electron Center of Mass Frame===
| |
| − | Relativistically, the x, y, and z components have the same magnitude, but opposite direction, in the conversion from the Moller electron's Center of Mass frame to the electron's Center of Mass frame.
| |
| − |
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>p^*_{2(x)}= -p^*_{1(x)}</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>p^*_{2(y)}= -p^*_{1(y)}</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>p^*_{2(z)}=-p^*_{1(z)}</math>
| |
| − | |}
| |
| − |
| |
| − | where previously it was shown
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>|p^*_{1}|\equiv |p^*_{2}|</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>E^*_{1}\equiv E^*_{2}</math>
| |
| − | |}
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | {| class="wikitable" align="center"
| |
| − | | style="background: gray" | <math>\theta_{1}^*=\pi-\theta_{1}^*</math>
| |
| − | |}
| |
| − |
| |
| − | ==Lorentz Transformation to Lab Frame==
| |
| − |
| |
| − | <center><math>
| |
| − | \begin{bmatrix}
| |
| − | E_{Lab} \\
| |
| − | P_{x(Lab)} \\
| |
| − | P_{y(Lab)} \\
| |
| − | P_{z(Lab)}
| |
| − | \end{bmatrix}=
| |
| − | \begin{bmatrix}
| |
| − | \gamma & 0 & 0 & \gamma \beta \\
| |
| − | 0 & 1 & 0 & 0 \\
| |
| − | 0 & 0 & 1 & 0 \\
| |
| − | \gamma \beta & 0 & 0 & \gamma
| |
| − | \end{bmatrix}
| |
| − | \cdot
| |
| − | \begin{bmatrix}
| |
| − | E_{CM} \\
| |
| − | P_{x(CM)} \\
| |
| − | P_{y(CM)} \\
| |
| − | P_{z(CM)}
| |
| − | \end{bmatrix}
| |
| − | </math></center>
| |
| − |
| |
| − | ==Weight==
| |
| − |
| |
| − |
| |
| − | [[Converting_to_barns| Using the theoretical differential cross section from previous]]
| |
| − | <center><math>\frac{d\sigma}{d\Omega}=\frac{ \alpha^2 }{4E^2}\frac{ (3+cos^2\theta)^2}{sin^4\theta}</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\alpha ^2=5.3279\times 10^{-5}</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>E\approx 53.013 MeV</math></center>
| |
| − |
| |
| − | We can take the Moller electron distribution of Theta in the Center of Mass frame, and multiply each given angle Theta by the expected differential cross section.
| |
| − |
| |
| − | <center>[[File:MolThetaCM_spread.png]]</center>
| |
| − |
| |
| − |
| |
| − | This causes the Moller Theta distribution in the Center of Mass frame to directly follow the theoretical differential cross section.
| |
| − |
| |
| − | <center>[[File:MolThetaCMWeightedTheory_LH2_11GeV.png]]</center>
| |
| − |
| |
| − |
| |
| − | The Lab frame distribution of Theta can also be weighted similarly. However, instead of having it be a differential cross section, [[Converting_to_barns|we can find the necessary number of particles.]]
| |
| − |
| |
| − | <center><math>\frac{d\sigma}{d\Omega} = \frac{\left(\frac{number\ of\ particles\ scattered/second}{d\Omega}\right)}{\left(\frac{number\ of\ incoming\ particles/second}{cm^2}\right)}=\frac{dN}{\mathcal L\, d\Omega} =differential\ scattering\ cross\ section</math></center>
| |
| − |
| |
| − |
| |
| − |
| |
| − | <center><math>\frac{d\sigma}{d\Omega} =\frac{dN}{\mathcal L\, d\Omega} =\frac{dN}{\Phi \rho \ell\, d\Omega}</math></center>
| |
| − |
| |
| − |
| |
| − | <center><math>\Rightarrow \int \frac{d\sigma}{d\Omega} \rho\ \ell \Phi d\Omega=\int dN</math></center>
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | Checking versus the given Luminosity for the experiment
| |
| − |
| |
| − |
| |
| − | <center><math>\mathcal L = \frac{1.33 \times 10^{35}}{cm^2\cdot s} \times \frac{10^{-24} cm^2}{barn}=1.33\times 10^{11} barns^{-1}s^{-1}</math></center>
| |
| − |
| |
| − |
| |
| − |
| |
| − | <center><math>\mathcal L \int \frac{d\sigma}{d\Omega} d\Omega = 1.33 \times 10^{11} barns^{-1} \times 2\pi\ \int \frac{d\sigma}{d\Omega}\ \sin{\theta} d\theta=N</math></center>
| |
| − |
| |
| − |
| |
| − | Limiting the range of Theta to within 5 to 40 degrees in the Lab frame:
| |
| − |
| |
| − | <center>[[File:DetectorRangeLab.png]]</center>
| |
| − |
| |
| − |
| |
| − |
| |
| − | We can find the corresponding angular range in the CM frame:
| |
| − | <center>[[File:DetectorRangeCM.png]]</center>
| |
| − |
| |
| − |
| |
| − |
| |
| − | Integrating the differential cross-section over the solid angle:
| |
| − | <center>[[File:DetectorIntegrationXSect_LH2_11GeV.png]]</center>
| |
| − | <pre>
| |
| − | IntegralDiffXSect->Integral()
| |
| − | </pre>
| |
| − |
| |
| − | 4.97493824519086629e+04 barns
| |
| − |
| |
| − |
| |
| − | Multiplying by the Luminosity, we find:
| |
| − |
| |
| − | <center><math>4.97\times 10^{4} barns \times\frac{1.35\times 10^{11}}{barns\cdot s}=6.71\times 10^{15}\frac{electrons}{s}</math></center>
| |
| − |
| |
| − | =Links=
| |
| − |
| |
| − | [[Run_in_GEMC#Cross_Section| Back]]
| |
