Context. Barium and S stars without technetium are red giants and are suspected of being members of binary systems due to their overabundances in heavy elements. These elements are produced by thes-process of nucleosynthesis, despite the stars not being evolved enough to be able to activate thes-process in their interiors. A companion formerly on the asymptotic giant branch (now a white dwarf) is supposed to be responsible for the barium- and S-star enrichment ins-process elements through mass transfer.Aims. This paper provides both long-period and revised orbits for barium and S stars, adding to previously published orbits. The sample of barium stars with strong anomalies (i.e., those classified as Ba3, Ba4, or Ba5 in the Warner scale) comprises all known stars of that kind, and in that sense forms a complete sample that allows us to investigate several orbital properties of these post-mass-transfer binaries in an unbiased way.Methods. Orbital elements are derived from radial velocities collected from a long-term radial-velocity monitoring campaign performed with the HERMES spectrograph mounted on theMercator1.2 m telescope. These new measurements were combined with older, CORAVEL measurements. With the aim of investigating possible correlations between orbital properties and abundances, we also collected a set of abundances for barium stars with orbital elements that is as homogeneous as possible. When unavailable in the literature, abundances were derived from high-resolution HERMES spectra.Results. We find orbital motion for all barium and extrinsic S stars monitored (except for the mild barium star HD 95345). We obtain the longest period known so far for a spectroscopic binary involving an S star, namely 57 Peg with a period of the order of 100−500 yr. We present the mass distribution for the barium stars, which ranges from 1 to 3M⊙, with a tail extending up to 5M⊙in the case of mild barium stars. This high-mass tail is mostly comprised of high-metallicity objects ([Fe/H] ≥ −0.1). The distribution of the companion masses was extracted from the barium-star mass distribution combined with the finding thatQ≡f(MBa,MWD)/sin3i=MWD3/(MBa+MWD)2is peaked at 0.057 ± 0.009 and 0.036 ± 0.027M⊙for strong and mild barium stars, respectively (f(MBa, MWD) is the mass function obtained from the orbital elements of spectroscopic binaries with one observable spectrum). Mass functions are compatible with WD companions whose masses range from 0.5 to 1M⊙. Strong barium stars have a tendency to be found in systems with shorter periods than mild barium stars, although this correlation is rather lose, with metallicity and WD mass also playing a role. Using the initial–final mass relationship established for field WDs, we derived the distribution of the mass ratioq′=MAGB, ini/MBa(whereMAGB, iniis the WD progenitor initial mass, i.e., the mass of the former primary component of the system) which is a proxy for the initial mass ratio (the less mass the barium star has accreted, the better the proxy). It appears that the distribution ofq′ is highly nonuniform, and significantly different for mild and strong barium stars, the latter being characterized by values mostly in excess of 1.4, whereas mild barium stars occupy the range 1−1.4.Conclusions. The orbital properties presented in this paper pave the way for a comparison with binary-evolution and nucleosynthesis models, which should account for the various significant correlations found between abundances and dynamical parameters (e.g. betweenMBaon one hand andMWD, [Fe/H], and [s/Fe] on the other hand, betweenq′ and [s/Fe], betweenPande, and betweenPand [s/Fe] altogether).