Using the same inputs in here I re-implemented the band structure and projected density of state calculations using Quantum Espresso, then I used the same command:
import pyprocar
pyprocar.bandsplot(code='qe',mode='plain', elimit=[-18,-8])
to plot the bands using PyProcar, I get:
Which is clearly different from the original one here. It seems to me the bands have been shifted by 2E_F.
Is there any explanation?