New function orient_corr

calculates g_m(r). Works, who knows if correctly.

Also some optimizations:
 - average is done with axis arg now in avg_hists
 - use numpy not zip, keep distances, positions, etc as arrays

	modified:   correlation.py