NOvA
(with focus on the mass hierarchy)
Jeff Hartnell
University of Sussex Solvay Workshop, Brussels 30th November 2017
NOvA (with focus on the mass hierarchy) Jeff Hartnell University - - PowerPoint PPT Presentation
NOvA (with focus on the mass hierarchy) Jeff Hartnell University of Sussex Solvay Workshop, Brussels 30 th November 2017 IntroducAon NOvA experiment and physics goals NuMI beam NOvA detectors Mass hierarchy via MSW maHer effect
(with focus on the mass hierarchy)
Jeff Hartnell
University of Sussex Solvay Workshop, Brussels 30th November 2017
IntroducAon
– NuMI beam – NOvA detectors
– Muon neutrino disappearance – NC analysis – Electron neutrino appearance
Jeff Hartnell, Solvay 2017 2
Jeff Hartnell, Solvay 2017 3
NOvA Overview
– measure beam composiAon – energy spectrum
– measure oscillaAons and search for new physics
Ash River
Ash River
810 km
The NOvA CollaboraAon
242 Collaborators 49 institutions 7 countries
Argonne, AtlanAco, Banaras Hindu University, Caltech, Cochin, InsAtute of Physics and Computer science of the Czech Academy of Sciences, Charles University, CincinnaA, Colorado State, Czech Technical University, Delhi, JINR, Fermilab, Goiás, IIT GuwahaA, Harvard, IIT Hyderabad, U. Hyderabad, Indiana, Iowa State, Jammu, Lebedev, Michigan State, Minnesota-Twin CiAes, Minnesota-Duluth, INR Moscow, Panjab, South Carolina, SD School
UCL,Virginia, Wichita State, William and Mary, Winona State
4 Jeff Hartnell, Solvay 2017
Physics Goals
Jeff Hartnell, Solvay 2017 5
Results from 3 different oscillation analyses
¨ Disappearance of
νµ CC events
¤ clear suppression as a
funcAon of energy
¤ 2016 analysis results
PRL 118.151802
sin2(2θ23)
32
events
¤ 2 GeV neutrinos
enhances maHer effects
¤ ±30% effect
¤ 2016 analysis results
in PRL 118.231801.
θ13, θ23, δCP , and Mass Hierarchy
∆m2
41, θ34, θ24
¨ Deficit of NC events?
¤ suppression of NCs could be evidence
neutrino
¤ Fit to 3+1model ¤ new!
Jeff Hartnell, Solvay 2017 6
Off-axis
On-axis
NuMI Beam
Jeff Hartnell, Solvay 2017 7
To APD 4 cm ⨯ 6 cm 1560 cm
A NO𝜉A cell
NO𝜉A detectors
Fiber pairs from 32 cells 32-pixel APD
Far detector:
14-kton, fine-grained, low-Z, highly-active tracking calorimeter → 344,000 channels
Near detector:
0.3-kton version of the same → 20,000 channels Extruded PVC cells filled with 11M liters of scintillator instrumented with 𝜇-shifting fiber and APDs
Jeff Hartnell, Solvay 2017 8
Long-baseline neutrino oscillations
𝜉𝜈 disappearance:
…to leading order
experimental data are consistent with unity (“maximal mixing”) Need a leap in precision on 𝜄23 (and m2 )
32
𝜉e appearance:
Daya Bay reactor experiment: sin2(2𝜄13) = 0.084 ± 0.005
…plus potentially large CPv and matter effect modifications!
Non-zero 𝜄 opens the long-baseline appearance channel, and…
Jeff Hartnell, Solvay 2017 9
StarAng with νμ
νµ ντ νe
L/E (km/GeV) 0.2 0.4 0.6 0.8 1
Oscillation Probability
1000 2000
How does the mass hierarchy come into play?
Jeff Hartnell, Solvay 2017 10
Δm2
31 and Δm2 32 differ by 3%
Small effect JUNO’s planned measurement involves this
MaHer Effect & Mass Hierarchy
maHer not anAmaHer
– electron density causes asymmetry (fake CPv!)
– different Feynman diagrams for νe and νe interacAons with electrons so different amplitudes
Jeff Hartnell, Solvay 2017 11
Arrows flip for antineutrinos
Jeff Hartnell, Solvay 2017 12
Long-baseline 𝜉𝜈→𝜉e
For fixed L/E = 0.4 km/MeV
A more quantitative sketch… At right: P(𝜉 ⎺𝜈→ 𝜉 ⎺e) vs. P(𝜉𝜈→𝜉e) plotted for a single neutrino energy and baseline
Jeff Hartnell, Solvay 2017 13
Long-baseline 𝜉𝜈→𝜉e
For fixed L/E = 0.4 km/MeV
A more quantitative sketch… At right: P(𝜉 ⎺𝜈→ 𝜉 ⎺e) vs. P(𝜉𝜈→𝜉e) plotted for a single neutrino energy and baseline Measure these probabilities (an example measurement
Also: Both probabilities ∝ sin2𝜄23
Jeff Hartnell, Solvay 2017 14
Non-maximal mixing scenario
then effect of octant is important
inverted% hierarchy normal% hierarchy Θ23%<%45o Θ23%>%45o
Jeff Hartnell, Solvay 2017 15
Effect of Increasing Energy
1 2 3 4 5 6 7 8
<P(!µ !e)> [%]
1 2 3 4 5 6 7 8
<P(!µ !e)> [%]
Normal Mass Hierarchy Inverted Mass Hierarchy " = 0 " = #/2 " = # " = 3#/2 L = 1300 km, <E> = 3.2 GeV sin 22$13 = 0.09𝜉𝜈→𝜉
For fixed L/E = 0.4 km/MeVA more quantitative sketch… 𝜉 ⎺𝜈→𝜉 ⎺ 𝜉𝜈→𝜉 ) ino
Increasing Energy 0.6 GeV 2 GeV 3 GeV T2K NOvA DUNE
[à bigger matter effect and hence bigger fake CP violation]
The measurements
Jeff Hartnell, Solvay 2017 16
νμ
e
νe
p μ p
1m 1mν
νµ CC νe CC NC
~5m ~2.5m
Long, straight track Shorter, wider, fuzzy shower Diffuse activity from nuclear recoil system
Event Types
Jeff Hartnell, Solvay 2017 17
Jeff Hartnell, Solvay 2017 18
(simulated 𝜉𝜈 CC event)
the Far and Near energy spectra
νμ Near Detector Data
Jeff Hartnell, Solvay 2017 19
–
μ ν
μ ν
–
μ
νμ Far Detector Data
Jeff Hartnell, Solvay 2017 20
Reconstructed neutrino energy (GeV)
1 2 3 4 5
Events / 0.25 GeV
5 10 15 20
Prediction, no systs.
σ 1- Prediction with systs. Backgrounds Data Normal Hierarchy
NOvA Preliminary
78 events observed in FD
– 473±30 with no oscilla5on – 82 at best oscillaAon fit – 3.9 beam BG + 2.7 cosmic
νμ Disappearance Result
Jeff Hartnell, Solvay 2017 21
23θ
2sin
0.3 0.4 0.5 0.6 0.7)
2eV
(10
32 2m ∆
2 2.5 3 3.5NOvA Preliminary
Normal Hierarchy, 90% CL NOvA 2016 T2K 2014 MINOS 2014
Best Fit (in NH):
32
sin2 θ23 = 0.40+0.03
−0.02(0.63+0.02 −0.03) No FC Correction
Maximal mixing excluded at 2.6σ
Driven by bins in oscillation dip (1-2 GeV). Forcing maximal mixing gives:
∆m2
32 = (2.46) × 10−3eV2
νμ
–𝜠𝒏𝟒𝟑
𝟑
Neutral Current Result
(NOvA’s first 2017 dataset result, presented at NuFact Sep/17)
Jeff Hartnell, Solvay 2017 22
NC Far Detector Data & Results
Jeff Hartnell, Solvay 2017 23 –
Observed 214 NC candidates Prediction 191.16 ± 13.82(stat.)±21.99 (syst.) No depletion of NC events observed NOvA sees no evidence for νs mixing
– θ
𝟐. 𝟐𝟘𝟏 ± 0.160 (𝑡𝑢𝑏𝑢. )−0.130
+0.080 (𝑡𝑧𝑡𝑢. )
θ
𝟐. 𝟐𝟘𝟏 ± 0.123 (𝑡𝑢𝑏𝑢. )−0.124
+0.143 (𝑡𝑧𝑡𝑢. )
θ
𝟐. 𝟐𝟖𝟘 ± 0.123 (𝑡𝑢𝑏𝑢. )−0.124
+0.142 (𝑡𝑧𝑡𝑢. )
θ
𝟐. 𝟐𝟖𝟕 ± 0.123 (𝑡𝑢𝑏𝑢. )−0.124
+0.142 (𝑡𝑧𝑡𝑢. )
No NC disappearance → R = ’s two degenerate best fit points θ23 Δ𝑛32
2
δ𝐷𝑄
–
“ ” Δ𝑛41
2
𝜾𝟑𝟓 𝜾𝟒𝟓
θ23 Δ𝑛32
2
δ𝐷𝑄
δ24
𝜄34
Jeff Hartnell, Solvay 2017 24
(simulated 𝜉e CC event)
in the Far Detector
as 𝜉e appearance
νe Near Detector Data
Jeff Hartnell, Solvay 2017 25
– NC, CC, beam νe each propagate differently – constrain beam νe using selected νµ CC spectrum – constrain νµ CC using Michel Electron distribuAon
beam νe up by 4% NC up by 17% νµ CC up by 10%
PredicAon
Jeff Hartnell, Solvay 2017 26
Total BG NC Beam νe νµ CC ντ CC Cosmics 8.2 3.7 3.1 0.7 0.1 0.5 NH, 3π/2, IH, π/2, 28.2 11.2 Signal events (±5% systematic uncertainty): Background by component (±10% systematic uncertainty):
¨ Extrapolate each component in
bins of energy and CVN output
¨ Expected event counts depend
CP
δ Total events expected
10 20 30 40 50 2 π π 2 π 3 π 2 POT equiv.
2010 × 6.05 NOvA FD =0.4-0.6
23θ
2sin NH IH
NOvA Simulation
νe Far Detector Data
Jeff Hartnell, Solvay 2017 27
events
Ø background 8.2 ± 0.8
>8σ electron neutrino appearance signal
CVN=0.991 E=1.63 GeV
Reconstructed neutrino energy (GeV) Events / 0.5 GeV Bin
5 10 15 20 1 2 3 1 2 3 1 2 3
0.75 < CVN < 0.87 0.87 < CVN < 0.95 0.95 < CVN < 1
NH
NOvA Preliminary
FD Data Best Fit Prediction Total Background Cosmic Background POT equiv.
2010 × 6.05
–
νμ
𝟏. 𝟓𝟏𝟓 δ π 𝜾𝟑𝟒 𝟏. 𝟕𝟑𝟒 , δ π
𝟒π 𝟑
σ
Joint νe + νμ Fit Contours
Jeff Hartnell, Solvay 2017 28
– Constrain sin2(2θ13)=0.085±0.05 – Constrain Δm2 and sin2θ23 with NOvA disappearance results
– best fit IH-NH, Δ𝜓2=0.47 – both octants & hierarchies allowed at 1σ – 3σ exclusion in IH, lower octant around 𝜺CP=π/2
δCP = 1.49π sin2(θ23) = 0.40
–
νμ
𝟏. 𝟓𝟏𝟓 δ π 𝜾𝟑𝟒 𝟏. 𝟕𝟑𝟒 , δ π
𝟒π 𝟑
σ
Contours
Jeff Hartnell, Solvay 2017 29
Antineutrino data will help resolve degeneracies, particularly for non-maximal mixing. Results planned for summer 2018
inverted% hierarchy normal% hierarchy Θ23%<%45o Θ23%>%45o
Future SensiAvity
Jeff Hartnell, Solvay 2017 30
Lower Octant Upper Octant
Year
2016 2018 2020 2022 2024
) σ Significance (
1 2 3 4 5 =0.022
13θ
2, sin
2eV
10 × =2.5
32 2m ∆ =0.625
23θ
2/2, sin π =3
CPδ Normal
and analysis improvements All projected beam intensiy
NOvA Simulation
µν +
eν NOvA joint
Hierarchy Octant CPV Year
2016 2018 2020 2022 2024
) σ Significance (
1 2 3 4 5 =0.022
13θ
2, sin
2eV
10 × =2.5
32 2m ∆ =0.403
23θ
2/2, sin π =3
CPδ Normal
and analysis improvements All projected beam intensiy
NOvA Simulation
µν +
eν NOvA joint
Hierarchy Octant CPV
Conclusions
With 6.05x1020 POT, NOvA finds:
– Maximal mixing excluded at 2.6σ
– Data prefers NH at low significance – IH, lower octant, 𝜺CP=π/2 region excluded at 3σ
– Neutral current events show no evidence of steriles – New νe and νμ results very soon
Jeff Hartnell, Solvay 2017 31
Backup slides
Jeff Hartnell, Solvay 2017 32
Future SensiAvity
Jeff Hartnell, Solvay 2017 33
Year
2016 2018 2020 2022 2024
) σ Significance (
1 2 3 4 5 =0.022
13θ
2, sin
2eV
10 × =2.5
32 2m ∆ =0.625
23θ
2/2, sin π =3
CPδ Normal
systematic uncertainty improvements 2016 analysis techniques with projected
µν +
eν NOvA joint
Hierarchy Octant CPV
NOvA Simulation
Year
2016 2018 2020 2022 2024
) σ Significance (
1 2 3 4 5 =0.022
13θ
2, sin
2eV
10 × =2.5
32 2m ∆ =0.403
23θ
2/2, sin π =3
CPδ Normal
systematic uncertainty improvements 2016 analysis techniques with projected
µν +
eν NOvA joint
Hierarchy Octant CPV
NOvA Simulation
Lower Octant Upper Octant
Jeff Hartnell, Solvay 2017 34
νμ à νe appearance probability
Jeff Hartnell, INSS 2016 34
[PDG, 2014]
MaHer Effect & Mass Hierarchy
maHer not anAmaHer
– electron density causes the asymmetry
– different Feynman diagrams for νe and νe interacAons with electrons...
Jeff Hartnell, Solvay 2017 35
Different Feynman Diagrams
neutrino interacAon with an electron
an5neutrino interacAon with an electron
Jeff Hartnell, Solvay 2017 36
+
Electron neutrinos and anAneutrinos are affected differently by interacAons with maHer à fake CP violaAon Why does the mass hierarchy affect oscillaAons involving electron (anA)neutrinos?
Jeff Hartnell, Solvay 2017 37
MaHer effect (neutrino case)
the mass eigenstates
– strength depends on electron neutrino content of each mass eigenstate
Jeff Hartnell, Solvay 2017 38
No matter effect matter effect for NH matter effect for IH Splitting changes size in presence of matter
AnAneutrino case
the mass eigenstates
– strength depends on electron neutrino content of each mass eigenstate
Jeff Hartnell, Solvay 2017 39
No matter effect matter effect for NH matter effect for IH Splitting changes size in presence of matter
Splivngs and mixing angles affected
the mass squared splivng in maHer (Δm2
M)
– e.g. simple 2-flavour case:
– Also see it in full 3-flavour equaAons (a few slides back)
Jeff Hartnell, Solvay 2017 40
Improved Event SelecAon
Jeff Hartnell, Solvay 2017 41
based on ideas from computer vision and deep learning
¨ Calibrated hit maps are
inputs to ConvoluAonal Visual Network (CVN)
¨ Series of image processing
transformaAons applied to extract abstract features
¨ Extracted features used as
inputs to a convenAonal neural network to classify the event
Improved Event SelecAon
Jeff Hartnell, Solvay 2017 42
based on ideas from computer vision and deep learning
¨ Calibrated hit maps are
inputs to ConvoluAonal Visual Network (CVN)
¨ Series of image processing
transformaAons applied to extract abstract features
¨ Extracted features used as
inputs to a convenAonal neural network to classify the event
Improved Event SelecAon
Jeff Hartnell, Solvay 2017 43
Improvement in sensitivity from CVN equivalent to 30% more exposure
based on ideas from computer vision and deep learning
¨ Calibrated hit maps are
inputs to ConvoluAonal Visual Network (CVN)
¨ Series of image processing
transformaAons applied to extract abstract features
¨ Extracted features used as
inputs to a convenAonal neural network to classify the event
[A. Aurisano et al., arXiv:1604.01444]
Jeff Hartnell, Solvay 2017 44
We consider multiple possible sources of systematic error
In each case:
through the extrapolation
as pull terms in the fit
quadrature) of the parameter measurement error is recorded
Systematic* Effect*on* sin2(θ23) Effect on Δm232 Normalisation ± 1.0% ± 0.21% Muon1E1scale ± 2.2% ± 0.81% Calibration ± 2.01% ± 0.21% Relative1E1scale ± 2.01% ± 0.91% Cross1sections1+1FSI ± 0.61% ± 0.51% Osc.1parameters ± 0.71% ± 1.51% Beam1backgrounds ± 0.91% ± 0.51% Scintillation1model ± 0.71% ± 0.11% All*systematics ± 3.4*% ± 2.4*% Stat.*Uncertainty ± 4.1*% ± 3.5*%
νμ Event SelecAon
Jeff Hartnell, Solvay 2017 45
νµCC events less than 5GeV
– Select events with long tracks – Suppress NC and cosmic backgrounds
muons
– track length – dE/dx along track – scaHering along track – track-only plane fracAon
well for muon variables
Muon ID
0.2 0.4 0.6 0.8 1Events
3 10 4 10 5 10 6 10 7 10Simulated selected events Simulated background Data
σ Full 1- POT
2010 × ND POT norm., 3.72
NOvA Preliminary
dE/dx Log-likelihood
3 − 2 − 1 − 1
Events
0.0 0.2 0.4 0.6
610 × Simulated Selected Events Simulated Background Data
σ Shape-only 1- POT
2010 × ND area norm., 3.72
NOvA Preliminary
Jeff Hartnell, Solvay 2017 46
Calorimetric energy (GeV)
1 2 3 4 5 Events 5 10 15 20
FD Data Best-fit prediction Background
POT-equiv.
2010 × 6.05
NOvA Preliminary
Hadronic energy (GeV)
1 2 3 4 5 Events 10 20 30
FD Data Best-fit prediction Background
POT-equiv.
2010 × 6.05
NOvA Preliminary
Length of main track (m)
5 10 15 20 25 Events 5 10 15 20
FD Data Best-fit prediction Background
POT-equiv.
2010 × 6.05
NOvA Preliminary
Our best fit oscillation prediction matches other distributions well
Goodness of fit
Jeff Hartnell, Solvay 2017 47
Events
5 10 15 20 Normal HierarchyNOvA Preliminary
POT-equiv. 20 10 × NOvA 6.05Prediction Data
Reconstructed neutrino energy (GeV)
1 2 3 4 5contribution
2χ
2 4 6 8 10 12χ Total
There is no significant pull in the oscillation fit from bins in the tail
23θ
2sin
0.3 0.4 0.5 0.6 0.7)
2eV
(10
32 2m ∆
2 2.5 3 3.5NOvA Preliminary
Normal HierarchyPOT-equiv.
2010 × NOvA 6.05 90% C.L. 0-5 GeV Analysis 90% C.L. 0-2.5 GeV Analysis
Jeff Hartnell, Solvay 2017 48
CP
δ ) σ Significance (
1 2 3 4 5 2 π π 2 π 3 π 2
NH lower octant lower octant n IH NH upper octant upper octant n IH POT equiv.
20
10 × 6.05 NOvA FD
ScaHering in a Nuclear Environment
Jeff Hartnell, Solvay 2017 49 Reco “q0” (=Ehad,vis) 103 Events
10 20 10 20 10 20 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 10 20 10 20 10 20 1.0 1.00.8 < |q|/GeV < 0.9 0.9 < |q|/GeV < 1 0.7 < |q|/GeV < 0.8 0.4 < |q|/GeV < 0.5 0.5 < |q|/GeV < 0.6 0.6 < |q|/GeV < 0.7 0.2 < |q|/GeV < 0.3 0.3 < |q|/GeV < 0.4 0.1 < |q|/GeV < 0.2
NOvA ND DataNOvA Preliminary
suggests unsimulated process between quasi- elasAc and delta producAon
Similar conclusions from MINERvA data reported in P.A. Rodrigues et al., PRL 116 (2016) 071802
Jeff Hartnell, Solvay 2017 50 Reco “q0” (=Ehad,vis) 103 Events
10 20 10 20 10 20 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 10 20 10 20 10 20 1.0 1.00.8 < |q|/GeV < 0.9 0.9 < |q|/GeV < 1 0.7 < |q|/GeV < 0.8 0.4 < |q|/GeV < 0.5 0.5 < |q|/GeV < 0.6 0.6 < |q|/GeV < 0.7 0.2 < |q|/GeV < 0.3 0.3 < |q|/GeV < 0.4 0.1 < |q|/GeV < 0.2
NOvA ND DatasystemaAcs
– hadronic energy scale – QE cross secAon modeling
pion producAon by 50% (P.A. Rodrigues et al,
arXiv:1601.01888.)
¨ Enable GENIE empirical Meson Exchange Current Model ¨ Reweight to match NOvA excess as a funcAon of 3-
momentum transfer
MEC model by S. Dytman, inspired by
ScaHering in a Nuclear Environment