The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows we find upper limits at 95% confidence level on an increase in VHE gamma-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges respectively. This puts upper limits on the flux increase during a radio GP of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.}

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the horizon or a few degrees below across an azimuthal range of about 80 degrees. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in the Earth crust or the surrounding ocean. In this paper we show how such air showers can be discriminated from the background of very inclined hadronic showers by using Monte Carlo simulations. Taking into account the orography of the site, the point source acceptance and the event rates expected have been calculated for a sample of generic neutrino fluxes from photo-hadronic interactions in AGNs. The analysis of about 30 hours of data taken towards the sealeads to a 90\% C.L. point source limit for tau neutrinos in the energy range from $1.0 \times 10^{15}$ eV to $3.0 \times 10^{18}$ eV of about $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 2.0 \times 10^{-4}$ GeV cm$^{-2}$ s$^{-1}$ for an assumed power-law neutrino spectrum with spectral index $\gamma$=-2. However, with 300 hours and in case of an optimistic neutrino flare model, limits of the level down to $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 8.4 \times 10^{-6}$ GeV cm$^{-2}$ s$^{-1}$ can be expected.

The dwarf spheroidal galaxy Ursa Major II (UMaII) is believed to be one of the most dark-matter dominated systems among the Milky Way satellites and represents a suitable target for indirect dark matter (DM) searches. The MAGIC telescopes carried out a deep observation campaign on UMaII between 2014 and 2016, collecting almost one hundred hours of good-quality data. This campaign enlarges the pool of DM targets observed at very high energy (E$\gtrsim$50GeV) in search for signatures of dark matter annihilation in the wide mass range between $\sim$100 GeV and $\sim$100 TeV. To this end, the data are analyzed with the full likelihood analysis, a method based on the exploitation of the spectral information of the recorded events for an optimal sensitivity to the explored dark matter models. We obtain constraints on the annihilation cross-section for different channels that are among the most robust and stringent achieved so far at the TeV mass scale from observations of dwarf satellite galaxies.

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 {\sigma} 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 time scale 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.

The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).

MAGIC (Major Atmospheric Gamma Imaging Cherenkov) is a system of two 17 m diameter, F/1.03 Imaging Atmospheric Cherenkov Telescopes (IACT). They are dedicated to the observation of gamma rays from galactic and extragalactic sources in the very high energy range (VHE, 30 GeV to 100 TeV). This submission contains links to the proceedings for the 35th International Cosmic Ray Conference (ICRC2017), held in Bexco, Busan, Korea from the 12th to the 17th of July, 2017.

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the MAGIC collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent~CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, found to worsen the statistical limits by about 36--42\%. Our constraints would have resulted much more competitive if the intrinsic pulse shape of the pulsar between 200GeV and 400GeV was understood in sufficient detail and allowed inclusion of events well below 400GeV.

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $\Gamma=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor ~1.7 (Reduced HV settings) with respect to standard settings (Nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15-30% and between 60 and 80% for Nominal HV, Reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from ~18%, under dark nights only, to up to ~40% in total with only moderate performance degradation.

The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.

It is widely believed that the bulk of the Galactic cosmic rays are accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between December 2014 and October 2016 we observed Cas A with the MAGIC telescopes, accumulating 158 hours of good-quality data. We derived the spectrum of the source from 100 GeV to 10 TeV. We also analysed $\sim$8 years of $Fermi$-LAT to obtain the spectral shape between 60 MeV and 500 GeV. The spectra measured by the LAT and MAGIC telescopes are compatible within the errors and show a clear turn off (4.6 $\sigma$) at the highest energies, which can be described with an exponential cut-off at $E_c = 3.5\left(^{+1.6}_{-1.0}\right)_{\textit{stat}} \left(^{+0.8}_{-0.9}\right)_{\textit{sys}}$ TeV. The gamma-ray emission from 60 MeV to 10 TeV can be attributed to a population of high-energy protons with spectral index $\sim$2.2 and energy cut-off at $\sim$10 TeV. This result indicates that Cas A is not contributing to the high energy ($\sim$PeV) cosmic-ray sea in a significant manner at the present moment. A one-zone leptonic model fails to reproduce by itself the multi-wavelength spectral energy distribution. Besides, if a non-negligible fraction of the flux seen by MAGIC is produced by leptons, the radiation should be emitted in a region with a low magnetic field (B$\lessapprox$100$\mu$G) like in the reverse shock.)

