How are you unable to make the pull distributions? Do you get an error, if so, what error? Over which entries do you want to plot the pulls - I only see one histogram? I suppose you want to plot the pull over each data point? Then this should help: Plot of Fit Residuals