Long Range Azimuthal Correlations of Multiple Gluons in Gluon Saturation Limit

We calculate the inclusive gluon correlation function for arbitrary number of gluons with full rapidity and transverse momentum dependence for the initial glasma state of the p-p, p-A and A-A collisions. The formula we derive via superdiagrams generates cumulants for any number of gluons. Higher order cumulants contain information on correlations between multiple gluons, and they are necessary for calculations of higher dimensional ridges as well as flow coefficients from multi-particle correlations.