In this work, the concentration of nine cannabinoids, six neutral cannabinoids (THC, CBD, CBC, CBG, CBN and CBDV) and three acidic cannabinoids (THCA CBGA and CBDA), was used to identify the Italian retailers of Cannabis sativa L. (hemp), reinforcing the idea that the practice of categorizing hemp samples only using THC and CBD is inadequate. A high-performance liquid chromatography/high-resolution mass spectrometry (HPLC-MS/MS) method was developed for screening and simultaneously analyzing the nine cannabinoids in 161 hemp samples sold by four retailers located in different Italian cities. The hemp samples dataset was analyzed by univariate and multivariate analysis with the aim to identify the hemp retailers without any other information on the hemp samples like Cannabis strains, seeds, soil and cultivation characteristics, geographical origin, product storage, etc. The univariate analysis highlighted that the hemp samples could not be differentiated by using any of the nine cannabinoids analyzed. To evaluate the real efficiency of the discrimination among the four hemp retailers a partial least squares discriminant analysis (PLS-DA) was applied. The PLS-DA results showed a very good discrimination between the four hemp retailers with an explained variance of 100% and low classification errors in both calibration (5%) and cross validation (6%). A total of 92% of the hemp samples were correctly classified by the cannabinoid variables in both fitting and cross validation. This work contributed to show that an analytical method coupled with multivariate analysis can be used as a powerful tool for forensic purposes.