Neutrino-Nucleus scattering in the SuSA model

The Super-Scaling Approach (SuSA) model, based on the analogies between electron and neutrino interactions with nuclei, is reviewed and its application to the description of neutrino-nucleus scattering is presented. The contribution of both one- and two-body relativistic currents is considered. A selection of results is presented where theoretical predictions are compared with cross section measurements from the main ongoing neutrino oscillation experiments.


Introduction
A huge experimental programme on accelerator-based neutrino oscillation experiments has been developed in recent years with the scope of improving our knowledge of neutrino properties and of potentially answering one of the fundamental open questions of modern physics, that is the origin of the matter/antimatter asymmetry in the universe, by measuring the weak CP-violating phase [1]. The success of this program partly relies on the control of systematic errors, largely due to uncertainties in the description of the neutrino interactions with the detector, typically made of medium-weight nuclei like carbon, oxygen or argon. With this motivation, intense theoretical activity has been carried out in parallel in order to provide an accurate description of neutrino-nucleus reactions in the GeV region, relevant for the experiments, for different processes ranging from quasielastic up to deep inelastic scattering, encompassing the excitation of nucleon resonances and the emission of two or more nucleons (see [2] for a comprehensive review of recent progresses and open challenges in the field).
In this article we review our work on the SuSA (Super-Scaling Approach) model, originally developed in Ref. [3] and subsequently refined and improved to the updated version SuSAv2 [4,5] (see also Refs. [6,7] for other recent reviews). The original SuSA, as will be explained in more detail in Sect. 2.2, is a phenomenological model which, while retaining the relativistic aspect of the Relativistic Fermi Gas (RFG), provides by construction -unlike the RFG -a good description of inclusive electron scattering data (e, e ′ ) in the quasielastic region. In the model initial and final state interaction effects, absent in the RFG, are directly extracted from (e, e ′ ) data. The updated version, SuSAv2, implements inputs from the relativistic mean field model (RMF), which provides a microscopic interpretation of the basic features of SuSA and also includes a more detailed description of the different spin and isospin channels. Moreover, the model has been extended to include two-body currents, able to excite two-particle-two-hole (2p2h) states; these, as will be shown in Sect. 4, play an important role in the analysis of neutrino-nucleus experimental data.
The scheme of this review is the following. In Sect. 2 we focus on the inclusive reaction, presenting the formalism, describing the model and testing it with electron scattering data. In Sect. 3 we discuss the semi-inclusive reaction an its relation to the inclusive one. In Sect. 4 we show some selected results and discuss the comparison with experimental data. Finally, in Sect. 5 we draw our conclusion and outline the future developments of our work.

.1 General formalism
In this review we mainly focus on the charged current (CC) inclusive process where a neutrino (antineutrino) of a given flavour l having four-momentum K µ = (ε, k) hits a nucleus A and a charged lepton l − (l + ) is detected in the final state with four-momentum K ′µ = (ε ′ , k ′ ) and scattering angle θ, while the residual system X is unobserved. The four-momentum transferred from the probe to the nucleus is Q µ = (ω, q), with ω = ε − ε ′ and q = k − k ′ . We shall assume the neutrino to be massless, while the outgoing lepton has finite mass m l . The corresponding cross section is obtained by contracting the leptonic and hadronic tensors W µν (q, ω) = n < A|J µ † (q, ω)|n >< n|J ν (q, ω)|A > δ(ω + EA − En), where χ = +1(−1) in the neutrino (antineutrino) case, |A > is the initial nuclear ground state and |n > are all the intermediate hadronic states that can be reached through the weak current operator J µ . The double differential cross section with respect to the momentum and scattering angle of the outgoing lepton can be expressed in terms of five response functions [3] (for the semi-inclusive reaction (ν l , l − N ), which will be briefly discussed in Sect. 3, the response functions are ten) where GF the Fermi weak constant, θc the Cabibbo angle and v0 = (ε+ε ′ ) 2 −q 2 . The indices C, L, T refer to the Coulomb, longitudinal and transverse components of the leptonic and hadronic currents with respect to q. The leptonic coefficients VK are related to the components of the tensor (2) and are given by where the variables δ 2 = m 2 l |Q 2 | , ρ = |Q 2 | q 2 and ρ ′ = q ε+ε ′ have been introduced. The response functions (6) embody the nuclear dynamics and are specific components of the hadronic tensor W µν , depending only upon ω and q .
The above decomposition into response functions is valid for all reaction channels -elastic, quasielastic, inelastic -each of which characterized by a different current operator J µ in the hadronic tensor (3). In the case of quasielastic scattering, corresponding to the interaction of the probe with a single nucleon, the weak current operator is where F1, F2, GA and GP are the Pauli, Dirac, axial and pseudoscalar weak form factors, respectively. In neutrino oscillation experiments the incident beam is not monochromatic. As a consequence, before comparing with experimental data, the cross section (4) must be folded with the normalized neutrino flux φ(ε) ( dε φ(ε) = 1):

