  We present evolutionary tracks of binary systems with high mass companionstars and stellar-through-intermediate mass BHs. Using Eggleton's stellarevolution code, we compute the luminosity produced by accretion from the donorduring its entire evolution. We compute also the evolution of the opticalspectrum of the binary system taking the disc contribution and irradiationeffects into account. The calculations presented here can be used to constrainthe properties of the donor stars in Ultraluminous X-ray Sources by comparingtheir position on the HR or color-magnitude diagrams with the evolutionarytracks of massive BH binaries. This approach may actually provide interestingclues also on the properties of the binary system itself, including the BHmass. We found that, on the basis of their position on the color-magnitudediagram, some of the candidate counterparts considered can be ruled out andmore stringent constraints can be applied to the donor masses.
Using Eggleton's stellar evolution code, we compute the luminosity produced by accretion from the donor during its entire evolution. We compute also the evolution of the optical spectrum of the binary system taking the disc contribution and irradiation effects into account.