We calculate differential distributions for diffractive dijets production in $ep\rightarrow e^{'}p\,jet\,jet$ reaction using off diagonal unintegrated gluon distributions (GTMDs). Different models are used. We concentrate on the contribution of exclusive $q\bar{q}$ dijets.
The results of our calculations are compared to H1 and ZEUS data. In general, except of one GTMD, our results are below the HERA data. This is in contrast with recent results where the normalization was adjusted to some selected distributions and no agreement with other observables was checked. We conclude that the calculated cross sections are only a small part of the measured ones which probably also contain processes with pomeron remnant, reggeon exchange, etc.
We also present azimuthal correlations between the sum and the difference of dijet transverse momenta. The cuts on transverse momenta of jets generate azimuthal correlations which can be misinterpreted as due to so-called elliptic GTMD.