The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab Nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for several nights during the period June 18-27, for more than 10 hours. One hour of observation was conducted simultaneously to a giant 22 GHz radio flare and a hint of signal at GeV energies seen by Fermi-LAT. The MAGIC observations did not show significant emission in any of the analysed time intervals. The derived flux upper limit, in the energy range 200--1250 GeV, is 4.8$\times 10^{-12}$ ph cm$^{-2}$ s$^{-1}$. We estimate the gamma-ray opacity during the flaring period, which along with our non-detection, points to an inefficient acceleration in the V404\,Cyg jets if VHE emitter is located further than $1\times 10^{10}$ cm from the compact object.

B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy {\gamma} rays in the inverse Compton process. We performed extensive observations with the MAGIC telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region, or bow shock nebula. No significant steady very high energy {\gamma}-ray emission was found. We derived a 95% confidence level upper limit of 3.0 x 10 -12 cm -2 s -1 on the average {\gamma}-ray emission from the binary system above 200 GeV. The upper limits obtained with MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below ~(2-10)% of the pulsar spin down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies.

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures to map the faint Cherenkov light emitted in extensive air showers onto their image sensors. Segmented reflectors fulfill these needs using mass produced and light weight mirror facets. However, as the overall image is the sum of the individual mirror facet images, alignment is important. Here we present a method to determine the mirror facet positions on a segmented reflector in a very direct way. Our method reconstructs the mirror facet positions from photographs and a laser distance meter measurement which goes from the center of the image sensor plane to the center of each mirror facet. We use our method to both align the mirror facet positions and to feed the measured positions into our IACT simulation. We demonstrate our implementation on the 4 m First Geiger-mode Avalanche Cherenkov Telescope (FACT).

The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

