use deltaR2: no need to call the extra sqrt here to get deltaR and waste cpu cycles.
Also, better compute eta and phi arrays for these jets in one loop so that they are not recomputed again in this nested loop