CFD Online Discussion Forums - PLM for Euler equations

I'm trying to add second order reconstruction based on the following Taylor expansions for the edge states from some advection notes:

$a^{n+1/2}_{i+1/2,L} = a^n_i + \frac{\Delta x}{2} \frac{\partial a}{\partial x}\Big| _i + \frac{\Delta t}{2} \frac{\partial a}{\partial t}\Big| _i$

$a^{n+1/2}_{i+1/2,R} = a^n_{i+1} + \frac{\Delta x}{2} \frac{\partial a}{\partial x}\Big| _{i+1} + \frac{\Delta t}{2} \frac{\partial a}{\partial t}\Big| _{i+1}$

My plan is to apply these two equations to each of my conserved quantities (rho, rho*u, E) and use the spatial derivative of the corresponding quantity in the Euler equations (rho*u, rho*u*u+p, u(E+p)) for the time derivative.
Like so

$(\rho u)^{n+1/2}_{i+1/2,L} = (\rho u)^n_i + \frac{\Delta x}{2} \frac{\partial (\rho u)}{\partial x}\Big| _i + \frac{\Delta t}{2} \Big(-\frac{\partial(\rho u^2 + p)}{\partial x}\Big)\Big| _i$

Then use something like centered differences for each of the spatial derivatives before converting back to the primitive variables for the Riemann solver.
However, I've never seen it done like this so I wonder if this is wrong or there's a better way to implement piecewise linear reconstruction without using eigenvalues/characteristics?

Edit: I'm using finite volume discretization.