  By using Monte-Carlo implementations of $k_T$-factorization formula with
running-coupling BK unintegrated gluon distributions for nucleus-nucleus
collisions, we compute higher order harmonic moments of the initial density
distribution for both RHIC(Au+Au@200GeV) and<EMAIL_ADDRESS>collisions. We
study their sensitivity to the size of the valence parton distribution in the
nucleon.
