(1) include the piezo stimulus in s.log and functions like make_df_multi_index(logs_stims)
(2) add a stimulus column with the actually presented stimulus (sound or piezo) in s.log
(3) s.stack() could have another argument: stack_trial: bool = False, if True it concatenates all the trials (and presented stimuli) for motion correction
(1) include the piezo stimulus in s.log and functions like make_df_multi_index(logs_stims)
(2) add a stimulus column with the actually presented stimulus (sound or piezo) in s.log
(3) s.stack() could have another argument: stack_trial: bool = False, if True it concatenates all the trials (and presented stimuli) for motion correction