update healtest to loop over pairs of points
update to use points less than specified max. use ring ordering scheme by default
add fake lines of sight along quasar pointing vectors
move correlation function estimation into standalone function