Probing the Weizsacker-Williams gluon Wigner distribution in pp collisions

We show that the diffractive forward production of two quarkonia, especially the $chi_{cJ}$ and $eta_c$ states, in proton-proton or proton-ion collisions can access the Weizsacker-Williams gluon Wigner distribution of the proton. We use the hybrid factorization approach where the collinear, double gluon PDF is applied for one of the protons and the $k_T$-dependent (Wigner or GTMD) distribution for the other. The production of quarkonia is treated in the NRQCD framework. A particularly concise formula is obtained for double $chi_{J=1}$ production.