ALICE à la recherche des particules étranges

1. Résumé

L’exercice proposé consiste à rechercher des particules étranges, produites dans des collisions au LHC et détectées par l’expérience ALICE. Cette recherche est basée sur la reconnaissance de leur désintégration typique dite V0, comme, par exemple, Kos → π+π-, Λ → p + π-, et Ξ- → Λ + π- (Λ → p + π-). -). L’identification des particules étranges est basée sur la topologie de leur désintégration et l’identification des produits de la désintégration; l’information des traces est utilisée pour calculer la masse invariante de la particule-mère et ainsi confirmer l’identité de cette particule.

Dans les sections suivantes l’expérience ALICE est brièvement présentée ainsi que ses objectifs de physique; ensuite la motivation de cette analyse. La méthodologie utilisée pour l’identification des particules étranges, ainsi que les outils sont présentés en détail; ensuite l’exercice est expliqué pas à pas, suivis par la présentation des résultats. La méthode pour rassembler et combiner tous les résultats est également présentée ainsi que quelques suggestions de discussion.


2. Introduction

ALICE (A Large Ion Collider Experiment), l’une des quatre grandes expériences du LHC au CERN, est destinée à l’étude de collisions d’ions lourds. Elle étudie également les collisions de protons, qui fournissent des données utilisées comme référence pour les collisions d'ions lourds. En plus, les donnés des collisions de protons sont utiles pour des études spécifiques de la physique proton-proton. Le détecteur ALICE a été conçu en vue des multiplicités de particules très élevées anticipées pour des collisions d’ions de plomb à l’énergie du LHC.


3. La physique de ALICE

Les quarks sont liés dans les protons et les neutrons par une force, dite interaction forte, par l’intermédiaire d’une échange de particules-messagères appelées gluons. L’interaction forte lie également les protons et les neutrons dans les noyaux atomiques.

Le quark, particule élémentaire qui constitue tous les hadrons, n’a jamais pu être observé isolé : les quarks ainsi que les gluons paraissent être prisonniers dans des particules composées, comme les protons et les neutrons. Le mécanisme à l’origine de ce «confinement » reste toujours inconnu.

Malgré le fait que la physique de l’interaction forte est bien comprise aujourd’hui, deux questions restent non résolues : l’origine du confinement et le mécanisme qui donne la masse aux quarks liés. La réponse à ces deux questions impliquerait la façon dont l’interaction forte modifie les propriétés du vide.

La théorie courante de l’interaction forte (appelée chromodynamique quantique) prévoit qu’à des températures et densités très élevées, les quarks et les gluons se libèrent des particules composées, mais peuvent donc exister libres.

Une telle transition devait avoir lieu quand la température dépasse une valeur critique, estimée à environ 100 000 fois plus élevée que la température au coeur du Soleil! De telles températures n’existent pus dans la nature depuis la naissance de l’Univers. En effet, selon la théorie du Big Bang la température de l’Univers était pendant quelques millionièmes de seconde au dessus de la valeur critique, et toute la matière de l’Univers était dans un état de plasma de quarks et de gluons.

Quand deux noyaux lourds s’approchent à une vitesse proche de celle de la lumière et entrent en collision, ces conditions extrêmes de température peuvent être recréées et libérer ainsi les quarks et les gluons. Quarks et gluons entrent en collisions en créant un environnement en équilibre thermique : le plasma de quarks et de gluons. Ce plasma se dilate et se refroidit à la température (1012 degrés) où les quarks et les gluons se regroupent pour former la matière ordinaire, à peine 10-23 seconds après le début de la collision. ALICE étudiera la formation et les propriétés de ce nouvel état de la matière.


4. Augmentation de la production de l’étrangeté comme signaturede la formation du plasma de quarks et de gluons

