We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich ( SZ ) power spectrum and cluster counts analyses , highlighting subtle dependencies on the parameterisation and data combination . In \Lambda CDM , with an X-ray mass bias corresponding the expected hydrostatic mass bias , i.e. , b \simeq 0.2 , our constraints from Planck SZ data are consistent with the latest results from SPT , DES-Y1 and KiDS+VIKING-450 . In \nu \Lambda CDM , without prior information on b , our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint , as well as ( 1 - b ) = 0.64 \pm 0.04 ( 68 % CL ) . For forecasts , we find that a combination of a mock cosmic variance limited SZ power spectrum with primary CMB and BAO can improve the uncertainty on the total neutrino mass by 14 % with respect to CMB combined with BAO . This requires masking the heaviest clusters and probing the small-scale SZ power spectrum up to multipoles of \ell _ { \mathrm { max } } = 10 ^ { 4 } . For a future competitive measurement of the total neutrino mass using CMB and the SZ power spectrum , but excluding BAO and lensing power spectrum , we find that a 1 % precision on the mass calibration of clusters is needed . Although SZ power spectrum-based measurements of the neutrino masses are challenging , we find that the SZ power spectrum can be used to tightly constrain intra-cluster medium properties .