  We study shock statistics in the scalar conservation law $\partial_t u +
\partial_x f(u)=0$, $x \in \R$, $t>0$, with a convex flux $f$ and spatially
random initial data. We show that the Markov property (in $x$) is preserved for
a large class of random initial data (Markov processes with downward jumps and
derivatives of L\'evy processes with downward jumps). The kinetics of shock
clustering is then described completely by an evolution equation for the
generator of the Markov process $u(x,t)$, $x\in \R$. We present four distinct
derivations for this evolution equation, and show that it takes the form of a
Lax pair. The Lax equation admits a spectral parameter as in Manakov (1976),
and has remarkable exact solutions for Burgers equation ($f(u)=u^2/2$). This
suggests the kinetic equations of shock clustering are completely integrable.
