Mass-loss rate determination for the massive binary V444 Cygni using 3-D Monte-Carlo simulations of line and polarization variability