The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. The telescopes are designed to measure Cherenkov light from air showers initiated by gamma rays in the energy regime from around 50 GeV to more than 50 TeV. The two telescopes were built in 2004 and 2009, respectively, with different cameras, triggers and readout systems. In the years 2011–2012 the MAGIC collaboration undertook a major upgrade to make the stereoscopic system uniform, improving its overall performance and easing its maintenance. In particular, the camera, the receivers and the trigger of the first telescope were replaced and the readout of the two telescopes was upgraded. This paper (Part I) describes the details of the upgrade as well as the basic performance parameters of MAGIC such as raw data treatment, linearity in the electronic chain and sources of noise. In Part II, we describe the physics performance of the upgraded system.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

Context. 1ES 1011+496 (z = 0.212) was discovered in very high-energy (VHE, E > 100 GeV) gamma rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to an incomplete characterization of the broadband spectral energy distribution (SED). Aims. We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron self-Compton (SSC) scenario is able to explain the observed broadband spectrum.Methods. We analyzed data in the range from VHE to radio data from 2011 and 2012 collected by MAGIC, Fermi-LAT, Swift, KVA, OVRO, and Metsahovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE.Results. The VHE spectrum was fit with a simple power law with a photon index of 3.69 +/- 0.22 and a flux above 150 GeV of (1.46 +/- 0.16) x 10(-11) ph cm(-2) s(-1). The source 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of 1.8 +/- 0.4 c, which is the highest speed statistically significant measured so far in a high-frequency-peaked BL Lac.Conclusions. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV, which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field with respect to the optical emission.

Aims. We investigate the extension of the very high-energy spectral tail of the Crab Pulsar at energies above 400 GeV.Methods. We analyzed similar to 320 h of good-quality Crab data obtained with the MAGIC telescope from February 2007 to April 2014.Results. We report the most energetic pulsed emission ever detected from the Crab Pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with those measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10GeV by the Large Area Telescope (LAT) on board the Fermi satellite. When making a joint fit of the LAT and MAGIC data above 10 GeV the photon indices of the spectra differ by 0.5 +/- 0.1.Conclusions. Using data from the MAGIC telescopes we measured the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least 5x10(6). These results strongly suggest IC scattering off low-energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse gamma-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E greater than or similar to 100GeV) for a total of 253 h from 2009 to 2014. The active nuclei of NGC1275, the central dominant galaxy of the cluster, and IC310, lying at about 0.6 degrees from the centre, have been detected as point-like VHE gamma-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC1275 spectrum, which is described well by a power law with a photon index Gamma = 3.6 +/- 0.2(stat) +/- 0.2(syst) between 90GeV and 1200 GeV. We do not detect any diffuse gamma-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes alpha. For alpha = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be less than or similar to 1-2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to less than or similar to 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B-0, that depend on the CR distribution. For alpha = 2.2, B-0 greater than or similar to 5-8 mu G, which is below the similar to 25 mu G inferred from Faraday rotation measurements, whereas for alpha less than or similar to 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our gamma-ray flux upper limits independently of the magnetic field strength.

The BL Lac object 1ES 1011+496 was discovered at very high energy (VHE, E > 100GeV) γ -rays by Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) in Spring 2007. Before that the source was little studied in different wavelengths. Therefore, a multiwavelength (MWL) campaign was organized in Spring 2008. Along MAGIC, the MWL campaign included the Metsähovi Radio Observatory, Bell and Kungliga Vetenskapsakademien (KVA) optical telescopes and the Swift and AGILE satellites. MAGIC observations span from 2008 March to May for a total of 27.9 h, of which 19.4 h remained after quality cuts. The light curve showed no significant variability yielding an integral flux above 200 GeV of (1.3 ± 0.3) × 10−11 photons cm−2 s−1. The differential VHE spectrum could be described with a power-law function with a spectral index of 3.3 ± 0.4. Both results were similar to those obtained during the discovery. Swift X-ray Telescope observations revealed an X-ray flare, characterized by a harder-when-brighter trend, as is typical for high synchrotron peak BL Lac objects (HBL). Strong optical variability was found during the campaign, but no conclusion on the connection between the optical and VHE γ -ray bands could be drawn. The contemporaneous spectral energy distribution shows a synchrotron-dominated source, unlike concluded in previous work based on non-simultaneous data, and is well described by a standard one-zone synchrotron self-Compton model. We also performed a study on the source classification. While the optical and X-ray data taken during our campaign show typical characteristics of an HBL, we suggest, based on archival data, that 1ES 1011+496 is actually a borderline case between intermediate and high synchrotron peak frequency BL Lac objects.

We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi -LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV—the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi -LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

In February-March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z=0.212) in flaring state at very-high energy (VHE, E>100GeV). The flux reached a level more than 10 times higher than any previously recorded flaring state of the source. We present the description of the characteristics of the flare presenting the light curve and the spectral parameters of the night-wise spectra and the average spectrum of the whole period. From these data we aim at detecting the imprint of the Extragalactic Background Light (EBL) in the VHE spectrum of the source, in order to constrain its intensity in the optical band. For this we implement the method developed by the H.E.S.S. collaboration in which the intrinsic energy spectrum of the source is modeled with a simple function, and the EBL-induced optical depth is calculated using a template EBL model. The likelihood of the observed spectrum is then maximized, including a normalization factor for the EBL opacity among the free parameters. From the data collected differential energy spectra was produced for all nights of the observed period. Evaluating the changes in the fit parameters we conclude that the spectral shape for most of the nights were compatible, regardless of the flux level, which enabled us to produce an average spectrum from which the EBL imprint could be constrained. The likelihood ratio test shows that the model with an EBL density 1.07(-0.20,+0.24)_stat+sys, relative to the one in the tested EBL template (Dominguez et al.2011), is preferred at the 4.6 sigma level to the no-EBL hypothesis, with the assumption that the intrinsic source spectrum can be modeled as a log-parabola. This would translate into a constraint of the EBL density in the wavelength range [0.24 um,4.25 um], with a peak value at 1.4 um of F=12.27_-2.29^+2.75 nW m^-2 sr^-1, including systematics.

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ~50GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220GeV is (0.66+/-0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the sigma of a 2-dimensional Gaussian distribution, at those energies is < 0.07 degree, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11-18% in flux normalization and +/-0.15 for the energy spectrum power-law slope.

The BL Lac object 1ES 1011+ 496 was discovered at very high energy (VHE, E > 100GeV) gamma-rays by Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) in Spring 2007. Before that the source was little studied in different wavelengths. Therefore, a multiwave-length (MWL) campaign was organized in Spring 2008. Along MAGIC, the MWL campaign included the Mets " ahovi Radio Observatory, Bell and Kungliga Vetenskapsakademien (KVA) optical telescopes and the Swift and AGILE satellites. MAGIC observations span from 2008 March to May for a total of 27.9 h, of which 19.4 h remained after quality cuts. The light curve showed no significant variability yielding an integral flux above 200GeVof (1.3 +/- 0.3) x10(-11) photons cm(-2) s(-1). The differential VHE spectrum could be described with a power-law function with a spectral index of 3.3 +/- 0.4. Both results were similar to those obtained during the discovery. Swift X-ray Telescope observations revealed an X-ray flare, characterized by a harder-when-brighter trend, as is typical for high synchrotron peak BL Lac objects (HBL). Strong optical variability was found during the campaign, but no conclusion on the connection between the optical and VHE gamma-ray bands could be drawn. The contemporaneous spectral energy distribution shows a synchrotron-dominated source, unlike concluded in previous work based on non-simultaneous data, and is well described by a standard one-zone synchrotron self-Compton model. We also performed a study on the source classification. While the optical and X-ray data taken during our campaign show typical characteristics of an HBL, we suggest, based on archival data, that 1ES 1011+496 is actually a borderline case between intermediate and high synchrotron peak frequency BL Lac objects.

The MAGIC stereoscopic system collected 69 hours of Crab Nebula data between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereoscopic software provided an unprecedented precision of spectral and night-by-night light curve determination at gamma rays. We derived a differential spectrum with a single instrument from 50 GeV up to almost 30 TeV with 5 bins per energy decade. At low energies, MAGIC results, combined with Fermi-LAT data, show a flat and broad Inverse Compton peak. The overall fit to the data between 1 GeV and 30 TeV is not well described by a log-parabola function. We find that a modified log-parabola function with an exponent of 2.5 instead of 2 provides a good description of the data (χred2=35/26). Using systematic uncertainties of the MAGIC and Fermi-LAT measurements we determine the position of the Inverse Compton peak to be at (53±3stat+31syst−13syst) GeV, which is the most precise estimation up to date and is dominated by the systematic effects. There is no hint of the integral flux variability on daily scales at energies above 300 GeV when systematic uncertainties are included in the flux measurement. We consider three state-of-the-art theoretical models to describe the overall spectral energy distribution of the Crab Nebula. The constant B-field model cannot satisfactorily reproduce the VHE spectral measurements presented in this work, having particular difficulty reproducing the broadness of the observed IC peak. Most probably this implies that the assumption of the homogeneity of the magnetic field inside the nebula is incorrect. On the other hand, the time-dependent 1D spectral model provides a good fit of the new VHE results when considering a 80 μG magnetic field. However, it fails to match the data when including the morphology of the nebula at lower wavelengths.

