Investigation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si1.svg"> <mml:mrow> <mml:msubsup> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi>c</mml:mi> <mml:mi>s</mml:mi> <mml:mn>0</mml:mn> </mml:mrow> <mml:mo>*</mml:mo> </mml:msubsup>

We have used the parton and hadron cascade model PACIAE together with the Dynamically Constrained Phase-space Coalescence model (DCPC) to study the T c s 0 * ( 2870 ) 0 production in pp collision at s = 7 TeV, following the LHCb observation of T c s 0 * ( 2870 ) 0 in the B − → D − D 0 K S 0 decays in pp collisions at s = 7, 8, and 13 TeV [PRL 134(2025)101901]. The final hadronic states of the pp collisions at s = 7 TeV are first simulated by the PACIAE model. Four sets of T c s 0 * ( 2870 ) 0 candidates are then recombined by the DCPC model using the constituent meson pair of D 0 K S 0 , D + K − , D * + K * − , and D * 0 K ¯ * 0 based on the above simulated final hadronic states, respectively. We calculate their rapidity distributions, transverse momentum spectra, and angular distribution between the two component mesons, as well as angular distribution between D component meson and T c s 0 * ( 2870 ) 0 . Our results show that the yields of four T c s 0 * ( 2870 ) 0 candidates follow the magnitude order of D 0 K S 0 > D + K − > D * + K * − ∼ D * 0 K ¯ * 0 . Similar ordering behavior is also observed in the aforementioned distributions.