HAGAR is a system of seven Non-imaging Atmospheric Cherenkov Telescopes located at Hanle in the Ladakh region of the Indian Himalayas at an altitude of 4270 meters {\it amsl}. Since 2008, we have observed the Crab Nebula to assess the performance of the HAGAR telescopes. We describe the analysis technique for the estimation of $\gamma$-ray signal amidst cosmic ray background. The consolidated results spanning nine years of the Crab nebula observations show long term performance of the HAGAR telescopes. Based on about 219 hours of data, we report the detection of $\gamma$-rays from the Crab Nebula at a significance level of about 20$\sigma$, corresponding to a time averaged flux of (1.64$\pm$0.09) $\times10^{-10}$ photons cm$^{-2}$ sec$^{-1}$ above 230 GeV. Also, we perform a detailed study of possible systematic effects in our analysis method on data taken with the HAGAR telescopes.

Context:The HAGAR Telescope Array at Hanle, Ladakh has been regularly monitoring the nearby blazar Mkn 421 for the past 7yrs. Aims: Blazars show flux variability in all timescales across the electromagnetic spectrum. While there is abundant literature characterizing the short term flares from different blazars, comparatively little work has been done to study the long term variability. We aim to study the long term temporal and spectral variability in the radiation from Mkn 421 during 2009-2015. Methods: We quantify the variability and lognormality from the radio to the VHE bands, and compute the correlations between the various wavebands using the z-transformed discrete correlation function. We construct the Spectral Energy Distribution (SED) contemporaneous with HAGAR observation seasons and fit it with a one zone synchrotron self Compton model to study the spectral variability. Results: The flux is found to be highly variable across all time scales. The variability is energy dependant, and is maximum in the X-ray and Very High Energy (VHE) bands. A strong correlation is found between the Fermi-LAT (gamma) and radio bands, and between Fermi-LAT and optical, but none between Fermi-LAT and X-ray. Lognormality in the flux distribution is clearly detected. This is the third blazar, following BL~Lac and PKS~2155$+$304 to show this behavior. The SED can be well fit by a one zone SSC model, and variations in the flux states can be attributed mainly due to changes in the particle distribution. A strong correlation is seen between the break energy $\gamma_b$ of the particle spectrum and the total bolometric luminosity.

Since the detection of very high energy (VHE) $\gamma$-rays from Mrk 501, its broad band emission of radiation was mostly and quite effectively modeled using one zone emission scenario. However, broadband spectral and flux variability studies enabled by the multiwavelength campaigns carried out during the recent years have revealed rather complex behavior of Mrk 501. The observed emission from Mrk 501 could be due to a complex superposition of multiple emission zones. Moreover new evidences of detection of very hard intrinsic $\gamma$-ray spectra obtained from {\it Fermi}--LAT observations have challenged the theories about origin of VHE $\gamma$-rays. Our studies based on {\it Fermi}--LAT data indicate the existence of two separate components in the spectrum, one for low energy $\gamma$-rays and the other for high energy $\gamma$-rays. Using multiwaveband data from several ground and space based instruments, in addition to HAGAR data, the spectral energy distribution of Mrk~501 is obtained for various flux states observed during 2011. In the present work, this observed broadband spectral energy distribution is reproduced with a leptonic, multi-zone Synchrotron Self-Compton model.

Context: The nearby (z=0.031) TeV blazar Mrk421 was reported to be in a high state of flux activity since November, 2009. Aims: To investigate possible changes in the physical parameters of Mrk421 during its high state of activity using multiwavelength data. Methods: We have observed this source in bright state using High Altitude GAmma Ray (HAGAR) telescope array at energies above 250 GeV during February 13 - 19, 2010. Optical, X-ray and gamma-ray archival data are also used to obtain the SEDs and light curves. Results: Mrk421 was found to undergo one of its brightest flaring episodes on February 17, 2010 by various observations in X-rays and gamma-rays. HAGAR observations during February 13 - 19, 2010 at the energies above 250 GeV show an enhancement in the flux level, with a maximum flux of ~ 7 Crab units being detected on February 17, 2010. We present the spectral energy distributions during this flaring episode and investigate the correlation of the variability in X-ray and gamma-ray bands. Conclusions: Our multiwavelength study suggests that the flare detected during February 16 and 17, 2010 could arise due to a passing shock in the jet.

