Implement trial-masked robust detrending 

https://www.sciencedirect.com/science/article/pii/S0165027021000157

It's almost already doable with the current code. 

What is missing is a function that takes event times and the max and min of your wanted epoch (so what you want to mask) and output appropriate weights as a `(n_channels * n_times)` matrix.

