I posted the exact same question on odForce, but unfortunately got no answer, I waited a few days, so I thought maybe someone can help me here.

As you know, we can use dopfield expression to get the data in dops, one of them is impulse when an impact happens, there are usually a lot of subimpacts (specially in case of fractured objects) that each one has a different impulse, I want to know how I can access the max impulse for each piece?

You can find the max impulse with CHOP's; Import the impacts to SOP's then fetch the impact points in CHOP's, find the max with a math CHOP (OP tab, Combine Channels to Maximum).

Thanks Alejandro,

But CHOPs is very slow, specially when you have hundreds of pieces and each piece has a lot of impacts, so I was wondering if there's a fast way to do that, either with a python code or maybe a simple expression.