A search for new particles in an extended Higgs sector, characterized by a massive resonance decaying to a pair of scalar bosons , which themselves decay to pairs of b quarks, is presented. The analysis is restricted to the mass ranges from 25 to 100 GeV and from 1 to 3 TeV. For these mass ranges, the decay products of each are expected to merge into a single large radius jet. Jet substructure and flavor identification techniques are used to identify these jets from background jets. The search is based on LHC proton-proton collision data at \sqrt{s}=13~\mathrm{TeV}, collected with the CMS detector in -, corresponding to an integrated luminosity of 138~\mathrm{fb}^{-1}. Model-specific limits at confidence level on the production cross section times branching fraction as a function of mass for are set, where both the and branching fractions are assumed to be . These are the first such limits on this process, and range between 1 and 30 fb for the considered mass ranges.