But this type of Lagrangians presents a huge drawback: we don't know how to quantize them in an exact way.* Things get messy with interactions. We only can talk, with some rigor of asymptotic states: states long before or long after any interactions resemble those of the free fields. Therefore, the real photon is the excitation of thequantumthe quantum electromagnetic potential that in the limit $t\to \pm\infty$ tends to the free photon as defined above.
This answer is only a sketch, a complete answer would require me to write a book on the topic. If you want to know more, I encourage you to read any book on QFT, like Peskin & Schroeder, Weinberg, Srednicki, etc.
* In an interacting theory, the classical equations of motion are non-linear, and can't be solved using a Fourier expansion that produces creation and annihilation operators. In the path integral formulation, we only know how to solve Gaussian integrals (i.e. free fields). To solve the path integrals for interacting fields we still need approximate methods like perturbative expansions or lattice QFT. According to Peskin & Schroeder:
No exactly solvable interacting field theories are known in more than two spacetime dimensions, and even there the solvable models involve special symmetries and considerable technical complication.