On the accuracy of explicitly correlated methods to generate potential energy surfaces for scattering calculations and clustering: application to the HCl–He complex