The radiation mechanism of very high energy $\gamma$-ray emission from blazars and crucial parameters like magnetic field, and size of the emitting region are not well understood yet. To understand the above mentioned properties of blazars, we observed five nearby TeV $\gamma$-ray emitting blazars (Mrk421, Mrk501, 1ES2344+514, 1ES1218+304 and 3C454.3) and one radio galaxy (M87) using the High Altitude GAmma Ray (HAGAR) telescope. HAGAR is an array of seven telescopes located at Hanle, India to detect Cherenkov light caused by extensive air showers initiated by $\gamma$-rays. Mrk421 was observed to undergo one of its brightest flaring episodes on 2010 February 17, and detected by various experiments in X-rays and $\gamma$- rays. HAGAR observations of this source during 2010 February 13 - 19, in the energies above 250 GeV show an enhancement in the flux level, with a flux of 6-7 Crab units being detected on 2010 February 17. We present the spectral energy distribution of the source during this flaring episode. In addition to this, the analysis procedure to extract $\gamma$-ray signal from HAGAR data is discussed and preliminary results on all the AGNs are presented.

We report our observations on 4 AGNs viz, Mkn421, Mkn501, 1ES1426+428 and ON231 belonging to a sub-class called blazars. The observations were carried out using the Pachmarhi Array of Cherenkov Telescopes and span about 6 years period from 2000 to 2005. We discuss our methods of analysis adopted to extract the gamma ray signal from cosmic ray background. We present our results on the emission of TeV gamma rays from these objects. Also, we report on the status of the new high altitude experiment, a 7 telescope array at Hanle in the Himalayas at an altitude of about 4200 m above mean sea level for detecting celestial gamma rays.

The Pachmarhi Array of Cerenkov Telescopes consists of a distributed array of 25 telescopes that are used to sample the atmospheric Cerenkov Photon showers. Each telescope consists of 7 parabolic mirrors each viewed by a single photo-multiplier tube. Reconstruction of photon showers are carried out using fast timing information on the arrival of pulses at each PMT. The shower front is fitted to a plane and the direction of arrival of primary particle initiating the shower is obtained. The error in the determination of the arrival direction of the primary has been estimated using the {\it split} array method. It is found to be $\sim 2.4^\prime ~$ for primaries of energy $ > 3 \~TeV$. The dependence of the angular resolution on the separation between the telescopes and the number of detectors are also obtained from the data.

Pachmarhi Array of \v Cerenkov Telescopes (PACT) consists of a 5$\times$5 array of \v Cerenkov telescopes deployed over an area of 100 $m$ $\times$ 80 $m$, in the form of a rectangular matrix. The experiment is based on atmospheric \v Cerenkov technique using wavefront sampling technique. Each telescope consists of 7 parabolic mirrors mounted para-axially on an equatorial mount. At the focus of each mirror a fast phototube is mounted. In this experiment a large number of parameters have to be measured and recorded from each of the 175 phototubes in the shortest possible time. Further, the counting rates from each phototube as well as the analog sum of the 7 phototubes from each telescope (royal sum) need to be monitored at regular intervals during the run. In view of the complexity of the system, the entire array is divided into four smaller sectors each of which is handled by an independent field signal processing centre (FSPC) housed in a control room that collects, processes and records information from nearby six telescopes that belong to that sector. The distributed data acquisition system (DDAS) developed for the purpose consists of stand-alone sector data acquisition system (SDAS) in each of the four FSPC's and a master data acquisition system (MDAS). MDAS running in the master signal processing centre (MSPC) records data from each of the 25 telescopes. The data acquisition and monitoring PCs (SDAS and MDAS) are networked through LAN. The entire real time software for DDAS is developed in C under $linux$ environment. Most of the hardware in DDAS are designed and fabricated in-house. The design features and the performance of the entire system along with some other auxiliary systems to facilitate the entire observations will be presented.

