Difference between revisions of "Verification of Relativistic Components"
Jump to navigation
Jump to search
![FnlMom backwards.png](/./images/d/df/FnlMom_backwards.png)
![MolMom reverse.png](/./images/6/69/MolMom_reverse.png)
![FnlTheta reverse.png](/./images/3/3f/FnlTheta_reverse.png)
![MolTheta reverse.png](/./images/1/13/MolTheta_reverse.png)
![FnlMomCM reverse.png](/./images/1/10/FnlMomCM_reverse.png)
![MolMomCM reverse.png](/./images/3/3e/MolMomCM_reverse.png)
![FnlThetaCM reverse.png](/./images/f/f7/FnlThetaCM_reverse.png)
![MolThetaCM reverse.png](/./images/e/ec/MolThetaCM_reverse.png)
(→Setup) |
|||
Line 20: | Line 20: | ||
The differences in these plots come from fact that the "Forward" plots come from a GEANT simulation in which a small amount of energy was lost traveling through the target material. | The differences in these plots come from fact that the "Forward" plots come from a GEANT simulation in which a small amount of energy was lost traveling through the target material. | ||
− | + | {| border="1" | |
− | + | |- | |
+ | ! KE<sub>i</sub> | ||
+ | ! Px<sub>i</sub> | ||
+ | ! Py<sub>i</sub> | ||
+ | ! Pz<sub>i</sub> | ||
+ | ! x<sub>i</sub> | ||
+ | ! y<sub>i</sub> | ||
+ | ! z<sub>1</sub> | ||
+ | ! KE<sub>f</sub> | ||
+ | ! Px<sub>f</sub> | ||
+ | ! Py<sub>f</sub> | ||
+ | ! Pz<sub>f</sub> | ||
+ | ! x<sub>f</sub> | ||
+ | ! y<sub>f</sub> | ||
+ | ! z<sub>f</sub> | ||
+ | ! KE<sub>m</sub> | ||
+ | ! Px<sub>m</sub> | ||
+ | ! Py<sub>m</sub> | ||
+ | ! Pz<sub>m</sub> | ||
+ | ! x<sub>m</sub> | ||
+ | ! y<sub>m</sub> | ||
+ | ! z<sub>m</sub> | ||
+ | |- | ||
+ | | 11000 | ||
+ | | 0 | ||
+ | | 0 | ||
+ | | 11000.5 | ||
+ | | 0 | ||
+ | | 0 | ||
+ | | -510 | ||
+ | | 10999.1 | ||
+ | | 0.433025 | ||
+ | | -0.858867 | ||
+ | | 10999.6 | ||
+ | | 0 | ||
+ | | 0 | ||
+ | | -509.276 | ||
+ | | 0.905324 | ||
+ | | -0.433025 | ||
+ | | 0.858867 | ||
+ | | 0.905366 | ||
+ | | 0 | ||
+ | | 0 | ||
+ | | -509.276 | ||
+ | |} | ||
=Theta in the CM Frame= | =Theta in the CM Frame= | ||
<center>[[File:FnlThetaCM_reverse.png]][[File:MolThetaCM_reverse.png]]</center> | <center>[[File:FnlThetaCM_reverse.png]][[File:MolThetaCM_reverse.png]]</center> |
Revision as of 04:16, 2 February 2016
Setup
Verifying that the kinematics are valid, the Moller momentum from GEANT simulations are extracted using
awk '{print $16, $17, $18}' MollerScattering_NH3_Large.dat > Just_Mol_Mom.dat
Then using the kinematics equations displayed above, used to reverse calculate the scattered electron's momentum and angles. The root histograms of an incoming 11GeV electron (in blue) are overlaid with the relativistic calculations starting with the Moller electron (Reverse Moller in red) and show good agreement.
Final Momentum in the Lab Frame
![FnlMom backwards.png](/./images/d/df/FnlMom_backwards.png)
![MolMom reverse.png](/./images/6/69/MolMom_reverse.png)
Final Thetas in the Lab Frame
![FnlTheta reverse.png](/./images/3/3f/FnlTheta_reverse.png)
![MolTheta reverse.png](/./images/1/13/MolTheta_reverse.png)
Momentum in the CM Frame
![FnlMomCM reverse.png](/./images/1/10/FnlMomCM_reverse.png)
![MolMomCM reverse.png](/./images/3/3e/MolMomCM_reverse.png)
The differences in these plots come from fact that the "Forward" plots come from a GEANT simulation in which a small amount of energy was lost traveling through the target material.
KEi | Pxi | Pyi | Pzi | xi | yi | z1 | KEf | Pxf | Pyf | Pzf | xf | yf | zf | KEm | Pxm | Pym | Pzm | xm | ym | zm |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
11000 | 0 | 0 | 11000.5 | 0 | 0 | -510 | 10999.1 | 0.433025 | -0.858867 | 10999.6 | 0 | 0 | -509.276 | 0.905324 | -0.433025 | 0.858867 | 0.905366 | 0 | 0 | -509.276 |
Theta in the CM Frame
![FnlThetaCM reverse.png](/./images/f/f7/FnlThetaCM_reverse.png)
![MolThetaCM reverse.png](/./images/e/ec/MolThetaCM_reverse.png)