Here we present a detailed analysis of the properties and evolution of different dwarf galaxies, candidate to host the coalescence of black hole binary systems (BHB) generating GW150914-like events. By adopting a novel theoretical framework coupling the binary population synthesis code SeBa with the Galaxy formation model GAMESH, we can investigate the detailed evolution of these objects in a well resolved cosmological volume of 4 cMpc, having a Milky Way-like (MW) galaxy forming at its center. We identify three classes of interesting candidate galaxies: MW progenitors, dwarf satellites and dwarf galaxies evolving in isolation. We find that: (i) despite differences in individual histories and specific environments the candidates reduce to only nine representative galaxies; (ii) among them, ∼ 44% merges into the MW halo progenitors by the redshift of the expected signal, while the remaining dwarfs are found as isolated or as satellites of the MW and their evolution is strongly shaped by both peculiar dynamical history and environmental feedback; (iii) a stringent condition for the environments where GW150914-like binaries can form comes from a combination of the accretion history of their DM halos and the radiative feedback in the high redshift universe; (iv) by comparing with the observed catalogues from DGS and ALLSMOG surveys we find two observed dwarfs respecting the properties predicted by our model. We finally note how the present analysis opens the possibility to build future strategies for host galaxy identification. the LIGO detectors: GW170104 (SNR = 13, Abbott et al. 2017a) and GW170814 (SNR = 18, Abbott et al. 2017b) 1 .The masses inferred for the components of the BHB generating GW150914 are M1 ∼ 36 M⊙ and M2 ∼ 29 M⊙; BHBs associated with GW170104 (and GW170814) are massive as well, with inferred component masses of M1 ∼ 31 M⊙ (∼ 30 M⊙), M2 ∼ 19 M⊙ (∼ 25 M⊙), respectively 2 .These massive BH binary systems are believed 1 Note that during the second run the lightest black hole binary so far observed was also detected (GW170608, Abbott et al. 2017c) with components in the range inferred for low-mass X-ray binaries: M 1 ∼ 12 M ⊙ and M 2 ∼ 7 M ⊙ , respectively. 2 During the referee process of our paper the LIGO-Virgo Collaboration announced the detection of four coalescing massive black