[time_mat{b},D,~,~,~,n_trial(b)] = psth_onset(ca,rocmat(:,bv_list(b)),time_wd,null_bv,fr,'clean_bv_time',[],'windowend','all','nbout',nbo,'subtract_wd',[],'min_len',0 ...