neutrinos are produced in states of given flavors  electron , muon or tau  that are mixtures of mass eigenstates .
as different masses propagate with different phases , the mixture at some distance from a neutrino source may be different from the initial one , resulting in a disappearance of the initial flavor and an appearance of other flavors in the neutrino beam .
quantitatively , the phase of the state of mass @xmath0 and momentum @xmath1 at the time @xmath2 from the production is @xmath3 at the distance @xmath4 from the source and , therefore , the disappearance and appearance probabilities depend on the squared mass differences , @xmath5 , ratio @xmath6 , @xmath7 being the neutrino energy , and mixing angles , @xmath8 , describing how mass eigenstates mix to form states of definite flavor .
this phenomenon , called the neutrino oscillations , has been well established in the past two decades and is now considered one of the major discoveries in particle physics , celebrated by two recent nobel prizes  @xcite .
the values of the oscillation parameters  two squared mass differences and three mixing angles  are being measured with an increasing precision , thanks to a global effort that produced a number of experiments currently taking data or planned in the near future .    within the next two decades
, the precision is expected to be sufficient to begin testing the unitarity of the mixing matrix  @xcite , discover the neutrino - mass hierarchy  @xcite and unambiguously determine the value of the dirac phase , @xmath9 , violating the charge - parity ( cp ) symmetry of neutrino mixing  @xcite .
these results are going to have profound consequences for possible extensions of the standard model , building models of the neutrino masses and our understanding of the matter - antimatter asymmetry in the universe .    in this topical review ,
we discuss systematic uncertainties in ongoing and future long - baseline ( @xmath103001300 km ) neutrino - oscillation experiments using conventional beams  @xcite .
such neutrino beams are tertiary products , originating from the decay of mesons
predominantly pions  produced in interactions of the primary proton beam with a target . as the resulting neutrinos are not monoenergetic , to extract the oscillation parameters from the collected event distribution , their energies have to be reconstructed on an event - by - event basis from the measured kinematics of particles in the final state . depending on the beamline
, the relevant neutrino energies extend from a few hundred mev to a few gev , at which the dominant interaction mechanisms change from quasielastic to resonant and nonresonant meson production  @xcite .
because the oscillation parameters are extracted from the energy dependence of event distributions , their accurate determination requires an accurate reconstruction of neutrino energy .
while the reconstruction would involve smallest uncertainties for scattering off free protons and deuteron  @xcite , owing to the low cross sections involved , it is necessary to employ nuclear targets as detector materials to ensure high statistics of collected events . as
a consequence , however , description of nuclear effects turns out to be one of the largest sources of systematic uncertainties in the oscillation analysis of modern long - baseline experiments .
additionally , when energy reconstruction requires the kinematics of produced hadrons to be measured , detector effects may play an important role and uncertainties of the detector response contribute to the systematic uncertainties of the extracted oscillation parameters  @xcite .    over the past decade
, the extensive neutrino - scattering program has yielded a wealth of experimental cross sections for carbon or hydrocarbon , ch  @xcite , and much fewer results for other targets , such as water@xcite , argon @xcite , iron @xcite and lead @xcite .
it is important to note that there seem to remain puzzling tensions between different measurements  such as @xcite or @xcite  which have attracted a sizable attention of the theoretical community .
a number of approaches developed to describe nuclear response to other probes has been extended to neutrino interactions , and the understanding of nuclear effects relevant to neutrino - oscillation experiments has clearly improved  @xcite ; see also recent reviews  @xcite .
however , many problems still await a quantitative explanation . while for detector targets containing hydrogen the pion - production events on free protons can be separated out  @xcite , current and future neutrino - oscillation experiments are going to collect data predominantly for interactions with the targets of atomic numbers @xmath11 ranging from 12 to 40 , for which accurate and complete nuclear models permitting fully trustable data analysis are not available yet .
the main difficulty of developing theoretical models useful to oscillation experiments stems from the flux average over polychromatic beams  @xcite . as any data bin may receive contributions from a range of neutrino energies and different interaction mechanisms , understanding the source of discrepancies between theoretical results and experimental data
is a formidable task .
for example , resonant pion production followed by pion absorption in medium leads to the same final state as single - nucleon knockout . at the kinematics of the long - baseline oscillation experiments
, nuclear models must be able to describe relativistic products of interaction .
additionally , experiments employing tracking detectors require predictions of exclusive cross sections for different hadronic final states .
this topical review is organized as follows .
after general remarks regarding systematic uncertainties of measurements of cp violation in , we present our procedure of the oscillation analysis in . in sections
[ sec : nucleareffects ] and [ sec : detectoreffects ] we illustrate the importance of nuclear and detector effects for an accurate reconstruction of neutrino energy and an unbiased extraction of oscillation parameters , discussing selected results from our previous publications  @xcite .
finally in we summarize this review .
in recent and ongoing oscillation experiments , it has been a common practice to use the spread between different theoretical descriptions of nuclear effects , implementations of the same model in monte carlo generators , detector response parametrizations etc . as an estimate of the associated systematic uncertainties .
while this method seems very effective , a word of caution is in order .
non - negligible correlations between the ingredients used to find the uncertainties  such as common assumptions of nuclear models or fine - tuning of event generators to the same data  may lead to sizable underestimates , making this procedure insufficient and unreliable .    for the next 20 years
, the long - baseline neutrino - oscillation program is most likely going to employ conventional beams  @xcite .
this technique of producing neutrino beams has been used for a few decades : an intense proton beam is impinged on a target typically made of beryllium or graphite in order to produce mesons , mainly pions .
the target is embedded within a horn  an electromagnet producing a toroidal magnetic field  that focuses secondary particles of a  selected charge and defocuses those of the opposite charge .
the horn polarization determines whether the resulting beam will be made of neutrinos or antineutrinos .
it is important to point out that the beams and their contaminations by wrong - sign mesons in the neutrino and antineutrino modes are not related in a meaningful way .
therefore , while using the same beamline , target and horn , the neutrino or antineutrino runs of an experiment could be considered two separate experiments .    an accurate determination of the neutrino flux  its flavor composition and energy spectrum  exclusively from the proton - beam parameters , meson - production data for thick targets and horn configuration is a challenging task .
state - of - the - art methods of controlling systematic uncertainties of a neutrino beam have been developed for the minos  @xcite and minerva  @xcite experiments , largely thanks to using two horns of adjustable positions at the numi beamline  @xcite .
it is not clear whether such a feature  difficult to realize for off - axis beams  is going to be available for future oscillation experiments .
therefore , the beam normalization uncertainty of about 5% achieved in the minos and minerva experiment , can be expected also in future oscillation experiments employing conventional beams .
note , however , that recently proposed new concepts of neutrino beams , such as @xmath12storm  @xcite , would allow the neutrino flux to be determined with an accuracy of the order of 1% .    to extract the probability of oscillation between the neutrino flavors @xmath13 and @xmath14 , experiments collect event distributions @xmath15 with respect to a set of observables @xmath16 , such as the charged lepton s energy and cosine of its production angle , @xmath17 in the above equation , the normalization factor @xmath18 depends on the beam power , data collecting time , fiducial mass etc . the flux expected in the detector , @xmath19 , and the oscillation probability @xmath20 are both functions of the true neutrino energy @xmath7 .
the differential cross section @xmath21 describes likelihood for a neutrino of the flavor @xmath14 and energy @xmath7 to produce an event of kinematics @xmath16 .
the detection efficiency is denoted as @xmath22 .
systematic uncertainties in neutrino - oscillation experiments can be effectively reduced by using an unoscillated event distribution measured in the near detector system to predict the distribution expected in the far detector .
this method of exploiting the cancelation of correlated uncertainties  giving the best results when the near and far detectors are functionally identical  has been used with great success in the reactor experiments daya bay  @xcite , reno  @xcite and double chooz  @xcite to measure the mixing angle @xmath23 with remarkably high precision . barring the differences in the backgrounds and in the geometric acceptances of the beam
, the ratio of the energy - unfolded event distributions in a disappearance experiment , @xmath24 , has in this case the simple interpretation @xmath25 where the far - to - near flux ratio , @xmath26 , reduces to the ratio of the squared distances from the neutrino source , @xmath27 and the oscillation probability in the near detector is @xmath28 .    unlike in the reactor experiments , in long - baseline studies ,
the relevant cross sections are currently known with an accuracy of 1020% and the procedure of energy unfolding in much more involved . additionally , in the near and far sites , different detectors are usually employed and
the differences in backgrounds and beam acceptances are non - negligible . for a disappearance measurement
, minos can serve as a benchmark of how well the near - to - far comparison does reduce systematic uncertainties  @xcite .    in @xmath29 appearance measurements ,
the situation is even more involved .
the final and initial flavors of neutrino are different , which leads to a dependence of the near - to - far event distributions ratio on the cross sections ratio @xmath30 . as the beams are designed to minimize the @xmath31 contamination in order to reduce the background in the far detector , direct measurements of the @xmath31 cross sections suffer from the statistics lower typically by 23 orders of magnitude than in the @xmath32 case and from larger uncertainties coming from both the flux and detector response  @xcite .
for example , in the recent t2k analysis  @xcite , the total systematic uncertainty is estimated to be 16.6% , with the flux and detector - response contributions of 12.9 and 8.4% , respectively .
there are no obvious methods to reduce those uncertainties and measurements for electron antineutrinos are going to be even more challenging .
however , it is important to note that all available data are consistent with lepton universality .    to emphasize the importance of reduction of systematic uncertainties in appearance measurements ,
let us make a simple estimate of those needed for an experiment to probe cp violation at 3@xmath33 confidence level for 75% of the possible @xmath9 values , the goal set by the p5 advisory panel  @xcite .
we consider the cp asymmetry defined as  @xcite @xmath34 the energy average @xmath35 is taken over the full width at half the first oscillation maximum . in vacuum @xmath36 is proportional to @xmath37 and , therefore , for @xmath9 values in the vicinity of 0 and @xmath38 , the uncertainties of @xmath36 are very similar to those of @xmath9 .
baseline dependence of the cp asymmetry   for different @xmath9 values.,scaledwidth=80.0% ]    shows the cp asymmetry @xmath36 as a function of the baseline for the @xmath9 values of @xmath39 , 0 , and @xmath40 and the remaining oscillation parameters from a recent global fit  @xcite . for the 1500-km baseline ,
matter effects introduce the asymmetry of @xmath1047% and @xmath9 can modify it by @xmath41 . for 75% of the @xmath9 values ,
the cp asymmetry can become as small as 5% , which at 3@xmath33
confidence level translates into @xmath101.5% accuracy . for equal contributions of statistical and systematic uncertainties ,
the latter are required not to exceed @xmath101% .
fortunately , if the maximal cp violation is realized in nature , as suggested by recent results of the t2k  @xcite and nova  @xcite experiments , less stringent conditions are going to be sufficient for the @xmath9 discovery .
discussing the role of nuclear and detector effects in the oscillation analysis , in sections [ sec : nucleareffects ] and [ sec : detectoreffects ] we present the results obtained using the software package globes  @xcite . simulating data according to  , we keep all the oscillation parameters fixed , setting their values to @xmath42    unless specified otherwise , we perform the @xmath43 disappearance analysis for a setup similar to that of the t2k experiment  @xcite with a narrow - band off - axis beam  @xcite , considering carbon as the detector material .
the main motivation for selecting the @xmath44 channel is that it can effectively be treated as two - flavor oscillations with the probability  @xcite @xmath45 where @xmath46 which allows for a clear interpretation of the obtained results .
as detailed in , the experimental setup consists of the near and far detectors simulated following  @xcite .
the exposure and beam intensity are chosen to ensure that the statistics of unoscillated events is high enough to constrain the analyzed systematic uncertainty .
our treatment of systematic uncertainties , the @xmath47 implementation and the method of determination of the confidence regions is the same as in @xcite .
we assume 20% systematic uncertainties for the shape ( bin - to - bin uncorrelated ) and for the overall normalization ( bin - to - bin correlated ) . adding a gaussian prior for each systematic error , the pull method is used and the final @xmath47 profile is obtained after the minimization over the nuisance parameters .
@lrr & baseline & fiducial mass + far detector & 295.0  km & 22.5  kt + near detector & 1.0  km & 1.0  kt +    as neutral - current background is expected not to play an important role in the disappearance analysis , we neglect it  unless specified otherwise  and consider the following mechanisms of charged - current ( cc ) neutrino interaction : quasielastic ( qe ) involving a single nucleon or more nucleons ( @xmath48 ) in the final state , resonance excitation , and deep - inelastic scattering . in all cases , meson - production events are obtained using the monte carlo generator genie 2.8.0  @xcite .
qe events are simulated using either the relativistic fermi gas ( rfg ) model from genie or the realistic spectral function ( sf ) approach  @xcite implemented in the @xmath49 package of additional modules  @xcite . for the @xmath48 contribution
, we employ the empirical procedure developed by dytman  @xcite or the effective approach of increasing the value of the axial mass @xmath50 in the qe cross sections .    performing the kinematic energy reconstruction
, we use the well - known formula  @xcite @xmath51 setting the invariant hadronic mass @xmath52 to the nucleon mass @xmath53 for mesonless events and fixing the separation energy @xmath54 to 34 mev .
the same value of @xmath54 is added for every detected nucleon in the calorimetric method , @xmath55 in which the neutrino energy is reconstructed summing the charged - lepton s energy @xmath56 , the kinetic energies of the knocked - out nucleons @xmath57 and the total energy of any other particle produced @xmath58 .    in an idealized scenario in which nuclear and detector effects would not play a significant role ,
the energy reconstruction would be perfect . in a realistic case
, however , there is a non - negligible probability of energy misreconstruction @xcite . in our calculations ,
both nuclear and detector effects are encoded in the migration matrices , @xmath59 , the elements of which describe the probability that an event of type @xmath60 with a true energy in the bin @xmath61 is reconstructed with an energy in the bin @xmath62 .
the simulated ( true and fitted ) event distributions are then obtained as @xmath63 where @xmath60 runs over the four types of interactions considered , @xmath62 and @xmath61 label the energy bins and @xmath64 stands for the true number of @xmath60 events in the bin @xmath61 .
as the oscillation parameters are extracted from the energy distribution of collected events , a precise reconstruction of neutrino energy is a prerequisite for a precise oscillation analysis .
it is important to note that the kinematic reconstruction method , valid for particular final states , relies on a selection of event subsample and requires an accurate simulation of both the signal events and the background contribution  @xcite . analyzing various aspects of qe interactions with one or more nucleons in the final state  @xcite , in this section
we discuss how uncertainties in the description of nuclear effects may affect the outcome of @xmath32 disappearance experiments employing the kinematic method of energy reconstruction .
while modern neutrino - scattering measurements have been predominantly performed for the carbon target , the future oscillation experiments are going to require detailed knowledge of the cross sections for oxygen and argon . to demonstrate the subtleties of extrapolating the results between nuclei
, we are going to neglect differences between the migration matrices for carbon ( @xmath65 ) and oxygen ( @xmath66 ) in the @xmath32 disappearance analysis  @xcite .
we consider a t2k - like experimental setup with water cherenkov detectors , and model nuclear effects relying exclusively on genie .
in particular , in this analysis qe scattering is described within the rfg model of bodek and ritchie  @xcite , with a high - momentum tail  inspired by the effect of nucleon - nucleon correlations  added to the nucleon momentum distribution given by the step function .
the globes analysis for 5 years of collecting data with the 750 kw beam  @xcite yields @xmath10880 cc qe events , @xmath10480 pionless events resulting from other mechanisms of interaction and 275 neutral - current background events .
impact of the nuclear target dependence of the cross sections in genie on the oscillation analysis . left : event distributions for oxygen ( solid line ) and carbon ( dashed line ) expected for the flux  @xcite .
right : 1 , 2 and 3@xmath33 confidence regions in the @xmath67 plane obtained when the simulated data for oxygen are fitted using the migration matrices for carbon ( solid lines ) and oxygen ( shaded areas ) .
the triangle and circle show the best fit point and true values of the oscillation parameters , respectively.,title="fig:",scaledwidth=39.0% ] impact of the nuclear target dependence of the cross sections in genie on the oscillation analysis . left : event distributions for oxygen ( solid line ) and carbon ( dashed line ) expected for the flux  @xcite .
right : 1 , 2 and 3@xmath33 confidence regions in the @xmath67 plane obtained when the simulated data for oxygen are fitted using the migration matrices for carbon ( solid lines ) and oxygen ( shaded areas ) .
the triangle and circle show the best fit point and true values of the oscillation parameters , respectively.,title="fig:",scaledwidth=39.0% ]    the @xmath68 distributions of events obtained using the oxygen and carbon migration matrices are shown in the left panel of .
the carbon result turns out to be somewhat shifted toward lower energies with respect to the oxygen one , which translates into a noticeable effect for the oscillation analysis .
the right panel of shows the 1 , 2 and 3@xmath33 confidence regions in the @xmath69 plane for the true event rates simulated using the oxygen migration matrices .
the shaded areas represent the analysis in which nuclear effects are accurately taken into account by using the oxygen migration matrices to calculate the fitted rates . on the other hand
, the solid lines correspond to the fitted rates for the carbon migration matrices and show the effect of an inaccurate description of nuclear effects .
the extracted and true values of the oscillation parameters turn out to differ at a @xmath101@xmath33 confidence level .
we observe a non - negligible difference between the energy reconstruction for the oxygen and carbon targets described by means of the rfg model in genie .
should the differences in their shell structures and density distributions be accounted for within a more sophisticated model , the effect can be expected to be even larger .
note also that carbon and oxygen are considered rather similar nuclei compared to argon .
the result of shows the importance of developing an accurate theoretical description of nuclear effects for the targets employed in oscillation experiments and of their extensive testing against scattering data .
it also suggests that employing the same target in the near and far detectors can be an effective way to reduce systematic uncertainties of an oscillation experiment .
as neutrino scattering data are collected using polychromatic beams , every data bin receives contributions from a range of neutrino energies .
moreover , the same detected final state may be produced by different interaction mechanisms .
therefore , an accurate description of experimental cross sections is a challenging task , requiring a theoretical model to cover broad kinematic region and describe multiple scattering processes  @xcite . with different probe s energies and reaction mechanisms intertwined
, it may be very difficult to identify the origin of discrepancies between theoretical and experimental results for neutrinos .    on the other hand , when the interaction process involves only one nucleon
with the remaining @xmath70 ones acting as spectators  neutrino and electron scattering are subject to the same nuclear effects .
the nuclear cross section in this regime is a convolution of the elementary cross section ( different for electrons and neutrinos ) with the hole and particle spectral functions ( common for electrons and neutrinos ) , describing the ground - state properties of the target and propagation of the struck nucleon , respectively  @xcite . as a consequence ,
models of _ nuclear effects _ in neutrino interactions can be validated by systematic comparisons to electron scattering data at the kinematics of interest .
we want to emphasize that it is highly improbable that a theoretical approach unable to reproduce electron scattering data would be able to describe nuclear effects in neutrino interactions .
electron scattering cross sections
double differential and for monoenergetic beams  allow for much easier understanding of discrepancies between theoretical results and experimental data than neutrino ones .
importance of an accurate description of nuclear effects on the example of the double differential c@xmath71 cross sections . left : comparison of the rfg ( dashed lines ) and sf  @xcite ( solid lines ) calculations with the data for scattering angle @xmath72 and beam energies 1.35 , 1.50 and 1.65 gev  @xcite .
right : comparison of the sf results obtained using genie  @xcite and genie@xmath73 t  @xcite with the experimental data for scattering angle @xmath74 and beam energy 961 mev @xcite.,title="fig : " ] importance of an accurate description of nuclear effects on the example of the double differential c@xmath71 cross sections . left : comparison of the rfg ( dashed lines ) and sf  @xcite ( solid lines ) calculations with the data for scattering angle @xmath72 and beam energies 1.35 , 1.50 and 1.65 gev  @xcite .
right : comparison of the sf results obtained using genie  @xcite and genie@xmath73 t  @xcite with the experimental data for scattering angle @xmath74 and beam energy 961 mev @xcite.,title="fig : " ]    as an illustrative example , in the left panel of we compare the qe predictions of the sf calculations  @xcite and the rfg model of smith and moniz  @xcite with the data  @xcite .
whereas the rfg results can not reproduce heights , shapes and positions of the qe peaks , the sf calculations turn out to be in very good agreement with the experimental points , although they do not involve any adjustable parameters .    using the realistic ground - state description  @xcite ,
the sf calculations take into account both the shell structure and correlations between nucleons .
the latter cause partial depletion of the shells  lowering the qe peaks  and give rise to deeply bound nucleons of high momenta , producing the tails at low values of the final electron energy @xmath75 and rendering the qe peaks asymmetric .
in addition , the sf calculations  @xcite include the effect of final - state interactions , essential to reproduce the positions of the qe peaks accurately .
all these effects are neglected in the rfg model , which treat the nucleus as a fragment of noninteracting nuclear matter of uniform density in a constant potential .
the sf approach is available in genie 2.8.0 , however , comparisons with electron scattering data have revealed some issues with the energy and momentum conservations in its implementation that affect both the position and width of the qe peak , as shown in the right panel of for the experimental cross sections  @xcite . on the other hand ,
the sf implementation in the @xmath49 package  @xcite of additional modules to genie is in good agreement with electron scattering data and , therefore , from now on we are going to rely on it to generate qe events .    to explore how the difference between the rfg model and the sf approach in the qe interaction channel influences the oscillation analysis  @xcite ,
let us consider a @xmath32 disappearance experiment similar to t2k , employing carbon as the target material .
assuming 5 years of data taking at the beam power of 750 kw in the globes analysis leads to @xmath10650 ( @xmath10730 ) cc qe events in the sf approach ( rfg model ) , together with @xmath10410 cc events of qe topology coming from other interaction mechanisms and @xmath10250 neutral - current background events .
the difference between the cc qe event numbers in the two considered approaches , resulting from the difference between the corresponding cross sections presented in the left panel of , contributes to a non - negligible effect for the oscillation analysis .
impact of differences between the quasielastic @xmath76 cross sections in the sf approach and the rfg model on the oscillation analysis . left : mesonless cross sections employed in the analysis .
right : confidence regions in the @xmath67 plane obtained when the data simulated within the sf approach are fitted using the migration matrices from the rfg ( solid lines ) and sf ( shaded areas ) calculations . ]
the right panel of shows the 1 , 2 and 3@xmath33 confidence regions in the @xmath69 plane for the true event rates from the sf approach .
the shaded areas are obtained for the fitted rates from the same approach , while the solid lines correspond to those from the rfg model .
we observe a shift between the true values of the oscillation parameters  represented by the dot  and the extracted best - fit point  marked by the triangle  at a @xmath101@xmath33 confidence level , both in the mixing angle and the squared mass difference .
as we analyze the difference between the rfg model and the sf approach only for the cc qe predictions , the presented result should be regarded as a lower limit of the influence of nuclear - model uncertainties . for resonance excitation
, similar differences can be expected based on the available theoretical results  @xcite . in the deep - inelastic region ,
the current modeling of nuclear effects has been shown inaccurate by the recent results from the minerva experiment  @xcite .    from now on , we rely on the sf approach to simulate cc qe events and employ a generalization of the standard formula   for the kinematic energy reconstruction  @xcite : in the case of mesonless events , we employ the invariant hadronic mass @xmath77 , regardless of the number of knocked - out nucleons , but when at least one meson is detected in the final state , we set @xmath52 to the @xmath78 resonance mass @xmath79 gev .
to discuss the effect of uncertainties of the @xmath48 cross sections on the oscillation analysis , we consider a t2k - like disappearance experiment running in the neutrino and antineutrino mode with the same flux  @xcite and compare the results obtained using two data - driven phenomenological methods to account for @xmath48 processes  @xcite . in the genie@xmath73 t approach ,
the @xmath48 estimate from genie 2.8.0  @xcite is added to the qe contributions obtained using the @xmath49 package  @xcite . in the effective approach ,
the modifications of the qe cross sections due to @xmath48 mechanisms are accounted for by applying in the sf calculations an effective value of the axial mass @xmath80 gev , as suggested by a number of experimental results for nuclear targets ranging from carbon to iron  @xcite .
the rationale for considering these two methods is that the spread of the cross sections they predict is a good representation of the spread between the available experimental data .
shows that the effective calculations are in good agreement with the cc qe results from nomad  @xcite and minerva  @xcite , for both neutrinos and antineutrinos .
they also reproduce the energy dependence of the neutrino and antineutrino data from miniboone  @xcite , but not their absolute normalization . to better illustrate this feature
, we also present the miniboone cross sections divided by a factor of 1.2 , consistent with the ratio of the detected to predicted events @xmath81 in the first miniboone analysis  @xcite .
the genie@xmath73 t approach is in very good agreement with the miniboone data for neutrinos , which were used to determine the strength of the @xmath48 contribution in genie  @xcite , but this is not the case for antineutrinos . owing to their large uncertainties
, the t2k data  @xcite can not discriminate between the two calculations .
comparison of the total cc inclusive @xmath32 cross sections in the left panel of shows that the two considered methods yield the results in good agreement with the nomad data  @xcite , collected in the region dominated by deep inelastic scattering .
while the flux - averaged on - axis measurement from t2k  @xcite does not distinguish the two approaches , the off - axis one  @xcite shows a distinct preference for the effective calculations . on the other hand ,
the sciboone point  @xcite clearly favors the genie@xmath73 t calculations .
note that the t2k  @xcite and sciboone  @xcite measurements do not involve energy unfolding or selection of particular final states and , therefore , they can be expected to involve the lowest uncertainties .
the puzzling difference between the t2k and sciboone data has important consequences for @xmath32 disappearance studies , as illustrated in the right panel of by the confidence regions obtained for the true event rates from the genie@xmath73 t calculations and the expected number of @xmath106000 unoscillated events .
the shaded areas represent the results for the genie@xmath73 t fitted rates and the solid lines correspond to the fitted rates from the effective calculations .
the high value of @xmath47 per degree of freedom in the best - fit point , 66.0/26 , clearly indicates that the differences between the two considered approaches are too large to be neglected in a precise oscillation analysis .
same as figure [ fig:2p2h_nu ] but for muon antineutrinos.,title="fig:",scaledwidth=53.0% ] same as figure [ fig:2p2h_nu ] but for muon antineutrinos.,title="fig:",scaledwidth=39.0% ]    for completeness , in we present the results for muon antineutrinos .
owing to the large differences between the @xmath82 cc qe cross sections , the two considered approaches yield very different predictions for the inclusive cc @xmath82 cross section , as shown in the left panel of . however , there are no experimental data available to discriminate between them .
the observed differences translate into a large effect for the oscillation analysis . using the fitted rates from the effective calculations and the true event rates from genie@xmath73 t
leads to a severe discrepancy between the extracted and true values of the oscillation parameters , see the right panel of . because treating the normalization of the qe event sample  with any number of nucleons  as arbitrary yields very similar results , the observed behavior can be traced back to the different shapes the two approaches predict for the migration matrices and for the energy dependence of the cc qe cross section .    in final remarks of this section
we note that the significance of nuclear effects , analyzed here for the carbon nucleus ( @xmath65 ) , can be expected to increase further for the argon target @xmath83 .
although the considered experimental setup and our procedure of data analysis do not closely follow those of any existing experiment , the discussed results point toward the crucial role that an accurate description of nuclear effects plays in monte carlo simulations of precise oscillation experiments .
for the success of precise oscillation studies it is , therefore , of great importance to continue efforts aiming to provide new cross - section measurements , improve existing monte carlo generators , develop more accurate nuclear models and determine their uncertainties in comparisons with available data .
in the oscillation analysis of long - baseline experiments , neutrino energy is inferred from the kinematics of particles produced in the interaction . as a consequence , the reconstructed energy may be affected by finite detection capabilities  energy resolutions , efficiencies and thresholds  and uncertainties they are known with .
we explore this issue within a realistic scenario , in which particles are detected according to their efficiencies and thresholds , and measured energies and angles are smeared by finite detector resolutions . to analyze how uncertainties of these detector effects may affect the oscillation analysis , we obtain simulated event distributions in the far detector within the realistic scenario and analyze them partly neglecting the detector effects implemented in the migration matrices .
accounting for the effect of finite detector resolution , we smear observables according to the normal distributions centered at their true values . for muons ,
this procedure is applied both to the momentum and production angle , using the realistic parameters  @xcite @xmath84 the energy resolutions of for @xmath85 s producing electromagnetic showers and other hadrons are set to @xmath86 respectively .
their values , as well as those of the energies appearing in the above equation , are expressed in units of gev .    in our calculations ,
the detection thresholds correspond to the measured kinetic energy of 20 mev for mesons and 40 mev for protons . for the sake of simplicity ,
the efficiencies are treated as energy - independent and set to 60% for @xmath85 s , 80% for other mesons , and 50% for protons .
owing to difficulties of their reconstruction in neutrino events , we assume that neutrons always escape detection .    in the context of the @xmath32 disappearance analysis  @xcite , we analyze the t2k - like setup described in , following the same treatment of systematic uncertainties and applying the same method of @xmath47 calculations .
our calculations are performed for the expected number of unoscillated events @xmath104900 .    instead of varying individual parameters involved in the detector description ,
we vary the migration matrices making a linear combinations of the those calculated with and without detector effects . in this manner
, we estimate the general sensitivity of the oscillation analysis to the detector performance , which can be expected to characterize a broad class of experiments rather than a particular ( highly idealized ) setup .
impact of detector - related uncertainties on the oscillation analysis employing the kinematic energy reconstruction . left : comparison of the event distributions for the detector performance overestimated by 10% ( dashed line ) and 30% ( solid line ) with an accurate estimate of detector effects ( shaded histogram ) . right :
@xmath87 confidence regions in the @xmath67 plane obtained when the data simulated with detector effects are fitted using the migration matrices neglecting them at the 10% , 20% , 30% level ( lines ) and accounting for them accurately ( shaded area ) .
, scaledwidth=39.0% ]    as shown in , the kinematic method of energy reconstruction  based only on the muon s energy and its production angle  turns out to be largely insensitive to detector effects and even the uncertainties at the level of 30% do not significantly affect the outcome of the oscillation analysis .
it is a consequence of very precise reconstruction of muon kinematics achieved in modern experiments and assumed in this analysis , see  .
same as figure [ fig : deteff_kin ] but for the calorimetric energy reconstruction.,scaledwidth=39.0% ]    on the other hand , we observe that detector effects play an important role in the calorimetric reconstruction method , in which the neutrino energy is inferred from the energies deposited in detector by all interaction products .
shows that in order to avoid a significant bias in the extracted values of the oscillation parameters , the detector response has to be determined with an accuracy of at least 10% .
this behavior can be traced back to much larger uncertainties in the hadron - energy determination , see , entering the calorimetric method of neutrino energy reconstruction .
while the results of figures [ fig : deteff_kin ] and [ fig : deteff_cal ] are obtained for the narrow beam peaked at @xmath100.6  gev  @xcite , the conclusions on influence of detector effects seem to be valid in more general case , as suggested by the findings for a  wide - band beam with the peak at @xmath101.6 gev and a non - negligible contribution of energies above 3 gev  @xcite .    owing to the very complicated oscillation probability ,
an accurate neutrino - energy reconstruction is of particular importance for appearance experiments , observing electron ( anti)neutrinos from oscillations of muon ( anti)neutrinos . in this context
, we analyze the role of missing energy  @xcite  carried away from event by undetected particles  on the @xmath9 sensitivity for an experiment similar to dune  @xcite .    in the considered setup ,
a wide - band neutrino beam is produced in interactions of the initial 1.08-mw proton beam with the target material and aimed at the far detector of fiducial mass 40 kton , located 1300 km from the source .
we assume 6 years of collecting data , 3 in @xmath12 mode and 3 in @xmath88 mode . for the signal ,
we consider 2% uncertainties of normalization ( bin - to - bin correlated ) and shape ( bin - to - bin uncorrelated ) .
for the background , only a global normalization uncertainty of 5% is taken into account  @xcite .
the simulated event distributions are obtained accounting for all detector effects  resolutions , efficiencies , and thresholds . on the other hand , extracting the oscillation parameters , we construct the migration matrices as linear combinations of those calculated with and without the shift resulting from the missing energy . in this way , for the purpose of this analysis , we single out the role of missing energy from other detector effects .
impact of an underestimation of the missing energy on the oscillation analysis employing the calorimetric energy reconstruction . left : reconstructed energy distributions for deep - inelastic scattering of electron neutrino at true energy 2.95 gev , obtained accounting for 0% , 50% and 100% of the shift due to the missing energy .
right : @xmath87 confidence regions in @xmath89 plane obtained when the simulated data are fitted using the migration matrices accounting for 90% , 80% and 70% ( lines ) of the missing energy and when all the missing energy is accounted for ( shaded area ) .
the dot shows the assumed true values of the oscillation parameters.,scaledwidth=39.0% ]    its influence on the reconstructed - energy distribution is illustrated in the left panel of for the example of deep - inelastic scattering of electron neutrino .
it clearly appears that owing to undetected particles in the final state , the maximum of the distribution is shifted to an energy lower than the true neutrino energy .
as the composition of particles in the final state differs for neutrinos and antineutrinos and depends on the interaction channel and the probe s energy , this is also the case for the energy carried away by undetected particles .
this can be expected to have important consequences for the oscillation analysis .
as presented in the right panel of , even a 20% underestimation of the missing energy introduces a sizable bias in the extracted @xmath9 value .
should the missing energy be underestimated by 30% , the analysis would exclude the true value of @xmath9 at a confidence level between 2 and 3@xmath33 .
this result illustrates the importance of an accurate determination of detector response in test - beam exposures and the relevance of a realistic simulation of nuclear effects in neutrino interactions , including intranuclear cascade .
we would like to emphasize that while our results are meant to point out the importance of various aspects of detector effects , much more detailed studies are necessary to draw truly quantitative conclusions for specific experiments .
in particular , our treatment of the missing - energy uncertainty  assumed to be equal for neutrinos and antineutrinos and independent of the energy and interaction channel  may be regarded as simplistic .
however , as more realistic sensitivity estimates would require an accurate knowledge of the detector response and inclusion of nuclear - model uncertainties , out of necessity , we leave them for future investigations within experimental collaborations .
thanks to tremendous progress in experimental neutrino physics over the last two decades , systematic uncertainties in oscillation studies have been greatly reduced . as a consequence ,
however , those arising from description of nuclear effects and related to determination of detector response have become very important in current appearance measurements .
this is going to be even more so in future experiments aiming to unambiguously discover the lepton - sector contribution to violation of charge - particle symmetry .    in this topical review , we have discussed a few examples illustrating the relevance of an accurate description of quasielastic interactions with any number of knocked - out nucleons and argued that for pion production a similar effect is expected .
we pointed out that comparisons to electron scattering data provide opportunity to test nuclear models and estimate their uncertainties ; they also allow for much clearer interpretation of discrepancies than neutrino data .
discussing detector effects , we have drawn attention to their larger importance for experiments employing the calorimetric method of energy reconstruction and emphasized the key role of an accurate estimate of the missing energy for precise oscillation studies .    as a final remark
we reiterate that future experiments are going to require coordinated efforts toward improving theoretical descriptions of nuclear effects , implementing them in monte carlo generators and determining detector response in extensive test - beam exposures .
however , the discoveries coming in the next two decades have the potential to overshadow those of the past two decades .
this work is supported by the national science foundation under grant phy-1352106 .
999 kajita t 2016 030501 mcdonald a b 2016 030502 kim s - b 2015 93 an f ( juno collaboration ) 2016 030401 patterson r b 2015 _ ann .
nucl . part .
sci . _ * 65 * 177 acciarri r ( dune collaboration ) 2015 long - baseline neutrino facility ( lbnf ) and deep underground neutrino experiment ( dune ) .
conceptual design report .
volume 2 : the physics program for dune at lbnf ( arxiv:1512.06148 ) abe k ( hyper - kamiokande proto - collaboration ) 2015 _ prog .
phys . _ * 2015 * 053c02 abe k ( t2k collaboration ) 2011 a * 659 * 106 ayres d s ( nova collaboration ) 2007 the nova technical design report , fermilab - design-2007 - 01 acciarri r ( dune collaboration ) 2016 long - baseline neutrino facility ( lbnf ) and deep underground neutrino experiment ( dune ) .
conceptual design report .
volume 1 : the lbnf and dune projects ( arxiv:1601.05471 ) formaggio j a and zeller g p 2012 1307 moreno o , donnelly t w , van orden j w and ford w p 2015 d * 92 * 053006 ankowski a m , benhar o , coloma p , huber p , jen c - m , mariani c , meloni d , and vagnoni e 2015 d * 92 * 073014 ankowski a m , coloma p , huber p , mariani c and vagnoni e 2015 d * 92 * 091301(r ) wu q ( nomad collaboration ) 2008 b * 660 * 19 lyubushkin v v ( nomad collaboration ) 2009 _ eur .
phys . j. _
c * 63 * 355 aguilar - arvalo a a ( miniboone collaboration ) 2010 d * 81 * 092005 aguilar - arvalo a a ( miniboone collaboration ) 2010 d * 82 * 092005 nakajima y ( sciboone collaboration ) 2011 d * 83 * 012005 aguilar - arvalo a a ( miniboone collaboration ) 2011 d * 83 * 052007 aguilar - arvalo a a ( miniboone collaboration ) 2011 d * 83 * 052009 fields l ( minerva collaboration ) 2013 022501 fiorentini g a ( minerva collaboration ) 2013 022502 abe k ( t2k collaboration ) 2013 d * 87 * 092003 aguilar - arvalo a a ( miniboone collaboration ) 2013 d * 88 * 032001 abe k ( t2k collaboration ) 2014 d * 90 * 052010 tice b g ( minerva collaboration ) 2014 231801 higuera a ( minerva collaboration ) 2014 261802 aguilar - arvalo a a ( miniboone collaboration ) 2015 d * 91 * 012004 walton t ( minerva collaboration ) 2015 d * 91 * 071301(r ) abe k ( t2k collaboration ) 2015 d * 91 * 112002 eberly b ( minerva collaboration ) 2015 d * 92 * 092008 le t ( minerva collaboration ) 2015 b * 749 * 130 abe k ( t2k collaboration ) 2015 d * 92 * 112003 rodrigues p a ( minerva collaboration ) 2016
071802 mousseau j ( minerva collaboration ) 2016 d * 93 * 071101 abe k ( t2k collaboration ) 2016 d * 93 * , 112012 marshall c m ( minerva collaboration ) 2016 d * 94 * 012002 abe k ( t2k collaboration ) 2014 d * 90 * 072012 abe k ( t2k collaboration ) 2016 first measurement of the muon neutrino charged current single pion production cross section on water with the t2k near detector ( arxiv:1605.07964 ) anderson c ( argoneut collaboration ) 2012 161802 acciarri r ( argoneut collaboration ) 2014 d * 89 * 112003 acciarri r ( argoneut collaboration ) 2014 d * 90 * 012008 acciarri r ( argoneut collaboration ) 2014 261801 tzanov m ( nutev collaboration ) 2006 d * 74 * 012008 adamson p ( minos collaboration ) 2010 d * 81 * 072002 abe k ( t2k collaboration ) 2016 d * 93 * 072002 adamson p ( minos collaboration ) 2016 measurement of single @xmath85 production by coherent neutral - current @xmath12 fe interactions in the minos near detector ( arxiv:1608.05702 ) benhar o and meloni d 2009 d * 80 * 073003 leitner t and mosel u 2010 c * 81 * 064614 martini m , ericson m , and chanfray g 2012 d * 85 * 093012 nieves j , sanchez f , ruiz simo i , and vicente vacas m j 2012 d * 85 * 113008 meloni d and martini m 2012 b * 716 * 186 lalakulich o , and mosel u , and gallmeister k 2012 c * 86 * 054606 erratum : 2014 c * 90 * 029902(e ) martini m , ericson m , and chanfray g 2012 d * 87 * 013009 coloma p and huber p 2013 221802 coloma p , huber p , jen c - m and mariani c 2014 d * 89 * 073015 mosel u , and lalakulich o , and gallmeister k 2014 151802 jen c - m , ankowski a m , benhar o , furmanski a p , kalousis l n , and mariani c 2014 d * 90 * 093004 ankowski a m , benhar o and sakuda m 2015 d * 91 * 033005 ericson m , garzelli m v , giunti c and martini m 2016 d * 93 * 073008 ankowski a m , benhar o , mariani c and vagnoni e 2016 d * 93 *
113004 alvarez - ruso l , hayato y and nieves j 2014 075015 benhar o , huber p , mariani c and meloni d 2015 neutrino - nucleus interactions and the determination of oscillation parameters ( arxiv:1501.06448 ) lu x - g , coplowe d , shah r , barr g , wark d , and weber a 2015 d * 92 * 051302 benhar o , coletti p and meloni d 2010 105 132301 michael d g ( minos collaboration ) 2008 a * 596 * 190 aliaga l ( minerva collaboration ) 2014 a * 743 * 130 abramov a g , galyaev n a , garkusha v i , hylen j , novoskoltsev f n , ryabov a d and zarucheisky v g 2002 a * 485 * 209 adey d ( @xmath12storm collaboration ) 2013 @xmath12storm  neutrinos from stored muons : proposal to the fermilab pac , fermilab - proposal-1028 ( arxiv:1308.6822 ) an f p ( daya bay collaboration ) 2015 111802 choi j h ( reno collaboration ) 2016 211801 abe y ( doublechooz collaboration ) 2014 jhep*10*(2014)086 michael d g
( minos collaboration ) 2006 191801 abe k ( t2k collaboration ) 2014 d * 89 * 092003 abe k ( t2k collaboration ) 2014 241803 abe k ( t2k collaboration ) 2015 d * 91 * 112010 wolcott j ( minerva collaboration ) 2016 081802 ritz s ( hepap subcommittee ) 2014 building for discovery : strategic plan for u.s .
particle physics in the global context ( http://science.energy.gov/~/media/hep/hepap/pdf/may-2014/final_p5_report_interactive_060214.pdf ) gonzalez - garcia m c , maltoni m , salvado j and schwetz t 2012 jhep*12*(2012)123 abe k ( t2k collaboration ) 2016 181801 adamson p ( nova collaboration ) 2016 d * 93 * 051104 huber p , lindner m and winter w 2005 195 huber p , kopp j , lindner m , rolinec m and winter w 2007 432 coloma p , huber p , kopp j and winter w 2013 d * 87 * 033004 huber p , lindner m , schwetz t and winter w 2009 jhep*11*(2009)044 nunokawa h , parke s j and zukanovich funchal r 2005 d * 72 * 013009 andreopoulos c 2010 a * 614 * 87 benhar o , fabrocini a , fantoni s and sick i 1994 a * 579 * 493 katori t 2015 _ aip conf .
* 1663 * 030001 ahn m h ( k2k collaboration ) 2003 041801 bodek a and ritchie j l 1981 d * 23 * 1070 benhar o , day d and sick i 2008 80 189 smith r a and moniz e j 1972 b * 43 * 605 erratum : 1975 b * 101 * 547 baran d t 1988 400 sealock r m 1989 1350 benhar o and meloni d 2007 a * 789 * 379 gran r ( k2k collaboration ) 2006 d * 74 * 052002 aguilar - arvalo a a ( miniboone collaboration ) 2008 032301 mariani c ( k2k collaboration ) 2008 _ aip conf .
* 981 * 247 adamson p ( minos collaboration ) 2015 d * 91 * 012005