The most recent experimental data for all measured production and decay channels of the bottomoniumlike states Z b ð10610Þ and Z b ð10650Þ are analyzed simultaneously using solutions of the Lippmann-Schwinger equations which respect constraints from unitarity and analyticity. The interaction potential in the open-bottom channels B ðÃÞBà þ c:c: contains short-range interactions as well as a one-pion exchange. It is found that the long-range interaction does not affect the line shapes as long as only S waves are considered. Meanwhile, the line shapes can be visibly modified once D waves, mediated by the strong tensor forces from the pion exchange potentials, are included. However, in the fit they get balanced largely by a momentum dependent contact term that appears to be needed also to render the results for the line shapes independent of the cutoff. The resulting line shapes are found to be insensitive to various higherorder interactions included to verify stability of the results. Both Z b states are found to be described by the poles located on the unphysical Riemann sheets in the vicinity of the corresponding thresholds. In particular, the Z b ð10610Þ state is associated with a virtual state residing just below the BB à =BB à threshold while the Z b ð10650Þ state most likely is a shallow state located just above the B ÃBà threshold.