Le diagnostic et l'étude des propriétés du plasma de quarks et de gluons (QGP) peuvent être réalisés à l'aide de quarks qui n'existent pas dans la matière qui nous entoure. Une des signatures expérimentales s'appuie sur l'idée d'augmentation de la production de l'étrangeté. Il s'agissait d'une des premières observables du plasma de quarks et de gluons, proposée en 1980.

Contrairement aux quarks up et down, les quarks étranges ne sont pas apportés dans la réaction par les noyaux en collision. Par conséquent, tout quark ou antiquark étranges observés par les expériences ont été créés à partir de de l'énergie cinétique des noyaux en collision. Comme la masse des quarks étranges est de l’ordre de grandeur de la température à laquelle les protons, les neutrons et les autres hadrons se dissolvent en quarks, l'abondance des quarks étranges est sensible aux conditions, la structure et la dynamique de la phase de la matière déconfinée. Ainsi une production élevée de quarks étranges signale que les conditions de déconfinement ont été atteintes.

Dans la pratique, l'augmentation de l'étrangeté peut être observée en comptant le nombre de particules étranges, c'est-à-dire les particules contenant au moins un quark étrange, et calculer le rapport particules étranges-particules non-étranges. Si ce rapport est supérieur à celui donné par les modèles théoriques qui ne prévoient pas la création de QGP, l’augmentation est établie.

Pour le cas de collisions d’ions de plomb, le nombre de particules étranges est normalisé au nombre de nucléons participant à l’interaction.


5. Particules étranges

Les particules étranges sont des hadrons contenant au moins un quark étrange. Cela se caractérise par le nombre quantique de l’« étrangeté ». Le méson étrange neutre le le plus léger est le Kos (ds) et le baryon étrange neutre le plus léger est le Λ (uds), de la famille des hypérons.

On étudiera leur désintégration, par exemple, Kos → π+π-, Λ → p + π-. Dans ces désintégrations le nombre quantique de l'étrangeté n'est pas conservé, puisque les produits de décomposition sont uniquement composés de quarks up et down. Par conséquent, il ne s’agit pas de désintégrations fortes (qui en outre seraient très rapides, avec un τ = 10-23 s) mais de désintégrations faibles, où l'étrangeté peut être conservée (ΔS = 0) ou modifiée d’une unité (ΔS = 1). Pour ces désintégrations la vie moyenne τ est entre 10-8 et 10-10 s. Pour les particules avec des vitesses proches de celle de la lumière, cela signifie que la particule se désintègre à une distance (en moyenne) de quelques cm du point de production (du point de l'interaction).


6. Comment identifier les particules étranges

Le but de l'exercice est la recherche de particules étranges produits par des collisions au LHC et enregistrées par l'expérience ALICE.

Comme mentionné dans la section précédente, les particules étranges ne vivent pas longtemps; elles se désintègrent peu après leur production. Cependant, elles vivent assez longtemps pour voyager sur quelques cm de distance du point d'interaction (IP), où elles ont été produites (vertex primaire). Leur recherche est donc basée sur l'identification de leurs produits de désintégration, qui doivent provenir d'un vertex secondaire commun.

Les particules étranges neutres, comme les Kos et les Λ , se désintègrent en donnant un motif caractéristique, appelé V0. La particule-mère disparaît à quelques cm du point d'interaction et deux particules de charges opposées apparaissent à sa place ; leur trajectoire est incurvée dans des directions opposées par le champ magnétique du solénoïde d’ALICE.

Dans ce qui suit les traces rouges indiquent des particules de charge positive; les traces vertes indiquent des particules de charge négative.

Les désintégrations que nous chercherons sont:

Kos → π+π-

Λ → p + π-

anti Λ→ p- + π+

Nous voyons que pour un état final de deux pions le motif de la décroissance est quasi-symétrique, alors que dans l'état final d’un pion et un proton, le rayon de courbure du proton est plus important que celui du pion : en raison de sa masse plus élevée, le proton transporte la plupart de l'impulsion initiale.

