We study, within the Monte-Carlo Glauber model, centrality dependence of the midrapidity charged multiplicity density dN_{ch}/deta and of the anisotropy coefficients varepsilon _{2,3} in Pb + Pb collisions at sqrt{s}=5.02 TeV and in Xe + Xe collisions at sqrt{s}=5.44 TeV. Calculations are performed for versions with and without nucleon meson cloud. The fraction of the binary collisions, alpha , has been fitted to the data on dN_{ch}/deta in Pb + Pb collisions. We obtain alpha approx 0.09(0.13) with (without) meson cloud. The effect of meson cloud on the dN_{ch}/deta is relatively small. For Xe + Xe collisions for 0–5% centrality bin we obtain dN_{ch}/deta approx 1149 and 1134 with and without meson cloud, respectively. We obtain varepsilon _2(mathrm {Xe})/varepsilon _2(mathrm {Pb})sim 1.45 for most central collisions, and varepsilon _2(mathrm {Xe})/varepsilon _2(mathrm {Pb}) close to unity at . We find a noticeable increase of the eccentricity in Xe + Xe collisions at small centralities due to the prolate shape of the Xe nucleus. The triangularity in Xe + Xe collisions is bigger than in Pb + Pb collisions at . We obtain varepsilon _3(text{ Xe })/varepsilon _3(text{ Pb })sim 1.3 at .