Hemodynamic response estimation from fNIRS signal through a modeling approach exploiting the reference channel