Imagine I have a scattering region (denoted as sample).
Scattering matrix and transfer matrix gives the same information about scattering. The scattering matrix tells us how incoming modes are scattered into outgoing modes, whereas the transfer matrix tells us how the modes on the left are related to the right (or vice versa).
$$ \begin{pmatrix} \Psi_R^{(+)} \\ \Psi_L^{(-)} \end{pmatrix} = S \begin{pmatrix} \Psi_L^{(+)} \\ \Psi_R^{(-)} \end{pmatrix} $$ Here, $S = \begin{pmatrix} S_{11} & S_{12}\\ S_{21} & S_{22} \end{pmatrix}$
The transfer matrix is written in the following way $$ \begin{pmatrix} \Psi_L^{(+)} \\ \Psi_L^{(-)} \end{pmatrix}= M \begin{pmatrix} \Psi_R^{(+)} \\ \Psi_R^{(-)} \end{pmatrix} $$ Where, $M = \begin{pmatrix} M_{11} & M_{12}\\ M_{21} & M_{22} \end{pmatrix}$. The elements of $M$ and $S$ are related (see http://assets.press.princeton.edu/chapters/s8695.pdf for example).
My question is how to compose transfer matrices and scattering matrices when I have multiple scattering regions. Imagine I have the following series of scattering regions each defined by the scattering matrix $S$ and transfer matrix $M$.
In the transfer matrix approach I just need to multiply a string of transfer matrices to compose them, i.e., the transfer matrix of the whole system is $M_{\text{eff}} = M_1 M_2 M_3 M_4$.
How to compose scattering matrices to get a matrix ($S_{\text{eff}}$) that has two incoming and two outgoing modes?