add histograms
add fig_bmass.py