Abstract-Recently, there has been significant research effort towards achieving high data rates in the millimeter wave bands by employing large antenna systems. These systems are considered to have only a fraction of the RF chains compared to the total number of antennas and employ analog phase shifters to steer the transmit and receive beams in addition to the conventional beamforming/combining invoked in the baseband domain. This scheme, which is popularly known as hybrid beamforming, has been extensively studied in the literature. To the best of our knowledge, all the existing schemes focus on obtaining the beamforming/combining matrices that maximize the system capacity computed using a Gaussian input alphabet. However, this choice of matrices may be suboptimal for practical systems, since they employ a finite input alphabet, such as QAM/PSK constellations. Hence, in this paper, we consider a hybrid beamforming/combining system operating with a finite input alphabet and optimize the analog as well as digital beamforming/combining matrices by maximizing the mutual information (MI). This is achieved by an iterative gradient ascent algorithm that exploits the relationship between the minimum mean-squared error and the MI. Furthermore, an iterative algorithm is proposed for designing a codebook for the analog and digital beamforming/combining matrices based on a vector quantization approach. Our simulation results demonstrate that the proposed gradient ascent algorithm achieves an ergodic rate improvement of up to 0.4 bits per channel use (bpcu) compared to the Gaussian input scenario. Furthermore, the gain in the ergodic rate achieved by employing the vector quantization based codebook is about 0.5 bpcu compared to the Gaussian input scenario.Index Terms-Hybrid processing, mm-wave communication, mutual information, beamforming and combining.