To measure spin-dependent parton distribution functions in the production of W bosons at the Relativistic Heavy Ion Collider, an accurate model for distributions of charged leptons from the W boson decay is needed. We present single-spin lepton-level cross sections of order $alpha_S$ for this process, as well as resummed cross sections, which include multiple parton radiation effects. We also present a program RhicBos for the numerical analysis of single-spin and double-spin cross sections in the Drell-Yan process, W and Z boson production.