timepoints = np.append(timestarts[:, np.newaxis], timeends[:, np.newaxis], axis=1)