Calculates probability distribution functions for the spherical coordinates of a set of vectors between atoms in a gromacs trajectory relative to a reference vector