The detailed knowledge of excited state proton transfer mechanisms in complex environments is of paramount importance in chemistry. However, the definition of an effective reaction coordinate and the understanding of the driving force of the reaction can be difficult from both the experimental and the theoretical points of view. Here we analyzed by theoretical approaches the mechanism and the driving forces of the excited state proton transfer reaction occurring between the 7-hydroxy-4-(trifluoromethyl)coumarin photoacid and the 1-methylimidazole base molecules in toluene solution. In particular, we compared the intrinsic and the dynamical reaction pathways, obtained by integrating the reaction coordinate, and by performing ab initio simulations of molecular dynamics, respectively. Time-dependent density functional theory and polarizable solvation continuum models were adopted to define the excited state potential energy surface. Results were analyzed by means of the D(CT) electronic density based index. Our findings suggest that the reaction coordinate is mainly composed of several intra- and intermolecular modes of the reactants. An analysis of both the intrinsic coordinate and the dynamical results shows that the charge transfer induced by electronic excitation of the coumarin molecule is the main proton transfer driving force. With regards to the methodological validation, the combination of ab initio molecular dynamics with time-dependent density functional theory appears to be feasible and reliable to study excited state proton transfer reactions in the condensed phase.