The high redshift Epoch of Reionization can be probed with several emission lines in order to obtain a complete picture of the reionization process. In this study we produce simulations in large sky fields of CII, Lyman alpha and 21 cm lines emitted at redshifts of z≈8 to z≈6. We were able to predict the intensity and spatial fluctuations of these lines using observational data, galaxy properties from catalogs made with the Millennium and Millennium II simulations and an altered version of the Simfast21 code especially adapted to simulate these lines. In particular we found that there are good prospects for the detection of Lyman alpha emission and CII emission with expected average intensities of 1×10¹ Jy sr⁻¹ and 3.29×10² Jy sr⁻¹, respectively. The combined observation of these lines has the potential to establish a complete picture of the characteristics, spatial distribution and relative distribution of the gas in its several phases. Taking advantage of the intensity spatial fluctuations that can be detected with the intensity mapping technique we also consider the potential of cross correlation between the 21 cm line with the CII line or the Lyman alpha line as a way to remove foregrounds.