How to Plot DOS using pymatgen and python?

Hii users
I want to plot DOS similar to figure added below using pymatgen and python. Could anyone please help me how to do this? Also what is mean by get_spd_dos and get_element_dos? I am confused is in spd s orbital is the sum of all elements valence s orbital and in element dos is it total dos of individual element or valence orbital dos of individual element?
Screenshot 2024-02-19 172850