Adding an argument `include_tmax` for `.get_data()`

We ran into a failure today in https://github.com/vferat/pycrostates/pull/14 because raw.get_data() was returning one less sample than expected, namely, the tmax sample.

from pathlib import Path

import mne
print (mne.__version__)

directory = Path(mne.datasets.sample.data_path()) / 'MEG' / 'sample'
fname = directory / 'sample_audvis_filt-0-40_raw.fif' 
raw = mne.io.read_raw_fif(fname, preload=True, verbose=False)

print (raw.get_data().shape)
print (raw.get_data(tmin=0, tmax=raw.times[-1]).shape)

OUT:

0.24.0
(376, 41700)
(376, 41699)

Do you consider this to be a bug in the handling of tmax, or should an additional argument be added, include_tmax, like for raw.crop()?

The relevant entry point is Raw.get_data(), with raw.crop() as the comparison named in the issue. Reproduce the two calls shown against the sample FIF data, then resolve whether the intended behavior is an include_tmax option or a correction to the existing handling; done means the selected behavior produces the expected sample shape.