Simulation of UV-vis spectra over the whole absorption window requires not only an accurate approach for the estimation of electronic transitions' energies, but also to take properly into account solvent and band-broadening effects. Although the first, solvent, can be easily introduced in quantum-chemical calculations at both implicit and explicit levels, the second, band broadening, is more troublesome to evaluate. To this end, in this contribution, we propose a protocol aimed to correctly simulate the whole absorption UV-vis spectra of organic molecules based on time-dependent density functional theory, coupled with implicit and explicit solvent models, and on the use of a fitting procedure with respect to the experimental data in order to define the best band broadening. This protocol is applied to the simulation of absorption spectra of 10H-phenothiazine and three related molecules, N,N-diphenylamine, iminodibenzyl and 10H-phenothiazine-5-oxide, which appear as impurities during its industrial synthesis. The obtained results show not only that the main peak positions are reproduced with an error not exceeding 10nm (in the 200-400nm range) but also that the overall shape of the UV-vis spectra can be correctly simulated.