Usually, in surface hopping simulations, the excitation pulse is supposed to be instantaneous. Thus, dynamics start in the excited state without any explicit pulse description. Nevertheless, the instantaneous pulse hypothesis is not adequate in excitation by thermal (incoherent) light, like solar irradiation. Thermal light is better represented by a continuous field acting on the molecule, slowing transporting electronic density to the excited state. The mixed quantum-classical pulsed ensemble (MQC-PE) method makes it possible to simulate the action of such thermal sources with minimal changes in conventional DC-FSSH. It simulates the incoherent light as an ensemble of pulses, following the Chenu-Brumer approach. In MQC-PE, the initial state L from where the trajectory should start is sampled with the probability distribution given by the black body radiance at temperature T [see inequality (31)]. The dynamics is carried out using conventional DC-FSSH. After the simulations are finished, ensemble averages are computed by translating the initial time of each trajectory to span the whole illumination window and weighting each trajectory by the field-induced initial-state population [eq (10) of reference [25]].