Weinberg is not the best person to teach you Noether Theorem. He usually uses the so called Levy and Gell-Mann procedure of the 1960. The procedure is useful for giving a quick way of finding the Noether currents and their divergences even in cases where the symmetry is not exact. So, the procedure is about finding the currents not proving the Noether theorem. It works fine for internal symmetries but it get very messy in the more complicated space-time symmetries such as conformal and super-Poincare symmetries. Below, I will explain the Levy & Gell-Mann method. But before I do that, I need to say some thing about the First Noether Theorem. When reading textbooks or listening to your professor, you often here the sentence “continuous symmetry implies conservation law”. This is, however, not what Noether proved in her first theorem. The correct statement of the theorem is “Global continuous symmetry implies, and is implied by, an identity, we now call The Noether Identity”. Mathematically, the first theorem says: giving a global group of transformations:
<br />
\varphi_{r}(x) \rightarrow\varphi_{r}(x) + \Gamma_{r}{}^{a}( \varphi , \partial \varphi) \ \epsilon_{a},<br />
where \epsilon_{a} are the infinitesimal constant parameters of the group and \Gamma's are Lie algebra-valued functions, then
<br />
\{ \frac{\delta S[ \varphi ]}{\delta \epsilon_{a}} = 0 \} \Longleftrightarrow \{ \frac{ \delta \mathcal{L}}{\delta \varphi_{s}} \ \Gamma_{s}{}^{a} (\varphi) + \partial^{\mu}\mathcal{J}_{\mu}^{a}=0 \}.<br />
where
<br />
\mathcal{J}^{a}{}_{\mu} = \frac{ \partial \mathcal{L}}{ \partial ( \partial^{\mu} \varphi_{r} )} \ \Gamma^{a}{}_{r}( \varphi ) + X^{a}{}_{\mu}(x) \mathcal{L},<br />
is the Noether current,X^{a}{}_{\mu} is the transformation matrix: \delta x^{\mu} = X^{\mu a}(x) \ \epsilon_{a},
<br />
\frac{\delta \mathcal{L}}{\delta \varphi_{r}} \equiv \frac{\partial \mathcal{L}}{\partial \varphi_{r}} - \partial_{\mu} \left( \frac{\partial \mathcal{L}}{\partial \left( \partial_{\mu}\varphi_{r}\right)}\right),<br />
is the Euler derivative, and
S[ \varphi ] = \int_{\Omega}d^{n}x \ \mathcal{L}( x, \varphi, \partial \varphi ),
is the action integral over simply contractible domain \Omega \subseteq \mathbb{R}^{n}. As you can see the statement “symmetry implies conserved current” is just the on-shell statement of the Noether identity, i.e. when the fields satisfy the Euler-Lagrange equations of motion
\frac{\delta \mathcal{L}}{\delta \varphi_{r}}=0.
The proof is easily done by calculating the variational derivative of the action integral. This is made up of the sum of the variation of the Lagrangian and of the variation in the region of integration;
\frac{ \delta S[\varphi]}{\delta \epsilon_{a}} = \int d^{n}x \ \frac{\delta \mathcal{L}}{\delta \epsilon_{a}} + \int \ \mathcal{L} \ \frac{ \delta }{ \delta \epsilon_{a}}( d^{n}x ),
where, from the Jacobian of the transformation,
<br />
\frac{\delta}{\delta \epsilon_{a}}\ ( d^{n}x ) \approx d^{n}x \ \partial_{\mu}( \frac{ \delta x^{\mu}}{\delta \epsilon_{a}}) = d^{n}x \ \partial_{\mu}X^{\mu a}, \ \ (J)<br />
and
<br />
\frac{\delta\mathcal{L}}{\delta \epsilon_{a}}= \frac{ \partial \mathcal{L}}{ \partial x^{\mu}}\ \frac{\delta x^{\mu}}{\delta \epsilon_{a}} + \frac{ \partial \mathcal{L}}{\partial \varphi_{r}}\ \frac{\delta \varphi_{r}}{\delta \epsilon_{a}}+ \frac{\partial \mathcal{L}}{ \partial ( \partial_{\mu} \varphi_{r})} \ \partial_{\mu}( \frac{\delta \varphi_{r}}{ \delta \epsilon_{a}} ),<br />
or
<br />
\frac{ \delta \mathcal{L}}{ \delta \epsilon_{a}} = X^{\mu a}\ \partial_{\mu}\mathcal{L} + \{ \frac{ \partial \mathcal{L}}{ \partial \varphi_{r}} - \partial_{\mu}( \frac{\partial \mathcal{L}}{ \partial ( \partial_{\mu} \varphi_{r} ) } )\} \Gamma^{a}{}_{r} + \partial_{\mu} \left( \frac{\partial \mathcal{L}}{\partial ( \partial_{\mu} \varphi_{r} ) } \ \Gamma^{a}{}_{r} \right)<br />
Putting this and eq(J) in the variation derivative of the action, we find
<br />
\frac{\delta S[ \varphi ]}{\delta \epsilon_{a}} = \int_{\Omega} d^{n}x \ \left[ \frac{\delta\mathcal{L}}{\delta\varphi_{r}} \ \Gamma^{a}{}_{r} + \partial_{\mu} \left( X^{\mu a} \ \mathcal{L} + \frac{\partial \mathcal{L}}{\partial ( \partial_{\mu} \varphi_{r} ) }\ \Gamma^{a}{}_{r}\right) \right] .<br />
The statement of Noether first theorem follows from the above equation, because \Omega is simply contractible.
Ok, let us try to explain the Levy and Gell-Mann method of generating Noether currents of the global symmetries of the Lagrangian \mathcal{L}(\varphi, \partial \varphi ). It is assumed that for constant \epsilon the Lagrangian is invariant, that is
\frac{\delta \mathcal{L}}{\delta \epsilon} = 0.
Gell-Mann and Levy then considered the transformation
<br />
\varphi \rightarrow \varphi + \epsilon (x) \ F( \varphi ),<br />
and studied the change in \mathcal{L}. From the fact that \epsilon (x) is infinitesimal, so that \epsilon^{2}(x) \approx 0, it follows that both ( \delta \mathcal{L}/ \delta \epsilon ) and ( \delta \mathcal{L} / \delta \partial_{\mu}\epsilon) are independent of \epsilon. So, they could define the object
<br />
J^{\mu}(x) \equiv \frac{\delta \mathcal{L}}{ \delta ( \partial_{\mu}\epsilon )},<br />
and then prove that
\frac{\delta\mathcal{L}}{\delta \epsilon} = \partial_{\mu} \left( \frac{ \delta \mathcal{L}}{ \delta ( \partial_{\mu}\epsilon )} \right) \equiv \partial_{\mu}J^{\mu},
when the field satisfies the Euler-Lagrange equation. In general, one can show that
<br />
\frac{\delta \mathcal{L}}{\delta \epsilon(x)} - \partial_{\mu}\left( \frac{\delta \mathcal{L}}{\delta ( \partial_{\mu}\epsilon )}\right) = \left[ \frac{\partial \mathcal{L}}{\partial \varphi} - \partial_{\mu} \left( \frac{ \partial \mathcal{L}}{\partial ( \partial_{\mu} \varphi )}\right) \right] \ F( \varphi )<br />
Since for constant \epsilon the transformation above is a global symmetry of the Lagrangian, \delta \mathcal{L}/ \delta \epsilon = 0, thus the “current” constructed for non-constant \epsilon is the conserved current of the global symmetry. If the global transformation does not leave \mathcal{L} invariant, then the divergence of the associated current is given directly by
<br />
\partial_{\mu}J^{\mu} = \frac{\delta \mathcal{L}}{ \delta \epsilon (x)}.<br />
So, when you make \epsilon \rightarrow \epsilon (x) and calculate the change in \mathcal{L}, the Noether current of the global transformation will show up as the coefficient of \partial_{\mu}\epsilon, while \partial_{\mu}J^{\mu} appears with \epsilon (x),
<br />
\delta \mathcal{L} = J^{\mu}\ \partial_{\mu}\epsilon + \epsilon \ \partial_{\mu}J^{\mu}.<br />
Sam