Difference between revisions of "Verification of Relativistic Components"

From New IAC Wiki
Jump to navigation Jump to search
Line 67: Line 67:
 
|}
 
|}
  
The particles start out with p=11000.511 GeV and have on average ~11000 GeV by the time they interact
+
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>

Revision as of 04:51, 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.pngMolMom reverse.png

Final Thetas in the Lab Frame

FnlTheta reverse.pngMolTheta reverse.png


Momentum in the CM Frame

FnlMomCM reverse.pngMolMomCM 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 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

Screen Shot 2016-02-01 at 9.30.23 PM.png


Screen Shot 2016-02-01 at 9.29.40 PM.png

Theta in the CM Frame

FnlThetaCM reverse.pngMolThetaCM reverse.png