In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008−57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise.
High mass X-ray binaries are among the brightest X-ray sources in the Milky Way, as well as in nearby Galaxies. Thanks to their highly variable emissions and complex phenomenology, they have attracted the interest of the high energy astrophysical community since the dawn of X-ray Astronomy. In more recent years, they have challenged our comprehension of physical processes in many more energy bands, ranging from the infrared to very high energies. In this review, we provide a broad but concise summary of the physical processes dominating the emission from high mass X-ray binaries across virtually the whole electromagnetic spectrum. These comprise the interaction of stellar winds with the high gravitational and magnetic fields of compact objects, the behaviour of matter under extreme magnetic and gravity conditions, and the perturbation of the massive star evolutionary processes by presence in a binary system. We highlight the role of the INTEGRAL mission in the discovery of many of the most interesting objects in the high mass X-ray binary class and its contribution in reviving the interest for these sources over the past two decades. We show how the INTEGRAL discoveries have not only contributed to significantly increase the number of high mass X-ray binaries known, thus advancing our understanding of the population as a whole, but also have opened new windows of investigation that stimulated the multi-wavelength approach nowadays common in most astrophysical research fields. We conclude the review by providing an overview of future facilities being planned from the X-ray to the very high energy domain that will hopefully help us in finding an answer to the many questions left open after more than 18 years of INTEGRAL scientific observations. ; The INTEGRALteams in the participating countries acknowledge the continuous support from their space agencies and funding organizations: the Italian Space Agency ASI (via different agreements including the latest one, 2019-35HH, and the ASIINAF agreement 2017-14-H.0), the French Centre national d'études spatiales (CNES), the Russian Foundation for Basic Research (KP, 19-02-00790), the Russian Science Foundation (ST, VD, AL; 19-12-00423), the Spanish State Research Agency (via different grants including ESP2017-85691-P, ESP2017-87676-C5-1-R and Unidad de Excelencia María de Maeztu – CAB MDM-2017-0737). IN is partially supported by the Spanish Government under grant PGC2018-093741-B-C21/C22 (MICIU/AEI/FEDER, UE). LD acknowledges grant 50 OG 1902.
