estimateAB function

Helper function to find region around the main peak of a distribution