After the calculation we analyze the wave packet to get the energy dependent reaction probability. To calculate this quantity we use the method taken from Ref. . The energy dependent reaction probability is defined as
Eq. (13) can also be derived rigorously rather than with the hand-waving argument presented here. However, the proof is elaborate, albeit straightforward. Therefore, we will only outline it here. First, take a plane wave along the SF z-axis, expand it into spherical Bessel functions and rotate it to the BF-frame. Then, project the plane wave onto a specific (J,M) state and normalize this state. Next, write the spherical Bessel functions in their asymptotic expression and take the inner product with the initial wave packet. As in Ref.  one gets two integrals, one of which is negligible. Neglegting this integral¸ one obtains Eq. (13), apart from a complex factor which disappears in calculating Pr(E).
By taking the surface not at rs, but at Rs we can use this method also to calculate Pnr(E), the non-reaction probability. This last quantity is interesting, because we can use it to assess the convergence of the calculation. If the results are converged, Pr(E) + Pnr(E) are equal to 1.