Scaling and SuSAv2
The phenomenon of scaling is sometimes found to occur in different fields, including solid state, atomic, molecular, nuclear and hadronic physics, when an interacting probe scatters from composite many-body systems. Under some circumstances it is found that the response of the complex system no longer depends on two independent variables, but only on a particular combination of those, called the scaling variable. This phenomenon is very well known in high energy physics where the inelastic nucleon response functions are shown to depend only on the Bjorken variable x [8]. A similar phenomenon is observed in lepton-nucleus scattering, where the nuclear response functions are found to scale with a single scaling variable, denoted as y. This indicates that the probe (lepton) interacts with the nucleus' constituents, in this case with the nucleons in the nucleus. The phenomenon of y-scaling emerges from the analysis of quasielastic (QE) (e, e ′ ) reactions and has been studied in detail in [9,10,11]. The scaling function is defined as the QE (e, e ′ ) differential cross section divided by a single-nucleon cross section averaged over the Fermi gas. For high enough values of the momentum transfer, q, the scaling function does only depend on a single variable, y, given as a particular combination of the energy (ω) and momentum (q) transferred in the process. The scaling variable y is (up to a sign) the minimum value of the missing momentum allowed by kinematics [9,10,11]. In the QE domain the basic mechanism in (e, e ′ ) reactions on nuclei corresponds to elastic scattering from individual nucleons in the nuclear medium. This implies that the inclusive (e, e ′ ) cross section is mainly constructed from the exclusive (e, e ′ N ) process, including the contribution of all nucleons in the target and integrating over all (unobserved) ejected nucleon variables (see also the following section). This approach constitutes the basis of the Impulse Approximation (IA). Thus, the double differential (e, e ′ ) inclusive cross section is given as the sum of two response functions corresponding to the longitudinal, L, and transverse, T , channels, where σMott is the Mott cross section and the vL,T are kinematical factors that involve leptonic variables (see [6,12] for explicit expressions). In terms of the scaling functions the nuclear responses are RL,T (q, ω) = 1 kF fL,T (q, ω)GL,T (q, ω) , where kF is the Fermi momentum and GL,T are defined as the responses of a moving nucleon and include relativistic corrections arising from the presence of the medium. Their explicit expressions can be found in [4,6,13]. The terms fL,T are the scaling functions that show a mild dependence upon the momentum transfer q (first-kind scaling) and very weak dependence on the nuclear system considered (second-kind scaling). The occurrence of both types of scaling is denoted as superscaling.
Hence to the extent that at some kinematics the above f functions are the same for all nuclei and do only depend on a single variable, denoted as the scaling variable ψ(q, ω), but not separately on q, one says that superscaling occurs. The superscaling variable ψ is given in terms of q and ω as where the dimensionless Fermi kinetic energy ξF = T F m N , energy transfer λ = ω 2m N , momentum transfer κ = q 2m N and squared four-momentum transfer τ = κ 2 − λ 2 have been introduced. Here the y-variable introduced above is given approximately by as discussed in [11]. A phenomenological energy shift E shif t , fitted for each nucleus to the experimental position of the quasielastic peak [14], is included in the definition of the scaling variable ψ ′ ≡ ψ(q, ω − E shif t ).
Scaling and superscaling properties of electron-nucleus interactions have been analyzed in detail in a series of previous works [9,10,11,14,15,16,17] . The importance of this phenomenon to test the validity of any nuclear model aiming to describe electron scattering reactions has been clearly proven. The model, denoted as the Superscaling Approach (SuSA), is entirely based on the phenomenology, making use of a unique, universal, scaling function extracted from the analysis of the longitudinal electron scattering data (see Fig. 1). Notice that the behavior and properties of the experimental superscaling function constitute a strong constraint for any theoretical model describing QE electron scattering. Not only should the superscaling behavior be fulfilled, but also the specific shape of the longitudinal scaling function, f exp L , must be reproduced. The SuSA model assumes the longitudinal phenomenological scaling function to be valid also in describing the transverse channel, i.e., fL = fT .   [18]. The RFG scaling function is also shown as reference. Right panel: RMF longitudinal scaling functions for (e, e ′ ) at different q values compared with (e, e ′ ) longitudinal scaling data from [18]. In the inner smaller plot a reduced-χ 2 analyses shows a minimum at q = 650 MeV/c. Figure adapted from [4].
In recent years we have explored in detail the extension of SuSA to neutrino-nucleus scattering. Due to the complexity of the weak process, with an increased number of nuclear responses compared with the pure electromagnetic one, and the particular role played by the axial term in the weak current, we have developed an improved version of the superscaling model, called SuSAv2 [4]. Contrary to the original SuSA [3,13,19] where a universal scaling function based on electron scattering data is used, the new SuSAv2 model incorporates relativistic mean field (RMF) effects [16,17] in the longitudinal and transverse nuclear responses, as well as in the isovector and isoscalar channels that is of great importance in order to describe charged-current (CC) neutrino reactions that are purely isovector [20].
The origin of the SuSAv2 approach is based on the capability of the RMF to describe properly the scaling behavior of the electron scattering data. As shown in previous works [16,17], RMF is one of the few microscopic models capable of reproducing the asymmetric shape of the phenomenological scaling function with a long tail extended to high values of the transfer energy (large values of ψ ′ ). Moreover, RMF produces an enhancement in the transverse scaling function, a genuine dynamical relativistic effect linked to the lower components in the wave functions, that is supported by the analysis of data. The RMF framework to finite nuclei has proven to successfully reproduce the scaling behavior shown by data at low to intermediate q values (see Fig. 1). However, the model clearly fails at higher momentum transfers where Final State Interactions (FSI) are expected to be weaker. This is due to the RMF strong energy-independent scalar and vector potentials used in the final state that lead to too much asymmetry in the scaling functions and shift the QE peak to very high transfer energies, in clear disagreement with data. To remedy this shortfall of the RMF model, the SuSAv2 incorporates both the RMF scaling functions at low-to-intermediate q values and the Relativistic Plane Wave Impulse Approximation (RPWIA) ones at higher q by using a q-dependent blending function that smoothly connects the two regimes (see [5,21] for details). A similar solution to this drawback of the RMF model has been taken in the recent Energy-Dependent RMF (ED-RMF) approach [22,23] where RMF potentials are multiplied by a blending function inspired by the SuSAv2 one that scales them down as the kinetic energy of the scattered nucleon increases, also preventing non-orthogonality issues. This model predicts both lepton and nucleon kinematics, showing a similar agreement on electron and neutrino data with SuSAv2.
In summary, the SuSAv2 model and, for extension, the ED-RMF one reproduce the experimental longitudinal scaling data, gives rise to an enhancement in the electromagnetic transverse channel, i.e., f (e,e ′ ) T > f (e,e ′ ) L , takes into account the differences in the isoscalar/isovector scaling functions, of crucial interest for neutrino scattering processes, and finally avoids the problems of the RMF model in the region of high momentum transfer, where FSI effects are negligible. One of the basic merits of SuSAv2 is the translation of sophisticated and demanding microscopic calculations into relatively straightforward parametrizations and, hence, easing its implementation in the MonteCarlo simulations employed in the analysis of neutrino oscillation experiments.
To conclude, it is important to point out that SuSAv2 is not restricted to the QE kinematics domain. On the contrary, it can be extended to the inelastic region by generalizing the superscaling variable (11) to the excitation of any inelastic state having invariant mass W . One can thus define [24] and evaluate the inelastic hadronic tensor as the integral between the limits Wmin, Wmax imposed by the kinematics, where f (ψW ) is the SuSAv2 superscaling function previously described. This procedure is based on the assumption that the nuclear effects, encoded in the function f , are the same in the QE and inelastic regions and requires the knowledge of the elementary tensor w µν inel across the full inelastic spectrum.  [27,28] (solid lines) and GRV98 PDFs [29,30] (dashed lines) at different kinematics (incident electron beam and scattering angle) in terms of the energy transferred to the nucleus (ω). Experimental data taken from [31,32]. The y-axis represents d 2 σ/dΩ/dω in nb/GeV/sr. The value of q at the QE peak (q QE ) is shown as reference. The SuSAv2 inelastic model has been applied to the electron scattering case using phenomenological fits of the single-nucleon inelastic structure functions w1 and w2 extracted from e-p and e-d data [25,26]. In Fig. 2, a comparison of the SuSAv2-inelastic model with 12 C(e, e ′ ) data using different inelastic structure functions shows a preference for the Bosted-Christy parametrization [25,26]. In general, the SuSAv2 model (QE+inelastic) provides a very good description of data for very different kinematics once Meson Exchange Currents (MEC) are also incorporated, denoted as SuSAv2-MEC (see Sect. 2.4 and [5,33,34]). The inclusion of the full inelastic spectrum in the SuSAv2 model for weak interactions is still in progress.
An alternative approach to the study of the resonant pion production was taken in Refs. [3,35,36], where a scaling function to be used in the ∆ resonance region, different from the quasielastic one, was extracted from electron scattering data and multiplied by the appropriate weak N → ∆ transition form factors to get the neutrino-nucleus cross section in this region. This method provides a phenomenological description valid at transferred energies below the ∆ peak, while at higher ω it fails due to the opening of other inelastic channels. Some results corresponding to this method, referred to as SuSA-∆ approach, will be shown in Sect. 4.

