In the coming decade the Gaia satellite will measure the positions and velocities of an unprecedented number of stars in our Galaxy, with unprecedented precision. Among many firsts, this revolutionary new data set will include full six-dimensional phase space information for millions of stars in the Galactic halo, including stars in many tidal streams. These streams, the products of hierarchical accretion, can be used to infer the Galactic gravitational potential thanks to the common origin of the stars in each one. We present a method for doing so by maximizing the information content (i.e. clumpiness) of the action space of the stream stars. This statistical approach eliminates the need to assign stars to particular streams. Using a toy m...