over several picoseconds and a wide range of probe wavelengths around the LSPR peak are followed by data fitting using the two-temperature model. The constraints, imposed by the second-order data, lead us to modify the model and force us to include the contribution of nonthermalized electrons in the early stages of the dynamics.