Decay12_AfterBoost_New

p = Plane(y = 0)
p_{3} = Vector((0, 0, 0), (-1, -1, 5))
p_{1} = Vector((0, 0, 0), (-3, -3, -6))
p_{2}: (-p_{1})
Notθ_{1} = Angle(Vector[(0, 0, 1)], p_{1})
θ_{2} = Angle(p_{2}, Vector[UnitVector[(0, 0, 1)]])
A = (-10.619236277471119, -3.697387906902769, 3.422626759704169)
p_{a} = Vector((0, 0, 0), A)
B = (2.47978, -2.72225, 4.0805396093259025)
p_{b} = Vector((0, 0, 0), B)
u: p_{a} ⊗ p_{b}
productionplane = OrthogonalPlane((0, 0, 0), Vector[p_{a} ⊗ p_{b}])
decayplane = OrthogonalPlane((0, 0), Vector[p_{1} ⊗ p_{3}])
ϕ_{1} = Angle(decayplane, productionplane)
a = Length(p_{1})
b = Length(p_{2})
c = Length(p_{3})
mp_{3}: (-p_{3})
θ_{1} = Angle(p_{1}, mp_{3})