An effective method to compute vibrationally resolved optical spectra of large molecules in solution including Duschinsky, temperature and Herzberg-Teller effects