Relativistic model for CC MEC and 2p2h responses
Multinucleon knockout processes give a non-negligible contribution to the inclusive neutrino cross section for the intermediate energies involved in the experiments [5,37,38,39,40]. In Ref. [41] we developed a model of two-particle two-hole excitations of the RFG induced by weak meson-exchange currents for inclusive CC neutrino scattering. The model is fully relativistic and includes the diagrams of Fig. 3, involving one-pion exchange and ∆ excitation, taken from the pion production model of [42].
The 2p2h matrix element of MEC depends on the momenta, spin and isospin coordinates ( ) of the two holes, h1, h2, and the two particles, p ′ 1 , p ′ 2 . It is the sum of four contributions corresponding in Fig. 3 to the seagull (diagrams a,b), pion in flight (c), pion-pole (d,e) and ∆(1232) excitation (f,g,h,i). Their explicit expressions are given in Ref. [41]. The inclusive hadronic tensor in the 2p-2h channel is computed by integration over all the 2p-2h excitations of the RFG. Momentum conservation enforces where j µ (1 ′ , 2 ′ , 1, 2)A is the antisymetrized MEC matrix element The factor 1/4 in Eq. (17) accounts for the antisymmetry of the two-body wave function in isospin formalism, to avoid double counting in the number of final 2p-2h states. Due to azimuthal symmetry around the z axis -in the q direction-we fix the azimuthal angle of particle 1' φ ′ 1 = 0, and multiply by a factor 2π. The energy delta-function enables integrating over p ′ 1 . Then Eq. (16) is reduced to a seven dimensions integral that is computed numerically [43,44]. The Dirac matrix elements of the currents are also computed numerically. The 2p-2h inclusive cross section requires one to compute the five weak response functions, R CC,CL,LL,T,T ′ , for (νµ, µ − ). All of these responses were computed and analyzed in Ref. [41]. The five response functions have been parametrized in the kinematic range 100 < q < 2000 [5]; the parametrization is convenient because in neutrino scattering there is an additional integration over the incident neutrino flux. This parametrization has been implemented in the Monte Carlo event generator GENIE [45,34].
Using the parametrization of [5] does not allow one to modify the internal parameters of the MEC model. This is why we have developed approximations to the 2p2h responses in order to speed up the calculation. Specifically, in the so-called frozen nucleon approximation we set the momenta of the two holes h1 = h2 = 0, allowing us to integrate over the initial states analytically [46]. In the modified convolution approximation (MCA) [47] we write the 2p2h responses as a convolution of two 1p-1h responses multiplied by the elementary hadronic tensor w µν (p ′ 1 , p ′ 2 ; h1, h2) evaluated for convenient averaged values of h1 and h2. The resulting MCA responses are a good approximation to the exact result with 4D integration only over the momentum and energy comunicated to one of the nucleons.