Pachmarhi Array of \v{C}erenkov Telescopes (PACT) has been designed to search for celestial TeV $\gamma-$rays using the wavefront sampling technique. PACT, located at Pachmarhi, (latitude 22$^\circ$ 28$^\prime$ N, longitude 76$^\circ$ 26$^\prime$ E, altitude 1075 m) consists of 25 telescopes deployed over an area of 80 m $\times$ 100m. Each telescopes consists of 7 parabolic reflectors, each viewed by a fast phototube behind a 3$^\circ$ mask at the focus. The density and the arrival time of the photons at the PMT are recorded for each shower. The energy threshold and collection area of the array are estimated, from Monte Carlo simulations, to be $\sim$ 900 GeV and 10$^5$ m$^2$ respectively. The accuracy in determination of arrival angle of a shower is estimated to be about 0.1$^\circ$ in the near vertical direction. About 99% of the off-axis hadronic events could be rejected from directional information alone. Further, at least 75% of the on-axis hadronic events could be rejected using species sensitive parameters derived from timing and density measurements. These cuts on data to reject background would retain $\sim$ 44% of the $\gamma-$ray signal. The sensitivity of the array for a 5$\sigma$ detection of $\gamma-$ray signal at a threshold energy of 1 TeV has been estimated to be $\sim$ 4.1 $\times$ 10$^{-12}$ photons cm$^{-2}$ s$^{-1}$ for an on source exposure of 50 hours. The PACT set-up has been fully commissioned and is collecting data. The details of the system parameters and sensitivity will be presented.

Recently an Indian Astronomical Observatory has been set up at Hanle (32$^\circ$ 46$^\prime$ 46$^{\prime\prime}$ N, 78$^\circ$ 57$^\prime$ 51$^{\prime\prime}$ E, 4515m amsl) situated in the high altitude cold desert in the Himalayas. The Observatory has 2-m aperture optical-infrared telescope, recently built by the Indian Institute of Astrophysics. We have carried out systematic simulations for this observation level to study the nature of \v{C}erenkov light pool generated by gamma ray and proton primaries incident vertically at the top of the atmosphere. The differences in the shape of the lateral distributions of \v{C}erenkov light with respect to that at lower altitudes is striking. This arises primarily due to the proximity of the shower maximum to the observation site. The limited lateral spread of the \v{C}erenkov light pool and near 90% atmospheric transmission at this high altitude location makes it an ideal site for a gamma ray observatory. This results in a decrease in the gamma ray energy threshold by a factor of 2.9 compared to that at sea-level. Several parameters based on density and timing information of \v{C}erenkov photons, including local and medium range photon density fluctuations as well as photon arrival time jitter could be efficiently used to discriminate gamma rays from more abundant cosmic rays at tens of GeV energies.

Monte Carlo calculations have been done to study the possibilities of using different spatial distributions, local fluctuations and the number of Cerenkov photons in gamma ray showers to increase the sensitivity of the PACT experiment. After matching the trigger rates in data with the predictions of simulations, the energy threshold of the array for gamma rays is found to be about 900 GeV. It is shown that a lateral distribution parameter $\beta$ can increase the signal to noise ratio in the experiment. A typical season of 15 nights should be able to conclusively establish any source with Crab level fluxes. Two other possible parameters for increase of sensitivity are also discussed. Preliminary results on the application of the $\beta$ parameter to the data from PACT experiment are discussed.

The main inference from the experiments of the '80s that the time-averaged energy spectra of pulsars had to steepen in the GeV-TeV energy region has been reinforced in the '90s from experiments with higher sensitivities. However, results from several experiments from the past and the more sensitive experiments at present can be reconciled by invoking a possibly different component arising in the TeV region. The results of the preliminary analysis of the data being collected with the PACT array will be presented.

