Unfolding of differential energy spectra in the MAGIC experiment

The paper describes the different methods, used in the MAGIC experiment, to unfold experimental energy distributions of cosmic ray particles (gamma-rays). Questions and problems related to the unfolding are discussed. Various procedures are proposed which can help to make the unfolding robust and reliable. The different methods and procedures are implemented in the MAGIC software and are used in most of the analyses.Comment: Submitted to NIM

Implementation of the Random Forest Method for the Imaging Atmospheric Cherenkov Telescope MAGIC

The paper describes an application of the tree classification method Random Forest (RF), as used in the analysis of data from the ground-based gamma telescope MAGIC. In such telescopes, cosmic gamma-rays are observed and have to be discriminated against a dominating background of hadronic cosmic-ray particles. We describe the application of RF for this gamma/hadron separation. The RF method often shows superior performance in comparison with traditional semi-empirical techniques. Critical issues of the method and its implementation are discussed. An application of the RF method for estimation of a continuous parameter from related variables, rather than discrete classes, is also discussed.Comment: 16 pages, 8 figure

Probing quantum gravity using photons from a flare of the active galactic nucleus Markarian 501 observed by the MAGIC telescope

We analyze the timing of photons observed by the MAGIC telescope during a flare of the active galactic nucleus Mkn 501 for a possible correlation with energy, as suggested by some models of quantum gravity (QG), which predict a vacuum refractive index \simeq 1 + (E/M_{QGn})^n, n = 1,2. Parametrizing the delay between gamma-rays of different energies as \Delta t =\pm\tau_l E or \Delta t =\pm\tau_q E^2, we find \tau_l=(0.030\pm0.012) s/GeV at the 2.5-sigma level, and \tau_q=(3.71\pm2.57)x10^{-6} s/GeV^2, respectively. We use these results to establish lower limits M_{QG1} > 0.21x10^{18} GeV and M_{QG2} > 0.26x10^{11} GeV at the 95% C.L. Monte Carlo studies confirm the MAGIC sensitivity to propagation effects at these levels. Thermal plasma effects in the source are negligible, but we cannot exclude the importance of some other source effect.Comment: 12 pages, 3 figures, Phys. Lett. B, reflects published versio

The Great Markarian 421 Flare of 2010 February: Multiwavelength Variability and Correlation Studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of 2010 February, when an extraordinary flare reaching a level of ∼27 Crab Units above 1 TeV was measured in very high energy (VHE) γ-rays with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory. This is the highest flux state for Mrk 421 ever observed in VHE γ-rays. Data are analyzed from a coordinated campaign across multiple instruments, including VHE γ-ray (VERITAS, Major Atmospheric Gamma-ray Imaging Cherenkov), high-energy γ-ray (Fermi-LAT), X-ray (Swift, Rossi X-ray Timing Experiment, MAXI), optical (including the GASP-WEBT collaboration and polarization data), and radio (Metsahovi, Owens Valley Radio Observatory, University of Michigan Radio Astronomy Observatory). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare "decline" epochs. The main flare statistics allow 2 minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of ∼25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ ⪆ 33) and the size of the emission region (δ-1RB≲ 3.8 × 1013cm) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10 minute binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship, from linear to quadratic to lack of correlation to anticorrelation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain with the classic single-zone synchrotron self-Compton model.</p