Dihadron Fragmentation Functions from the NJL-jet model and their QCD Evolution

We present results for dihadron fragmentation functions from the NJL-jet model evolved from the model scale to a typical experimental scale of $4 \mathrm{GeV}^2$. The numerical method used in this evolution is discussed in detail. The effect of evolution on the shapes of the dihadron fragmentation functions is discussed for a doubly favored process ($u \to \pi^+\pi^-$), as well as a singly favored ($u \to \pi^+K^-$) process involving light quarks. Finally, we explore the production of $K^+ K^-$ pairs from an initial $u$, $d$ or $s$ quark.