Pachmarhi Array of \v{C}erenkov Telescopes(PACT), consisting of a distributed array of 25 telescopes is used to sample the atmospheric \v{C}erenkov Photon showers. The shower front is fitted to a plane and the direction of arrival of primary particle is obtained. The accuracy in the estimation of the arrival direction of showers has been estimated to be $\sim 0^{\circ}.1~$ using `split' array method. The angular resolution is expected to be even better when a spherical front is used for direction reconstruction or correction for the curvature of the front is applied. This is the best angular resolution among all the currently operating atmospheric \v{C}erenkov telescopes in the world.

The Crab nebula has proved to be the nearest to a standard candle in VHE $\gamma-$ ray astronomy. Results on the gamma ray emission from the nebula at various energies have come in the last decade mostly from imaging telescopes. The aim of the new Pachmarhi Atmospheric Cerenkov Telescope (PACT) array has been to use the temporal and spatial distribution of Cerenkov photons in distinguishing between proton and gamma ray showers. The array, with timing information from 175 mirrors, is ideally suited for precise estimation of the arrival direction. Preliminary results from the recent data taken on Crab has shown a 12$\sigma$ signal with the flux in good agreement with those of the other experiments.

Pachmarhi Array of \v{C}erenkov Telescopes (PACT), based on wavefront sampling technique, has been used for detecting TeV gamma rays from galactic and extra-galactic $\gamma $-ray sources. The Blazar, Mkn 421 was one such extra-galactic source observed during the winter nights of 2000 and 2001. We have carried out a preliminary analysis of the data taken during the nights of January, 2000 and 2001. Results show a significant gamma ray signal from this source during both these periods above a threshold energy of 900 GeV. The source was contemporaneously observed by CAT imaging telescope during the first episode of January 2000 while HEGRA CT1 was observing the source during the second episode. Both these observations have detected variable $\gamma $-ray emission this source and they reported that it was flaring during both these periods. The light curve in the TeV gamma ray range derived from the first PACT observations during both these episodes is in agreement with that reported by other experiments. The analysis procedure and the preliminary results will be presented and discussed.

We study the propagation of Cerenkov photons generated by Very High Energy gamma-rays and hadrons in the atmosphere. The photon production height distributions are estimated from semi-empirical methods and compared with those derived by standard simulation techniques. Incident spectra at various observation altitudes are then derived after applying wavelength dependent corrections due to photon attenuation in the atmosphere during the propagation of photons from the height of production to the height of observation. These are generated both for gamma- and hadron primaries of various energies. The derived production height distributions agree very well with those generated by the simulation package `CORSIKA' at all energies and for both gamma-$ray and proton primaries. The incident photon spectra are found to be both altitude dependent and primary energy dependent. The peak of the incident spectrum shifts towards the shorter wavelength with increasing altitude of observation for a given primary. Also the peak of the photon spectrum shifts towards the shorter wavelength with increasing energy of the primary at given altitude. The fraction of the UV component in the incident Cerenkov spectrum is estimated both for gamma-ray and hadronic primaries at various observation altitudes and energies. Hadron generated Cerenkov spectra are marginally richer in UV light and the difference increases slightly at higher altitudes. The fraction of the UV to the visible light in the Cerenkov spectrum could be a useful parameter to separate gamma-rays from cosmic ray background only if one can measure this fraction very accurately.

Very High Energy (VHE) gamma -ray Astronomy observations are planned to be carried out at Pachmarhi in the central Indian state of Madhya Pradesh using the well known atmospheric Cerenkov technique. An array of 25 Cerenkov telescopes is currently under construction. Using this array it is proposed to sample the Cerenkov light pool at various distances from the shower core in order to estimate the lateral distribution parameters of the shower. In this paper we discuss the scientific motivation of this concept of enriching the gamma -ray signal as compared to the standard imaging technique. The current status of the detector development and the expected results will be presented.