One of the crucial observations from the electron-proton DIS data at HERA is the predominance of gluons in the hadronic wavefunction and the burgeoning of gluon densities at
high energies or small Bjorken-x. It is conjectured that repulsive gluon recombination and screening effects rival against bremsstrahlung to attenuate this resulting growth of high energy
cross-sections, in conformity with constraints imposed by unitarity. As an outcome of these competing effects, a semi-hard saturation scale, Q_{s}(x) dynamically emerges which allows
the use of effective field theory techniques to systematically compute the otherwise non-linear, many-body gluodynamics of this novel regime of QCD. This physical intuition is concretized
in a Wilsonian renormalization group (RG) framework known as the Color Glass Condensate (CGC) effective field theory (EFT).
In this talk, I will discuss the basics of this formalism that describes the initial state of nucleonic matter in heavy ion collisions. Starting with the classical version originally
theorized for an infinite nucleus I will discuss the generalization for realistic nuclei. This will be followed by the computational tools that reflect the power of the CGC EFT to extend
ab-initio techniques for phenomenological analyses of collision data. In the last part of this talk, I will show the application of these theoretical tools towards a first computation of
the differential cross-section for inclusive photon+dijet production at small x.