Finally, the critical addition in this work compared with the previous JAM19 analysis [1]} is the inclusion of unidentified charged hadron data, along with charged pions and kaons, in the SIDIS off deuterium targets from the COMPASS Collaboration [2]}, [3]} at CERN.
Since the SIDIS data
\(\mathrm {d}\sigma _{\mbox{\tiny SIDIS}}^{h^\pm }/\mathrm {d}Q^2\mathrm {d}x_{\mbox{\tiny \rm Bj}}\mathrm {d}z_h\) 
are differential in \(x_{\mbox{\tiny \rm Bj}}\)  and \(z_h\) , they combine information on both PDFs and FFs, which appear in the description of SIA, Drell-Yan, and DIS data.
Furthermore, as illustrated in Fig. REF , the SIDIS data have significant overlap in \(x_{\mbox{\tiny \rm Bj}}\)  and \(z_h\)  with the \(x_{\mbox{\tiny \rm Bj}}\)  and \(x_{\mbox{\tiny \rm F}}\)  range of inclusive DIS and Drell-Yan data, respectively, and the \(z_h\)  range of SIA data, so that the combined analysis constitutes a genuine test of their universality.
For the COMPASS SIDIS data we use the same kinematic cuts on \(W^2\)  and \(Q^2\)  as for inclusive DIS, and restrict the fragmentation variable to \(0.2 < z_h < 0.8\)  in order to exclude data from the target fragmentation region and avoid large-\(z\)  threshold corrections.
