bw_ranges functions are ready to get a GRanges object to reduce the analysis to only bins or loci that overlap a second set of loci. This has been used only for debugging, but I think it is an interesting feature to have, which would also involve a nice refactor.
Raise the use of selection to the plotting functions and include it also on bed_ functions.