Calculate alignment angles of different subsystems

[spflueger]

Once different subsystem have been identified between two topologies, the alignment angles for a Rotation that connects one with the other have to be computed from the four momentum information. So a function of the form

def calculate_alignment_angles(topology: Topology, reference: Topology) -> Tuple[float, float, float]:
    pass

has to be implemented. The return value are the three euler angles that "align" topology with reference.

See p.6, just before Eq. (37) of Marangotto, “Helicity Amplitudes for Generic Multi-Body Particle Decays Featuring Multiple Decay Chains”, 2020, 10.1155/2020/6674595.

[redeboer]

Closed by #245 and #307