Validation vs electron scattering
The model introduced in the previous sections allows one to describe inclusive lepton-nucleus scattering in the kinematic region including the quasielastic reaction, the excitation of 2p2h states and the inelastic spectrum. Validation against the large amount of existing high-precision electron scattering data is a necessary test to be performed before using nuclear models in the analysis of neutrino oscillation experiments. Such tests have been successfully carried out for the SuSAv2 model in a very wide range of kinematics from low/intermediate energies up to the highly inelastic regime. Only at very low energy and momentum transfers (q <300 MeV/c and ω <50 MeV) does the model fail to reproduce the data, as do all models based on the impulse approximation, which is clearly not appropriate in the low-energy regime. SuSAv2 predictions have been extensively compared with the available (e, e ′ ) data on carbon and oxygen in [21,33] and with the more recent JLab data on argon and titanium in [48]. Here we just show a few representative examples, choosing the kinematics of particular interest for ongoing neutrino experiments.
In Fig. 4 (top left panel) the double differential 12 C(e, e ′ ) cross section is shown as a function of the energy transfer ω, for an incident electron beam of energy E and scattering angle θ, corresponding to the momentum transfer q QE at the quasielastic peak. The separate QE, 2p2h and inelastic contributions are shown. The agreement with the data is excellent across the full spectrum. Note that 2p2h excitations are essential in order to describe the "dip" region between the QE and ∆ peaks.
An important feature of the SuSAv2 model, based on the RMF theory, is its capability to reproduce not only cross sections but also the separated longitudinal and transverse responses, as shown in the top right panel of Fig. 4, where RL (solid) and RT (dashed) are plotted versus ω for q =570 MeV/c, together with the separate QE, 2p2h and inelastic contributions. Note that while the longitudinal response is almost purely quasielastic, with a small 2p2h contributions arising from relativistic corrections, the transverse response receives contributions from all the three processes, which sizeably overlap at these kinematics.
In Fig. 4 (bottom panels) the same kind of comparison is shown for the oxygen (left) and argon (right) nuclei, showing that the superscaling approach allows for a consistent description of different nuclei.
These results give us confidence on the reliability of the model for the application to neutrino scattering. Before showing the comparison with neutrino data, in the next section we address the study of semi-inclusive electron scattering.

