A quark model analysis of the Sivers function