The functions in operators.py (eg ufl.grad) are nice user-friendly functions. There isn't a corresponding nice function for ufl.differentiation.ReferenceGrad.
I propose that we add a function reference_grad to operators.py (and import it into init.py)