We consider the Laplacian with drift in $\mathbb R^n$ defined by $\Delta_\nu = \sum_{i=1}^n(\frac{\partial^2}{\partial x_i^2} + 2 \nu_i\frac{\partial }{\partial{x_i}})$ where $\nu=(\nu_1,\ldots,\nu_n)\in \mathbb R^n\setminus\{0\}$. The operator $\Delta_\nu$ is selfadjoint with respect to the measure $d\mu_\nu(x)=e^{2\langle\nu,x\rangle}dx$. This measure is not doubling but it is locally doubling in $\mathbb R^n$. We define, for every $M>0$ and $k \in \mathbb N$, the operators $$ W^k_{\nu,M,*}(f) = \sup_{t>0}\left|A^k_{\nu,M,t}(f)\right|,\hspace{5mm}g_{\nu,M}^k(f) = \left(\int_0^\infty\left|A^k_{\nu,M,t}(f)\right|^2\frac{dt}{t}\right)^{\frac{1}{2}},\,k\geq 1, $$ the $\rho$-variation operator $$ V_\rho\left( \{A^k_{\nu,M,t}\}_{t>0}\right)(f)= \sup_{0
2, $$ and, if $\{t_j\}_{j\in \mathbb N}$ is a decreasing sequence in $(0,\infty)$, the oscillation operator $$ O(\{A_{\nu,M,t}^k\}_{t>0},\{t_j\}_{j\in \mathbb N})(f)=\Big(\sum_{j\in \mathbb N}\;\;\sup_{t_{j+1}\leq \varepsilon <\varepsilon '\leq t_j}|A^k_{\nu,M,\varepsilon}(f)-A^k_{\nu,M,\varepsilon '}(f)|^2 \Big)^{1/2}. $$ where $A^k_{\nu,M,t}=t^k\partial^k_t(I-t\Delta_\nu)^{-M}$, $t>0$. We denote by $T_{\nu,M}^k$ any of the above operators. We analyze the boundedness of $T^k_{\nu,M}$ on $L^p(\mathbb R^n,\mu_\nu)$ into itself, for every $1