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

First bounds on the high-energy emission from isolated Wolf-Rayet binary systems

High-energy gamma-ray emission is theoretically expected to arise in tight binary star systems (with high mass loss and high velocity winds), although the evidence of this relationship has proven to be elusive so far. Here we present the first bounds on this putative emission from isolated Wolf-Rayet (WR) star binaries, WR 147 and WR 146, obtained from observations with the MAGIC telescope.Comment: (Authors are the MAGIC Collaboration.) Manuscript in press at The Astrophysical Journal Letter

Systematic search for VHE gamma-ray emission from X-ray bright high-frequency BL Lac objects

All but three (M87, BL Lac and 3C 279) extragalactic sources detected so far at very high energy (VHE) gamma-rays belong to the class of high-frequency peaked BL Lac (HBL) objects. This suggested to us a systematic scan of candidate sources with the MAGIC telescope, based on the compilation of X-ray blazars by Donato et al. (2001). The observations took place from December 2004 to March 2006 and cover sources on the northern sky visible under small zenith distances zd < 30 degrees at culmination. The sensitivity of the search was planned for detecting X-ray bright F(1 keV) > 2 uJy) sources emitting at least the same energy flux at 200 GeV as at 1 keV. In order to avoid strong gamma-ray attenuation close to the energy threshold, the redshift of the sources was constrained to values z<0.3. Of the fourteen sources observed, 1ES 1218+304 and 1ES 2344+514 have been detected in addition to the known bright TeV blazars Mrk 421 and Mrk 501. A marginal excess of 3.5 sigma from the position of 1ES 1011+496 was observed and has been confirmed as a source of VHE gamma-rays by a second MAGIC observation triggered by a high optical state (Albert et al. 2007). For the remaining sources, we present here the 99% confidence level upper limits on the integral flux above ~200 GeV. We characterize the sample of HBLs (including all HBLs detected at VHE so far) by looking for correlations between their multi-frequency spectral indices determined from simultaneous optical, archival X-ray, and radio luminosities, finding that the VHE emitting HBLs do not seem to constitute a unique subclass. The absorption corrected gamma-ray luminosities at 200 GeV of the HBLs are generally not higher than their X-ray luminosities at 1 keV.Comment: 15 pages, 7 figures, 5 tables, submitted to ApJ (revised version

Fermi LAT Observations of LS I +61 303: First detection of an orbital modulation in GeV Gamma Rays

This Letter presents the first results from the observations of LSI +61 303 using Large Area Telescope data from the Fermi Gamma-Ray Space Telescope between 2008 August and 2009 March. Our results indicate variability that is consistent with the binary period, with the emission being modulated at 26.6 +/- 0.5 days. This constitutes the first detection of orbital periodicity in high-energy gamma rays (20 MeV-100 GeV, HE). The light curve is characterized by a broad peak after periastron, as well as a smaller peak just before apastron. The spectrum is best represented by a power law with an exponential cutoff, yielding an overall flux above 100 MeV of 0.82 +/- 0.03(stat) +/- 0.07(syst) 10^{-6} ph cm^{-2} s^{-1}, with a cutoff at 6.3 +/- 1.1(stat) +/- 0.4(syst) GeV and photon index Gamma = 2.21 +/- 0.04(stat) +/- 0.06(syst). There is no significant spectral change with orbital phase. The phase of maximum emission, close to periastron, hints at inverse Compton scattering as the main radiation mechanism. However, previous very high-energy gamma ray (>100 GeV, VHE) observations by MAGIC and VERITAS show peak emission close to apastron. This and the energy cutoff seen with Fermi suggest the link between HE and VHE gamma rays is nontrivial.Comment: 7 pages, 5 figures, accepted for publication in ApJ Letters 21 July 200

Fermi Observations of the Very Hard Gamma-ray Blazar PG 1553+113

We report the observations of PG 1553+113 during the first ~200 days of Fermi Gamma-ray Space Telescope science operations, from 4 August 2008 to 22 February 2009 (MJD 54682.7-54884.2). This is the first detailed study of PG 1553+113 in the GeV gamma-ray regime and it allows us to fill a gap of three decades in energy in its spectral energy distribution. We find PG 1553+113 to be a steady source with a hard spectrum that is best fit by a simple power-law in the Fermi energy band. We combine the Fermi data with archival radio, optical, X-ray and very high energy (VHE) gamma-ray data to model its broadband spectral energy distribution and find that a simple, one-zone synchrotron self-Compton model provides a reasonable fit. PG 1553+113 has the softest VHE spectrum of all sources detected in that regime and, out of those with significant detections across the Fermi energy bandpass so far, the hardest spectrum in that energy regime. Thus, it has the largest spectral break of any gamma-ray source studied to date, which could be due to the absorption of the intrinsic gamma-ray spectrum by the extragalactic background light (EBL). Assuming this to be the case, we selected a model with a low level of EBL and used it to absorb the power-law spectrum from PG 1553+113 measured with Fermi (200 MeV - 157 GeV) to find the redshift which gave the best fit to the measured VHE data (90 GeV - 1.1 TeV) for this parameterisation of the EBL. We show that this redshift can be considered an upper limit on the distance to PG 1553+113.Comment: Accepted for publication in the Astrophysical Journal (28 pages, 5 figures

Observations of Milky Way Dwarf Spheroidal galaxies with the Fermi-LAT detector and constraints on Dark Matter models

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)Comment: 25 pages, 4 figures, accepted to ApJ, Corresponding authors: J. Cohen-Tanugi, C. Farnier, T.E. Jeltema, E. Nuss, and S. Profum

MAGIC upper limits on the very high energy emission from GRBs

The fast repositioning system of the MAGIC Telescope has allowed during its first data cycle, between 2005 and the beginning of year 2006, observing nine different GRBs as possible sources of very high energy gammas. These observations were triggered by alerts from Swift, HETE-II, and Integral; they started as fast as possible after the alerts and lasted for several minutes, with an energy threshold varying between 80 and 200 GeV, depending upon the zenith angle of the burst. No evidence for gamma signals was found, and upper limits for the flux were derived for all events, using the standard analysis chain of MAGIC. For the bursts with measured redshift, the upper limits are compatible with a power law extrapolation, when the intrinsic fluxes are evaluated taking into account the attenuation due to the scattering in the Metagalactic Radiation Field (MRF).Comment: 25 pages, 9 figures, final version accepted by ApJ. Changet title to "MAGIC upped limits on the VERY high energy emission from GRBs", re-organized chapter with description of observation, removed non necessaries figures, added plot of effective area depending on zenith angle, added an appendix explaining the upper limit calculation, added some reference

Open and Hidden Charm Production in 920 GeV Proton-Nucleus Collisions

The HERA-B collaboration has studied the production of charmonium and open charm states in collisions of 920 GeV protons with wire targets of different materials. The acceptance of the HERA-B spectrometer covers negative values of xF up to xF=-0.3 and a broad range in transverse momentum from 0.0 to 4.8 GeV/c. The studies presented in this paper include J/psi differential distributions and the suppression of J/psi production in nuclear media. Furthermore, production cross sections and cross section ratios for open charm mesons are discussed.Comment: 5 pages, 9 figures, to be published in the proceedings of the 6th International Conference on Hyperons, Charm & Beauty Hadrons (BEACH04), Chicago, IL, June 27 - July 3, 200

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