We develop an Efficient Wake procedure for computing the distribution of hadrons originating from jet wakes in heavy ion collisions — the hydrodynamic response of a droplet of quark-gluon plasma to the energy and momentum deposited in it by high-energy partons propagating through it. The procedure employs the linearity of linearized hydrodynamics and takes account of the effects of both longitudinal expansion and transverse radial flow on the hydrodynamic evolution of the wakes and on the resulting particle production at the freezeout hypersurface.
It makes repeated use of template solutions to linearized hydrodynamics in a Bjorken flow background with no transverse flow, templates that need only be computed once, and uses suitable rotations and boosts to map fluctuations from these templates to fluctuations at a point on the freezeout hypersurface in a way that incorporates the effects of the radial flow. We benchmark this procedure by comparing its results to results obtained from full $(3+1)$-dimensional nonlinear hydrodynamics calculations, find reasonable agreement, and find that our Efficient Wake procedure yields a {\it much} better description of the distribution of hadrons originating from jet wakes than does the older oversimplified procedure employed in the Hybrid Model. And, the Efficient Wake procedure {\it is} computationally efficient: it is at least tens of thousands of times faster than full nonlinear hydrodynamics calculations. Hence, we anticipate that when our new procedure is implemented in Monte Carlo analyses of jets in heavy ion collisions, for example in the Hybrid Model, it will greatly improve the description of the soft component of many jet and jet substructure observables as compared to experimental data.