Ne21 = AtomicMassTable.GetElement(10, 21)
H2 = AtomicMassTable.GetElement(1, 2)
Ne20_mass = Ne20[3]
H1_mass = H1[3]
Ne21_mass = Ne21[3]
H2_mass = H2[3]
Kr = KINEMATICS()
Kr.m[0] = Ne20_mass
Kr.m[1] = H2_mass
Kr.m[2] = Ne21_mass
Kr.m[3] = H1_mass
Kr.K0 = 100
Kr.Eex2 = 0
Kr.Eex3 = 0
lab2 = []
K2 = []
def ReturnLines(theta=30):
degree = math.pi / 180
Kr.calculate(theta * degree, 0)
MagneticFieldB = 2
tr2 = HELIOS.trajectory(Kr.V2, Kr.thetalab2 * degree, 10, Kr.m[2], MagneticFieldB)
tr3 = HELIOS.trajectory(Kr.V3, Kr.thetalab3 * degree, 1, Kr.m[3], MagneticFieldB)
return tr2, tr3
