import pyfar as pf
import numpy as np
n_samples = 64
delay_samples = n_samples // 2 + 1/2
ir = pf.signals.impulse(n_samples)
ir = pf.dsp.linear_phase(ir, delay_samples, unit='samples')
start_samples = pf.dsp.find_impulse_response_delay(ir)
ax = pf.plot.time(ir, unit='ms', label='impulse response')
ax.axvline(
    start_samples/ir.sampling_rate*1e3,
    color='k', linestyle='-.', label='start sample')
ax.legend()