Diffuse interstellar bands (DIBs) are broad and weak absorption features that appear in spectra recorded towards stars that lie behind interstellar (IS) clouds. They are known as a long-standing challenge in astronomical spectroscopy because of their unknown carriers. However, they are also a potential IS tracers, in particular to probe distant clouds because they are not easily saturated. We devised automated DIB-fitting methods appropriate for cool star spectra and multiple IS components. The data were fitted with a combination of a synthetic stellar spectrum, a synthetic telluric transmission, and empirical DIB profiles. The initial number of DIB components and their radial velocity were guided by HI 21 cm emission spectrum, or, when available in the spectral range, IS neutral sodium (NaI) absorption lines. For NaI, radial velocities of IS NaI lines and DIBs were maintained linked during a global simultaneous fit. Using the tools, we have analyzed Gaia-ESO Survey (GES) spectra of stars that probe between 2-10 kpc long line-of-sight (LOS) in few different regions of the Milky Way (fields). For all fields, the DIB strength and target extinction are well correlated. We show that, toward distance-distributed target stars, DIBs can be used to locate dense clouds in the Galaxy (i.e., Galactic arm) and probe a LOS kinematical structure. In addition, we have devised the DIB-DIB global fit to study the relation between DIBs and can be a potential tool to trace interstellar environmental conditions.