Computing Kinematic Arrays

To generate arrays of kinematic quantities over all center-of-mass angles, use kinematics_table_at_beam_energy().

from reaction_kinematics import Reaction

rxn = Reaction("p", "3H", "n", "3He")
# Proton + Tritium Reaction

data = rxn.kinematics_table_at_beam_energy(1.2)

print(data["theta4"])
print(data["e3"])

This will return a dictionary containing the following:

  • coscm : cos(θ_CM)
  • theta_cm: CM angle (rad)
  • theta3 : Ejectile lab angle (rad)
  • theta4 : Recoil lab angle (rad)
  • e3 : Ejectile energy (MeV)
  • e4 : Recoil energy (MeV)
  • v3 : Ejectile velocity (c)
  • v4 : Recoil velocity (c)