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.
- ID de l'article: 2501.00147
- Titre: Estimation des incertitudes systématiques du détecteur pour le détecteur lointain de T2K
- Auteur: Michael Reh (représentant la collaboration T2K)
- Classification: hep-ex (physique expérimentale des hautes énergies)
- Date de publication/Conférence: 25e Atelier international sur les neutrinos provenant d'accélérateurs (25th International Workshop on Neutrinos from Accelerators)
- Lien de l'article: https://arxiv.org/abs/2501.00147
T2K est une expérience de long parcours mesurant les paramètres d'oscillation des neutrinos liés à la disparition de νμ(ν̄μ) et à l'apparition de νe(ν̄e). T2K utilise Super-Kamiokande (SK) comme détecteur lointain, et les erreurs systématiques du détecteur SK constituent l'une des principales sources d'incertitude systématique dans l'analyse d'oscillation de T2K. La procédure d'estimation des erreurs du détecteur quantifie les incertitudes en ajustant les simulations Monte Carlo des neutrinos atmosphériques de SK aux données dans le cadre d'une chaîne de Markov Monte Carlo (MCMC). Cette procédure a été améliorée en un ajustement MCMC à 540 dimensions, incluant le traitement des erreurs systématiques pour les nouveaux échantillons d'analyse νeCC1π± et NCπ0 améliorés.
L'expérience T2K nécessite une mesure précise des paramètres d'oscillation des neutrinos, mais les erreurs systématiques du détecteur lointain Super-Kamiokande sont devenues l'un des principaux facteurs limitant la précision des mesures. Les erreurs de modélisation du détecteur et les incertitudes de reconstruction des événements affectent directement la précision de l'extraction des paramètres d'oscillation des neutrinos.
- Signification physique majeure: La mesure précise des paramètres d'oscillation des neutrinos est cruciale pour comprendre la physique des particules fondamentales et l'évolution cosmique
- Domination des erreurs systématiques: Les erreurs systématiques du détecteur SK constituent la principale source d'incertitude dans l'analyse T2K actuelle
- Exigences de précision: Avec l'amélioration de la précision statistique, le contrôle des incertitudes systématiques devient plus critique
Les méthodes traditionnelles d'estimation des erreurs systématiques du détecteur dépendent souvent d'hypothèses de modèles physiques spécifiques, manquent de généralité et de robustesse suffisantes, et ont du mal à s'adapter aux changements des conditions du détecteur et aux améliorations des choix de modélisation.
Développer une méthode d'estimation des erreurs systématiques du détecteur indépendante du modèle et robuste, capable de quantifier avec précision les incertitudes de modélisation du détecteur SK et de les transférer efficacement à l'analyse d'oscillation du faisceau T2K.
- Établissement d'un cadre d'estimation des erreurs systématiques du détecteur basé sur MCMC, utilisant les données de neutrinos atmosphériques comme échantillon de contrôle
- Proposition d'une méthode de paramétrisation générale, utilisant des paramètres de translation (β) et d'élargissement (α) pour décrire les décalages systématiques des variables d'identification des particules
- Développement d'un programme d'ajustement MCMC à 540 dimensions, supportant les nouveaux échantillons d'analyse νeCC1π± et NCπ0
- Réalisation d'une amélioration significative de l'efficacité computationnelle, réduisant le temps de convergence de l'analyse du niveau mensuel au niveau hebdomadaire
- Construction d'une chaîne complète de propagation des erreurs, de l'ajustement des neutrinos atmosphériques à la matrice de covariance de l'analyse du faisceau T2K
Entrée: Données de neutrinos atmosphériques de SK et événements de simulation Monte Carlo
Sortie: Matrice de covariance des erreurs systématiques du détecteur pour différents intervalles d'énergie et types d'échantillons dans l'analyse du faisceau T2K
Contraintes: Maintenir la rationalité physique et assurer la signification statistique de l'estimation des erreurs
T2K adopte une méthode de paramétrisation générale pour décrire les erreurs de modélisation du détecteur:
Où:
- L: valeur nominale de la variable d'identification des particules (PID)
- α: paramètre d'élargissement, contrôlant la variation de la largeur de la distribution
- β: paramètre de translation, contrôlant le décalage de la position centrale de la distribution
Estimation des paramètres utilisant le cadre MCMC:
Étapes:
- MCMC propose des combinaisons de paramètres {α, β}
- Modification des valeurs PID des événements MC selon la formule
- Construction d'histogrammes des événements MC modifiés par bin
- Calcul de la fonction de vraisemblance de forme entre les données et le MC modifié
- Répétition du processus pour obtenir la distribution a posteriori de {α, β}
Transfert des résultats d'ajustement des neutrinos atmosphériques à l'analyse du faisceau T2K via un programme MC jouet:
- Extraction aléatoire de paramètres {α, β} de la distribution a posteriori MCMC
- Modification des valeurs PID des événements MC du faisceau T2K
- Observation de la migration des événements entre différents échantillons d'analyse
- Statistiques des variations du spectre énergétique des neutrinos
- Construction de la matrice de covariance
Cette méthode ne dépend pas d'hypothèses de modèles physiques spécifiques et peut capturer diverses erreurs potentielles de modélisation du détecteur, notamment:
- Imprécision de la modélisation de l'intensité de la diffusion lumineuse
- Écarts dans les caractéristiques de réponse des photomultiplicateurs
- Variations des propriétés optiques de l'eau
Séparation des incertitudes basée sur la topologie réelle des événements au-dessus du seuil Cherenkov, garantissant que les erreurs systématiques affectent correctement les différents échantillons de signal et de bruit de fond.
L'ajustement MCMC amélioré à 540 dimensions peut décrire plus finement les caractéristiques des erreurs systématiques pour différents types d'événements et régions cinématiques.
- Données de neutrinos atmosphériques SK: Données d'observation de neutrinos atmosphériques d'environ 30 ans, couvrant la gamme d'énergie et les types d'événements pertinents pour T2K
- Données du faisceau T2K: Ensemble de données Run 1-4, incluant les modes de faisceau νμ et ν̄μ
- Simulation Monte Carlo: Nouveau produit MC utilisant des modèles d'interaction de neutrinos mis à jour
- Fonction de vraisemblance de forme: Évaluation de la cohérence entre la distribution des données et celle du MC
- Convergence de la distribution a posteriori: Assurance de la fiabilité de l'estimation des paramètres via des diagnostics MCMC
- Rationalité physique de la matrice de covariance: Vérification de la rationalité des corrélations d'erreur
- Échantillon d'électron monoanneau (Single ring e-like)
- Échantillon de muon monoanneau (Single ring μ-like)
- Échantillon νeCC1π+ (nouveau)
- Échantillon NCπ0 amélioré (nouveau)
Quantification réussie des décalages systématiques des variables PID clés telles que l'identification électron/muon via les paramètres {α, β}, comme le montre la figure 2 illustrant l'effet de différentes combinaisons de paramètres sur la forme de la distribution.
La figure 3 présente les spectres énergétiques prédits a posteriori pour différents échantillons d'analyse T2K, les erreurs bin reflétant l'incertitude 1σ provenant de la matrice d'erreur du détecteur SK. Des différences significatives dans les erreurs bin sont observées entre l'analyse 2022 et 2023, suite à la suppression de certaines erreurs systématiques "temporaires".
Réalisation d'améliorations de performance significatives:
- Temps de convergence: Réduction de O(mois) à O(semaines)
- Dimension des paramètres: Extension à un ajustement MCMC à 540 dimensions
- Couverture des échantillons: Support de l'estimation des erreurs systématiques pour les nouveaux échantillons d'analyse
Vérification de la rationalité de l'application des résultats d'ajustement des neutrinos atmosphériques à l'analyse du faisceau T2K par comparaison du chevauchement entre neutrinos atmosphériques et neutrinos du faisceau en termes de gamme cinématique, topologie des événements et conditions du détecteur.
Les expériences d'oscillation des neutrinos de long parcours font face universellement aux défis des erreurs systématiques du détecteur, avec des expériences telles que NOνA et DUNE développant des méthodes d'estimation d'erreur correspondantes.
La modélisation précise des détecteurs à eau Cherenkov a toujours été un sujet important dans ce domaine, impliquant plusieurs aspects tels que les propriétés optiques, la réponse des photomultiplicateurs et les algorithmes de reconstruction des événements.
Par rapport à d'autres méthodes, le cadre proposé dans cet article présente:
- Généralité plus forte: Méthode de paramétrisation indépendante du modèle
- Précision plus élevée: Description fine de l'espace de paramètres à 540 dimensions
- Robustesse meilleure: Capacité à s'adapter aux changements des conditions du détecteur
- Établissement réussi d'un cadre d'estimation des erreurs systématiques du détecteur SK basé sur MCMC
- Réalisation d'un transfert efficace des erreurs des données de neutrinos atmosphériques à l'analyse du faisceau T2K
- Fourniture d'un moyen important de contrôle des erreurs systématiques pour les analyses T2K à venir
- Limitations des hypothèses: Hypothèse que les caractéristiques des erreurs systématiques des neutrinos atmosphériques et du faisceau sont similaires
- Dépendance du modèle: Dépendance toujours présente de la précision fondamentale de la simulation MC
- Complexité computationnelle: L'ajustement MCMC à 540 dimensions nécessite toujours des ressources computationnelles considérables
- Optimisation supplémentaire de l'algorithme MCMC pour améliorer l'efficacité de convergence
- Exploration d'une classification des événements plus fine et d'une modélisation des erreurs systématiques
- Combinaison avec d'autres échantillons de contrôle pour vérifier la précision de l'estimation des erreurs
- Forte innovativité de la méthode: Proposition d'une méthode de paramétrisation des erreurs systématiques indépendante du modèle
- Conception expérimentale rationnelle: L'utilisation de neutrinos atmosphériques comme échantillon de contrôle possède une forte motivation physique
- Implémentation technique complète: L'ajustement MCMC à 540 dimensions démontre une capacité technique de haut niveau
- Valeur pratique élevée: Service direct à l'amélioration de la précision de l'analyse physique de T2K
- Analyse théorique limitée: Manque de discussion approfondie sur les fondements théoriques de la méthode de paramétrisation
- Moyens de vérification uniques: Dépendance principalement des données de neutrinos atmosphériques, manque de vérification indépendante
- Quantification de l'incertitude: Discussion insuffisante des incertitudes introduites par la méthode elle-même
- Contribution académique: Fourniture d'une méthode importante de contrôle des erreurs systématiques pour les expériences de physique des neutrinos
- Valeur pratique: Impact direct sur la précision des résultats physiques de l'expérience T2K
- Promotion technologique: La méthode peut être généralisée à d'autres expériences de physique des particules similaires
Cette méthode est particulièrement adaptée à:
- Expériences d'oscillation des neutrinos de long parcours
- Grands détecteurs à eau Cherenkov
- Mesures de physique des particules nécessitant un contrôle des erreurs systématiques haute précision
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)
Évaluation générale: Cet article est un travail de physique expérimentale très technique qui fournit une contribution méthodologique importante au contrôle des erreurs systématiques de l'expérience T2K. Bien qu'il présente certaines insuffisances dans l'analyse théorique, sa valeur pratique et son innovativité technique en font un travail important dans le domaine de la physique des neutrinos.