up to date, need to remake histograms in for loop, change to angular diam dist, and add estimator