Dear all,
I am using Kwant to study a system with both spins and electron and hole, so the hopping matrix is 4X4. Now I want to know how to separate the spins, electron and hole. I study this through the "2.6. Superconductors: orbital degrees of freedom, conservation laws and symmetries"
smatrix = kwant.smatrix(syst, energy)
data.append(smatrix.submatrix((0, 0), (0, 0)).shape[0] -smatrix.transmission((0, 0), (0, 0)) +smatrix.transmission((0, 1), (0, 0)))

For my system, if I have 3 leads, and the hopping and onsite energy is set in this order:
(e↑,0,0,0)- spin up electron,first row of the 4X4matrix
(0,e↓,0,0)- spin down electron,second row of the 4X4matrix
(0,0, h↑ ,0)- spin up hole,third row of the 4X4matrix
(0,0, 0, h ↓ )- spin down hole,fourth row of the 4X4matrix

I want to obtain the transmissions from 0 →2.
smatrix = kwant.smatrix(syst, energy)
The transmission from h↑  to e↑  is: smatrix.transmission((2, 1), (0, 2))
The transmission from h ↓   to e↑  is: smatrix.transmission((2, 1), (0, 3))
Is my understanding correct?