A new theoretical approach is presented for the general treatment of nonadiabatic hybrid dynamics (mixing classical and quantum approach) and applied to the postionization of rare-gas trimers. There was an important disagreement between trajectory surface hopping (TSH) or mean field (MF) approaches and the experimental results; noteworthy, with the new method qualitative and almost quantitative agreement is found for the fragmentation ratios of ionic monomers and dimers. For the first time in the theory as in the experiment, the dimers prevail for argon while monomers strongly dominate for the heavier rare gases, krypton and xenon. A new compromise between MF and TSH approaches is proposed and the new method is found quite robust with results not too sensitive to various possible implementations.