We have presented wide field deep NIR survey of eight regions containing BRCs to study the initial configuration of the YSOs and their nature of interaction with parent molecular clouds. The completeness limit for our observations in K band is ~18 magnitude. For our sample of BRC having distance around 2-3 kpc, we can detect upto sub-solar mass protostars. In addition to mid-infrared images and near-infrared extinction maps, we have identified and classified numerous YSOs, analyzed their spacings, and performed basic spatial distribution measurements and analyses. We identified and classified 100-300 YSOs in different regions with no more than 10-50 expected to be residual contaminators such as broad-line AGNs. Of those, 60 % are classified as IR excess stars, 20% are classified as CTTS, and 15%, 5% are classified as Class II and Class 0 protostars, respectively. We have also presented a purely algorithmic method to isolate local density enhancements in point source distributions from a more diffuse, poorly sampled, and potentially varying density background that uses no smoothing. This method was applied to our survey, extracting eight cluster cores of 10 or more YSO members. Of the identified YSOs, ~60% are members of one of these cores. We have demonstrated that protostars are found in regions of marginally higher surface densities than the more evolved pre-main-sequence stars with disks. The mean of the median spacings of YSOs (regardless of class) in the eight cores isolated here is ~0.06 pc. There suggestion by other investigations (e.g., Teixeira et al. 2006) that Jeans fragmentation is a worthwhile starting point for understanding primordial structure in star-forming regions. In that light, we infer the natal cloud properties from the mean YSO spacing and an assumed temperature of 20K, arriving at properties that are similar to those reported as probable cluster-forming clumps in infrared dark clouds by Rathborne et al. (2006).