A first working implementation of the pair_correlation_function method. Serious refactoring and design improvements needed