This can also get simplified by using `deltaR2` instead: just pass as argument the square of the `dR`  The *_DR parameters could be replaced by some *_DR2 ones, and as such the same simplification could be applied also (I think) in all other cases in which deltaR is used in this producer.