We discuss the algorithmic problem of
minimal coupling gauge fields of the Yang–Mills type
to Quantum Gravity in the approach known as Causal Dynamical Triangulations (CDT) as a
step towards studying, ultimately, systems of gravity coupled with bosonic and fermionic matter.
We first describe the algorithm for general dimensions and gauge groups
and then focus on the results obtained from simulations of 2d CDT coupled
to gauge fields with U(1) and SU(2) gauge groups,
where we studied both observables related to gravity and gauge fields,
and compared them with analogous simulations in the static flat case.