Semi-inclusive scattering
A pivotal difference between electron and neutrino scattering experiments is that the neutrino energy is known only as a broad distribution, while the energy in electron beams is typically defined with high precision. Hence, in electron scattering experiments by, for instance, detecting one nucleon in coincidence with the scattered electron, the hadronic final state is completely determined if one restricts one's attention to excitation energies of the residual system below the two-nucleon emission threshold. On the contrary, in the neutrino case, even in scenarios in which one or a few hadrons are detected in coincidence with the scattered lepton, there is no way of fully determining the hadronic final state because the beam energy and consequently the excitation energy of the residual system is unknown. This means that one needs to 'integrate' for all possible final states compatible with the given kinematics. Formally, this integration is equivalent to the integration over the neutrino flux (or flux folding procedure) introduced in Sect. 2.1.
To obtain a somewhat deeper understanding of how the inclusive cross sections and scaling emerge, in this section we briefly discuss semi-inclusive (coincidence) electron scattering, focusing on the reaction (e, e ′ N ) in which the scattering electron and a nucleon are assumed to be detected in the final state. Then, in addition to the electron kinematical variables introduced above, we have an outgoing nucleon with 4-momentum P µ N = (EN , pN ) involving 3-momentum pN and polar and azimuthal angles θN and φN , respectively, together with energy EN = p 2 N + m 2 N . No other particles are assumed to be detected, although, depending on the specific kinematics, they must be present (see below). The magnitude of the nucleon's 3-momentum is given by pN = |pN |. Apart from the detected nucleon, the final state contains an undetected hadronic system having missing 4-momentum (EB, pB), namely, a total energy of EB and a missing 3-momentum pB ≡ pm. In the following we shall assume that the detected nucleon is a proton. One then has The undetected hadronic system has invariant mass MB (M 0 B at threshold with MB ≥ M 0 B ) and total energy which defines the kinetic energy of the unobserved final-state system, TB. From (20) one has where M 0 A is the target ground-state mass and TN = EN − mN is the kinetic energy of the detected nucleon. This leads to an expression for the so-called missing-energy, where Es = M 0 B +mN −M 0 A is the separation energy and the (typically very small) recoil kinetic energy difference has been neglected. Defining the excitation energy of the residual system E = MB − M 0 B , we have The magnitude of the missing-momentum pm is given by Depending on the value of the missing-energy, the residual system may be the daughter nucleus in its ground state (this defines the threshold for the reaction to become possible); or it may be in a discrete excited state (lying below the threshold where a second nucleon can be ejected), and, although they de-excite by γ-decay, that process is slow on the nuclear timescale and thus these states may be treated effectively as stationary states. Then, at a well-defined threshold a second nucleon must be emitted (this is not optional: there are no nuclear states involving one nucleon and a residual bound nucleus above this point). As Em continues to increase more and more particles enter in the final state in addition to the one special nucleon that is assumed to be detected. At even larger missing energy (roughly 140 MeV) pion production becomes possible (still with the lepton and one nucleon assumed to be detected).
From (22 -24) it is clear that for fixed values of the observable parameters the values of Em and pm are determined for each value of E. In Fig. 5 trajectories are shown for selected kinematics, namely, E l = 3800 MeV, θ l = 7 deg, TN = 140 MeV, φN = 180 deg, representing a "typical" situation. What is varied here is the polar angle for the detected proton, θN . As stated above, as one goes along a given trajectory the electron energy E that determines where on the trajectory one finds oneself must vary. The lower boundary defines the threshold for the semi-inclusive reaction to  occur. In effect, each event where an electron and a proton are detected in coincidence corresponds to a specific trajectory and point on that trajectory. One sees a striking pattern to the behavior one should expect when going along a given trajectory. The strength in the Rome spectral function [53,54], which should provide a good starting point for the characteristics to be expected in semi-inclusive reactions, is extremely localized. One sees the largest concentration of strength where the p-shells are located (at around Em = 20 MeV) with less where the broad s-shell is located (at around Em = 50 MeV); at still larger values of Em (and pm) the spectral function does have some strength, although it is spread over a wide region in the Em − pm plane and is too small to be seen in this representation. Furthermore, we note that pion production cannot occur until one reaches Em ∼ mπ and that it is not appreciable until Em ∼ m∆ − mN ∼ 300 MeV.
We notice that in the case of the CCQE interaction, studied in Ref. [52], one gets a Em − pm trajectory plot that is almost identical to the one in Fig. 5. The reason is that the muon mass does not play a strong role in these trajectories, especially for the kinematics shown here for which E l is much larger than the muon mass.
In passing we note that, while other choices of two variables to replace Em and pm can of course be made, the generic behavior seen here strongly suggests that the present choice is a good one and that other choices may not reflect the highly localized nature of the nuclear response.   In Fig. 6 we represent the six-fold differential cross section for fixed electron and proton kinematics as a function of the missing-energy. The model used was PWIA together with the Rome spectral function [53,54], as described, e.g., in Refs. [51,55]. Equation (22) tells us that for fixed nucleon and final lepton energies, if one neglects the kinetic energy of the residual system TB, which under typical QE conditions is always very small, then one finds a one-to-one linear relation between the initial lepton energy and the missing energy. Thus, in Fig. 6 the beam energy is shown as a second x-axis. By varying the angle θN , the cross section changes its magnitude and its shape. One observes two prominent peaks corresponding to the p-shells, a wide bump for the s-shell and a background that extends up to high missing-energies. Clearly, as expected, the p-shell strength is largest, the s-shell strength is smaller and the high-Em strength is completely negligible, being down by several orders of magnitude. By examining these results in the light of the trajectories shown in Fig. 5 we see that the general behavior we expect to occur is borne out. For example, the trajectories for θN = 60 and 80 degrees both pass though the p-shell region near its peak. However, one trajectory intersects the s-shell region more than the other one does and this results in relatively different amounts from the s-shell compared with the p-shell. Or, large values of θN correspond to very small cross sections, as they should, since neither the p-nor s-shell regions are crossed. We do have some knowledge about this generic behavior of the distribution of strength from inclusive electron scattering, (e, e ′ ). Inclusive scattering corresponds to performing integrals over specific regions in the Em − pm plane [6,9,56]. In Fig. 7 we show the typical situation for kinematics at values of q and ω where one is below the QE peak (y < 0, left panel) or above the peak (y > 0, right panel). In each case the inclusive QE cross section is obtained by integrating the semi-inclusive (e, e ′ N ) cross sections over the shaded regions. Accordingly one can see that the two classes of cross section are intrinsically related. Clearly, if one had complete knowledge of the semi-inclusive response for a wide range of kinematics then the integrations could be performed to yield the inclusive response. Unfortunately this is not the case. Note that having a model for the inclusive cross section (the total hadronic cross section) does not mean that such a model will be valid for the semi-inclusive cross section (which constitutes the integrand of the former). An example of this is the RFG model which is not unreasonable for the total (inclusive) QE cross section, but is poor for the semi-inclusive response.
An example is that of the model used for the semi-inclusive response discussed above. If employed for the inclusive cross section one finds for 16 O that somewhat over 50% of the inclusive cross section stems from the p-shells, about 25% comes from the s-shell region and the rest comes from a broad region at higher missing-energy. This strength at higher missing energy is partially responsible for the asymmetry found in the scaling functions. Note that inclusive scattering at high momentum transfers and hence the scaling functions involve broad integrals, whereas semi-inclusive scattering at fixed final-state electron and nucleon momenta involves a trajectory in the Em − pm plane as the beam energy is varied. In detailed analyses one finds much more strength in that case coming from the valence knockout region, with much less arising from the high-Em region. This has consequences for CCν reactions as discussed in Ref. [52].

