We present results for the azimuthal anisotropy of charged hadron distributions in A+A, p+A, d+A, and 3 He+A collisions within the IP-Glasma+music model. Obtained anisotropies are due to the fluid dynamic response of the system to the fluctuating initial geometry of the interaction region. While the elliptic and triangular anisotropies in peripheral Pb+Pb collisions at √ s = 2.76 TeV are well described by the model, the same quantities in √ s = 5.02 TeV p+Pb collisions underestimate the experimental data. This disagreement can be due to neglected initial state correlations or the lack of a detailed description of the fluctuating spatial structure of the proton, or both. We further present predictions for azimuthal anisotropies in p+Au, d+Au, and 3 He+Au collisions at √ s = 200 GeV. For d+Au and 3 He+Au collisions we expect the detailed substructure of the nucleon to become less important.