code: def get_site_dos(self, site):
          
          site_dos = functools.reduce(add_densities, self.pdos[site].values())
          return Dos(self.efermi, self.energies, site_dos)
docs: Get the total Dos for a site (all orbitals).
  
          Args:
              site: Site in Structure associated with CompleteDos.
  
          Returns:
              Dos containing summed orbital densities for site.