The file "TightBindingModel_rs.dat" contains an array with the real-space tight-binding parameters and staggered orbital order in the following format:
x,y,z,orb1,orb2,hop 
Here x,y,z is the hopping vector, orb1 and orb2 are integer numbers describing combined sublattice/orbital index ({1,2,3,4} --> {A,d_xz; A,d_yz; B,d_xz; B,d_yz}) and hop is the value set for that particular hoppings integral.
To obtain the band structure, a Fourier transformation must be performed.

Remaining .dat files contain the data loaded in the Python script "LoadData_and_Plot.py" which reproduces theory results from the main text.