Towards a first estimate of the gluon Sivers function from AN data in pp collisions at RHIC