In QFT we define the scattering matrix from the scattering amplitude as
$$ S_{fi} = \lim\limits_{t\rightarrow \infty}\lim\limits_{t_0\rightarrow -\infty }\left\langle f\left|U(t,t_0)\right|i \right\rangle = \delta_{fi} +(-i)(2\pi)^4\delta^4(p^\mu_f-p^\mu_i)\underbrace{\mathcal{M}_{fi}}_{\text{Scatt. Amp.}} + \dots \; , $$ such that once we wish to obtain the probability of a scattering event, we need to evaluate $|M_{fi}|^2$. In the case where we want to obtain the full probability, we should sum over all of these individual amplitudes.
I understand $S_{fi}$ as the unitary transformation matrix that relates the initial and final states, while $M_{fi}$ is a probability density of a certain event occuring. What is the physical interpretation of $S_{fi}$. Is it also a probability density?
In Sakurai, Modern Quantum Mechanics, a so-called transition matrix $T_{fi}$ is defined from $\langle f| V(r) | i\rangle$, how is this related, why is it useful, and how does one interpret it?