Nous rechercherons également des désintégrations en cascade de particules étranges chargées, tels que le Ξ-; cette particule se désintègre en π- et Λ; ; ensuite le Λ se désintègre en π- et proton; le pion initial s’appelle bachelor (célibataire) représenté en violet.

Ξ-→π-Λ→ π- p + π-

La recherche de V0s est basée sur la topologie de la désintégration et l'identification des produits de désintégration ; une confirmation supplémentaire de l'identité de la particule est le calcul de sa masse; cela se fait en utilisant les informations (masse et impulsion) des produits de décomposition comme décrit dans la section suivante.


7. Le calcul de la masse (invariante)

Nous considérons la désintégration ci-dessus du kaon neutre en deux pions chargés. Soit E, p et m l'énergie totale, l'impulsion et la masse de la particule - mère (K0)
Soit E1, p1 et m1 totale, l'impulsion et la masse de la particule-fille numéro 1 (π+); et E2, p2 et m2 l'énergie totale, l'impulsion et la masse de la particule-fille numéro 2 (π-).

Conservation de l'énergie E = E1 + E2 (1)
Conservation de l'impulsion p = p1 + p2 (2)
De la relativité: E2 = p2 + m2 (3)

Cela s'applique, bien entendu, également pour les particules-filles:

E12 = p12 + m12 (4)
E22 = p22 + m22 (5)

où p1=|p1| et p2=|p2| sont les longueurs des vecteurs p1 et p2.

Des équations ci-dessus, nous trouvons que:

m2 = E2 - p2 = (E1+E2) 2 -(p1+p2) 2
= E12 + E22 + 2E1E2 - p1.p2 - p2.p2 - 2p1.p2 (6)

Le produit de deux vecteurs est la somme des produits des composantes x, y et z de deux vecteurs:

p1.p2 = p1x p2x + p1y p2x + p1z p2z (7)
p1.p1 = p1x2 + p1y2 + p1z2 = p12 (8)
p1.p1 = p1x2 + p1y2 + p1z2 = p12 (9)
m2 = E12 + E22 + 2E1E2 - p12 - p22 - 2p1.p2
= m12 + m12 + 2E1E2 - 2p1.p2

Nous pouvons donc calculer la masse de la particule initiale et la masse et les composants de l'impulsion des particules-filles.

Les masses des particules-filles m1 et m2 sont connues: différents détecteurs d'ALICE identifient les particules.

Les impulsions des particules-filles p1, p2 se déduisent de la mesure du rayon de courbure de leur trajectoire due au champ magnétique connu. Dans l'exercice, nous utilisons les trois composantes du vecteur d’impulsion de chaque trace associée à la désintégration V0, comme dans les équations ci-dessus.

Le calcul de la masse invariante donne généralement des distributions comme indiqué ci-dessous. La distribution sur la gauche est la masse calculée pour les paires de pion-proton ; le pic correspond à des Λ et le continuum est le bruit de fond de combinaisons aléatoires de pions et de protons qui apparaissent venir du même vertex secondaire ou qui ont été mal identifiés; la distribution à droite est la masse calculée pour les paires de pions négatifs et positifs ; le pic correspond à Kos.


8. Les outils et comment les utiliser

L'exercice est effectuée dans le cadre de ROOT, en utilisant une version simplifiée du programme de visualisation utilisé par ALICE : dans une fenêtre de terminal qui est déjà ouverte sur votre ordinateur (de sorte que vous êtes dans le répertoire approprié) vous tapez “root masterclass.C”. Une petite fenêtre apparait, comme dans l'image.


Cela offre trois options : démonstration, étudiant pour l'analyse des événements et enseignant pour la collecte et fusion des résultats.

La démonstration donne des exemples de désintégrations de Kos, Λ, anti-Λ et Ξ-.

L’option «étudiant » pour l'analyse des événements et de la recherche visuelle pour V0s ouvre une fenêtre, comme le montre la figure suivante.

