A stochastic model has been developed to study the hierarchical fragmentation process of young massive clusters in external galaxies considering close binary components along with individual ones. Stellar masses for individual ones have been generated from truncated Pareto distribution and stellar masses for close binary components have been generated from a truncated Bi-variate Gumbel Exponential distribution. The above distribution is identified by fitting the observed bi-variate distribution of masses of eclipsing binary stars computed from the light curves catalogued in the package Binary Maker 3.0. The resulting mass spectra computed at different projected distances, show signature of mass segregation. Degree of mass segregation becomes reduced due to the inclusion of binary fraction. This might be due to the reduction of massive stars and inclusion of less massive stars rather than inclusion of single massive stars and the effect of line of sight length projected to an observer. © 2016, Springer Science+Business Media Dordrecht.