  Cross sections for the concurrent production of pairs of quarkonia (J/psi,
Upsilon) and/or gauge bosons (W, Z) from double-parton scatterings (DPS) in
high-energy proton-nucleus and nucleus-nucleus collisions at the LHC are
calculated. The estimates are based on next-to-leading-order perturbative QCD
predictions, including nuclear modifications of the parton densities, for the
corresponding single-scattering cross sections. Expected event rates for
J/psi+J/psi, J/psi+Upsilon, J/psi+W, J/psi+Z, Upsilon+Upsilon, Upsilon+W,
Upsilon+Z, and same-sign W+W production in their (di)leptonic decay modes,
after typical acceptance and efficiency losses, are given for pPb and PbPb
collisions.
