I would appreciate if someone could walk me through this derivation.
I have a transfer function in the frequency domain, which has two poles
$$\tilde{H}(\omega) = \Big(\frac{1}{1 + i \omega \tau_1}\Big)\Big(\frac{1}{1 + i \omega \tau_2}\Big)$$
For a single pole, the inverse Fourier transform is: $$\mathcal{F}^{-1}\Big\{\Big(\frac{\tau}{1 + i \omega \tau}\Big)\Big\} = \Theta(t) ~ e^{\displaystyle -\frac{t}{\tau}} $$
where $\Theta(t)$ is the Heaviside function since this is a causal filter. (See this answer.) I also plugged the above $\tilde{H}$ into Wolfram Alpha and got the expected result:
$$\mathcal{F}^{-1}\Big\{\tilde{H}(\omega)\Big\} = \frac{1}{\tau_2 - \tau_1} \Bigg( ~ e^{\displaystyle -\frac{t - t_d}{\tau_2}} - e^{\displaystyle -\frac{t - t_d}{\tau_1}} \Bigg)\Theta(t - t_d)$$
where $t_d$ is the time after which the causal filter is effective. I'm not sure if it's a reversal of integration limits, integration by parts, or simply some of the terms cancelling after integration and plugging in the limits that might be responsible for the subtraction.
And I'm also not 100% sure how to work in the $t_d$ term.