Polarimetry has proven to be a useful observational tool in the study of Be stars. In recent years, Monte Carlo radiative transfer methods have often been used in modeling polarimetric measurements, since they allow accurate predictions of flux and polarization via numerical simulation of the interaction of light with matter, even in systems with complex configurations. We have extended these two complementary techniques to the investigation of binary Be stars by developing a Monte Carlo code that models the flux and polarization phase curves of a binary-disk system. We here present the latest results matching geometrical model configurations with our own spectropolarimetric observations of the Algol binary β Lyrae.