For the first time an excess of photons above an energy threshold of 730 GeV from the giant radio galaxy M87 has been measured at a significance level above 4 sigma. The data have been taken during the years 1998 and 1999 with the HEGRA stereoscopic system of 5 imaging atmospheric Cherenkov telescopes. The excess of 107.4 +- 26.8 events above 730 GeV corresponds to an integral flux of 3.3% of the Crab flux or N_gamma(E > 730 GeV) = (0.96 +- 0.23) * 10^(-12) phot.cm^(-2)s^(-1). M87 is located at the center of the Virgo cluster of galaxies at a relatively small redshift of z = 0.00436 and is a promising candidate among the class of giant radio galaxies for the emission of TeV gamma-radiation. The detection of TeV gamma-rays from M87 -- if confirmed -- would establish a new class of extragalactic source in this energy regime since all other AGN detected to date at TeV energies are BL Lac type objects.

The unidentified TeV source in Cygnus is now confirmed by follow-up observations from 2002 with the HEGRA stereoscopic system of Cherenkov Telescopes. Using all data (1999 to 2002) we confirm this new source as steady in flux over the four years of data taking, extended with radius 6.2 arcmin (+-1.2 arcmin (stat) +-0.9 arcmin (sys)) and exhibiting a hard spectrum with photon index -1.9. It is located in the direction of the dense OB stellar association, Cygnus OB2. Its integral flux above energies E>1 TeV amounts to \~5% of the Crab assuming a Gaussian profile for the intrinsic source morphology. There is no obvious counterpart at radio, optical nor X-ray energies, leaving TeVJ2032+4130 presently unidentified. Observational parameters of this source are updated here and some astrophysical discussion is provided. Also included are upper limits for a number of other interesting sources in the FoV, including the famous microquasar Cygnus X-3.

TeV gamma-rays from the BL Lac object 1ES1959+650 have been measured during the years 2000 and 2001 with a significance of 5.2 sigma at a value of 5.3% of the Crab flux and in May 2002 during strong outbursts with > 23 sigma at a flux level of up to 2.2 Crab, making 1ES1959+650 the TeV Blazar with the third best event statistics. The deep observation of 197.4 h has been performed with the HEGRA stereoscopic system of 5 imaging atmospheric Cherenkov telescopes (IACT system). 1ES1959+650 is located at a redshift of z = 0.047, providing an intermediate distance between the nearby Blazars Mkn 421 and Mkn 501, and the much more distant object H1426+428. This makes 1ES1959+650 an important member of the class of TeV Blazars in view of the absorption of TeV photons by the diffuse extragalactic background radiation (DEBRA). The differential energy spectrum of 1ES1959+650 during the flares can be fitted by a power law with a spectral index of 2.83 +- 0.14_stat +- 0.08_sys or by a power law with an exponential cut-off at (4.2^(+0.8)_(-0.6)_stat +- 0.9_sys) TeV and a spectral index of 1.83 +- 0.15_stat +- 0.08_sys. The low state differential energy spectrum obtained with lower statistics can be described by a pure power law with a spectral index of 3.18 +- 0.17_stat +- 0.08_sys.

This preprint encompasses thirteen individual contributions of the HEGRA collaboration to the 28th International Cosmic Ray conference in Tsukuba, Japan. The papers report on results on M87, M31, H1426+428, Mkn421, 1ES1959+650, Scan of the TeV sky, Performance of the system, GeV sources, 54 AGN, and SN1006

A search for potential point sources of very high energy gamma rays has been carried out on the data taken simultaneously by the HEGRA AIROBICC and Scintillator arrays from August 1994 to March 2000. The list of sought sources includes supernova remnants, pulsars, AGNs and binary systems. The energy threshold is around 15 TeV. For the Crab Nebula, a modest excess of 2.5 standard deviations above the cosmic ray background has been observed. Flux upper limits (at 90% c.l.) of around 1.3 times the flux of the Crab Nebula are obtained, on average, for the candidate sources. A different search procedure has been used for an all-sky search which yields absolute flux upper limits between 4 and 9 crabs depending on declination, in the band from 0 to 60 degrees.

