• BL Lacertae is a blazar at the redshift of z = 0.069, eponym of the BL Lac blazar type. It is also a prototype of the low-frequency-peaked BL Lac (LBL) subclass. It was first detected in sub-TeV gamma-ray range by MAGIC in 2005. In 2015, MAGIC observations of BL Lacertae were triggered by the Fermi-LAT analysis report in the MAGIC group, and were performed during 10 individual nights between 15th and 28th June for a total of 8.6 h. The measured gamma-ray flux varied from 40\% to 10\% of the Crab nebula flux above 200 GeV in the nights from 15th to 17th June. In particular, a fast variability was found during the nights of 15th and 17th June. % with a halving time scale of (33+/-9) minutes, statistical errors only. We also performed multi-wavelength (MWL) observations in the radio, optical, UV, X-ray and gamma-ray bands, and %the MWL light curves indicate that the very high energy gamma-ray flare was a so-called orphan flare, %i.e. no simultaneous activity in other wavebands was observed. the MWL light curves indicate that no apparent simultaneous activity in other wavebands accompanying the very high energy gamma-ray flare in June 2015 like an another occurrence of an orphan very high energy flare. In this proceedings we will present the results of the campaign and discuss their implications on our understanding of the object.
  • The extragalactic background light (EBL) is the radiation accumulated through the history of the Universe in the wavelength range from the ultraviolet to the far infrared. Local foregrounds make the direct measurement of the diffuse EBL notoriously difficult, while robust lower limits have been obtained by adding up the contributions of all the discrete sources resolved in deep infrared and optical galaxy observations. Gamma-ray astronomy has emerged in the past few years as a powerful tool for the study of the EBL: very-high-energy (VHE) photons traversing cosmological distances can interact with EBL photons to produce e$^+$e$^-$ pairs, resulting in an energy-dependent depletion of the gamma-ray flux of distant sources that can be used to set constraints on the EBL density. The study of the EBL is one of the key scientific programs currently carried out by the MAGIC collaboration. We present here the results of the analysis of 32 VHE spectra of 12 blazars in the redshift range 0.03 - 0.94, obtained with over 300 hours of observations with the MAGIC telescopes between 2010 and 2016. A combined likelihood maximization approach is used to evaluate the density and spectrum of the EBL most consistent with the MAGIC observations. The results are compatible with state-of-the-art EBL models, and constrain the EBL density to be roughly within $\simeq 20\%$ of the nominal value in such models. The study reveals no anomalies in gamma-ray propagation in the large optical depth regime - contrary to some claims based on meta-analyses of published VHE spectra.
  • The wealth of recent data from imaging air Cherenkov telescopes (IACTs), ultra-high energy cosmic-ray experiments and neutrino telescopes have fuelled a renewed interest in hadronic emission models for gamma-loud blazars. We explore physically plausible solutions for a lepto-hadronic interpretation of the stationary emission from high-frequency peaked BL Lac objects (HBLs). The modelled spectral energy distributions are then searched for specific signatures at very high energies that could help to distinguish the hadronic origin of the emission from a standard leptonic scenario. By introducing a few basic constraints on parameters of the model, such as assuming the co-acceleration of electrons and protons, we significantly reduced the number of free parameters. We then systematically explored the parameter space of the size of the emission region and its magnetic field for two bright gamma-loud HBLs, PKS 2155-304 and Mrk 421. For all solutions close to equipartition between the energy densities of protons and of the magnetic field, and with acceptable jet power and light-crossing timescales, we inspected the spectral hardening in the multi-TeV domain from proton-photon induced cascades and muon-synchrotron emission inside the source. Very-high-energy spectra simulated with the available instrument functions from the future Cherenkov Telescope Array (CTA) were evaluated for detectable features as a function of exposure time, source redshift, and flux level. Over a large range of model parameters, the spectral hardening due to internal synchrotron-pair cascades, the "cascade bump", should be detectable for acceptable exposure times with the future CTA for a few nearby and bright HBLs.
  • Imaging atmospheric Cherenkov telescopes (IACTs) such as the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes endeavor to reach the lowest possible energy threshold. In doing so the trigger system is a key element. Reducing the trigger threshold is hampered by the rapid increase of accidental triggers generated by ambient light, the so-called Night Sky Background (NSB). In this paper we present a topological trigger, dubbed Topo-trigger, which rejects events on the basis of their relative orientation in the telescope cameras. We have simulated and tested the trigger selection algorithm in the MAGIC telescopes. The algorithm was tested using MonteCarlo simulations and shows a rejection of 85% of the accidental stereo triggers while preserving 99 % of the gamma rays. A full implementation of this trigger system would achieve an increase in collection area between 10 and 20% at the energy threshold. The analysis energy threshold of the instrument is expected to decrease by ~8 %. The selection algorithm was tested on real MAGIC data taken with the current trigger configuration and no $\gamma$-like events were found to be lost.
  • The Cherenkov Telescope Array (CTA) is the next generation ground-based very high energy gamma-ray observatory. The Large-Sized Telescope (LST) of CTA targets 20 GeV -- 1 TeV gamma rays and has 1855 photomultiplier tubes (PMTs) installed in the focal plane camera. With the 23 m mirror dish, the night sky background (NSB) rate amounts to several hundreds MHz per pixel. In order to record clean images of gamma-ray showers with minimal NSB contamination, a fast sampling of the signal waveform is required so that the signal integration time can be as short as the Cherenkov light flash duration (a few ns). We have developed a readout board which samples waveforms of seven PMTs per board at a GHz rate. Since a GHz FADC has a high power consumption, leading to large heat dissipation, we adopted the analog memory ASIC "DRS4". The sampler has 1024 capacitors per channel and can sample the waveform at a GHz rate. Four channels of a chip are cascaded to obtain deeper sampling depth with 4096 capacitors. After a trigger is generated in a mezzanine on the board, the waveform stored in the capacitor array is subsequently digitized with a low speed (33 MHz) ADC and transferred via the FPGA-based Gigabit Ethernet to a data acquisition system. Both a low power consumption (2.64 W per channel) and high speed sampling with a bandwidth of $>$300 MHz have been achieved. In addition, in order to increase the dynamic range of the readout we adopted a two gain system achieving from 0.2 up to 2000 photoelectrons in total. We finalized the board design for the first LST and proceeded to mass production. Performance of produced boards are being checked with a series of quality control (QC) tests. We report the readout board specifications and QC results.
  • The Cherenkov Telescope Array (CTA) is the the next generation facility of imaging atmospheric Cherenkov telescopes; two sites will cover both hemispheres. CTA will reach unprecedented sensitivity, energy and angular resolution in very-high-energy gamma-ray astronomy. Each CTA array will include four Large Size Telescopes (LSTs), designed to cover the low-energy range of the CTA sensitivity ($\sim$20 GeV to 200 GeV). In the baseline LST design, the focal-plane camera will be instrumented with 265 photodetector clusters; each will include seven photomultiplier tubes (PMTs), with an entrance window of 1.5 inches in diameter. The PMT design is based on mature and reliable technology. Recently, silicon photomultipliers (SiPMs) are emerging as a competitor. Currently, SiPMs have advantages (e.g. lower operating voltage and tolerance to high illumination levels) and disadvantages (e.g. higher capacitance and cross talk rates), but this technology is still young and rapidly evolving. SiPM technology has a strong potential to become superior to the PMT one in terms of photon detection efficiency and price per square mm of detector area. While the advantage of SiPMs has been proven for high-density, small size cameras, it is yet to be demonstrated for large area cameras such as the one of the LST. We are working to develop a SiPM-based module for the LST camera, in view of a possible camera upgrade. We will describe the solutions we are exploring in order to balance a competitive performance with a minimal impact on the overall LST camera design.
  • MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located on the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a major upgrade. The main subsystems upgraded were the MAGIC-I camera and its trigger system and 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 about 50 GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectra above 220 GeV 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 energies of a few hundred GeV is below 0.07degree, while the energy resolution is around 16%. We investigate 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 slope of the energy spectrum.
  • QSO B0218+357 is a blazar located at a cosmological redshift of z=0.944. It is gravitationally lensed by a spiral galaxy at a redshift of z=0.68. The blazar and its lens are well studied in the radio through X-ray bands, and several blazar outbursts were detected by Fermi-LAT at energies above 100 MeV. Strong gravitational lensing was invoked to explain the two components appar- ent in the radio and GeV light curves, separated by 10-12 days. In July 2014 another outburst was observed by Fermi-LAT, triggering follow-up observations with the MAGIC telescopes at energies above 100 GeV. The observations were scheduled at the expected time of arrival of the component delayed by the strong gravitational field of the lens, resulting in a firm detection of QSO B0218+357. Using the combined Fermi-LAT and MAGIC data sets, we report on variability of this unique blazar, the most distant among all currently known very high energy sources.
  • So far the current ground-based Imaging Atmospheric Cherenkov Telescopes (IACTs) have energy thresholds in the best case in the range of ~30 to 50 GeV (H.E.S.S. II and MAGIC telescopes). Lowest energy gamma-ray showers produce low light intensity images and cannot be efficiently separated from dominating images from hadronic background. A cost effective way of improving the telescope performance at lower energies is to use novel photosensors with superior photon detection efficiency (PDE). Currently the best commercially available superbialkali photomultipliers (PMTs) have a PDE of about 30-33%, whereas the silicon photomultipliers (SiPMs, also known as MPPC, GAPD) from some manufacturers show a photon detection efficiency of about 40-45%. Using these devices can lower the energy threshold of the instrument and may improve the background rejection due to intrinsic properties of SiPMs such as a superb single photoelectron resolution. Compared to PMTs, SiPMs are more compact, fast in response, operate at low voltage, and are insensitive to magnetic fields. SiPMs can be operated at high background illumination, which would allow to operate the IACT also during partial moonlight, dusk and dawn, hence increasing the instrument duty cycle. We are testing the SiPMs for Cherenkov telescopes such as MAGIC and CTA. Here we present an overview of our setup and first measurements, which we perform in two independent laboratories, in Munich, Germany and in Barcelona, Spain.
  • The MAGIC telescopes are two Imaging Atmospheric Cherenkov Telescopes (IACTs) located on the Canary island of La Palma. With 17m diameter mirror dishes and ultra-fast electronics, they provide an energy threshold as low as 50 GeV for observations at low zenith angles. The first MAGIC telescope was taken in operation in 2004 whereas the second one joined in 2009. In 2011 we started a major upgrade program to improve and to unify the stereoscopic system of the two similar but at that time different telescopes. Here we report on the upgrade of the readout electronics and digital trigger of the two telescopes, the upgrade of the camera of the MAGIC I telescope as well as the commissioning of the system after this major upgrade.
  • In this contribution we describe the design of the Data AcQuisition (DAQ) and online analysis software of the MAGIC telescopes after the 2012 upgrade. Although the final stereo trigger requires coincidence between the two telescopes, the actual data acquisition is performed independently, producing two separate data streams. Events are first readout and built from the front-end electronics and then stored in the DAQs' internal ring buffer for further processing: pre-calibration and signal extraction. The pixel signals, previously used only for data quality monitoring, are now also sent "on-the-fly" to the centralized online analysis program MOLA, which acts as a single client for the two DAQ data streams, and uses this information to provide preliminary high level analysis results. The integrated DAQ and online analysis programs allows an immediate feedback in case of a rapid gamma-ray flare of the pointed astrophysical source.
  • MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located at the Canary Island of La Palma, designed to observe gamma rays with energies above 50 GeV. Recently it has undergone an upgrade of the camera, digital trigger and readout systems. The upgrade has led to an improvement in the performance of the telescopes, especially in the lower energy range. We evaluate the performance of the upgraded MAGIC telescopes using Monte Carlo simulations and a large sample of Crab Nebula data. We study differential and integral sensitivity of the system, its angular resolution as well as its energy resolution.
  • Recently the readout of the MAGIC telescopes has been upgraded to a new system based on the Domino Ring Sampler version 4 chip. We present the analysis techniques and the signal extraction performance studies of this system. We study the behaviour of the baseline, the noise, the cross-talk, the linearity and the time resolution. We investigate also the optimal signal extraction. In addition we show some of the analysis techniques specific to the readout based on the Domino Ring Sampler version 2 chip, previously used in the MAGIC II telescope.
  • Very high energy (VHE, E >100 GeV) gamma-rays are absorbed via interaction with low-energy photons from the extragalactic background light (EBL) if the involved photon energies are above the threshold for electron-positron pair creation. The VHE gamma-ray absorption, which is energy dependent and increases strongly with redshift, distorts the VHE spectra observed from distant objects. The observed energy spectra of the AGNs carry, therefore, an imprint of the EBL. The detection of VHE gamma-ray spectra of distant sources (z = 0.11 - 0.54) by current generation Imaging Atmospheric Cherenkov Telescopes (IACTs) enabled to set strong upper limits on the EBL density, using certain basic assumptions about blazar physics. In this paper it is studied how the improved sensitivity of the Cherenkov Telescope Array (CTA) and its enlarged energy coverage will enlarge our knowledge about the EBL and its sources. CTA will deliver a large sample of AGN at different redshifts with detailed measured spectra. In addition, it will provide the exciting opportunity to use gamma ray bursts (GRBs) as probes for the EBL density at high redshifts.
  • The unfolding program TRUEE is a software package for the numerical solution of inverse problems. The algorithm was first applied in the FORTRAN77 program RUN. RUN is an event-based unfolding algorithm which makes use of the Tikhonov regularization. It has been tested and compared to different unfolding applications and stood out with notably stable results and reliable error estimation. TRUEE is a conversion of RUN to C++, which works within the powerful ROOT framework. The program has been extended for more user-friendliness and delivers unfolding results which are identical to RUN. Beside the simplicity of the installation of the software and the generation of graphics, there are new functions, which facilitate the choice of unfolding parameters and observables for the user. In this paper, we introduce the new unfolding program and present its performance by applying it to two exemplary data sets from astroparticle physics, taken with the MAGIC telescopes and the IceCube neutrino detector, respectively.
  • The large majority of extragalactic very high energy (VHE; E>100 GeV) sources belongs to the class of active galactic nuclei (AGN), in particular the BL Lac sub-class. AGNs are characterized by an extremely bright and compact emission region, powered by a super-massive black hole (SMBH) and an accretion disk, and relativistic outflows (jets) detected all across the electro-magnetic spectrum. In BL Lac sources the jet axis is oriented close to the line of sight, giving rise to a relativistic boosting of the emission. In radio galaxies, on the other hand, the jet makes a larger angle to the line of sight allowing to resolve the central core and the jet in great details. The giant radio galaxy M 87 with its proximity (1 6Mpc) and its very massive black hole ((3-6) x 10^9 M_solar) provides a unique laboratory to investigate VHE emission in such objects and thereby probe particle acceleration to relativistic energies near SMBH and in jets. M 87 has been established as a VHE emitter since 2005. The VHE emission displays strong variability on time-scales as short as a day. It has been subject of a large joint VHE and multi-wavelength (MWL) monitoring campaign in 2008, where a rise in the 43 GHz VLBA radio emission of the innermost region (core) was found to coincide with a flaring activity at VHE. This had been interpreted as a strong indication that the VHE emission is produced in the direct vicinity of the SMBH black hole. In 2010 again a flare at VHE was detected triggering further MWL observations with the VLBA, Chandra, and other instruments. At the same time M 87 was also observed with the Fermi-LAT telescope at GeV energies and the European VLBI Network (EVN). In this contribution preliminary results from the campaign will be presented.
  • Very high energy (VHE, energy $E \gtrsim 100$\,GeV) \gamma-rays from cosmological sources are attenuated due to the interaction with photons of the extragalactic background light (EBL) in the ultraviolet to infrared wavelength band. The EBL, thus, leaves an imprint on the observed energy spectra of these objects. In the last four years, the number of extragalactic VHE sources discovered with imaging atmospheric Cherenkov telescopes (IACTs), such as MAGIC, H.E.S.S., and VERITAS, has doubled. Furthermore, the measurements of the \emph{Fermi} satellite brought new insights into the intrinsic spectra of the sources at GeV energies. In this paper, upper limits on the EBL intensity are derived by considering the most extensive VHE source sample ever used in this context. This is accomplished by constructing a large number of generic EBL shapes and combining spectral informations from \emph{Fermi} and IACTs together with minimal assumptions about the source physics at high and very high \gamma-ray energies. The evolution of the EBL with redshift is accounted for and the possibility of the formation of an electromagnetic cascade and the implications on the upper limits are explored. The EBL density at $z=0$ is constrained over a broad wavelength range between 0.4 and 100\,\mu m. At optical wavelengths, the EBL density is constrained below 24\,nW\,m$^{-2}$\,sr$^{-1}$ and below 5\,nW\,m$^{-2}$\,sr$^{-1}$ between 8\,\mu m and 31\,\mu m.
  • MAGIC has been performing optically triggered Target of Opportunity observations of flaring blazars since the beginning of its scientific operations. The alerts of flaring blazars originate from Tuorla Blazar Monitoring Programme, which started the optical monitoring of candidtate TeV blazars in 2002 and has now collected up to eight years of data on more than 60 blazars. These ToO observations have resulted in the discovery of five new VHE gamma-ray emitting blazars (S5 0716+714, 1ES 1011+496, Mrk 180, ON 325 and B3 2247+381). In addition part of the discovery of BL Lac and the discovery of 3C 279 was made during a high optical state. In this contribution we present a detailed analysis of the optical light curves which are then compared to MAGIC observations of the same sources. We aim to answer the question: "Is there a connection between optical and VHE -ray high states in blazars or have we just been lucky?"
  • The Crab Pulsar Wind Nebula is the best studied source of $\gamma$-ray astrophysics. The contribution of the various soft radiation fields to the Inverse Compton component of its high energy emission, the strenght of the internal magnetic field and the maximum energies reached by primary electrons are however still matter of study. The MAGIC stereoscopic system recorded almost 50 hours of Crab Nebula data in the last two years, between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereo software provided an unprecedented differential energy spectrum spanning three decades in energy, from 50 GeV up to 45 TeV. At low energies, the MAGIC results, combined with the Fermi/LAT data, yield a precise measurement of the Inverse Compton peak. In addition, we present light curves of the Crab Nebula at different time scales, including a measurement simultaneous to one of the Crab Nebula flares recently detected by both Fermi/LAT and AGILE. Using the MAGIC spectrum together with multiwavelength data, we discuss the implications for the modeling of the Crab Nebula.
  • In the spring of 2008 MAGIC organised multi-wavelength (MWL) observations of the blazar 1ES 1011+496. 1ES 1011+496 is a high-frequency peaked BL Lac object discovered at VHE gamma-rays by MAGIC in spring 2007 during an optical outburst reported by the Tuorla Blazar Monitoring Programme. MAGIC re-observed the source during the 2008 MWL campaign which also included the Mets\"ahovi, KVA, Swift and AGILE telescopes. This was the first MWL campaign on this source that also included VHE coverage. MAGIC observed 1ES 1011+496 from March 4th to May 24th 2008 for a total of 27.9 hours, of which 20 h remained after quality cuts. The observations resulted in a detection of the source a ~7 sigma significance level with a mean flux and spectral index similar to those during the discovery. Here we will present the results of the MAGIC observations of the source in combination with contemporaneous observations at other wavelengths (radio, optical, X-rays, high energy gamma-rays) and discuss their implications on the modelling of the spectral energy distribution.
  • The MAGIC telescope observed the region around the distant blazar 3C 66A for 54.2 hr in 2007 August-December. The observations resulted in the discovery of a gamma-ray source centered at celestial coordinates R.A. = 2h23m12s and decl.= 43deg 0.7arcmin (MAGIC J0223+430), coinciding with the nearby radio galaxy 3C 66B. A possible association of the excess with the blazar 3C 66A is discussed. The energy spectrum of MAGIC J0223+430 follows a power law with a normalization of (1.7+-0.3_stat+-0.6_syst) x 10^-11 TeV^-1 cm-2 s-1 at 300 GeV and a photon index of 3.10+-0.31_stat+-0:2syst.
  • MAGIC has been performing optically triggered Target of Opportunity observations of flaring blazars since the beginning of its scientific operations. The alerts of flaring blazars originate from Tuorla blazar monitoring program, which started the optical monitoring of the candidate TeV blazars in 2002 and has now collected up to six years of data on 30 blazars. These Target of Opportunity observations have resulted in the discovery of three new VHE gamma-ray emitting blazars (Mrk180, 1ES1011+496 and S50716+714) and in addition the discoveries of BLLac and 3C279 were made during a high optical state. In this talk we present a detailed analysis of the optical light curves which are then compared to MAGIC observations of the same sources. We aim to answer the question "Is there a connection between optical and VHE gamma-ray high states in blazars or have we just been lucky?"
  • Very high energy gamma ray emission from the active galactic nucleus PG 1553+113 was observed during 2005 and 2006 by the MAGIC collaboration, for a total observation time of 18.8 hours. Here we present the results of follow up observations: more than 20 hours of good quality data collected by the MAGIC Telescope during the 2007 and 2008 campaigns. The obtained spectra are compared and combined with previous measurements, and corrected for absorption adopting different EBL models. Upper limits on the unknown source redshift are derived by assuming the absence of a break in the intrinsic spectrum, or alternatively by constraining the hardness of the intrinsic source spectrum.
  • We present the results of X-ray observations of the well-studied TeV blazar Mrk 421 with the Suzaku satellite in 2006 April 28. During the observation, Mrk 421 was undergoing a large flare and the X-ray flux was variable, decreasing by ~ 50 %, from 7.8x10^{-10} to 3.7x10^{-10} erg/s/cm^2 in about 6 hours, followed by an increase by ~ 35 %. Thanks to the broad bandpass coupled with high-sensitivity of Suzaku, we measured the evolution of the spectrum over the 0.4--60 keV band in data segments as short as ~1 ksec. The data show deviations from a simple power law model, but also a clear spectral variability. The time-resolved spectra are fitted by a synchrotron model, where the observed spectrum is due to a exponentially cutoff power law distribution of electrons radiating in uniform magnetic field; this model is preferred over a broken power law. As another scenario, we separate the spectrum into "steady" and "variable" components by subtracting the spectrum in the lowest-flux period from those of other data segments. In this context, the difference ("variable") spectra are all well described by a broken power law model with photon index Gamma ~ 1.6, breaking at energy epsilon_{brk} ~ 3 keV to another photon index Gamma ~ 2.1 above the break energy, differing from each other only by normalization, while the spectrum of the "steady" component is best described by the synchrotron model. We suggest the rapidly variable component is due to relatively localized shock (Fermi I) acceleration, while the slowly variable ("steady") component is due to the superposition of shocks located at larger distance along the jet, or due to other acceleration process, such as the stochastic acceleration on magnetic turbulence (Fermi II) in the more extended region.
  • The quasar 3C 279 is one of the best-studied flat spectrum radio quasars. It is located at a comparatively large redshift of z=0.536: E>100 GeV observations of such distant sources were until recently impossible both due to the expected steep energy spectrum and the expected attenuation of the gamma-rays by the extragalactic background light. Here we present results on the observation of 3C 279 with the MAGIC telescope in early 2006. We report the detection of a significant very high energy gamma-ray signal in the MAGIC energy range on the observation night of 2006 February 23.