Initial states of high energy heavy ion collisions are studied using a dipole model through the DIPSY event generator that dynamically includes saturation together with the fluctuations and correlations of the BFKL cascade. The eccentricities of first up to fourth moment are calculated at RHIC and LHC. Predictions are made for correlations and fluctuations in rapidity of the eccentricities, and conventional theoretical approximations are tested. A large set of initial state Au-Au, Cu-Au and Pb-Pb collisions have been generated and are published online.