Results
In Sect. 2.4, the validation of the SuSAv2-MEC model based on the RMF theory with (e, e ′ ) data has been proven as a solid benchmark to assess the validity of a nuclear model before its application to neutrino reactions. Here, we show the capability of the SuSAv2-MEC model to describe a wide range of kinematics of interest for neutrino oscillation experiments. In particular, we focus on the comparisons of the SuSAv2-MEC model with charged-current neutrino cross sections from different experiments, in particular T2K and MINERvA. Our analysis is mainly related to charged-current quasielastic-like events, also called CC0π, which are characterized by having no pions detected in the final state and are dominated by the QE and 2p2h channels. Finally, we extend the analysis to cross section measurements where pion production is also included. 0.975 < cosθ µ < 1.000 Figure 8: Comparison of the SuSAv2-MEC model with the T2K flux-integrated CCQE and 2p2h doubledifferential cross section for neutrino for scattering on 12 C (top panels) and 16 O (bottom panels) in units of 10 −39 cm 2 /GeV per nucleon (for carbon) and per nucleon target (for oxygen). The CC0π T2K data are from Ref. [57] and [58].
In Fig. 8 we show the comparison of the SuSAv2-MEC model with T2K CC0π double differential cross sections on 12 C (top panels) and 16 O (bottom panels). These data are compared with the QE and 2p2h MEC contributions showing an overall good agreement. Monte Carlo event generators also include other contributions that can mimic a CCQE-like event in the data analysis such as pionabsorption processes in the nucleus. Nevertheless, as it will be discussed later, these contributions are not particularly relevant at T2K kinematics. It is also important to note that the QE and 2p2h contributions in the SuSAv2-MEC model can be easily extrapolated from one nucleus to another by means of scaling rules, which are different for each nuclear regime. This is based on the assumption of 2nd-kind scaling, already mentioned in Sect. 2.2, and that has been proven in the SuSAv2-MEC approach for the QE and 2p2h channels [48,59]. It is worth commenting on the case of the most forward-angle region, associated with low kinematics, i.e., low energy and momentum transferred to the nucleus. Within this region, the SuSAv2-MEC model slightly overestimates the T2K data for C and O, being more noticeable in the latter case. This is mainly due to RMF scaling violations, related to low-energy nuclear effects and to different binding energies for each nucleus, which are not properly accounted for in the SuSAv2 approach. This limitation of the model has been solved in the recent ED-RMF approach [22,23] where the goodness of the RMF strong vector and scalar potentials at low-intermediate kinematics are present while also retaining the benefits of the SuSAv2 model description at larger kinematics, as described in Sect. 2.2. The differences introduced by this more accurate description of low-energy nuclear effects in the RMF (and ED-RMF) model can be observed in Fig. 9 where large discrepancies between oxygen and carbon predictions are observed between the RMF model and the SuSAv2-MEC one at very forward angles, being smaller as we move 0.98 < cosθ µ < 1.00 Figure 9: Comparison of double differential cross sections on 12 C (solid lines) and 16 O (dashed lines) at T2K kinematics within the SuSAv2 (1p1h) and RMF models. Results are displayed from forward to backward angles.
to more backward kinematics. At the same time, these differences are more prominent in the case of oxygen. This is connected with the assumptions made in the SuSAv2 approach which is based on the RMF analysis on 12 C but on the application of scaling rules to describe other nuclear targets instead of relying on particular RMF predictions for each nucleus. Therefore, at kinematics where scaling violations are present, these differences, while noticeable in the nucleus of reference ( 12 C), can be even more important for other nuclear targets. On the contrary, SuSAv2 assumes that scaling works well even at low kinematics, thus implying minor differences between 12 C and 16 O, mainly due to the scaling rule and differences in the energy shift introduced for each nucleus. This analysis is of relevance for T2K [60] where the accuracy of the C to O extrapolation plays an important role on the oscillation analysis and whose predictions have shown some discrepancies between C and O for very forward-going muons which may be explained by the RMF predictions. 10.00 < p || < 15.00 Figure 10: (Color online) The MINERvA "QE-like" and "CCQE" double differential cross sections forν µ (top panels) and ν µ (bottom panels) scattering on hydrocarbon versus the muon transverse momentum, in bins of the muon longitudinal momentum (in GeV/c). The curves represent the prediction of the SuSAv2+2p2h-MEC (blue) as well as the separate quasielastic (red) and 2p2h-MEC (orange) channels. MINERvA data and experimental fluxes are from Refs. [61] and [62].
The SuSAv2-MEC model is also compared in Fig. 10 with the MINERvA CCQE-like double differential (anti)neutrino measurements on hydrocarbon (CH) in terms of the transverse momentum of the outgoing muon (with respect to the antineutrino beam), in bins of the muon longitudinal momentum. Due to its relativistic nature, the SuSAv2-MEC model is well suited to describe these data [34,61] where the mean neutrino energy is around 3.5 GeV. An overall good agreement is reached without resorting to any tuning or additional parameters. According to Refs. [34,61], the MINERvA "QE-like" cross sections entail, besides pure quasielastic contributions, events that have post-FSI final states without mesons, prompt photons above nuclear de-excitation energies, heavy baryons, or protons above a kinetic energy threshold of 120 MeV, thus including zero-meson final states arising from resonant pion production followed by pion absorption in the nucleus and from multi-nucleon interactions. Apart from the "QE-like" points, the MINERvA "CCQE" signal is also shown, and corresponds to events initially generated in the GENIE event generator [63] as quasielastic (that is, no resonant or deep inelastic scatters, but including scatters from nucleons in correlated pairs with zero-meson final states), regardless of the final-state particles produced, thus including CCQE and 2p2h reactions. The difference between the two data sets is mainly related to pion production plus re-absorption and goes from ∼ 15% to ∼ 5% depending on the kinematics. The present SuSAv2-MEC results do not include processes corresponding to these pion re-absorption processes inside the nucleus and the comparison should be done with the "CCQE" data rather than with the "QE-like" ones. A more detailed analysis of these results together with a χ 2 test can be found in [34], where the SuSAv2 χ 2 shows its compatibility with data and with the MINERvA/GENIE predictions. It is also important to note that due to MINERvA's acceptance, the muon scattering angle is limited to θµ < 20 • as well as the muon kinematics (1.5 GeV < p || < 15 GeV, pT < 1.5 GeV) in both experimental and theoretical results. This implies important phase-space restrictions for large energy and momentum transfer to the nuclear target and makes the available phase space not so different from the T2K one, as shown in [64].
The previous T2K and MINERvA results are only related to the final lepton kinematics, the so-called inclusive measurements. At this point, it is worth mentioning that other models [38,65,66,67,68,69,70,71,39,72,73,74,75,76,77,78,79,80,81,82,83,84,85] have been also developed to address these CC inclusive neutrino interactions and, although similar agreement with data can be obtained, they are based on different assumptions about the nuclear properties and dynamics. Combined analyses of these models with more exclusive neutrino measurements where hadron kinematics and other nuclear effects can be analyzed in more detail would help to improve model selection for data analysis (see also the discussions in Sect. 3). The advantage of SuSAv2-MEC (and RMF) is that of being a fully relativistic model that has shown an overall good agreement with electron and neutrino scattering data and that can be extended without further approximations to the full-energy range of interest for present and future neutrino experiments. For this reason, the SuSAv2-MEC model (1p1h and 2p2h) has been recently implemented in the GENIE neutrino event generator [63] with the aim of improving the characterisation of the nuclear effects in neutrino cross section measurements and work is now in progress to implement this model together with the ED-RMF one in the NEUT event generator [86] for its application on the T2K oscillation analysis.
Based on these works, in Fig. 11 we study the T2K CC0π measurement of interactions with protons less than 500 MeV [87] in comparison with the SuSAv2-1p1h, SuSAv2-2p2h and pion-absorption predictions from GENIE (see [45] for details). At proton momentum below 500 MeV/c a clear dominance of the SuSAv2-1p1h channel is observed. Note that SuSAv2-GENIE's also shows an overestimation of data at very forward angles, as observed in Fig. 8 which is mainly due to low energy transfer scaling violations which are absent in the SuSAv2-model but present in the RMF theory. In order to analyze the low-energy nuclear effects and these scaling violations, we also show the comparison with the RMF model for the 1p1h channel together with the SuSAv2-2p2h and pionabsorption results from GENIE, noticing a remarkable improvement in the data comparison at very forward angles. Here, one of the drawbacks of the SuSAv2 model is that it only predicts lepton kinematics so its implementation in MC event generators has to rely on a factorization approach and on the information available in these generators to determine hadron kinematics (see [45] for details). This limitation will be addressed in the forthcoming implementation of the ED-RMF model in event generators which will provide full information about lepton and hadron kinematics in the final state together with a consistent description of the nuclear dynamics.
Due to aforementioned dominance of the 1p1h channel on CC0π interactions with low momentum protons observed in Fig. 11, it remains difficult to draw clear conclusions about the goodness of the 2p2h description. Nevertheless, this can be explored further using more semi-inclusive measurements with measurements of proton and muon kinematics as shown in Fig. 12, where the SuSAv2-MEC model in GENIE is compared with the so-called single transverse variables (STV), and in particular with the transverse momentum imbalance, δpT , defined in terms of the momentum imbalances between the outgoing muon and highest momentum proton in the plane transverse to the incoming neutrino (see [45,88] for details). These transverse kinematic imbalances allow one to study initial-  Figure 11: Comparisons of data and model predictions for differential CC0π muon-neutrino cross sections on 12 C in the T2K neutrino beam as a function of the muon kinematics when there are no protons with momenta above 500 MeV. Two 1p1h predictions are shown (one from RMF, the other from SuSAv2 implemented in GENIE), in addition to the SuSAv2 2p2h and pion absorption contributions from GENIE. The total contributions when using each of the two 1p1h models is also shown. Goodness of fit are calculated to be χ 2 RMF = 171.87 (59 bins) and χ 2 SuSA = 168.92 (60 bins), where the latter includes a single extra bin from -1.0 to -0.3 cos θ (not shown). The data points are taken from [87].  Figure 12: The regularised T2K measurement of CC0π muon-neutrino cross sections on 12 C at T2K kinematics as a function of the Single Transverse Variable [88] δp T compared to predictions from the GENIE-implemented SuSAv2 and Valencia 1p1h+2p2h models, each of which is added to GENIE's pion absorption prediction. The total contributions when using SuSAv2 and Valencia models is also displayed. A shape only comparison is also shown (right panel). Goodness of fit are calculated as follows. For δp T : χ 2 SuSA = 20.5, χ 2 V alencia = 27.1. The data points are taken from [87].
state nuclear dynamics but also to better isolate the 2p2h channel. In δpT the 1p1h channel is not expected to contribute significantly beyond the initial state nucleon momentum (∼230 MeV/c for carbon in a Fermi gas approach), thus implying that the high-δpT values will be dominated by 2p2h and other contributions, as observed in Fig. 12. The overestimation at high-δpT in the left panel may indicate that the 2p2h contribution is too strong. However, as discussed in [89], this overall over-prediction could potentially be improved by stronger nucleon FSI, which may improve the data agreement in the tail. On the other hand, this may also be explained by the approximations taken to produce semi-inclusive predictions from inclusive models, as described in [45]. It is expected that these drawbacks can be addressed in further works via the implementation of RMF models for the 1p1h channel and a full semi-inclusive 2p2h model. Apart from these limitations, it is worth noticing the almost perfect description of the shape of δpT in the right panel, improving the agreement reached by other descriptions also implemented in generators.
To conclude this section, we show the comparison of the SuSAv2-MEC model with measurements that also consider non-QE contributions, mainly pion production. This reaction channel is also of relevance for present and future neutrino oscillation experiments. Neutral current π 0 production constitutes an important background in the electron neutrino and antineutrino appearance analyses but also neutrino-induced pions emitted can mimic QE-like events if they are not properly detected or if they are absorbed due to FSI effects. Thus, the analysis and detection of these pion production events in coincidence with the final lepton and other hadrons is of paramount importance for the neutrino energy reconstruction.
In order to analyze the neutrino-induced pion production channel we make use of the phenomenological SuSA-∆ approach which extends the superscaling arguments observed in the QE regime to the ∆ resonance region (see [36] for details). In Fig. 13 (left panel) we show the νe-12 C inclusive differential cross section averaged with the T2K flux versus the reconstructed four-momentum transfer, Q 2 QE , and in Fig. 13 (right panel) the νµ total cross section on C8H8 target. Three contributions are shown in these plots, namely, the SuSAv2 QE and SuSAv2 2p2h-MEC channels together with the SuSA-∆ approach. Although an overall good agreement with data is observed in both panels, some underestimations are present at large kinematics, i.e., high Eν and Q 2 QE , which reveal the need for including higher resonances and deep inelastic scattering in the description. To address this, the SuSAv2-inelastic, which has been successfully applied for the analysis of the full inelastic spectrum in electron scattering (see Sect. 2.4 for details), will be soon extended to the analysis of the neutrino sector.  Figure 13: (Left panel) We show the CC-inclusive T2K flux-folded ν e -12 C Q 2 QE differential cross section per nucleon. (Right panel) The CC ν µ total cross section on C 8 H 8 is presented. Experimental data are from T2K [90] and SciBooNE [91]. Theoretical predictions for QE, non-QE (1π) and the 2p2h MEC are shown separately. Plots from Ref. [6].

