How to separate the data plot for spin up and spin down band structure into 2 different graphs, in Pymatgen?

I am a beginner user of the Pymatgen package. In order to process the data from VASP DFT calculation software, I use Pymatgen to visualize the output band structure. When I do spin polarised band calculation, I cannot separate the band structure of spin up and spin down into 2 different subplots or graphs. Can anyone tell me how can I deal with this? Because it looks messy when the package automatically plot it into the same graph.

Here I provide an example with my own python scripts to realize your purpose rather than using Pymatgen (You can save data firstly with Pymatgen and plot with python). I assume that you can perform correctly the band and dos calculations with the VASP code. The example I pick up is the monolayer FM NiBr$$_2$$ and the final result is the following: