PHYSICS @ INOBeyond 3-flavor osci%ation physics.. Monday 12 August 2013 Large fluctuations in ΔP...

17
PHYSICS @ INO SANDHYA CHOUBEY Harish- Chandra Research Institute, India NuFact 2013, August 19-24, Beijing, China 1

Transcript of PHYSICS @ INOBeyond 3-flavor osci%ation physics.. Monday 12 August 2013 Large fluctuations in ΔP...

  • PHYSICS @ INO

    SANDHYA CHOUBEYHarish-Chandra Research Institute, India

    NuFact 2013, August 19-24, Beijing, China

    1

  • Sandhya Choubey NuFact 2013, August 21

    The Unknowns...after θ13

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle θ23 …

    Beyond 3-flavor osci%ation physics..

    2

  • Sandhya Choubey NuFact 2013, August 21

    The Unknowns...after θ13

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle θ23 …

    Beyond 3-flavor osci%ation physics..

    3

    We wi% cover in this talk

    3 posters and a talk on details of INOby Kaur, Laksmi, Meghna

    using atmospheric neutrinos

    Plenary talk by Naba Mondal on Friday

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Neutrinos

    4

    such that this issue could be alleviated. It was shown that with this if one defined normal hierarchy as∆m2eff > 0 and inverted hierarchy as ∆m

    2eff < 0, then Pµµ is the same for normal and inverted hierarchies

    when θ13 = 0. Though this complication in the definition of the mass hierarchy does not make muchdifference to the final neutrino mass hierarchy sensitivity results when one marginalizes correctly over|∆m231| (for a very recent discussion see [28]), we will continue to use this definition for the normal andinverted hierarchy, especially since in this paper we will not marginalize our χ2 function over the oscilla-tion parameters, as discussed in the Introduction section.

    Energy[GeV]2 3 4 5 6 7 8 9 10 11

    (IH

    µ(N

    H)-P

    µµP

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6

    0.8L=7202 Km

    L=9110 Km

    L=11019 Km

    L=12611 Km

    !"cos-1 -0.8 -0.6 -0.4 -0.2 0

    (IH

    )µµ

    (NH

    )-PµµP

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6

    0.8=3 GeV!E

    =5 GeV!E

    =7 GeV!E

    =9 GeV!E

    Figure 1: The difference in the muon neutrino survival probability for normal and inverted hierarchy,∆Pµµ, for various neutrino energy and zenith angles. The left panel shows ∆Pµµ as a function of theneutrino energy for four different choices of the zenith angle corresponding to four different path lengthsfor the neutrinos in the earth matter. The right panel shows ∆Pµµ as a function of the neutrino zenithangle for four different choices of the neutrino energy.

    Parameter True value used in data

    ∆m221 7.5× 10−5 eV2

    sin2 θ12 0.3

    |∆m2eff | 2.4× 10−3 eV2

    δCP 0

    sin2 θ23 0.5

    sin2 2θ13 0.1

    Table 1: Benchmark true values of oscillation parameters used in the simulations.

    We show in Fig. 1 the difference between the survival probability predicted for the normal (PNHµµ )

    and inverted (P IHµµ ) hierarchy. This figure has been generated by calculating the full three generation

    6

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Large fluctuations in ΔP in both E as we% as Θ

    Ghosh, SC, arxiv:1306.1423

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Need good reconstruction in both E as we% as Θ - resolution

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    ΔP opposite for neutrino and antineutrino - charge id

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Muons

    5

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=11019 km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=12611 Km

    Figure 2: The difference in the predicted muon event rates for normal and inverted hierarchy, normalizedto the no-oscillation muon event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point of thebin. The solid black lines show the ∆Pµµ as in Fig. 1. The black dashed lines give the difference betweenthe predicted rate for normal and inverted hierarchy when only the Pµµ (disappearance) channel is takenand data binned in neutrino energy. The blue dot-dashed lines show the corresponding difference whenonly the Pµµ channel is taken and data binned in muon energy. The red dotted lines are obtained whenwe add the Peµ (appearance) channel and bin the data in muon energy. The pink dashed and cyan longdashed lines are obtained when we apply on the red dotted lines the muon energy resolution functionswith widths σEµ/Eµ = 0.15% and 2%, respectively.

    9

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=11019 km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )-1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=12611 Km

    Figure 2: The difference in the predicted muon event rates for normal and inverted hierarchy, normalizedto the no-oscillation muon event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point of thebin. The solid black lines show the ∆Pµµ as in Fig. 1. The black dashed lines give the difference betweenthe predicted rate for normal and inverted hierarchy when only the Pµµ (disappearance) channel is takenand data binned in neutrino energy. The blue dot-dashed lines show the corresponding difference whenonly the Pµµ channel is taken and data binned in muon energy. The red dotted lines are obtained whenwe add the Peµ (appearance) channel and bin the data in muon energy. The pink dashed and cyan longdashed lines are obtained when we apply on the red dotted lines the muon energy resolution functionswith widths σEµ/Eµ = 0.15% and 2%, respectively.

    9 Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to the ν-nucleon cross-section

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to presence of both νμ and νe in the atm flux

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to finite detector resolution

    Ghosh, SC, arxiv:1306.1423

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Muons

    6

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=11019 km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=12611 Km

    Figure 2: The difference in the predicted muon event rates for normal and inverted hierarchy, normalizedto the no-oscillation muon event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point of thebin. The solid black lines show the ∆Pµµ as in Fig. 1. The black dashed lines give the difference betweenthe predicted rate for normal and inverted hierarchy when only the Pµµ (disappearance) channel is takenand data binned in neutrino energy. The blue dot-dashed lines show the corresponding difference whenonly the Pµµ channel is taken and data binned in muon energy. The red dotted lines are obtained whenwe add the Peµ (appearance) channel and bin the data in muon energy. The pink dashed and cyan longdashed lines are obtained when we apply on the red dotted lines the muon energy resolution functionswith widths σEµ/Eµ = 0.15% and 2%, respectively.

    9

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=11019 km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )-1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=12611 Km

    Figure 2: The difference in the predicted muon event rates for normal and inverted hierarchy, normalizedto the no-oscillation muon event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point of thebin. The solid black lines show the ∆Pµµ as in Fig. 1. The black dashed lines give the difference betweenthe predicted rate for normal and inverted hierarchy when only the Pµµ (disappearance) channel is takenand data binned in neutrino energy. The blue dot-dashed lines show the corresponding difference whenonly the Pµµ channel is taken and data binned in muon energy. The red dotted lines are obtained whenwe add the Peµ (appearance) channel and bin the data in muon energy. The pink dashed and cyan longdashed lines are obtained when we apply on the red dotted lines the muon energy resolution functionswith widths σEµ/Eµ = 0.15% and 2%, respectively.

    9 Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to the ν-n cross-section

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to presence of both νμ and νe in the atm flux

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Smearing due to finite detector resolution

    Residual matter effects still present which gives MHGhosh, SC, arxiv:1306.1423Net sensitivity depends on the detector performance

    Two key ingredients: Resolution and Efficiency

  • Sandhya Choubey NuFact 2013, August 21

    MH @ INO (Only Muons)

    7

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =62"# =72"# =82"# =92"#=102"# =112"# =122"# =132"#=142"# =152"# =162"#

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =92"# =102"# =112"# =122"#=132"# =142"# =152"# =162"#=172"# =182"# =192"# =202"#

    Figure 6: Constant ∆χ2 contours in the reconstruction efficiency and energy resolution plane. The leftpanels shows the contours from the analysis which uses only the muon data from the experiments. Theright panels shows the contours for the analysis in which both the muon and hadron data were includedin the combined statistical analysis.

    present results for both 80 bins as well as 20 bins between 1 GeV and 11 GeV.

    The dependence of the mass hierarchy sensitivity to the muon energy resolution and muon recon-struction efficiency are shown in Figs. 5. In the left panel on Fig. 5 we show the ∆χ2 as a function ofthe muon energy resolution σEµ/Eµ, for different values of the muon reconstruction efficiency #. We varyσEµ/Eµ continuously between 20% and 2% and repeat this for # = 70% to 100%. For a given #, we seethat the ∆χ2 falls sharply as the σEµ/Eµ is increased from 2% to about 6%. Thereafter, the rate of fallof ∆χ2 reduces, and it falls steadily as σEµ/Eµ worsens. As expected, ∆χ

    2 is seen to increase with thereconstruction efficiency #. For comparison, we show by the black dot-dashed lines in the figure, the ∆χ2

    obtained in [28] using the full detector simulation results for the muon analysis.

    In the right panel of Fig. 5 we show the ∆χ2 as a function of the reconstruction efficiency #, fordifferent values of the muon energy resolution σEµ/Eµ. The ∆χ

    2 is seen to increase linearly with #. Thisis not surprising as in the way we have included it in our analysis, # linearly increases the statistics ofthe experiment. Since the mass hierarchy signal in an ICAL@INO like experiment for a 50× 10 kton-yrsexposure is still in the statistics dominated regime, the ∆χ2 grows linear with more data.

    Finally, in the left panel of Fig. 6, we show contours of constant ∆χ2 obtained in the reconstructionefficiency and energy resolution plane, from the analysis of the muon events in the detector. We show thecontours for ∆χ2 = 6 to 16. The figure shows that for energy resolution in the range of 2% to 3% andreconstruction efficiency above 90%, one would get a 4σ signal for the neutrino mass hierarchy from theanalysis of the muon data alone. We reiterate the the reconstruction efficiency can be compensated byadjusting the exposure of the experiment. On the other hand, a 3σ measurement of the mass hierarchyfrom the muon analysis alone seems to be extremely plausible for reasonable range of values for the muonenergy resolution and reconstruction efficiency.

    19

    10×50 kton-yr data

    Ghosh, SC, arxiv:1306.1423

    For σE/E95%Δχ2>14 (muons)For σE/E90%Δχ2>12 (muons)For σE/E80%Δχ2>8 (muons)

    INO μ angle resoln is goodWe take σcosΘ=0.025

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Neutrinos

    8

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µ!(NH) - Pµ!P eventsµ!Oscillated

    0=10"#=0.1, !/E!E#0=20"#=0.1, !/E!E#0=10"#=0.2, !/E!E#0=20"#=0.2, !/E!E#

    L=7043 Km

    !$cos-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1

    X(N

    H)-X

    (IH

    )

    -0.5

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    0.3 (IH)µ!

    (NH) - Pµ!P eventsµ!Oscillated

    0=10"#=0.1, !/E!E# 0=20"#=0.1, !/E!E# 0=10"#=0.2, !/E!E# 0=20"#=0.2, !/E!E#

    6.25 GeV

    Figure 4: The difference in the predicted hadron event rates for normal and inverted hierarchy, normalizedto the no-oscillation hadron event rate, shown as a function of energy. The left panels shows the energydependence of the hierarchy sensitivity, where the events are generated in the zenith angle bin marked bythe path-length traversed at the mid-point of the bin. The right panel shows the corresponding plot asa function of the neutrino zenith angle in the energy bin corresponding to neutrino energy of 6.25 GeV.The solid black lines show the ∆Pµµ as in Fig. 2. The red dotted lines show the difference in the neutrinoenergy spectrum when no detector resolutions are added. The pink dashed, cyan long-dashed, bluedot-dashed and green dot-dashed lines show how the hierarchy sensitivity in the the neutrino spectrumchanges as we include the detector resolutions in neutrino energy and neutrino zenith angle.

    about any drastic change in the spectral shape. Indeed the spectra for σEH/EH = 0.8/√EH + 0.2 and

    σEH/EH = 0.56/√EH + 0.02 look almost identical to each other. As in Fig. 2, the black solid lines in

    all the panels show the ∆Pµµ at the probability level. Likewise, the black dashed lines correspond todifference in the event spectrum between normal and inverted hierarchies in terms of the neutrino energywith only the disappearance channel. The red dotted lines in all the panels show the event spectrumin terms of the true muon energy, where we have included events from both the disappearance as wellas appearance channels but have not incorporated the muon energy resolutions. These lines again arethe same as the red dotted lines in Fig. 2. Therefore, the black solid lines, black dashed lines andthe red dotted lines and have been shown here again for comparison of the hadron spectrum with thecorresponding neutrino and the muon spectra.

    3.3 Earth Matter Effects in Neutrino Events

    If the data is classified in terms of the neutrino energy and zenith angle, then the event spectrum is givenas

    N ′thi (ν) = N∑

    i

    ∫ Emaxν

    Eminν

    dE′ν

    ∫ cosΘmaxν

    cosΘminν

    d cosΘ′ν REν RΘν

    (

    Eνi Ci ni(ν) + Eνi (1− Ci)ni(ν)

    )

    , (14)

    where all quantities are similar to those defined for the µ− events in Eq. (14) with the only difference thatthey now correspond to the neutrinos rather than to the muons. The charge identification efficiency is the

    15

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Detector resolutions *extremely crucial*Ghosh, SC, arxiv:1306.1423

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Neutrinos

    9

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µ!(NH) - Pµ!P eventsµ!Oscillated

    0=10"#=0.1, !/E!E#0=20"#=0.1, !/E!E#0=10"#=0.2, !/E!E#0=20"#=0.2, !/E!E#

    L=7043 Km

    !$cos-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1

    X(N

    H)-X

    (IH

    )

    -0.5

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    0.3 (IH)µ!

    (NH) - Pµ!P eventsµ!Oscillated

    0=10"#=0.1, !/E!E# 0=20"#=0.1, !/E!E# 0=10"#=0.2, !/E!E# 0=20"#=0.2, !/E!E#

    6.25 GeV

    Figure 4: The difference in the predicted hadron event rates for normal and inverted hierarchy, normalizedto the no-oscillation hadron event rate, shown as a function of energy. The left panels shows the energydependence of the hierarchy sensitivity, where the events are generated in the zenith angle bin marked bythe path-length traversed at the mid-point of the bin. The right panel shows the corresponding plot asa function of the neutrino zenith angle in the energy bin corresponding to neutrino energy of 6.25 GeV.The solid black lines show the ∆Pµµ as in Fig. 2. The red dotted lines show the difference in the neutrinoenergy spectrum when no detector resolutions are added. The pink dashed, cyan long-dashed, bluedot-dashed and green dot-dashed lines show how the hierarchy sensitivity in the the neutrino spectrumchanges as we include the detector resolutions in neutrino energy and neutrino zenith angle.

    about any drastic change in the spectral shape. Indeed the spectra for σEH/EH = 0.8/√EH + 0.2 and

    σEH/EH = 0.56/√EH + 0.02 look almost identical to each other. As in Fig. 2, the black solid lines in

    all the panels show the ∆Pµµ at the probability level. Likewise, the black dashed lines correspond todifference in the event spectrum between normal and inverted hierarchies in terms of the neutrino energywith only the disappearance channel. The red dotted lines in all the panels show the event spectrumin terms of the true muon energy, where we have included events from both the disappearance as wellas appearance channels but have not incorporated the muon energy resolutions. These lines again arethe same as the red dotted lines in Fig. 2. Therefore, the black solid lines, black dashed lines andthe red dotted lines and have been shown here again for comparison of the hadron spectrum with thecorresponding neutrino and the muon spectra.

    3.3 Earth Matter Effects in Neutrino Events

    If the data is classified in terms of the neutrino energy and zenith angle, then the event spectrum is givenas

    N ′thi (ν) = N∑

    i

    ∫ Emaxν

    Eminν

    dE′ν

    ∫ cosΘmaxν

    cosΘminν

    d cosΘ′ν REν RΘν

    (

    Eνi Ci ni(ν) + Eνi (1− Ci)ni(ν)

    )

    , (14)

    where all quantities are similar to those defined for the µ− events in Eq. (14) with the only difference thatthey now correspond to the neutrinos rather than to the muons. The charge identification efficiency is the

    15

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Net matter effects in neutrinos and muons is the same a+er putting detector resolutions - even for very optimistic neutrino resolutions

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=11019 km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP only disapp channelµ!Oscillated

    - only disapp channelµOscillated - eventsµOscillated

    =0.15µ/EµE

    "- events µSmeared =0.02µ/E

    µE"- events µSmeared

    L=12611 Km

    Figure 2: The difference in the predicted muon event rates for normal and inverted hierarchy, normalizedto the no-oscillation muon event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point of thebin. The solid black lines show the ∆Pµµ as in Fig. 1. The black dashed lines give the difference betweenthe predicted rate for normal and inverted hierarchy when only the Pµµ (disappearance) channel is takenand data binned in neutrino energy. The blue dot-dashed lines show the corresponding difference whenonly the Pµµ channel is taken and data binned in muon energy. The red dotted lines are obtained whenwe add the Peµ (appearance) channel and bin the data in muon energy. The pink dashed and cyan longdashed lines are obtained when we apply on the red dotted lines the muon energy resolution functionswith widths σEµ/Eµ = 0.15% and 2%, respectively.

    9

    Ghosh, SC, arxiv:1306.1423

    MuonsNeutrinos

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Neutrinos

    10

    !/E!E"0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24

    2 #$

    7

    8

    9

    10

    11

    12

    13

    14 energyµ! zenith angle, 20 µBinned in 80

    energyµ! zenith angle, 20 µ!Binned in 20

    Reconstruction efficiency 80%

    Charge id efficiency 99%

    Figure 7: The ∆χ2 corresponding to the mass hierarchy sensitivity as a function of the width of theneutrino resolution function. The black line shows the expected sensitivity when we use the binning inneutrino zenith angles with the neutrino zenith angle resolution function taken from the MONOLITHproposal [36]. The red dashed line shows the expected sensitivity when we use the binning in muon zenithangles with the muon zenith angle resolution width of 0.01 in cosΘµ. The black dot-dashed lines in thefigure show the ∆χ2 obtained from the muon analysis performed in [28], using the detector response tomuons obtained by the INO collaboration.

    just the muon data with 80% efficiency and energy resolution of 12%. If the hadron data is added to themuon data, then even with these modest detector response for the muons, we would get ∆χ2 ! 12− 13(cf. Table 3). If the muon energy resolution was improved to 2%, then we would have a more than 4σsensitivity to the mass hierarchy from the muon-plus-hadron analysis. Whereas, from the Fig. 4 we cansee that the sensitivity from the neutrino analysis can never match these numbers, even for extremelyoptimistic energy resolution of 10%. The reason can be traced to the event rates plots shown in Figs.2, 3 and 4. The effect of the detector resolutions on the neutrino spectrum make it comparable, andsometimes even worse than the muon spectrum. If we add the hadron spectral data to the muon spectralanalysis, we get an additional contribution to the hierarchy sensitivity which cannot be matched by theneutrino analysis.

    For comparison, in Fig. 7 we also show the ∆χ2 obtained when we bin the data in neutrino energyand muon zenith angle bins. We keep the reconstruction efficiency at 80% and charge identification effi-ciency at 99%. The muon zenith angle resolution function is taken with width 0.01. We see that the ∆χ2

    obtained in this case is even worse than what we had obtained with the neutrino zenith angle analysisusing the MONOLITH zenith angle resolution function.

    24

    using Θ-resoln ,om Monolith Δχ2=9 for 21% E-resolnwhich is the E-resoln for dets like INO/MIND- this is the same as you get ,om a *very* conservative muon analysis!!

    Ghosh, SC, arxiv:1306.1423

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Neutrinos

    11

    !/E!E"0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24

    2 #$

    7

    8

    9

    10

    11

    12

    13

    14 energyµ! zenith angle, 20 µBinned in 80

    energyµ! zenith angle, 20 µ!Binned in 20

    Reconstruction efficiency 80%

    Charge id efficiency 99%

    Figure 7: The ∆χ2 corresponding to the mass hierarchy sensitivity as a function of the width of theneutrino resolution function. The black line shows the expected sensitivity when we use the binning inneutrino zenith angles with the neutrino zenith angle resolution function taken from the MONOLITHproposal [36]. The red dashed line shows the expected sensitivity when we use the binning in muon zenithangles with the muon zenith angle resolution width of 0.01 in cosΘµ. The black dot-dashed lines in thefigure show the ∆χ2 obtained from the muon analysis performed in [28], using the detector response tomuons obtained by the INO collaboration.

    just the muon data with 80% efficiency and energy resolution of 12%. If the hadron data is added to themuon data, then even with these modest detector response for the muons, we would get ∆χ2 ! 12− 13(cf. Table 3). If the muon energy resolution was improved to 2%, then we would have a more than 4σsensitivity to the mass hierarchy from the muon-plus-hadron analysis. Whereas, from the Fig. 4 we cansee that the sensitivity from the neutrino analysis can never match these numbers, even for extremelyoptimistic energy resolution of 10%. The reason can be traced to the event rates plots shown in Figs.2, 3 and 4. The effect of the detector resolutions on the neutrino spectrum make it comparable, andsometimes even worse than the muon spectrum. If we add the hadron spectral data to the muon spectralanalysis, we get an additional contribution to the hierarchy sensitivity which cannot be matched by theneutrino analysis.

    For comparison, in Fig. 7 we also show the ∆χ2 obtained when we bin the data in neutrino energyand muon zenith angle bins. We keep the reconstruction efficiency at 80% and charge identification effi-ciency at 99%. The muon zenith angle resolution function is taken with width 0.01. We see that the ∆χ2

    obtained in this case is even worse than what we had obtained with the neutrino zenith angle analysisusing the MONOLITH zenith angle resolution function.

    24

    using Θ-resoln ,om Monolith Δχ2=9 for 21% E-resolnwhich is the E-resoln for dets like INO/MIND- this is the same as you get ,om a *very* conservative muon analysis!!

    Ghosh, SC, arxiv:1306.1423

    Using the hadron data to reconstruct the neutrino is not helping much

  • Sandhya Choubey NuFact 2013, August 21

    Matter Effects in Hadrons

    12

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013 Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Though significantly less than in muons

    Sandhya Choubey NuFact 2013, August 21

    The neutrino mass ordering (MH)…

    CP violation in the lepton sector…

    Octant of the mixing angle !23 …

    Beyond 3-flavor osci%ation physics..

    Monday 12 August 2013

    Detector resolution less important for hadrons

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.4

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2 (IH)µµ(NH)-PµµP events only disapp channelµ!Oscillated

    Oscillated hadrons events - eventsµOscillated

    +0.2 HE=0.8/H/EHE"Smeared hadrons with +0.02 HE=0.56/H/EHE"Smeared hadrons with

    L=7202 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.3

    -0.2

    -0.1

    0

    0.1

    0.2

    (IH)µµ(NH)-PµµP events only disapp channelµ!Oscillated

    Oscillated hadrons events - eventsµOscillated

    +0.2 HE=0.8/H/EHE"Smeared hadrons with +0.02 HE=0.56/H/EHE"Smeared hadrons with

    L=9110 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    (IH)µµ(NH)-PµµP events only disapp channelµ!Oscillated

    Oscillated hadrons events - eventsµOscillated

    +0.2 HE=0.8/H/EHE"Smeared hadrons with +0.02 HE=0.56/H/EHE"Smeared hadrons with

    L=11019 Km

    Energy[GeV]2 3 4 5 6 7 8 9 10

    X(N

    H)-X

    (IH

    )

    -1

    -0.8

    -0.6

    -0.4

    -0.2

    0

    0.2

    0.4

    0.6(IH)µµ(NH)-PµµP

    events only disapp channelµ!Oscillated

    Oscillated hadrons events - eventsµOscillated

    +0.2 HE=0.8/H/EHE"Smeared hadrons with +0.02 HE=0.56/H/EHE"Smeared hadrons with

    L=12611 Km

    Figure 3: The difference in the predicted hadron event rates for normal and inverted hierarchy, normalizedto the no-oscillation hadron event rate, shown as a function of energy. The four panels show the eventspectrum generated in four zenith angle bins, marked by the path-length traversed at the mid-point ofthe bin. The solid black lines, black dashed lines and red dotted lines show the things as in Fig. 2. Theblack dashed lines give the difference between the predicted rate for normal and inverted hierarchy whenonly the Pµµ (disappearance) channel is taken and data binned in neutrino energy. The blue dot-dashedlines show the corresponding difference when only the Pµµ channel is taken and data binned in muonenergy. The red dotted lines are obtained when we add the Peµ (appearance) channel and bin the datain muon energy. The pink dashed and cyan long dashed lines are obtained when we apply on the reddotted lines the muon energy resolution functions with widths σEµ/Eµ = 0.15% and 2%, respectively.

    14

    Ghosh, SC, arxiv:1306.1423Matter effects in hadrons present...wi% give MH sensitivity

    Hadrons are ta.ed with the muon direction

  • Sandhya Choubey NuFact 2013, August 21

    MH@INO(Muon+Hadron)

    13

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =62"# =72"# =82"# =92"#=102"# =112"# =122"# =132"#=142"# =152"# =162"#

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =92"# =102"# =112"# =122"#=132"# =142"# =152"# =162"#=172"# =182"# =192"# =202"#

    Figure 6: Constant ∆χ2 contours in the reconstruction efficiency and energy resolution plane. The leftpanels shows the contours from the analysis which uses only the muon data from the experiments. Theright panels shows the contours for the analysis in which both the muon and hadron data were includedin the combined statistical analysis.

    present results for both 80 bins as well as 20 bins between 1 GeV and 11 GeV.

    The dependence of the mass hierarchy sensitivity to the muon energy resolution and muon recon-struction efficiency are shown in Figs. 5. In the left panel on Fig. 5 we show the ∆χ2 as a function ofthe muon energy resolution σEµ/Eµ, for different values of the muon reconstruction efficiency #. We varyσEµ/Eµ continuously between 20% and 2% and repeat this for # = 70% to 100%. For a given #, we seethat the ∆χ2 falls sharply as the σEµ/Eµ is increased from 2% to about 6%. Thereafter, the rate of fallof ∆χ2 reduces, and it falls steadily as σEµ/Eµ worsens. As expected, ∆χ

    2 is seen to increase with thereconstruction efficiency #. For comparison, we show by the black dot-dashed lines in the figure, the ∆χ2

    obtained in [28] using the full detector simulation results for the muon analysis.

    In the right panel of Fig. 5 we show the ∆χ2 as a function of the reconstruction efficiency #, fordifferent values of the muon energy resolution σEµ/Eµ. The ∆χ

    2 is seen to increase linearly with #. Thisis not surprising as in the way we have included it in our analysis, # linearly increases the statistics ofthe experiment. Since the mass hierarchy signal in an ICAL@INO like experiment for a 50× 10 kton-yrsexposure is still in the statistics dominated regime, the ∆χ2 grows linear with more data.

    Finally, in the left panel of Fig. 6, we show contours of constant ∆χ2 obtained in the reconstructionefficiency and energy resolution plane, from the analysis of the muon events in the detector. We show thecontours for ∆χ2 = 6 to 16. The figure shows that for energy resolution in the range of 2% to 3% andreconstruction efficiency above 90%, one would get a 4σ signal for the neutrino mass hierarchy from theanalysis of the muon data alone. We reiterate the the reconstruction efficiency can be compensated byadjusting the exposure of the experiment. On the other hand, a 3σ measurement of the mass hierarchyfrom the muon analysis alone seems to be extremely plausible for reasonable range of values for the muonenergy resolution and reconstruction efficiency.

    19

    10 yrs data

    For σE/E95%Δχ2>20 (muons+hadrons)Δχ2>14 (muons)

    Ghosh, SC, arxiv:1306.1423

    For σE/E90%Δχ2>18 (muons+hadrons)Δχ2>12 (muons)

    For σE/E85%Δχ2>16 (muons+hadrons)Δχ2>10 (muons)

    Hadrons are ta.ed with the muon direction

  • Sandhya Choubey NuFact 2013, August 21

    MH@INO(Muon+Hadron)

    14

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =62"# =72"# =82"# =92"#=102"# =112"# =122"# =132"#=142"# =152"# =162"#

    Reconstruction efficiency0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

    µ/EµE

    !

    0.02

    0.04

    0.06

    0.08

    0.1

    0.12

    0.14

    0.16

    0.18

    0.2

    0.22 =92"# =102"# =112"# =122"#=132"# =142"# =152"# =162"#=172"# =182"# =192"# =202"#

    Figure 6: Constant ∆χ2 contours in the reconstruction efficiency and energy resolution plane. The leftpanels shows the contours from the analysis which uses only the muon data from the experiments. Theright panels shows the contours for the analysis in which both the muon and hadron data were includedin the combined statistical analysis.

    present results for both 80 bins as well as 20 bins between 1 GeV and 11 GeV.

    The dependence of the mass hierarchy sensitivity to the muon energy resolution and muon recon-struction efficiency are shown in Figs. 5. In the left panel on Fig. 5 we show the ∆χ2 as a function ofthe muon energy resolution σEµ/Eµ, for different values of the muon reconstruction efficiency #. We varyσEµ/Eµ continuously between 20% and 2% and repeat this for # = 70% to 100%. For a given #, we seethat the ∆χ2 falls sharply as the σEµ/Eµ is increased from 2% to about 6%. Thereafter, the rate of fallof ∆χ2 reduces, and it falls steadily as σEµ/Eµ worsens. As expected, ∆χ

    2 is seen to increase with thereconstruction efficiency #. For comparison, we show by the black dot-dashed lines in the figure, the ∆χ2

    obtained in [28] using the full detector simulation results for the muon analysis.

    In the right panel of Fig. 5 we show the ∆χ2 as a function of the reconstruction efficiency #, fordifferent values of the muon energy resolution σEµ/Eµ. The ∆χ

    2 is seen to increase linearly with #. Thisis not surprising as in the way we have included it in our analysis, # linearly increases the statistics ofthe experiment. Since the mass hierarchy signal in an ICAL@INO like experiment for a 50× 10 kton-yrsexposure is still in the statistics dominated regime, the ∆χ2 grows linear with more data.

    Finally, in the left panel of Fig. 6, we show contours of constant ∆χ2 obtained in the reconstructionefficiency and energy resolution plane, from the analysis of the muon events in the detector. We show thecontours for ∆χ2 = 6 to 16. The figure shows that for energy resolution in the range of 2% to 3% andreconstruction efficiency above 90%, one would get a 4σ signal for the neutrino mass hierarchy from theanalysis of the muon data alone. We reiterate the the reconstruction efficiency can be compensated byadjusting the exposure of the experiment. On the other hand, a 3σ measurement of the mass hierarchyfrom the muon analysis alone seems to be extremely plausible for reasonable range of values for the muonenergy resolution and reconstruction efficiency.

    19

    10 yrs data

    For σE/E95%Δχ2>20 (muons+hadrons)Δχ2>14 (muons)

    Ghosh, SC, arxiv:1306.1423

    For σE/E90%Δχ2>18 (muons+hadrons)Δχ2>12 (muons)

    For σE/E85%Δχ2>16 (muons+hadrons)Δχ2>10 (muons)

    Hadrons are ta.ed with the muon direction

    One can get Δχ2=9 with 5 years INO data

  • Sandhya Choubey NuFact 2013, August 21

    MH - μ+h vs ν Analysis

    15

    Muon+Hadron analysis gives Δχ2=16 for σE(μ)=6%E and σΘ(μ)=0.025 in cosΘ .... both are realistic for INO/MIND

    To get the same Δχ2=16 using the neutrino analysis, one needs σE(ν)=10%E and σΘ(ν)=7.3o....the energy resol requirement not possible at INO/MIND

    Better to use the hadron data separately in the χ2 than using it to reconstruct the neutrino energy and angle

  • Sandhya Choubey NuFact 2013, August 21

    Octant @ INO

    16

    (True)232sin

    0.35 0.4 0.45 0.5 0.55 0.6 0.65

    2

    0

    5

    10

    15

    20

    25

    30NORMAL HIERARCHY

    MarginalizedAll exptLBL+reactorINO

    All exptLBL+reactorINO For sin2θ232.44 (INO)Δχ2>9 (LBL+React)Δχ2>12.3 (INO+LBL+React)

    Sensitivity for IH is worse

    INO Co%ab 2013

  • Sandhya Choubey NuFact 2013, August 21

    Conclusions

    17

    With 5% Eμ-resoln and 90% efficiency, INO wi% have Δχ2=12 in 10 years ,om muon data alone

    If we tag the hadrons with the muons and add the hadron data separately into a combined analysis, then with 5% Eμ-resoln and 90% reconstruction efficiency INO could have Δχ2=9 in just 5 years

    Using the hadrons to reconstruct the neutrino energy and directions wi% always yield poorer results