Conclusions
In this paper we have summarized the basic ingredients that go into scaling analyses of inclusive electron scattering and charge-changing neutrino reactions with nuclei. Importantly, good agreement in the former case is viewed as a pre-requisite to being able to predict the latter -clearly if one fails to account for inclusive electron scattering then it is unreasonable to expect that inclusive CCν reactions will be adequately modeled.
In particular, we have emphasized the use of scaling of the first kind (independence of q at high energies) and of scaling of the second kind (independence of nuclear species). When both kinds of scaling are invoked we refer to it as Super-Scaling and accordingly our present focus has been placed on the so-called Super-Scaling Approach (SuSA) in its original form and in extensions of that original form. Such extensions have been introduced to account for the modest level of scaling violations seen in direct comparisons with inclusive electron scattering data.
In addition we have introduced the way 2p-2h MEC effects have been incorporated, showing that, while they are usually corrections to the cross sections being represented via the scaling functions obtained either phenomenologically or through use of specific models, they typically are required to get a successful picture of the response.
We have briefly summarized these basic ideas and then proceeded to show examples of the excellent agreement found with inclusive electron scattering cross sections measured for several light to mediumweight nuclei.
To place the discussions in context we have also provided a section outlining how semi-inclusive scattering and inclusive scattering are related, and specifically how the former probes particular regions of missing-energy and -momentum (characterized by so-called trajectories), while the latter involves integrations over the Em − pm plane of the semi-inclusive response. On the one hand, inclusive scattering, being a total hadronic cross section, is less dependent on the underlying details of nuclear structure. In fact, given that the kinematics of the single-nucleon knockout are handled relativistically and that there are sum rules that determine the integral of the inclusive response, then a single parameter can be used to fix the width of the response. For this reason many models yield the rough behaviour of the inclusive cross section, even the relativisitic Fermi gas model. On the other hand, the semi-inclusive cross section depends critically on how the strength is distributed as a function of Em and pm. Accordingly, models that fare reasonably well for inclusive scattering cannot be relied upon to properly represent the underlying nuclear structure needed for semi-inclusive reactions.
Finally, having discussed the foundations of Super-Scaling for electron scattering we show some selected results for CCν reactions. The agreements with existing data are found to be excellent, giving us confidence that at least for inclusive neutrino reactions the problem appears to be well in hand. Note, however, from the statements made in the previous paragraph, that this is not to be taken as proof that semi-inclusive CCν modeling should be expected to be robust. Indeed, different models yield quite different results for cross sections in which both a charged lepton and a nucleon are detected in coincidence (a trend in modern experiments), even when the corresponding inclusive cross sections do not differ significantly. The modeling of semi-inclusive reactions is part of our present projects.