Detection of cosmic rays in the SKALTA experiment Marek Bombara (P. J. Šafárik University Košice), Kysak, August 2011 1
Outline Cosmic rays Detection of cosmic rays SKALTA experiment SKALTA as a barometer Student work proposals 2
Cosmic rays Cosmic rays are natural source of particles with high energy. Cosmic rays consist of high-energy particles (primary) coming to Earth from outer space and a shower of secondary particles created in the earth atmosphere as a consequence of interaction of primary particle with the atmosphere. 3
Cosmic rays The composition of charged primary particles varies with energy, however, about 86% primary particles are protons, 11% alpha particles (helium nuclei), 1% heavier nuclei and 2% electrons. Neutral primary particles are composed of photons, neutrinos and antineutrinos. Energy spectrum of cosmic rays reaches values between 10 20 10 21 ev (for comparison most powerful accelerators on the Earth can accelerate particles up to ~10 12 ev). The spectrum shows deviations like knee and ankle whose origins are still not sufficiently explained. A big mystery remains the existence of particles with energy above 10 20 ev. Those particles should not exist in such a big amount as it is observed by cosmic ray experiments. Energy spectrum of primary cosmic rays 4
Detection of cosmic rays direct detection - by satellites, observing mostly primary cosmic rays dense network of ground detectors (e.g. Pierre Auger Observatory). There can be thousands of detectors distributed over big area (~ 1000 km 2 ), detecting whole shower of secondary cosmic rays. sparse network of ground detectors distributed over huge area (even millions km 2 ), which can not detect whole shower but can detect some very rare events occurring in Solar System, for example fission of the nucleus into two fragments both later creating shower in the Earth atmosphere - this can be studied by correlation between two working stations (in the distance ~ 1000 km). Each detecting method is thus focused on different type of physics. 5
Shower absorption in different materials LHC: ECAL CMS - PbWO4 Atmosphere as a calorimeter γ 6
SKALTA experiment The SKALTA (SlovaKiAn Large-area coincidence Time Array) experiment measures secondary cosmic rays originating from an interaction of primary cosmic ray particle with atmosphere. The working station is composed of three scintillation detectors each with the dimensions 60x60 cm. All three detectors are connected in coincidence. The detectors are arranged into a triangle with side length of 10 m. The area of the triangle defines a minimal size of the shower and therefore a minimal energy of original primary particle (> 1014 ev). Primary Cosmic Ray Atmosphere Secondary Cosmic Rays Shower Front Earth 10 m GPS Detectors Site Building ALTA Electronics PC Using the time difference among the signals from the detectors the point in the sky can be localized (up to a certain resolution) from which the original primary particle came from. By the measuring of the exact time by GPS (Global Positioning System) the data from other workstations (i.e. CZELTA in Czech Republic) can be compared and long distance correlations can be studied. 7
SKALTA as a barometer 4000 Pocet udalosti za den meranych experimentom SKALTA od 30.6.2010 do 17.1.2011 3500 SKALTA events per day huc Entries 147 Mean 116.4 RMS 55.71 3000 2500 2000 1500 1100 1000 Priemerny denny tlak vzduchu od 30.6.2010 do 17.1.2011 na kosickom letisku 1080 1060 1040 Average atmospheric pressure at Košice Airport per day htc Entries 147 Mean 115.4 RMS 56.32 1020 1000 980 960 940 920 900 Jul10 Aug10 Sep10 Okt10 Nov10 Dec10 Jan11 We can see anti-correlation by eye, even the distance between the locations (SKALTA and airport) is about 7 km. 8
SKALTA as a barometer We can quantitatively estimate the dependence of the cosmic ray flux on atmospheric pressure as: N = N 0 *e α*dp/p where N 0 is average number of detected events, dp = p - p 0, p is actual pressure and p 0 is average pressure at the place measured over very long period. By obtaining parameter α from the data we can measure (average for some time period) pressure p, if we know N. It means the SKALTA can be used as a barometer. 9
Student work proposals comparison of number of SKALTA events with sun activity (Exercise 1) comparison of number of SKALTA events with weather (Exercise 2) comparison of number of SKALTA events with numbers of Sunspots search in data for events incoming in very short time after each other (~seconds) - possible candidates for gamma ray bursts. Is SKALTA (or other similar stations) capable of supernova explosion detection? Czelta sites management http://czelta.utef.cvut.cz/graphs.php?site%5b%5d=praha_utef&site%5b%5d=sk_kosice_upjs&do=counts&from=31.05... To find source of correlation among stations. Example: Stav stanic - Dení!ek stanic - Download dat - Download meteodat - Nastavení - Po!ty událostí - Histogram - Sm"ry událostí - Odhla# se [guest]!"#$%&'()*"+$,! "#!$%&'!()!*%&!+,((!,+#--#,,!./0!1!/2&'!+,!*%&!+,((!,+#--#,,!./0 %&345&43 672%84#9:);,<=:+(>? 04&@2!8=8@A7!BC2!8BCDE%!3F@!&FB8FG!!"#$%&'()*+#H!IJKB27L&M=!F!&DECN=!7K!BC28L7!OF8LG4P@4!&F!(,*"-.H 1 of 2 20/6/11 12:47 10
Exercise 1: Is the cosmic ray measured by SKALTA (> 10 14 ev) affected by Sun? (Or alternatively, does it come from Sun?) Are the arrival directions of highenergy cosmic rays isotropic? Proposed analysis method: Compare numbers of events (showers of secondary cosmic rays) measured by SKALTA during the day and during the night. For example: measurement for day could be noon ± 3 hours (just to be sure the Sun was shining all that time) and night measurement could be midnight ± 3 hours. The outcome of the analysis should be a plot with (day and night) event number as a function of time (minimal time unit in the plot will be 24 hours). 11
Exercise 2: What is the relation between secondary cosmic ray flux and air temperature and density? Primary Cosmic Ray Atmosphere Shower Front T, ρ, p 10 m Detectors Secondary Cosmic Rays GPS Earth Site Building ALTA Proposed analysis method: Electronics PC The outcome of the analysis should be 3 plots with SKALTA event number, air temperature and air density as a function of time (minimal time unit in the plot will be 24 hours). Make a comparison among them and try to interpret the result. The measurement of the air temperature and air density were obtained from Hydrometeorological Institute of Slovak Academy of Sciences in Košice (In fact, air density was calculated by students from atmospheric pressure measurement.). 12
Data analysis: two ways 1.way: by web interface using for example MS Excel, or paper, pen and calculator... 2. way: by raw data using some programming language (for example C/C++, Python, Pascal, etc.). x 0 0 0 0 0 0 0.0 0 0 0 0 0 0 0.0 0.0 0.0 0.0 a 2011 03 24 15 05 45 645281133.9 2430 3164 4000 433 200 142 29.5 25.5 31.0 46.5 c 2011 03 24 15 05 51 673260.5 0 2322 3961 712 991 941 29.5 25.5 31.0 46.5 a 2011 03 24 15 05 55 854942857.6 1295 3119 4008 189 204 410 29.5 26.0 31.0 46.5 a 2011 03 24 15 05 56 657416869.7 1263 3118 4025 544 480 548 29.5 25.5 31.0 46.5 c 2011 03 24 15 06 51 674999.5 4095 2076 4017 776 1095 1006 29.5 25.5 31.0 46.5 a 2011 03 24 15 07 16 96301373.2 2140 3231 4007 576 213 125 29.5 25.5 30.5 46.5 a 2011 03 24 15 07 38 803003570.3 1634 3012 4000 1019 502 434 29.5 25.5 30.5 46.5 a 2011 03 24 15 07 44 893916939.4 2037 3279 4011 546 190 132 29.5 25.5 30.5 46.5 c 2011 03 24 15 07 51 674080.6 0 2082 3969 673 1062 1008 29.5 25.0 30.5 46.5 c 2011 03 24 15 08 51 670025.5 4095 2052 4012 793 1001 944 29.5 25.0 30.5 46.5 c 2011 03 24 15 09 51 670322.6 0 1908 3984 782 1077 977 29.5 25.0 30.5 46.5 13
Caveats during analysis (to correct for) Power cut? a (not c nor x) is correct... x 0 0 0 0 0 0 0.0 0 0 0 0 0 0 0.0 0.0 0.0 0.0 a 2011 03 24 15 05 45 645281133.9 2430 3164 4000 433 200 142 29.5 25.5 31.0 46.5 c 2011 03 24 15 05 51 673260.5 0 2322 3961 712 991 941 29.5 25.5 31.0 46.5 a 2011 03 24 15 05 55 854942857.6 1295 3119 4008 189 204 410 29.5 26.0 31.0 46.5 a 2011 03 24 15 05 56 657416869.7 1263 3118 4025 544 480 548 29.5 25.5 31.0 46.5 c 2011 03 24 15 06 51 674999.5 4095 2076 4017 776 1095 1006 29.5 25.5 31.0 46.5 a 2011 03 24 15 07 16 96301373.2 2140 3231 4007 576 213 125 29.5 25.5 30.5 46.5 a 2011 03 24 15 07 38 803003570.3 1634 3012 4000 1019 502 434 29.5 25.5 30.5 46.5 a 2011 03 24 15 07 44 893916939.4 2037 3279 4011 546 190 132 29.5 25.5 30.5 46.5 c 2011 03 24 15 07 51 674080.6 0 2082 3969 673 1062 1008 29.5 25.0 30.5 46.5 c 2011 03 24 15 08 51 670025.5 4095 2052 4012 793 1001 944 29.5 25.0 30.5 46.5 c 2011 03 24 15 09 51 670322.6 0 1908 3984 782 1077 977 29.5 25.0 30.5 46.5 14