  We investigate the inflation and reheating phenomenology in
scalar-Einstein-Gauss-Bonnet theory of gravity where a scalar field
non-minimally couples with the Gauss-Bonnet (GB) curvature term. Regarding the
inflationary phenomenology, we find -- (1) the inflation starts with a quasi
de-Sitter phase and has an exit at a finite e-fold, (2) the scalar and tensor
perturbations prove to be ghost free and do not suffer from gradient
instability, (3) the curvature perturbation amplitude as well as its tilt and
the tensor-to-scalar ratio turn out to be simultaneously compatible with the
recent Planck data for suitable values of the parameters. After the inflation
ends, the scalar field starts to decay to radiation with a constant decay
width. For our considered scalar potential and the GB coupling function, the
model results to an analytic power law solution of the Hubble parameter and a
logarithmic solution of the scalar field during the reheating era, where the
exponent of the Hubble parameter determines the effective EoS parameter
($w_\mathrm{eff}$) during the same. The stability of such reheating dynamics is
examined by dynamical analysis which ensures that $w_\mathrm{eff}$ can go
beyond unity and reach up-to the maximum value of $\mathrm{max}(w_\mathrm{eff})
= 1.56$. The scenario with $w_\mathrm{eff} > 1$ proves to be purely due to the
presence of the GB coupling function, which in turn may have important
consequences on enhancing the primordial gravitational waves' amplitude
observed today. The inflationary e-fold number gets further constrained by the
input of the reheating stage. We finally construct the complete forms of scalar
potential ($V(\phi)$) and the GB coupling ($\xi(\phi)$) function that smoothly
transits from inflation to reheating, and numerically solve the Hubble
parameter and the scalar field for such complete forms of $V(\phi)$ and
$\xi(\phi)$.
