How to calculate angular order parameters in solution?

Dear all,
Hello everyone, I would like to ask if there is a script or method for calculating the angular order parameters(tetrahedral order parameter) of water molecules in solution,or four-body order parameter ?