code: def get_site_spd_dos(self, site):
          
          spd_dos = dict()
          for orb, pdos in self.pdos[site].items():
              orbital_type = _get_orb_type(orb)
              if orbital_type in spd_dos:
                  spd_dos[orbital_type] = add_densities(spd_dos[orbital_type], pdos)
              else:
                  spd_dos[orbital_type] = pdos
          return {orb: Dos(self.efermi, self.energies, densities)
                  for orb, densities in spd_dos.items()}
docs: Get orbital projected Dos of a particular site
  
          Args:
              site: Site in Structure associated with CompleteDos.
  
          Returns:
              dict of {orbital: Dos}, e.g. {"s": Dos object, ...}