The nearby BL Lacertae (BL Lac) object Markarian~421 (Mkn~421) at a red shift z=0.031 was observed to undergo strong TeV Gamma-ray outbursts in the observational periods from December 1999 until May 2001. The time averaged energy spectrum is well described by a power law with an exponential cut-off energy at 3.6(+0.4-0.3)_{stat}(+0.9-0.8)_{sys} TeV. The energy spectra derived for different average flux levels ranging from 0.5 to 10 * 10^{-11}ph/cm**2/s follow a clear correlation of photon index and flux level. Generally, the energy spectra are harder for high flux levels. From January to April 2001 Mkn~421 showed rapid variability (doubling time as short as 20 minutes), accompanied with a spectral hardening with increasing flux level within individual nights. For two successive nights (MJD 51989-51991, March 21-23,2001), this correlation of spectral hardness and change in flux has been observed within a few hours. The cut-off energy for the Mkn~421 TeV spectrum remains within the errors constant for the different flux levels and differs by 2.6+/-0.6_stat+/-0.6_sys TeV from the value determined for Mkn~501. This indicates that the observed exponential cut-off in the energy spectrum of Mkn~421 is not solely caused by absorption of multi-TeV photons by pair-production processes with photons of the extragalactic near/mid infrared background radiation.

The detection of TeV gamma-rays from the blazar H 1426+428 at an integral flux level of (4 +- 2(stat) +- 1(syst)) 10^(-12) erg cm^(-2) s^(-1) above 1 TeV with the HEGRA imaging atmospheric Cherenkov telescope system is reported. H 1426+428 is located at a redshift of z=0.129, which makes it the most distant source detected in TeV gamma-rays so far. The TeV radiation is expected to be strongly absorbed by the diffuse extragalactic background radiation (DEBRA). The observed energy spectrum of TeV photons is in good agreement with an intrinsic power law spectrum of the source ~E^(-1.9) corrected for DEBRA absorption. Statistical errors as well as uncertainties about the intrinsic source spectrum, however, do not permit strong statements about the density of the DEBRA infrared photon field.

Using data from the HEGRA air shower array, taken in the period from April 1998 to March 2000, upper limits on the ratio I_gamma /I_CR of the diffuse photon flux I_gamma to the hadronic cosmic ray flux I_CR are determined for the energy region 20 TeV to 100 TeV. The analysis uses a gamma-hadron discrimination which is based on differences in the development of photon- and hadron-induced air showers after the shower maximum. A method which is sensitive only to the non-isotropic component of the diffuse photon flux yields an upper limit of I_gamma /I_CR (at 54 TeV) < 2.0*10^-3 (at the 90% confidence level) for a sky region near the inner galaxy (20 degrees < galactic longitude <60 degrees and |galactic latitude|<5 degrees). A method which is sensitive to both the isotropic and the non-isotropic component yields global upper limits of I_gamma /I_CR (at 31 TeV) <1.2*10^-2 and I_gamma /I_CR (at 53 TeV) < 1.4*10^-2 (at the 90% confidence level).

232 hours of data were accumulated from 1997 to 1999, using the HEGRA Stereoscopic Cherenkov Telescope System to observe the supernova remnant Cassiopeia A. TeV gamma ray emission was detected at the 5 sigma level, and a flux of (5.8 +- 1.2(stat) +- 1.2(syst)) 10^(-9) ph m^(-2) s^(-1) above 1 TeV was derived. The spectral distribution is consistent with a power law with a differential spectral index of -2.5 +- 0.4(stat) +- 0.1(syst) between 1 and 10 TeV. As this is the first report of the detection of a TeV gamma ray source on the "centi-Crab" scale, we present the analysis in some detail. Implications for the acceleration of cosmic rays depend on the details of the source modeling. We discuss some important aspects in this paper.

For distant extragalactic sources of gamma-rays in the PeV energy range, interactions of the gamma rays with intergalactic diffuse radiation fields will initiate a pair cascade. Depending on the magnetic fields in the vicinity of the source, the cascade can either result in an isotropic halo around an initially beamed source, or remain more or less collimated. Data recorded by the HEGRA system of imaging atmospheric Cherenkov telescopes are used to derive limits on the halo flux from the AGN Mrk 501. This is achieved by comparing the angular distribution of TeV gamma-rays during the 1997 burst phase -- where direct photons should dominate -- with the distribution during the 1998/99 quiescent state, where a steady-state halo contribution should be most pronounced. The results depend on the assumptions concerning the angular distribution of the halo; limits on the halo flux within 0.5 to 1 degr. from the source range between 0.1% and 1% of the peak burst flux.

