started cpp code to try grouping more rapidly. empiricalBeta better estimates nbins and draws guassian pdf on plot