I've learned that in canonical quantization you take a Lagrangian, transform to a Hamiltonian and then "put the hat on" the fields (make them an operator). Then you can derive the equations of motion of the Hamiltonian.
What is the reason that you cannot already put hats in the Dirac Lagrangian? Therefore write the Lagrangian with operators and go straight to Euler-Lagrange equations. Or is there any way this has been tried?
As mentioned in the comments you'd have to make an adjustment to get a real scalar from the operator expression. Can it be done right in the Lagrangian, without artificially putting hats on after a transformation?