Difference between revisions of "Verification of Relativistic Components"
| (4 intermediate revisions by the same user not shown) | |||
| Line 15: | Line 15: | ||
<center>[[File:FnlTheta_reverse.png]][[File:MolTheta_reverse.png]]</center> | <center>[[File:FnlTheta_reverse.png]][[File:MolTheta_reverse.png]]</center> | ||
| + | =Theta in the CM Frame= | ||
| + | <center>[[File:FnlThetaCM_reverse.png]][[File:MolThetaCM_reverse.png]]</center> | ||
=Momentum in the CM Frame= | =Momentum in the CM Frame= | ||
| Line 51: | Line 53: | ||
| 0 | | 0 | ||
| -510 | | -510 | ||
| − | | | + | | 10931.3 |
| − | | | + | | -7.51848 |
| − | | - | + | | -3.62463 |
| − | | | + | | 10931.8 |
| 0 | | 0 | ||
| 0 | | 0 | ||
| − | | -509. | + | | -509.361 |
| − | | | + | | 68.5938 |
| − | | | + | | 7.51848 |
| − | | | + | | 3.62463 |
| − | | | + | | 68.597 |
| 0 | | 0 | ||
| 0 | | 0 | ||
| − | | -509. | + | | -509.361 |
|} | |} | ||
| − | The particles start out with p=11000.511 | + | The particles start out with p=11000.511 MeV and have on average ~11000 MeV by the time they interact |
<center>[[File:Screen_Shot_2016-02-01_at_9.30.23_PM.png]]</center> | <center>[[File:Screen_Shot_2016-02-01_at_9.30.23_PM.png]]</center> | ||
| Line 74: | Line 76: | ||
<center>[[File:Screen_Shot_2016-02-01_at_9.29.40_PM.png]]</center> | <center>[[File:Screen_Shot_2016-02-01_at_9.29.40_PM.png]]</center> | ||
| − | = | + | |
| − | + | ||
| + | =Conclusion= | ||
| + | |||
| + | The energy loss must be dealt with by using both momentum components from a GEANT Moller Scattering data file. By using both components, they can be transformed to a CM Frame, then have the phi angle varied to cover 2π radians. A quick awk search shows that there are no Moller electrons with energy over 5500 MeV. This might be able to be fixed by switching the scattered electron with the Moller electron, allowing a larger range of energy to be covered. | ||
Latest revision as of 15:56, 4 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
Final Thetas in the Lab Frame
Theta in the CM Frame
Momentum in the CM Frame
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 | 10931.3 | -7.51848 | -3.62463 | 10931.8 | 0 | 0 | -509.361 | 68.5938 | 7.51848 | 3.62463 | 68.597 | 0 | 0 | -509.361 |
The particles start out with p=11000.511 MeV and have on average ~11000 MeV by the time they interact
Conclusion
The energy loss must be dealt with by using both momentum components from a GEANT Moller Scattering data file. By using both components, they can be transformed to a CM Frame, then have the phi angle varied to cover 2π radians. A quick awk search shows that there are no Moller electrons with energy over 5500 MeV. This might be able to be fixed by switching the scattered electron with the Moller electron, allowing a larger range of energy to be covered.