A bstract This work studies the self-interacting dark matter (SIDM) scenario in the general NMSSM and beyond, where the dark matter is a Majorana fermion and the force mediator is a scalar boson. An improved analytical expression for the dark matter (DM) self-interacting cross section which takes into account the Born level effects is proposed. Due to the large couplings and light mediator in SIDM scenario, the DM/mediator will go through multiple branchings if they are produced with high energy. Based on the Monte Carlo simulation of the showers in the DM sector, we obtain the multiplicities and the spectra of the DM/mediator from the Higgsino production and decay at the LHC for our benchmark points.