Radial distribution function

Add radial distribution function action.
C++ version works for any boxdimension.
Cython version works for 2D and 3D.
Depends on pele::distance and on mcpele::histogram.
Integer template parameters are handled like for the inverse power-law potential of pele.