Computation of the lineshape of electronic spectra in solution within polarizable continuum models: Accounting for both vibrational structure and inhomogeneous broadening