To evaluate the efficiency of our approach, we compared the results obtained using the CT-TSH method with the TSH method, with and without ODC decoherence correction. The comparisons included simulations of both nπ * and ππ * excitations. In addition, we tested four implementations of the CT-TSH algorithm, each of them with a distinct method for evaluating the nuclear density used in the computation of the quantum momentum. The nuclear density is represented as a sum of Gaussian functions, each Gaussian being associated with one trajectory. The four variants of the CT-TSH method concern the set of nuclear coordinates considered as arguments of the nuclear density and two different ways to determine the widths of the Gaussian functions.