Xgboost take k best predictions
After fitting the model you can use predict_proba( ) from the docs here

This returns a numpy array with the probability of each data example being of a given class. The three highest probabilities will be your best 3 predictions.