Once you have DensityOfStates object you can just ask it for the two types of spectra, with the relevant projection:
dos = DensityOfStates(configuration, ...)
projection_list = ProjectionList(elements=Sulphur)
ts = dos.tetrahedronSpectrum(projection_list)
gs = dos.gaussianSpectrum(projection_list)
You may be interested in the projection on individual atoms, then use ProjectionList(atoms=[12]) etc where 12 is the atom index (starting at 0) of the atom to project on; you can project on more than one at the time.
Tips on plotting are in the reference manual under the link above.