Data taken with the HEGRA system of imaging atmospheric Cherenkov telescopes during the 1997 flares of Mkn 501 are reanalyzed using an algorithm providing improved energy resolution. A resolution of 10% to 12% is obtained by accounting for the variation of the Cherenkov light yield with the height of the shower maximum in the atmosphere. The improved energy resolution is particularly relevant for the study of the high-energy cutoff in the spectrum, which might be caused by interactions with the intergalactic infrared background radiation. The reanalysis presented here confirms the results obtained in the previous analysis, but hints a steeper slope of the spectrum in the region around 20 TeV.

The HEGRA system of imaging atmospheric Cherenkov telescopes provides for specially selected classes of events an angular resolution of better than 3'. By comparing the measured angular distribution of TeV gamma rays from the Crab Nebula with the distribution expected on the basis of Monte Carlo simulations, and with measurements of gamma rays from the point source Mrk 501, we conclude that the rms size of the VHE gamma-ray emission region in the Crab Nebula is less than 1.5'.

For stereoscopic systems of imaging atmospheric Cherenkov telescopes (IACTs), a key parameter to optimize the sensitivity for VHE gamma-ray point sources is the intertelescope spacing. Using pairs of telescopes of the HEGRA IACT system, the sensitivity of two-telescope stereo IACT systems is studied as a function of the telescope spacing, ranging from 70 m to 140 m. Data taken during the 1997 outburst of Mrk 501 are used to evaluate both the detection rates before cuts, and the sensitivity for weak signals after cuts to optimize the significance of signals. Detection rates decrease by about 1/3 between the minimum and maximum spacing. The significance of signals is essentially independent of distance, in the range investigated.

The time profiles of Cherenkov images of cosmic-ray showers and of gamma-ray showers are investigated, using data gathered with the HEGRA system of imaging atmospheric Cherenkov telescopes during the 1997 outbursts of Mrk 501. Photon arrival times are shown to vary across the shower images. The dominant feature is a time gradient along the major axis of the images. The gradient varies with the distance between the telescope and the shower core, and is maximal for large distances. The time profiles of cosmic-ray showers and of gamma-ray showers differ in a characteristic fashion. The main features of the time profiles can be understood in terms of simple geometrical models. Use of the timing information towards improved shower reconstruction and cosmic-ray suppression is discussed.

A search for gamma-ray bursts (GRBs) above 20 TeV within the field of view (1 sr) of the HEGRA AIROBICC Cherenkov array (29N, 18W, 2200 m a.s.l.) has been performed using data taken between March 1992 and March 1993. The search is based on an all-sky survey using four time scales, 10 seconds, 1 minute, 4 minutes and 1 hour. No evidence for TeV-emission has been found for the data sample. Flux upper limits are given. A special analysis has been performed for GRBs detected by BATSE and WATCH. Two partially and two fully contained GRBs in our field of view (FOV) were studied. For GRB 920925c which was fully contained in our FOV, the most significant excess has a probability of 7.7 10**-8 (corresponding to 5.4 sigmas) of being caused by a background fluctuation. Correcting this probability with the appropriate trial factor, yields a 99.7% confidence level for this excess to be related to the GRB (corresponding to 2.7 sigmas). This result is discussed within the framework of the WATCH detection.

Using air showers induced by TeV gamma-rays from Mrk 501, the radial distribution of Cherenkov light is investigated. The shower geometry is reconstructed from the stereoscopic shower images obtained with the telescopes of the HEGRA IACT system. We observe a change in the shape of the light pool with shower energy, and with zenith angle. Data are well reproduced by Monte-Carlo air shower simulations.

First results concerning the performance characteristics of the HEGRA IACT array are given, based on stereoscopic observations of the Crab Nebula with four telescopes. The system provides a gamma-ray energy threshold around 0.5 TeV. The Crab signal demonstrates an angular resolution of about 0.1 deg. Shape cuts allow to suppress cosmic ray background by almost a factor 100, while maintaining 40% efficiency for gamma-rays. The Crab signal is essentially background free. For longer observation times of order 100 h, the system in its present form provides sensitivity to point sources at a level of 3% of the Crab flux. Performance is expected to improve further with the inclusion of the fifth telescope and the implementation of advanced algorithms for shower reconstruction.