We present the first detection of the nearby (z=0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS, helps to expand the small population of known TeV BL Lacs. The source was observed with the MAGIC telescopes between 2010 April and 2011 May, collecting 46 h of good quality data. These observations led to the detection of the source at 6.0 $\sigma$ confidence level, with a steady flux $\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm 2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential spectral photon index $\Gamma = 2.4 \pm 0.2_{\mathrm{stat}} \pm 0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the broad-band spectral energy distribution (SED) simultaneous with MAGIC observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests values for the SSC parameters that are quite common among known TeV BL Lacs except for a relatively low Doppler factor and slope of electron energy distribution. A thermal feature seen in the SED is well matched by a giant elliptical's template. This appears to be the signature of thermal emission from the host galaxy, which is clearly resolved in optical observations.

The high-energy non-thermal universe is dominated by power law-like spectra. Therefore results in high-energy astronomy are often reported as parameters of power law fits, or, in the case of a non-detection, as an upper limit assuming the underlying unseen spectrum behaves as a power law. In this paper I demonstrate a simple and powerful one-to-one relation of the integral upper limit in the two dimensional power law parameter space into the spectrum parameter space and use this method to unravel the so far convoluted question of the sensitivity of astroparticle telescopes.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Context. We present the results of a multi-year monitoring campaign of the Galactic Center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; $\geq$100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 hours of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 year observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio Arc.

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.

For studying variable sources like blazars, it is crucial to achieve unbiased monitoring, either with dedicated telescopes in pointing mode or survey instruments. At TeV energies, the High Altitude Water Cherenkov (HAWC) observatory monitors approximately two thirds of the sky every day. It uses the water Cherenkov technique, which provides an excellent duty cycle independent of weather and season. The First G-APD Cherenkov Telescope (FACT) monitors a small sample of sources with better sensitivity, using the imaging air Cherenkov technique. Thanks to its camera with silicon-based photosensors, FACT features an excellent detector performance and stability and extends its observations to times with strong moonlight, increasing the duty cycle compared to other imaging air Cherenkov telescopes. As FACT and HAWC have overlapping energy ranges, a joint study can exploit the longer daily coverage given that the observatories' locations are offset by 5.3 hours. Furthermore, the better sensitivity of FACT adds a finer resolution of features on hour-long time scales, while the continuous duty cycle of HAWC ensures evenly sampled long-term coverage. Thus, the two instruments complement each other to provide a more complete picture of blazar variability. In this presentation, the first joint study of light curves from the two instruments will be shown, correlating long-term measurements with daily sampling between air and water Cherenkov telescopes. The presented results focus on the study of the variability of the bright blazars Mrk 421 and Mrk 501 during the last two years featuring various flaring activities.

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC, are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray sources detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broad band emission can be modeled in the framework of a two zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment remains a challenge. Here we present a simple, yet extendable method, to align a segmented reflector using its Bokeh. Bokeh alignment does not need a star or good weather nights but can be done even during daytime. Bokeh alignment optimizes the facet orientations by comparing the segmented reflectors Bokeh to a predefined template. The optimal Bokeh template is highly constricted by the reflector's aperture and is easy accessible. The Bokeh is observed using the out of focus image of a near by point like light source in a distance of about 10 focal lengths. We introduce Bokeh alignment on segmented reflectors and demonstrate it on the First Geiger-mode Avalanche Cherenkov Telescope (FACT) on La Palma, Spain.

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and as they are composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment is a challenge. Here we present a computer vision based star tracking alignment method, which also works for limited or changing star light visibility. Our method normalizes the mirror facet reflection intensities to become independent of the reference star's intensity or the cloud coverage. Using two CCD cameras, our method records the mirror facet orientations asynchronously of the telescope drive system, and thus makes the method easy to integrate into existing telescopes. It can be combined with remote facet actuation, but does not require one to work. Furthermore, it can reconstruct all individual mirror facet point spread functions without moving any mirror. We present alignment results on the 4 meter First Geiger-mode Avalanche Cherenkov Telescope (FACT).

The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband emission shows a periodicity of $\sim$26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray observations. We present a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long timescale signature in the VHE emission. We focus on the search for super-orbital modulation of the VHE peak and on the search for correlations between TeV emission and optical determination of the extension of the circumstellar disk. A four-year campaign has been carried out by MAGIC. The source was observed during the orbital phases when the periodic VHE outbursts have occurred ($\phi$=0.55-0.75). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61$^{\circ}$303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past ($\phi$=0.75-1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases $\phi$=0.5--0.75 was found to show yearly variability consistent with the $\sim$4.5 years long-term modulation found in the radio band. This modulation of the TeV flux can be well described by a sine function with the best fit period of $1610\pm 58$ days. The complete dataset span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star but this may be affected by the strong, short timescale (as short as intra-day) variation displayed by the H$\alpha$ fluxes.

The BL Lac object 1ES 1011+496 was discovered at Very High Energy gamma-rays by MAGIC in spring 2007. Before that the source was little studied in different wavelengths. Therefore a multi-wavelength (MWL) campaign was organized in spring 2008. Along MAGIC, the MWL campaign included the Metsahovi radio observatory, Bell and KVA optical telescopes and the Swift and AGILE satellites. MAGIC observations span from March to May, 2008 for a total of 27.9 hours, of which 19.4 hours remained after quality cuts. The light curve showed no significant variability. The differential VHE spectrum could be described with a power-law function. Both results were similar to those obtained during the discovery. Swift XRT 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 SED shows a synchrotron dominated source, unlike concluded in previous work based on nonsimultaneous 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 Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 x 10^{-13} TeV cm^{-2} s^{-1} for the Geminga pulsar and 3.5 x 10^{-12} TeV cm^{-2} s^{-1} for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

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.

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.

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>100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6$^\circ$ 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 NGC 1275 spectrum, which is well described by a power law with a photon index of $3.6\pm0.2_{stat}\pm0.2_{syst}$ between 90 GeV and 1.2 TeV. We do not detect any diffuse gamma-ray emission from the cluster and set stringent constraints on its CR population. In order 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 below 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 <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\gtrsim5-8 \mu$G, which is below the 25 $\mu$G inferred from Faraday rotation measurements, whereas, for $\alpha\lesssim2.1$, the hadronic interpretation of the diffuse radio emission is in contrast with our gamma-ray flux upper limits independently of the magnetic field strength.

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.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, 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 (LAT), 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, Mets\"ahovi 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 LIDAR (LIght Detection And Ranging) 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 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) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton 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 On-Off problem, aka. Li-Ma problem, is a statistical problem where a measured rate is the sum of two parts. The first is due to a signal and the second due to a background, both of which are unknown. Mostly frequentist solutions are being used that are only adequate for high count numbers. When the events are rare such an approximation is not good enough. Indeed, in high-energy astrophysics this is often the rule rather than the exception. I will present a universal objective Bayesian solution that depends only on the initial three parameters of the On-Off problem: the number of events in the "on" region, the number of events in the "off" region, and their ratio-of-exposure. With a two-step approach it is possible to infer the signal's significance, strength, uncertainty or upper limit in a unified a way. The approach is valid without restrictions for any count number including zero and may be widely applied in particle physics, cosmic-ray physics and high-energy astrophysics. I apply the method to Gamma Ray Burst data.

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV {\gamma}-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the {\gamma}-ray spectrum at TeV energies. Aims. We aim to explore the very-high-energy domain to search for {\gamma}-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, triggered by optical and subsequent GeV {\gamma}-ray detec- tions. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of con- temporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp<~0.15 Le.

A recent report on the over saturation in SiPMs is puzzling. The measurements, using a variety of SiPMs, show an excess in signal far beyond the physical limit of the number of SiPM microcells without indication of an ultimate saturation. In this work I propose a solution to this problem. Different measurements and theoretical models of avalanche propagation indicate that multiple simultaneous primary avalanches produce an ever narrower and faster signal. This is because of a speed-up of effective avalanche propagation processes. It means that SiPMs, operated at their saturation regime, should become faster the more light they detect. Therefore, signal extraction methods that use the amplitude of the signal should see an over saturation effect. Measurements with a commercial SiPM illuminated with bright picosecond pulses in the saturation regime demonstrate that indeed the rising edge of the SiPM signal gets faster as the light pulses get brighter. A signal extractor based on the amplitude shows a nonlinear behavior in comparison to an integrating charge extractor. This supports the proposed solution for the over saturation effect. Furthermore I show that this effect can already be seen with a bandwidth of 300MHz, which means that it should be taken into account for fast sampling experiments.

The First Geiger-mode Avalanche photodiode (G-APD) Cherenkov Telescope (FACT) has been operating on the Canary island of La Palma since October 2011. Operations were automated so that the system can be operated remotely. Manual interaction is required only when the observation schedule is modified due to weather conditions or in case of unexpected events such as a mechanical failure. Automatic operations enabled high data taking efficiency, which resulted in up to two terabytes of FITS files being recorded nightly and transferred from La Palma to the FACT archive at ISDC in Switzerland. Since long term storage of hundreds of terabytes of observations data is costly, data compression is mandatory. This paper discusses the design choices that were made to increase the compression ratio and speed of writing of the data with respect to existing compression algorithms. Following a more detailed motivation, the FACT compression algorithm along with the associated I/O layer is discussed. Eventually, the performances of the algorithm is compared to other approaches.

The First G-APD Cherenkov Telescope (FACT) was built on the Canary Island of La Palma in October 2011 as a proof of principle for silicon based photosensors in Cherenkov Astronomy. The scientific goal of the project is to study the variability of active galatic nuclei (AGN) at TeV energies. Observing a small sample of TeV blazars whenever possible, an unbiased data sample is collected. This allows to study the variability of the selected objects on timescales from hours to years. Results from the first three years of monitoring will be presented. To provide quick flare alerts to the community and trigger multi-wavelength observations, a quick look analysis has been installed on-site providing results publicly online within the same night. In summer 2014, several flare alerts were issued. Results of the quick look analysis are summarized.