p(epsilon_a^1,epsilon_b^1) = delta(J, Tau)