We made R-band polarization measurements of 234 stars towards the direction of the MBM 33-39 cloud complex. The distance of the MBM 33-39 complex was determined as 120+/-10pc using polarization results and near-infrared photometry from the 2MASS survey. The magnetic field geometry of the individual clouds inferred from our polarimetric results reveals that the field lines are in general consistent with the global magnetic field geometry of the region obtained from previous studies. This implies that the clouds in the complex are permeated by the interstellar magnetic field. Multi-wavelength polarization measurements of a few stars projected on to the complex suggest that the size of the dust grains in these clouds is similar to those found in the normal interstellar medium of the Milky Way. We studied a possible formation scenario of the MBM 33-39 complex by combining the polarization results from our study with those from the literature and by identifying the distribution of ionized, atomic and molecular (dust) components of material in the region.