  We compute cross sections for two-jet production in deep inelastic scattering
(DIS), with one jet from initial state radiation (ISR) and the other from final
state radiation, with a summation of large logarithms up to
next-to-next-to-leading logarithmic (NNLL) accuracy. Use of the DIS event shape
1-jettiness ensures that events have two well-collimated jets. We calculate
distributions for three versions of 1-jettiness that have different sensitivity
to the transverse momentum of the ISR, and derive factorization theorems for
each of them using the soft collinear effective theory (SCET). The structure of
the transverse momentum dependence in the factorization theorems is different
for each 1-jettiness. We present numerical results for these three observables
with parameters for the HERA collider.