The number of known very high energy (VHE) blazars is similar to 50, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the gamma-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy gamma-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES 0033+595 with a statistical significance of 5.5 sigma. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parametrized with a power law with an integral flux above 150 GeV of (7.1 +/- 1.3) x 10(-12) photons cm(-2) s(-1) and a photon index of (3.8 +/- 0.7). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED, we used simultaneous opticalR-band data from the KVA telescope, archival X-ray data by Swift as well as INTEGRAL, and simultaneous high-energy (HE, 300 MeV-10 GeV) gamma-ray data from the Fermi Large Area Telescope (LAT) observatory. Using the empirical approach of Prandini et al. (2010) and the Fermi LAT and MAGIC spectra for this object, we estimate the redshift of this source to be 0.34 +/- 0.08 +/- 0.05. This is a relevant result because this source is possibly one of the 10 most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.

PG 1553+113 is a very high energy (VHE, E > 100 GeV) gamma-ray emitter classified as a BL Lac object. Its redshift is constrained by intergalactic absorption lines in the range 0.4 < z < 0.58. The MAGIC telescopes have monitored the source's activity since 2005. In early 2012, PG 1553+113 was found in a high state, and later, in April of the same year, the source reached its highest VHE flux state detected so far. Simultaneous observations carried out in X-rays during 2012 April show similar flaring behaviour. In contrast, the gamma-ray flux at E < 100 GeV observed by Fermi-LAT is compatible with steady emission. In this paper, a detailed study of the flaring state is presented. The VHE spectrum shows clear curvature, being well fitted either by a power law with an exponential cut-off or by a log-parabola. A simple power-law fit hypothesis for the observed shape of the PG 1553+113 VHE gamma-ray spectrum is rejected with a high significance (fit probability P = 2.6 x 10(-6)). The observed curvature is compatible with the extragalactic background light (EBL) imprint predicted by current generation EBL models assuming a redshift z similar to 0.4. New constraints on the redshift are derived from the VHE spectrum. These constraints are compatible with previous limits and suggest that the source is most likely located around the optical lower limit, z = 0.4, based on the detection of Ly alpha absorption. Finally, we find that the synchrotron self-Compton model gives a satisfactory description of the observed multiwavelength spectral energy distribution during the flare.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 2013 April 1 and August 10, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope, Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsahovi, and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a Light Detection and Ranging (LIDAR) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution (SED) between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) show evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton (SSC) model to five simultaneous broadband SEDs. We find that the SSC model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHEγ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations fromγ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded aγ-ray signal above 250 GeV of (3.7± 0.7) per cent of the Crab Nebula flux with a statistical significa nce of 9.9σ. The multiwavelength observations showed significant vari ability in essentially all energy bands, including a VHEγ-ray flare that lasted less than one night, which provided unp recedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of 2.97± 0.29 between∼150 GeV and 1 TeV and an integral flux of (9.3± 1.9) per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible w ith the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.

Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 10^4<γ< 6 × 10^5, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.

The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability timescale is estimated to be 6.4 ± 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad-line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.

Context: MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims: Our aim is to evaluate if the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods. We have observed MWC 656 with the MAGIC telescopes for $\sim$23 hours during two observation periods: between May and June 2012 and June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results: We have not detected the MWC 656 binary system at TeV energies with the MAGIC Telescopes in either of the two campaigns carried out. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of $\sim$5% of the Crab Nebula flux. The results obtained from the MAGIC observations do not support persistent emission of very high energy gamma rays from this system at a level of 2.4% the Crab flux.

Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. However, most of the gamma-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available gamma-ray instrumentation was easier to perform.Aims. Our goal is to characterize the source gamma-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity.Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign.Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density.Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

Aims. We perform an extensive characterization of the broadband emission of Mrk 421, as well as its temporal evolution, during the non-flaring (low) state. The high brightness and nearby location (z = 0.031) of Mrk 421 make it an excellent laboratory to study blazar emission. The goal is to learn about the physical processes responsible for the typical emission of Mrk 421, which might also be extended to other blazars that are located farther away and hence are more difficult to study.Methods. We performed a 4.5-month multi-instrument campaign on Mrk 421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. This extensive radio to very-high-energy (VHE; E > 100 GeV) gamma-ray dataset provides excellent temporal and energy coverage, which allows detailed studies of the evolution of the broadband spectral energy distribution.Results. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical /UV and X-rays extending over the duration of the campaign.Conclusions. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multiwavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk 421 during non-flaring activity. Such a temporally extended X-ray /VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.