def METHOD_NAME(self, lineshape = 'GAU', HWHM = 100.0):
    """
    Calculate density of states by computing the overlap of the two FCF spectra
    The two spectra are broadened using the chosen lineshape and Half-Width at Half-Maximum in cm-1
    """
    # Initialize spectrum
    self._spcinit()
    # Get stick spectra
    absstick = self._getstick(self.absrkf)
    emistick = self._getstick(self.emirkf)
    # Convolute with gaussians
    absspc = numpy.copy(self.spc)
    emispc = numpy.copy(self.spc)
    absspc = self.convolute(absspc, absstick, lineshape, HWHM)
    emispc = self.convolute(emispc, emistick, lineshape, HWHM)
    # Integrate
    # Calculate DOS
    self.spc[:, 1] = absspc[:, 1] * emispc[:, 1]
    METHOD_NAME = self.trapezoid(self.spc)
    return METHOD_NAME