2025-11-10T03:09:47.574386

Estimating detector systematic uncertainties for the T2K far detector

Reh
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.
academic

Estimation des incertitudes systématiques du détecteur pour le détecteur lointain de T2K

Informations fondamentales

  • 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

Résumé

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.

Contexte et motivation de la recherche

Problème fondamental

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.

Importance

  1. 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
  2. Domination des erreurs systématiques: Les erreurs systématiques du détecteur SK constituent la principale source d'incertitude dans l'analyse T2K actuelle
  3. Exigences de précision: Avec l'amélioration de la précision statistique, le contrôle des incertitudes systématiques devient plus critique

Limitations des méthodes existantes

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.

Motivation de la recherche

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.

Contributions principales

  1. É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
  2. 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
  3. Développement d'un programme d'ajustement MCMC à 540 dimensions, supportant les nouveaux échantillons d'analyse νeCC1π± et NCπ0
  4. 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
  5. 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

Explication détaillée de la méthode

Définition de la tâche

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

Architecture de la méthode fondamentale

1. Paramétrisation des erreurs systématiques

T2K adopte une méthode de paramétrisation générale pour décrire les erreurs de modélisation du détecteur:

L → α · L + β

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

2. Ajustement des neutrinos atmosphériques

Estimation des paramètres utilisant le cadre MCMC:

Étapes:

  1. MCMC propose des combinaisons de paramètres {α, β}
  2. Modification des valeurs PID des événements MC selon la formule
  3. Construction d'histogrammes des événements MC modifiés par bin
  4. Calcul de la fonction de vraisemblance de forme entre les données et le MC modifié
  5. Répétition du processus pour obtenir la distribution a posteriori de {α, β}

3. Construction de la matrice d'erreur

Transfert des résultats d'ajustement des neutrinos atmosphériques à l'analyse du faisceau T2K via un programme MC jouet:

  1. Extraction aléatoire de paramètres {α, β} de la distribution a posteriori MCMC
  2. Modification des valeurs PID des événements MC du faisceau T2K
  3. Observation de la migration des événements entre différents échantillons d'analyse
  4. Statistiques des variations du spectre énergétique des neutrinos
  5. Construction de la matrice de covariance

Points d'innovation technique

1. Indépendance du modèle

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

2. Séparation de la topologie des événements

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.

3. Espace de paramètres haute dimension

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.

Configuration expérimentale

Ensembles de données

  1. 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
  2. Données du faisceau T2K: Ensemble de données Run 1-4, incluant les modes de faisceau νμ et ν̄μ
  3. Simulation Monte Carlo: Nouveau produit MC utilisant des modèles d'interaction de neutrinos mis à jour

Indicateurs d'évaluation

  1. Fonction de vraisemblance de forme: Évaluation de la cohérence entre la distribution des données et celle du MC
  2. Convergence de la distribution a posteriori: Assurance de la fiabilité de l'estimation des paramètres via des diagnostics MCMC
  3. Rationalité physique de la matrice de covariance: Vérification de la rationalité des corrélations d'erreur

Échantillons d'analyse

  • Échantillon d'électron monoanneau (Single ring e-like)
  • Échantillon de muon monoanneau (Single ring μ-like)
  • Échantillon νeCC1π+ (nouveau)
  • Échantillon NCπ0 amélioré (nouveau)

Résultats expérimentaux

Résultats principaux

1. Quantification des décalages systématiques des variables PID

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.

2. Estimation de l'incertitude du spectre énergétique

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".

3. Amélioration de l'efficacité computationnelle

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 méthode

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.

Travaux connexes

Erreurs systématiques dans les expériences d'oscillation des neutrinos

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.

Techniques de modélisation du détecteur

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.

Avantages de cet article

Par rapport à d'autres méthodes, le cadre proposé dans cet article présente:

  1. Généralité plus forte: Méthode de paramétrisation indépendante du modèle
  2. Précision plus élevée: Description fine de l'espace de paramètres à 540 dimensions
  3. Robustesse meilleure: Capacité à s'adapter aux changements des conditions du détecteur

Conclusions et discussion

Conclusions principales

  1. Établissement réussi d'un cadre d'estimation des erreurs systématiques du détecteur SK basé sur MCMC
  2. Réalisation d'un transfert efficace des erreurs des données de neutrinos atmosphériques à l'analyse du faisceau T2K
  3. Fourniture d'un moyen important de contrôle des erreurs systématiques pour les analyses T2K à venir

Limitations

  1. Limitations des hypothèses: Hypothèse que les caractéristiques des erreurs systématiques des neutrinos atmosphériques et du faisceau sont similaires
  2. Dépendance du modèle: Dépendance toujours présente de la précision fondamentale de la simulation MC
  3. Complexité computationnelle: L'ajustement MCMC à 540 dimensions nécessite toujours des ressources computationnelles considérables

Directions futures

  1. Optimisation supplémentaire de l'algorithme MCMC pour améliorer l'efficacité de convergence
  2. Exploration d'une classification des événements plus fine et d'une modélisation des erreurs systématiques
  3. Combinaison avec d'autres échantillons de contrôle pour vérifier la précision de l'estimation des erreurs

Évaluation approfondie

Points forts

  1. Forte innovativité de la méthode: Proposition d'une méthode de paramétrisation des erreurs systématiques indépendante du modèle
  2. Conception expérimentale rationnelle: L'utilisation de neutrinos atmosphériques comme échantillon de contrôle possède une forte motivation physique
  3. Implémentation technique complète: L'ajustement MCMC à 540 dimensions démontre une capacité technique de haut niveau
  4. Valeur pratique élevée: Service direct à l'amélioration de la précision de l'analyse physique de T2K

Insuffisances

  1. Analyse théorique limitée: Manque de discussion approfondie sur les fondements théoriques de la méthode de paramétrisation
  2. Moyens de vérification uniques: Dépendance principalement des données de neutrinos atmosphériques, manque de vérification indépendante
  3. Quantification de l'incertitude: Discussion insuffisante des incertitudes introduites par la méthode elle-même

Influence

  1. Contribution académique: Fourniture d'une méthode importante de contrôle des erreurs systématiques pour les expériences de physique des neutrinos
  2. Valeur pratique: Impact direct sur la précision des résultats physiques de l'expérience T2K
  3. Promotion technologique: La méthode peut être généralisée à d'autres expériences de physique des particules similaires

Scénarios d'application

Cette méthode est particulièrement adaptée à:

  1. Expériences d'oscillation des neutrinos de long parcours
  2. Grands détecteurs à eau Cherenkov
  3. Mesures de physique des particules nécessitant un contrôle des erreurs systématiques haute précision

Références

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.