Monte Carlo Radiative transfer codes (MCRT) are a powerful tool to interpret the results of astronomical measurements. However, the majority of the MCRT codes are not able to calculate polarization maps. Our code is already capable of calculating the full Stokes vector (I,Q,U,V) for scattering of light by electrons and spherical grains. We are currently implementing the mechanisms of dichroic extinction and polarized dust emission and present first results.