How to do angle distribution between water dipole vectors and a normal vector?

Dear All

I wrote a script to generate water dipole vectors and normal vectors of a best fit plane. But I don't know how to continue to do angle distribution analysis with Forcite module, does anyone have an idea how to do that? Thank you.

Regards,

Jingjing Bu