Adding back the R-wave to ECG_Delineate?

[RichardBJ]

I find this routine remarkably great. (a couple of tiny little bugs for me); I use three outputs:
It uses the detected R peaks to seed the hunt for QRST...
i) Simple plot showing delineated events with peaks marked.
ii) Snazzy dictionary with Peak as the key and list of time stamps as value.
iii) sometimes the actually delineated record...
OK so the first thing I typically do is update the dict to include the original R-peaks, that is easy enough. Next, however, I really want the plot to also include the original R-wave peaks? I also typically want to drop the Q wave since these are rodents and they don't have much of a Q (although estimating Q is necessary for finding the other peaks!).

So I wondered about allowing EITHER a list input argument that lists the peaks you want to return and plot..... or even simpler, a Boolean stating whether to include the original R-wave in the outputs? Was happy to code this if required. I think that looks straightforward. Certainly the second option. What do you think?
Additionally, we do a lot of rodent work, so I am not sure what appetite there is for updating so that runs more contiguously rather than extra-package tweaking? [I could add a species argument to ECG delineate, default="human".... for example].