Comparison between the input reionization scenario and recovered model from the MCMC analysis using EoR power spectra as input to the phenomenological model considered in this work. Left to right columns represent C2RAY, 21CMFAST, and GRIZZLY simulation respectively. From top to bottom we show, respectively, the ratio of 21 cm signal to density power spectrum (${\mathrm{\Delta }}_{\delta T_{\text{b}}}^2/{\mathrm{\Delta }}_{\delta \delta }^2$) as a function of k-scales at different redshifts, the redshift evolution of ${\mathrm{\Delta }}_{\delta T_{\text{b}}}^2/{\mathrm{\Delta }}_{\delta \delta }^2$ for different scales and the corresponding reionization histories. The solid curves represent the inputs from simulations while the dotted curves stand for the MCMC best-fit prediction on the power spectrum model used in this work.