Tokai to Kamioka (T2K) is a long-baseline neutrino oscillation experiment that measures oscillation parameters related to both $ν_μ(\barν_μ)$ disappearance and $ν_e(\barν_e)$ appearance in a $ν_μ(\barν_μ)$ beam. T2K uses Super-Kamiokande (SK) as its far detector, and SK detector systematic errors are currently among the leading sources of systematic uncertainty in the T2K oscillation analysis. Therefore, accurate understanding of detector mis-modelling and event reconstruction uncertainties in SK is crucial to the extraction of the neutrino oscillation parameters. The detector error estimation procedure quantifies uncertainties by fitting SK atmospheric MC to data in a Markov Chain Monte Carlo (MCMC) framework. Uncertainties are separated based on true event topology above Cherenkov threshold in the SK atmospheric MC in order to capture how the uncertainties should affect the different signal and background samples when propagated to the T2K beam oscillation analysis. The procedure has been upgraded for the upcoming analysis cycle to include systematics targeting newly added $ν_eCC1Ï^\pm$ and $NCÏ^0$-enhanced analysis samples, among other changes, culminating in a 540 dimensional MCMC fit between atmospheric MC and data in SK. The far detector analysis chain and its integration into the broader T2K oscillation analysis will be discussed in these proceedings.
- Papier-ID: 2501.00147
- Titel: Estimating detector systematic uncertainties for the T2K far detector
- Autor: Michael Reh (im Namen der T2K-Kollaboration)
- Klassifizierung: hep-ex (Hochenergiephysik-Experiment)
- Veröffentlichungszeit/Konferenz: 25. Internationaler Workshop zu Neutrinos aus Beschleunigern (25th International Workshop on Neutrinos from Accelerators)
- Papierlink: https://arxiv.org/abs/2501.00147
T2K ist ein Langbaseline-Neutrinooszillationsexperiment zur Messung von Oszillationsparametern, die mit dem νμ(ν̄μ)-Verschwinden und νe(ν̄e)-Erscheinen verbunden sind. T2K nutzt Super-Kamiokande (SK) als Fernedetektor, wobei die Detektorsystematiken des SK-Detektors eine der Hauptquellen für systematische Unsicherheiten in der T2K-Oszillationsanalyse darstellen. Das Detektorfehler-Abschätzungsverfahren quantifiziert Unsicherheiten durch Anpassung von SK-Atmosphärenneutrino-Monte-Carlo-Simulationen an Daten im Rahmen einer Markov-Ketten-Monte-Carlo-(MCMC-)Analyse. Das Verfahren wurde auf eine 540-dimensionale MCMC-Anpassung erweitert, die die Behandlung von Systematiken für neu hinzugefügte νeCC1π±- und NCπ0-verstärkte Analysestichproben umfasst.
Das T2K-Experiment erfordert eine präzise Messung von Neutrinooszillationsparametern, aber die Detektorsystematiken des Super-Kamiokande-Ferndetektors sind zu einem der Hauptfaktoren geworden, die die Messgenauigkeit begrenzen. Fehler bei der Detektormodellierung und Unsicherheiten bei der Ereignisrekonstruktion beeinflussen direkt die Extraktionsgenauigkeit von Neutrinooszillationsparametern.
- Große physikalische Bedeutung: Die präzise Messung von Neutrinooszillationsparametern ist von großer Bedeutung für das Verständnis der fundamentalen Teilchenphysik und der kosmischen Entwicklung
- Systematiken dominieren: SK-Detektorsystematiken sind die Hauptquelle für Unsicherheiten in der aktuellen T2K-Analyse
- Präzisionsanforderungen: Mit verbesserter statistischer Genauigkeit wird die Kontrolle systematischer Unsicherheiten zunehmend kritisch
Traditionelle Methoden zur Abschätzung von Detektorsystematiken basieren häufig auf spezifischen physikalischen Modellannahmen, mangelt es an ausreichender Allgemeingültigkeit und Robustheit, und es ist schwierig, sich an Änderungen der Detektorbedingungen und Verbesserungen der Modellierungswahl anzupassen.
Entwicklung einer modellunabhängigen, robusten Methode zur Abschätzung von Detektorsystematiken, die die Modellierungsunsicherheiten des SK-Detektors genau quantifizieren kann und diese effektiv auf die T2K-Strahl-Oszillationsanalyse überträgt.
- Etablierung eines MCMC-basierten Rahmens zur Abschätzung von Detektorsystematiken, unter Verwendung von Atmosphärenneutrinodaten als Kontrollstichprobe
- Vorschlag einer universellen Parametrisierungsmethode, die Verschiebungs-(β) und Verbreiterungs-(α)-Parameter zur Beschreibung systematischer Versätze von Teilchenidentifikationsvariablen verwendet
- Entwicklung eines 540-dimensionalen MCMC-Anpassungsprogramms, das neu hinzugefügte νeCC1π±- und NCπ0-Analysestichproben unterstützt
- Realisierung einer signifikanten Verbesserung der Recheneffizienz, Verkürzung der Analyskonvergenzzeit von Monatsebene auf Wochenebene
- Aufbau einer vollständigen Fehlerübertragungskette, von der Atmosphärenneutrino-Anpassung zur Kovarianzmatrix der T2K-Strahl-Analyse
Eingabe: SK-Atmosphärenneutrinodaten und Monte-Carlo-Simulationsereignisse
Ausgabe: Detektorsystematik-Kovarianzmatrix für verschiedene Energiebereiche und Stichprobentypen in der T2K-Strahl-Analyse
Einschränkungen: Beibehaltung der physikalischen Plausibilität, Sicherstellung der statistischen Signifikanz der Fehlerabschätzung
T2K verwendet eine universelle Parametrisierungsmethode zur Beschreibung von Detektormodellierungsfehlern:
wobei:
- L: Nominalwert der Teilchenidentifikations-(PID-)Variable
- α: Verbreiterungsparameter, kontrolliert die Breitenänderung der Verteilung
- β: Verschiebungsparameter, kontrolliert die Verschiebung der Verteilungsmitte
Parameterabschätzung unter Verwendung des MCMC-Rahmens:
Schritte:
- MCMC schlägt {α, β}-Parameterkombinationen vor
- Modifizierung der PID-Werte von MC-Ereignissen nach der Formel
- Konstruktion von Histogrammen aus modifizierten MC-Ereignissen in Bins
- Berechnung der Formähnlichkeitsfunktion zwischen Daten und modifiziertem MC
- Wiederholung des Prozesses zur Erlangung der posterioren Verteilung von {α, β}
Übertragung der Atmosphärenneutrino-Anpassungsergebnisse auf die T2K-Strahl-Analyse durch ein Spielzeug-MC-Programm:
- Zufällige Stichprobenentnahme von {α, β}-Parametern aus der MCMC-posterioren Verteilung
- Modifizierung der PID-Werte von T2K-Strahl-MC-Ereignissen
- Beobachtung der Ereignismigration zwischen verschiedenen Analysestichproben
- Statistische Erfassung von Änderungen im Neutrinoenergiespektrum
- Konstruktion der Kovarianzmatrix
Diese Methode hängt nicht von spezifischen physikalischen Modellannahmen ab und kann verschiedene potenzielle Detektormodellierungsfehler erfassen, einschließlich:
- Ungenauigkeiten bei der Modellierung der Lichtstreuintensität
- Abweichungen in den Photomultiplier-Reaktionseigenschaften
- Änderungen der optischen Eigenschaften des Wassers
Trennung der Unsicherheiten basierend auf echter Ereignistopologie oberhalb der Cherenkov-Schwelle, um sicherzustellen, dass Systematiken die verschiedenen Signal- und Untergrundstichproben korrekt beeinflussen.
Die erweiterte 540-dimensionale MCMC-Anpassung ermöglicht eine feinere Beschreibung der systematischen Fehlermerkmale für verschiedene Ereignistypen und kinematische Bereiche.
- SK-Atmosphärenneutrinodaten: Etwa 30 Jahre Atmosphärenneutrino-Beobachtungsdaten, die den für T2K relevanten Energiebereich und Ereignistypen abdecken
- T2K-Strahldaten: Run-1-4-Datensätze, einschließlich νμ- und ν̄μ-Strahlmodi
- Monte-Carlo-Simulation: Neue MC-Produkte mit aktualisierten Neutrinointeraktionsmodellen
- Formähnlichkeitsfunktion: Bewertung der Konsistenz zwischen Daten- und MC-Verteilung
- Konvergenz der posterioren Verteilung: Sicherstellung der Zuverlässigkeit der Parameterabschätzung durch MCMC-Diagnose
- Physikalische Plausibilität der Kovarianzmatrix: Überprüfung der Angemessenheit von Fehlerkorrelationen
- Einfach-Ring-Elektron-Stichprobe (Single ring e-like)
- Einfach-Ring-Myon-Stichprobe (Single ring μ-like)
- νeCC1π+-Stichprobe (neu hinzugefügt)
- NCπ0-verstärkte Stichprobe (neu hinzugefügt)
Erfolgreiche Quantifizierung systematischer Versätze kritischer PID-Variablen wie Elektron-/Myonidentifikation durch {α, β}-Parameter, wie in Abbildung 2 die Auswirkungen verschiedener Parameterkombinationen auf die Verteilungsform zeigt.
Abbildung 3 zeigt die posterioren Vorhersagespektren für verschiedene T2K-Analysestichproben, wobei Bin-Fehler die 1σ-Unsicherheiten aus der SK-Detektorfehlermatrix widerspiegeln. Von der 2022er bis zur 2023er Analyse zeigen sich erhebliche Unterschiede in den Bin-Fehlern durch die Entfernung mehrerer "vorläufiger" Systematiken.
Realisierung signifikanter Leistungsverbesserungen:
- Konvergenzzeit: Verkürzung von O(Monate) auf O(Wochen)
- Parameterdimension: Erweiterung auf 540-dimensionale MCMC-Anpassung
- Stichprobenabdeckung: Unterstützung der Systematik-Abschätzung für neu hinzugefügte Analysestichproben
Durch Vergleich der Überlappung von Atmosphärenneutrinos und Strahnneutrinos in kinematischen Bereichen, Ereignistopologien und Detektorbedingungen wird die Angemessenheit der Anwendung von Atmosphärenneutrino-Anpassungsergebnissen auf die T2K-Strahl-Analyse validiert.
Langbaseline-Neutrinooszillationsexperimente stehen allgemein vor der Herausforderung von Detektorsystematiken, wobei Experimente wie NOνA und DUNE entsprechende Fehlerabschätzungsmethoden entwickeln.
Die präzise Modellierung von Wasser-Cherenkov-Detektoren ist ein wichtiges Thema in diesem Bereich und umfasst mehrere Aspekte wie optische Eigenschaften, PMT-Reaktion und Ereignisrekonstruktionsalgorithmen.
Im Vergleich zu anderen Methoden bietet der vorgeschlagene Rahmen:
- Stärkere Allgemeingültigkeit: Modellunabhängige Parametrisierungsmethode
- Höhere Präzision: Feine Beschreibung im 540-dimensionalen Parameterraum
- Bessere Robustheit: Fähigkeit, sich an Änderungen der Detektorbedingungen anzupassen
- Erfolgreiche Etablierung eines MCMC-basierten Rahmens zur Abschätzung von SK-Detektorsystematiken
- Realisierung einer effektiven Fehlerübertragung von Atmosphärenneutrinodaten zur T2K-Strahl-Analyse
- Bereitstellung eines wichtigen Mittels zur Kontrolle von Systematiken für bevorstehende T2K-Analysen
- Annahmebeschränkungen: Annahme, dass Atmosphärenneutrinos und Strahnneutrinos ähnliche Systematik-Merkmale aufweisen
- Modellabhängigkeit: Abhängigkeit von der grundlegenden Genauigkeit der MC-Simulation
- Rechenkomplexität: 540-dimensionale MCMC-Anpassung erfordert immer noch erhebliche Rechenressourcen
- Weitere Optimierung des MCMC-Algorithmus zur Verbesserung der Konvergenzeffizienz
- Erkundung feinerer Ereignisklassifizierung und Systematik-Modellierung
- Kombination mit anderen Kontrollstichproben zur Validierung der Genauigkeit der Fehlerabschätzung
- Starke Methodische Innovation: Vorschlag einer modellunabhängigen Parametrisierungsmethode für Systematiken
- Angemessenes Experimentaldesign: Verwendung von Atmosphärenneutrinodaten als Kontrollstichprobe hat starke physikalische Motivation
- Vollständige technische Implementierung: 540-dimensionale MCMC-Anpassung demonstriert hohes technisches Niveau
- Hoher praktischer Wert: Direkte Unterstützung der Genauigkeitsverbesserung der T2K-Physikanalyse
- Begrenzte theoretische Analyse: Mangelnde tiefgehende Diskussion der theoretischen Grundlagen der Parametrisierungsmethode
- Einseitige Validierungsmittel: Hauptsächlich abhängig von Atmosphärenneutrinodaten, mangelnde unabhängige Validierung
- Unzureichende Quantifizierung von Unsicherheiten: Unzureichende Diskussion der durch die Methode selbst eingeführten Unsicherheiten
- Akademischer Beitrag: Bereitstellung einer wichtigen Methode zur Kontrolle von Systematiken für Neutrinophysik-Experimente
- Praktischer Wert: Direkte Auswirkung auf die Genauigkeit der physikalischen Ergebnisse des T2K-Experiments
- Technische Verbreitung: Methode kann auf andere ähnliche Teilchenphysik-Experimente übertragen werden
Diese Methode ist besonders geeignet für:
- Langbaseline-Neutrinooszillationsexperimente
- Große Wasser-Cherenkov-Detektoren
- Teilchenphysik-Messungen, die hochpräzise Systematik-Kontrolle erfordern
1 K. A. et al., The T2K experiment, NIM-A 659, 106 (2011)
2 S. F. et al., The Super-Kamiokande detector, NIM-A 501, 418 (2003)
3 K. A. et al., Measurements of neutrino oscillation parameters from the T2K experiment using 3.6×10²¹ protons on target, Eur. Phys. J. C 83, 782 (2023)
Gesamtbewertung: Dies ist ein technisch anspruchsvolles experimentelles Physik-Papier, das einen wichtigen methodischen Beitrag zur Systematik-Kontrolle des T2K-Experiments leistet. Obwohl es in der theoretischen Analyse gewisse Mängel aufweist, machen sein praktischer Wert und seine technische Innovativität es zu einer wichtigen Arbeit im Bereich der Neutrinophysik.