  This thesis studies gluon saturation in hadronic matter at high energy by
calculating next-to-leading order (NLO) corrections to inclusive and
diffractive deep inelastic scattering cross sections in the Color Glass
Condensate (CGC) effective field theory. We demonstrate that the large soft
gluon logarithm is correctly factorized into the Balitsky-Kovchegov (BK)
renormalization group equation by accurately connecting the NLO scattering
kinematics to the rapidity scale of the dipole amplitude in the scattering.
This brings the perturbative expansion under control and enables us to do
precision comparisons between theory and data. We fit the initial condition of
the BK evolution equation to HERA inclusive deep inelastic scattering data by
combining the NLO accuracy inclusive cross sections with beyond leading order
BK evolution prescriptions. This results in the state-of-the-art accuracy
comparison between CGC theory and HERA data, and determination of the dipole
amplitude initial shape which is a necessary input for all NLO CGC
phenomenology. In this introductory part, the effect of the NLO BK equation on
the fits is assessed, and an alternative form for the NLO loop correction to
the inclusive cross sections is derived which enables the consistent setting of
the dipole amplitude rapidity scale in the NLO corrections. Diffraction is
studied in this thesis in the CGC formalism, and we calculate the tree-level $q
\bar q g$ NLO contribution to the diffractive deep inelastic scattering
structure functions where the $q \bar q g$ Fock state scatters off the target
and becomes the diffractively produced system. This contribution has previously
been known in the literature only in leading $\log(Q^2)$ accuracy valid at
large $Q^2$, and only for the structure function $F^D_T$. The $q \bar q g$
contribution to both structure functions $F^D_T$ and $F^D_L$ are presented in
full NLO accuracy.
