including minor updates ... code used for results presented here: https://indico.cern.ch/event/303563/contribution/10/material/slides/0.pdf