Question for calculating angular autocorrelation function of OH ion by pymatgen

Dear all,

Based on AIMD by vasp, I want to obtain the angular autocorrelation function of OH in Li2OHCl structure, as described in following figure, but i don’t know whether this function can be realized by pyamtgen? Or any other software can realize this function?

Please give me some suggestions.

Thank you

image