Monte Carlo Glauber modeling of initial state anisotropies in ultrarelativistic heavy-ion collisions