Spontaneous Raman spectroscopy is a formidable tool to probe molecular vibrations. Under electronic resonance conditions, the cross section can be selectively enhanced enabling structural sensitivity to specific chromophores and reaction centers. The addition of an ultrashort, broadband femtosecond pulse to the excitation field allows for coherent stimulation of diverse molecular vibrations. Within such a scheme, vibrational spectra are engraved onto a highly directional field, and can be heterodyne detected overwhelming fluorescence and other incoherent signals. At variance with spontaneous resonance Raman, however, interpreting the spectral information is not straightforward, due to the manifold of field interactions concurring to the third order nonlinear response. Taking as an example vibrational spectra of heme proteins excited in the Soret band, we introduce a general approach to extract the stimulated Raman excitation profiles from complex spectral lineshapes. Specifically, by a quantum treatment of the matter through density matrix description of the third order nonlinear polarization, we identify the contributions which generate the Raman bands, by taking into account for the cross section of each process.