A common problem on this site is to compute the Fourier transform of the Heaviside function that is $0$ on the negative reals and $1$ on the positive reals. A standard technique is to consider the approximations $H_\alpha(x)$ defined to be $0$ on the negative reals and $e^{-\alpha x}$ on the positive reals. We then compute the Fourier transform and obtain (ignoring a constant factor) $$\widehat{H_\alpha}(w) = \frac{1}{\alpha + iw}.$$
One can show that taking the limit as $\alpha$ goes to zero gives the correct answer.
I am uncomfortable with one part of this argument. We are claiming that since $H_\alpha\rightarrow H$ as distributions, $\widehat{H_\alpha}\rightarrow \widehat H$. On one hand, this seems true because $\widehat H_\alpha(\phi)=H(\hat{\phi})$, so just staring at the definitions seems to prove it. On the other hand, I am having trouble making sense of the following computation.
Let's consider approximating $H$ a different way, using the characteristic functions of the interval $[0,M]$. We should obtain $$\widehat H(w) = \lim_{M\rightarrow \infty} \int_0^M e^{-iwt} \ dt = \lim_{M\rightarrow \infty} \frac{1}{\sqrt{2\pi}} \frac{1-e^{iwM}}{iw}.$$
But this limit does not appear to be $\hat H$. In fact, it doesn't seem to exist.
What is going on here? If my reasoning is correct, how can one show the latter limit gives the right answer?