  The high mass X-ray binaries (HMXBs) provide an exciting framework to
investigate the evolution of massive stars and the processes behind binary
evolution. HMXBs have shown to be good tracers of recent star formation in
galaxies and might be important feedback sources at early stages of the
Universe. Furthermore, HMXBs are likely the progenitors of gravitational wave
sources (BH--BH or BH--NS binaries that may merge producing gravitational
waves). In this work, we investigate the nature and properties of HMXB
population in star-forming galaxies. We combine the results from the population
synthesis model MOBSE (Giacobbo et al. 2018) together with galaxy catalogs from
EAGLE simulation (Schaye et al. 2015). Therefore, this method describes the
HMXBs within their host galaxies in a self-consistent way. We compute the X-ray
luminosity function (XLF) of HMXBs in star-forming galaxies, showing that this
methodology matches the main features of the observed XLF.