La colonne de gauche offre un nombre d'options: Instructions, Navigation des événements, « finder » de V0 et cascades, calculatrice, sélection de ce qui est affiché (traces, géométrie du détecteur,...). En outre il y a l’animation d’un événement et « Encyclopédie », avec une brève description du détecteur d'ALICE et ses composantes principales ainsi que les motifs de la désintégration V0.

La visualisation de l'événement montre trois vues du détecteur ALICE (vue tridimensionnelle, projection rφ et projection rz). Vous pouvez sélectionner les informations affichées pour chaque événement. Si vous cliquez sur la case correspondante, vous voyez tous les clusters et les traces de l'événement ; si vous cliquez sur les cases du finder V0 (et cascade), les V0s (et cascades) sont mis en évidence, s’ils existent. Lorsqu’un V0 est trouvé, le reste des traces et des clusters de l'événement peuvent être supprimés et seules les traces associées au V0 sont indiquées. La convention de couleur est que les traces positives de V0 sont rouges, les traces négatives sont vertes (et les « bachelors », dans le cas des cascades, violettes).


En cliquant sur chaque trace, les valeurs des composantes de l’impulsion et la masse des particules, (celle avec la probabilité maximale, selon les algorithmes d'identification de particules) apparaissent dans une petite fenêtre (figure suivante, à droite). Cette information peut être copiée sur la calculatrice, qui calcule ensuite la masse invariante de la particule-mère, à l'aide de la formule expliquée dans la section précédente.

Le programme comprend quatre histogrammes de masse invariante (pour Kos, Λ, anti-Λ et Ξ-). Après avoir inspecté chaque désintégration V0, vous pouvez identifier la particule-mère par les produits de désintégration et la valeur de la masse invariante (un tableau de référence de masses de certaines particules est donnée dans le cadre de la calculatrice, voir figure). Ensuite, vous appuyez sur le bouton correspondant (c'est un kaon; c'est une Lambda etc..). De cette façon, vous ajoutez une entrée à l'histogramme correspondant. Les histogrammes de la masse invariante peuvent être affichés en cliquant sur le bouton de masse invariante, au-dessus de la visualisation de l'événement. Pour mettre à jour leur contenu, vous devez cliquer à l'intérieur de chaque histogramme.


Calcul et histogrammes de rapidité ont été également appliqués ; ceux-ci ont été prévus pour les étudiants de l'Université et sont au-delà de la portée de cet exercice.


9. L'exercice - analyser les événements et trouver les hadrons étranges

La partie analyse consiste de l'identification et le comptage de particules étranges dans un échantillon donné, généralement contenant 30 événements. Lors du démarrage de l'exercice, vous devez cliquer sur « student » et sélectionner l'échantillon d'événements que vous allez analyser. Actuellement, il y a 6 échantillons de différents événements avec des données de collisions des protons à 7 TeV d’énergie au centre de masse.

Quand vous inspectez un événement, vous devez initialement cliquer sur les boutons pour visualiser les clusters et les traces ; vous pouvez observer la complexité des événements et le nombre élevé de traces produites par les collisions à l'intérieur des détecteurs. La plupart de ces traces sont des pions.

En cliquant sur « V0 » et « Cascades » les traces des désintégrations de V0 - si il y en a - et des cascades - si il y en a - apparaissent en surbrillance. En cliquant sur chaque trace, vous obtenez l'information sur la trace – la charge, les trois composantes du vecteur de l'impulsion et la masse de la particule plus probable associée à la trace. Cela a été trouvé par l'information fournie par les différents détecteurs utilisés pour l'identification des particules. Par les produits de décomposition, vous pouvez déjà deviner ce qu'est la particule-mère ; pour le confirmer, vous calculez la masse invariante, comme expliqué dans la section 7 et comparez sa valeur avec les valeurs données au tableau de votre calculatrice.

Si la masse est 497 MeV ±13 MeV (dans l’intervalle [484, 510] MeV) c’est un Kos;
Si la masse est 1115 MeV ±5 MeV (dans l’intervalle [1110, 1120] MeV) et les particules-filles sont un proton et un pion négatif, alors c’est un Λ.
Si la masse est 1115 MeV ±5 MeV (dans l’intervalle [1110, 1120] MeV) et les particules-filles sont un antiproton et un pion positif, alors c’est un anti-Λ.
Pour une désintégration en cascade, si la masse calculée par les 3 traces est 1321 ±10 MeV (dans l’intervalle [1311, 1331] MeV ) alors c’est un Ξ-.

Selon le résultat, vous cliquez sur le bouton « C'est un Kaon, Lambda, etc ». De cette façon, cette entrée est ajoutée dans l'histogramme de masse invariante correspondant.

Il peut arriver que la masse calculée ne corresponde pas à l'une des valeurs ci-dessus. Il s'agit de bruit de fond: les traces apparaissent venir d' un vertex secondaire, qui dans ce cas a été mal identifié. Pour cet exercice, nous ne comptons pas ces V0.

10. Présentation des résultats


Le tableau présenté ici résume les résultats. La colonne "real data" contient le nombre de Kos, Λ, anti-Λ et Ξ- que vous avez trouvé (à condition que vous n’avez pas oublié d’appuyer sur le bouton « C'est un Kaon, Lambda etc. »).


Vous pouvez également regarder les histogrammes de la masse invariante et vérifier le nombre d'entrées pour chaque type de particule. Lorsque vous avez analysé tous les événements de votre échantillon de données, enregistrez les résultats sur un fichier en suivant les instructions à l'intérieur du programme d'analyse.

11. Collection de tous les résultats

En sélectionnant l'option «Teacher» dans le menu initial de MasterClass, vous pouvez recueillir tous les résultats. Sous « Teacher Controls», vous sélectionnez l'option «Get Files » et obtenez, un à un, les fichiers avec les résultats de l'analyse de chaque échantillon de données. Évidemment, vous devez transférer d'abord les fichiers avec les résultats sur l’ordinateur le l’« enseignant »! Puis, sous « Résultats », vous pouvez consulter le tableau avec la statistique complète.


12. Analyse de grande statistique

La visualisation de l'événement est un outil puissant qui aide à vérifier la qualité des données et leur reconstruction et donne une idée à quoi ressemblent les événements. Toutefois, dans la vraie vie, l’analyse des données ne se fait pas visuellement - ça serait beaucoup trop long et fastidieux. Pour analyser les millions d'événements que nous recueillons quotidiennement au LHC nous exécutons des programmes, et c'est ce que vous ferez ici, afin de rechercher des V0s dans un plus grand échantillon d'événements.

Modifiez le répertoire sur votre fenêtre de terminal (en donnant la commande: cd MasterClass_extended) et écrivez root MasterClassExtended.C. Dans l'espace « put your name here», donnez une combinaison de caractères qui formeront le nom du fichier de résultats. Choisissez un échantillon de données à analyser (actuellement il y a 6 échantillons de données de collisions proton-proton à 7 TeV avec 2000 événements chaque échantillon) ; puis choisissez « student » pour procéder à l'analyse.

Sous «Analysis tools», vous pouvez analyser 100 événements ou 2000 événements et calculer la masse invariante de paires de particules, comme π+π-. Vous pouvez voir que la masse invariante est une distribution continue – c'est parce que les paires de pions combinées sont aléatoires, n’ayant pas d’origine commune (vertex secondaire) et peuvent donner n'importe quelle valeur de masse. Il s'agit de bruit de fond. Lorsque vous continuez à la sélection de V0, seulement paires de traces provenant d'un vertex secondaire commun sont considérés ; leur masse invariante est calculée selon les données des traces et la masse des produits identifiés de désintégration. Vous pouvez sélectionner Kos ou Λ (y compris les anti-Λ). Chaque fois que l'analyse de tous les événements de l'échantillon est fini (observez le terminal derrière le menu) cliquez sur l'écran avec les distributions de masses invariantes pour afficher l'histogramme correspondant.

Afin de trouver le nombre de particules d'un certain type, par exemple Kos, vous devez trouver le nombre d'événements dans le pic après soustraction du bruit de fond. Afin de décrire le bruit de fond par une courbe (polynôme de second degré), vous choisissez d'abord la région du fit à l'aide du curseur et cliquez sur «Fit background». Lorsque vous cliquez sur l'écran, la fonction ajustée est superposée sur l'histogramme et vous pouvez vérifier visuellement si l'ajustement est raisonnable. Vous décrivez ensuite le signal par une distribution gaussienne, en sélectionnant d'abord la région du pic. Pour la soustraction du bruit de fond, les coefficients de la polynôme de deuxième degré sont utilisés ; en cliquant sur l'histogramme, vous obtenez le nombre total d'événements dans le pic, le nombre d'événements correspondant au bruit de fond et ceux qui sont le signal, ainsi que la valeur moyenne de la gaussienne (la masse) et sa largeur, sigma.

Afin de sauver les histogrammes, vous devez d'abord sélectionner le répertoire où les histogrammes seront sauvegardés (valeur par défaut est «Teacher») et ensuite le sous-répertoire : K0s (pour Kaons), Lambda (pour Λ et anti-Λ ensemble). En cliquant sur 1,2,3 ou 4 vous enregistrez l'histogramme affiché sur la partie supérieure gauche, supérieure droite, inférieure gauche ou inférieure droite de l'écran. Il pourrait être difficile d'ajuster les données qui correspondent à 2000 événements – il sera plus facile avec la statistique plus haute une fois que vous avez ajouté tous les résultats, consultez la section suivante.

13. Collection de tous les résultats de l’ analyse de grande statistique

Choisissez « Teacher mode » dans le menu « large-scale analysis ». L’option «Get Files» fonctionne comme suit : en cliquant au dessus de « Teacher », le répertoire par défaut, vous pouvez choisir les sous-répertoires de Kaons (K0s), Λ et anti-Λ (Lambda). Une fois que vous avez sélectionné le répertoire (« Teacher » peut être modifié pour un répertoire de votre choix, où vous avez enregistré tous les histogrammes) et le sous-répertoire, cliquez sur 1,2,3 ou 4 – tous les fichiers d'histogrammes sur ce sous-répertoire seront ajoutés (l’ information sur le nombre et la taille des fichiers est affichée dans la fenêtre de terminal derrière le menu; lorsque vous cliquez sur l'écran, l'histogramme fusionné sera affiché dans la partie supérieure gauche (1), supérieure droite (2), inférieure gauche (3) ou inférieure droite (4) de l'écran. Afin de trouver le nombre total de particules d’une certaine espèce suivez la procédure pour ajuster le bruit de fond et le signal décrit dans la section précédente.

14. Calcul du rendement des particules

À l'aide d'information qui vous sera donnée lorsque vous avez terminé l'analyse (par exemple l'efficacité pour chaque type de particule) calculez le rendement (nombre de particules produites par interaction) pour chaque type de V0.

15. Aperçu de l’exercice

Analyse visuelle

  1. Étudiez les exemples de désintégrations de V0 et découvrez comment utiliser les outils (V0 finder, calculatrice, histogrammes).
  2. Analyser visuellement un échantillon de 30 événements V0.
  3. Ajouter les résultats de l'analyse de tous les groupes.

Analyse de grande statistique

  1. Analyser 2000 événements / observer la distribution de masse du bruit de fond.
  2. Analyser 2000 événements en cherchant les Kos / adjuster le bruit de fond/ adjuster le pic
  3. Analyser 2000 événements en cherchant les Λ et anti-Λ / adjuster le bruit de fond/ adjuster le pic.
  4. Ajouter les résultats de l'analyse de tous les groupes.