Abstract. We obtain a model-independent expression for the Lamb shift in muonic hydrogen. This expres- The latter are controlled by the chiral theory, which allows for their model-independent determination. In this paper we give the missing piece for their complete expression including the pion and Delta particles. Out of this analysis, and the experimental measurement of the Lamb shift in muonic hydrogen, we determine the electromagnetic proton radius: r p = 0.8412 (15) fm. This number is at 6.8σ variance with respect to the CODATA value. The accuracy of our result is limited by uncomputed terms of OThis parametric control of the uncertainties allows us to obtain a model-independent estimate of the error, which is dominated by hadronic effects.The recent measurement [1,2] of the Lamb shift in muonic hydrogen, E(2P 3/2 ) − E(2S 1/2 ),and the associated determination of the root mean square electric radius of the proton: r p = 0.84087(39) fm has led to a lot of controversy. The reason is that this number is 7.1σ away from the CODATA value, r p = 0.8775(51) fm [3]. This last number is an average of determinations coming from hydrogen spectroscopy and electron-proton scattering. It should be mentioned though that the latter have been recently been challenged in refs. [4,5], and its inclusion would certainly diminish this tension. Leaving this aside, in order to asses the significance of the discrepancy, it is of fundamental importance to perform the computation (in particular of the errors) in a model-independent way. In this letter we revisit the theoretical derivation of the Lamb shift in muonic hydrogen with this aim in mind. In this respect, the use of effective field theories is specially useful. They help organizing the computation by providing with power counting rules that asses the importance of the different contributions. This becomes increasingly necessary as higher-order effects are included. Even more important, these power a e-mail: peset@ifae.es b e-mail: pineda@ifae.es counting rules allow to parametrically control the size of the uncalculated terms and, thus, give an educated estimate of the error. This discussion specially applies to the muonic hydrogen, as its dynamics is characterized by several scales:By considering ratios between them, the main expansion parameters are obtained:This approach to the problem has been followed in [6][7][8] (see [9] for a review of these computations) with a combined use of Heavy Baryon Chiral Perturbation Theory (HBChPT) [10] (see also [11]), Non-Relativistic QED (NRQED) [12] and, specially, potential NRQED (pNRQED) [13][14][15]. Particularly relevant for us is ref. [7], which contains detailed information on the application of pNRQED to the muonic hydrogen. We refer to it for details (and to [16] where a more detailed account of the hadronic computation presented here is given).Since pNRQED describes degrees of freedom with E ∼ m r α 2 , any other degree of freedom with larger energy is integrated out. This implies treating the proton