Cosmic-Ray Astronomy at the Highest Energies with Ten Years of Data of the Pierre Auger Observatory
Total Page:16
File Type:pdf, Size:1020Kb
Universit´e Pierre et Marie Curie Ecole´ doctorale des Sciences de la Terre et de l’environnement et Physique de l’Univers, Paris Laboratoiredephysiquenucl´eaire et des hautes ´energies Cosmic-ray astronomy at the highest energies with ten years of data of the Pierre Auger Observatory Par Lorenzo Caccianiga Th`ese de doctorat de physique Dirig´ee par: Piera Luisa Ghia Co-encadr´ee par: Julien Aublin Pr´esent´ee et soutenue publiquement le 14 septembre 2015 Devant un jury compos´ede: Piera Luisa Ghia - Directrice de th`ese Miguel Mostafa - Rapporteur Elisa Resconi - Rapporteur Corinne Berat - Examinateur Delphine Hardin - Examinateur Cyril Lachaud - Examinateur ii Abstract Identifying the sources of the ultra-high energy cosmic rays (UHECRs, above ∼ 1018 eV), the most energetic particles known in the universe, would be an important leap forward for both the astrophysics and particle physics knowledge. However, developing a cosmic- ray astronomy is arduous because magnetic fields, that permeate our Galaxy and the extra-Galactic space, deflect cosmic rays that may lose the directional information on their sources. This problem can be reduced by studying the highest energy end of the cosmic ray spectrum. Indeed, magnetic field deflections are inversely proportional to the cosmic ray energy. Moreover, above ∼ 4·1019 eV, cosmic rays interact with cosmic photon backgrounds, losing energy. This means that the sources of the highest energy cosmic rays observed on Earth can be located only in the nearby universe (∼ 200 Mpc or less). The largest detector ever built for detecting cosmic rays at such high energies is the Pierre Auger Observatory, in Argentina. It combines a 3000 km2 surface array of water Cerenkovˇ detectors with fluorescence telescopes to measure extensive air showers initiated by the UHECRs. This thesis was developed inside the Auger Collaboration and was devoted to study the highest energy events observed by Auger, starting from the selection and reconstruction up to the analysis of their distribution in the sky. Moreover, since the composition at these energies is unknown, we developed a method to select proton-like events, since high Z cosmic rays are too much deflected by magnetic fields to be used for cosmic-ray astronomy. In the first chapter, we will give a description of the current general knowledge of cosmic rays: their measurements at Earth, their propagation and consideration on possible acceleration sites. In chapter 2, the principle of detection of UHECRs will be summarized, describing the properties of the extensive air showers (EAS) they produce as they interact with the atmosphere and how they can be detected. Chapter 3 will focus on the Pierre Auger Observatory, describing its technical features and the key results it obtained in 10 years of operation. Chapter 4 will detail the reconstruction of events recorded by the surface detector of the Pierre Auger Observatory, that will be the ones used in this thesis work. Chapter 5 will describe the first part of the work: the building of the dataset of the highest energy events recorded by Auger. A thorough study of all the selected events will be described, together with the improvement to the event selection and reconstruction that were motivated by the findings of this study. In chapter 6 we will present the anisotropy tests to which this dataset has been subjected and the cross-checks made in this thesis work. In chapter 7 we will describe a multi-variate analysis method to separate light from heavy events in the dataset. Finally, in chapter 8, we will perform anisotropy tests on the proton-enriched samples obtained this way. i ii Contents 1 Cosmic Rays 1 1.1 Cosmic rays at the Earth: energy spectrum, mass composition and arrival directions . 2 1.1.1 Spectrum . 2 1.1.2 Composition . 4 1.1.3 Arrival directions . 4 1.2 Cosmic ray propagation . 5 1.2.1 Magnetic field deflections . 6 1.2.2 Interactions with background radiations . 8 1.3 UHECR origin . 10 1.3.1 Acceleration mechanism . 10 1.3.2 Candidate sources . 12 1.3.3 Alternative models . 15 2 UHECR Extensive Air Showers: Physics and Detection 17 2.1 Extensive Air Showers . 18 2.1.1 Electromagnetic showers . 19 2.1.2 Hadronic showers . 20 2.1.3 Shower simulations and hadronic interaction models . 21 2.2 Principles of Extensive Air Shower Detection . 24 2.2.1 Arrays of particle detectors . 24 2.2.2 Light detectors . 25 2.2.3 Other techniques . 26 3 The Pierre Auger Observatory 29 3.1 The site . 30 3.2 The Surface Detector . 31 3.2.1 Cerenkovˇ stations . 31 3.2.2 Data communication and acquisition . 32 3.2.3 Maintenance . 33 3.2.4 Aperture and exposure . 33 3.3 Fluorescence Detector . 34 3.3.1 The telescopes . 34 3.3.2 FD event reconstruction . 34 3.3.3 Atmospheric monitoring . 37 iii 3.3.4 Hybrid exposure calculation . 37 3.4 Enhancements . 38 3.4.1 Low-energy extensions . 38 3.4.2 Radio and Microwave detection . 39 3.4.3 Planned upgrades . 39 3.5 Highlight of the Pierre Auger Observatory results . 40 3.5.1 Energy Spectrum . 40 3.5.2 Mass composition: depth of shower maximum from FD . 42 3.5.3 Mass-sensitive observables from SD and hadronic models . 45 3.5.4 Arrival directions . 48 3.5.5 Search for neutral particles . 50 4 Reconstruction of surface detector events 57 4.1 SD Triggers . 58 4.1.1 Station trigger: T1 & T2 . 58 4.1.2 T3: array trigger . 59 4.2 Signal treatment . 60 4.2.1 The Vertical Equivalent Muon (VEM) . 61 4.2.2 Conversion of the FADC trace into VEM units . 61 4.2.3 Saturation recovery . 61 4.3 Selection and quality . 62 4.3.1 Station selection and T4 physics trigger . 63 4.3.2 T5: fiducial trigger . 64 4.4 Reconstruction . 65 4.4.1 Time fit: first estimate of the arrival direction . 65 4.4.2 Lateral distribution function and global fit . 65 4.4.3 Uncertainty on the arrival direction . 68 4.4.4 Uncertainty on S1000 ........................... 69 4.5 From S1000 to primary energy . 70 4.5.1 From S1000 to S38: the constant intensity cut . 71 4.5.2 From S38 to energy: the \golden" hybrid events . 72 4.5.3 From S38 to energy . 73 5 Auger dataset of extremely high energy cosmic rays 75 5.1 Studies on the selection criteria . 76 5.1.1 PMT selection . 76 5.1.2 Station selection: misidentified non-triggered stations . 78 5.1.3 Event selection . 81 5.2 Detector aging and its impact on event reconstruction . 82 5.2.1 Shape of the signal and effect on the trigger . 82 5.2.2 Effects on the event rate . 84 5.2.3 Effects on the highest energy dataset . 85 5.3 Energy and angular resolutions . 86 5.3.1 Angular resolution . 86 5.3.2 Energy resolution . 87 5.4 The final dataset and consistency checks . 90 iv CONTENTS 6 Anisotropy studies 99 6.1 Context of the study . 100 6.2 Analysis strategy . 102 6.3 Autocorrelation . 104 6.4 Cross correlation with astrophysical catalogues . 105 6.4.1 Catalogues . 105 6.4.2 Cross-correlation with flux-limited samples . 106 6.4.3 Cross-Correlation with Bright AGNs . 108 6.5 Single-source correlation (CenA) . 111 6.6 Search for correlation with neutrino \catalogues" . 113 6.6.1 Analysis strategy . 113 6.6.2 Magnetic fields deflections . 114 6.6.3 Preliminary results . 116 6.7 Conclusions . 118 7 Multi-variate analysis for mass discrimination 121 7.1 The sample of simulated events . 122 7.2 Mass-sensitive observables . 123 7.2.1 Longitudinal development . 123 7.2.2 Lateral development . 126 7.2.3 Muon content . 129 7.3 Multi-variate analysis . 135 7.3.1 Methods training and fine-tuning . 139 7.4 Conclusions . 142 8 An application of MVA: correlations test with candidate proton events145 8.1 The proton-enhanced datasets . 145 8.1.1 Cross-check with the Universality method . 147 8.2 Correlation tests . 149 8.3 Conclusions . 150 Appendices A Multi Variate analysis parameters and results plot 157 A.1 Fine-tuned parameters . 157 Bibliography 159 v vi Chapter 1 Cosmic Rays Omnia qua propter debent per.