added a/b/position_angle estimator, given g1g2 estimates, for input to GLMoments() or RunFDNT()
(#22)