If $f$ is a modular form of weight $k$, it is well known that $$ D(f)=f' -\tfrac k{12}E_2f $$ is modular of weight $k+2$. Here $E_2$ is the Eisenstein series. I wanted to ask if there is an extension of this fact for Hilbert modular forms.
When looking up Hilbert modular forms (Shimura 1975), I was only able to find the following differential operators: for a Hilbert modular form $f$ on $\mathbb H^n$, we define $$ D_{j,t}(f)=\left(\frac{t}{2i y_j}+\frac{\partial}{\partial z_j}\right)f $$ where $z_j=x_j+iy_j$ is the $j$-th coordinate. I don't think this is the same operator (is it?)
I am really looking for a many variable version of the operator $D(f)=f' -\frac k{12}E_2f$. I don't really have a background in analysis or number theory so this question may be really easy for experts.