Cosmic rays have valuable information about universe surroundings us. Finding energy, mass and arrival direction of primary cosmic ray particle are the most important aspects of extensive air shower studies. In order to determine these parameters, arrival direction and core position of extensive air showers should be determined at first. In this article, a new method has been introduced that utilizes arrival time information of secondary particles of extensive air showers for finding their core location and correcting plane wave front approximation so that calculate the arrival direction of extensive air showers. This method does not need number sensitive detectors-detectors which are sensitive to the number of crossing particles- and consequently there is no need for lateral distribution models. This model has been developed by analysis of simulated events generated by CORSIKA package.

The most important goal of studying an extensive air shower is to find the energy, mass and arrival direction of its primary cosmic ray. In order to find these parameters, the core position and arrival direction of the shower should be determined. In this paper, a new method for finding core location has been introduced that utilizes trigger time information of particle detectors. We have also developed a simple technique to reconstruct the arrival direction. Our method is not based upon density-sensitive detectors which are sensitive to the number of crossing particles and is also independent of lateral distribution models. This model has been developed and examined by the analysis of simulated shower events generated by the CORSIKA package.

A new technique has been devised for the analysis of extensive air shower data in observing the effect of the moon on this data. In this technique the number of EAS events with arrival directions falling in error circles centered about the moving moon is compared to the mean number of events falling in error circles with centers randomly chosen in the sky. For any assumed angular radius of the error circle the deficit in EAS event count in the direction of moon i.e., Nsky-Nmoon which is a moon-related effect is interpreted as the shadow of the moon. A simple theoretical model has been developed to relate Nsky to the angular radius of the error circle and has been applied to the counts from the moon's direction in order to extract the physical parameters of the shadow of the moon. The technique and the theoretical model has been used on 1.7 *10^5 EAS events recorded at Alborz observatory.

The study of the anisotropy of the arrival directions is an essential tool to investigate the origin and propagation of cosmic rays primaries. A simple way of recording many cosmic rays is to record coincidences between a number of detectors. We have monitored multi-TeV cosmic rays by a small array of water cherenkov detectors in Tehran(35 43 N, 51 20 E, 1200m a.s.l). More than 1.1*10^6 extensive air shower events were recorded. In addition to the Compton- Getting effect due to the motion of the earth in the Galaxy, an anisotropy has been observed which is due to a unidirectional anisotropy of cosmic ray flow along the Galactic arms.

Ultra-High-Energy (UHE) ($E>100 $TeV) Extensive Air Showers (EASs) have been monitored for a period of five years (1997-2003), using a small array of scintillation detectors in Tehran, Iran. The data have been analyzed to take in to account of the dependence of source counts on zenith angle. Because of varying thickness of the overlaying atmosphere, the shower count rate is extremely dependent on zenith angle. During a calendar year different sources come in the field of view of the array at varying zenith angles and have different effective observation time equivalent to zenith in a day. High energy gamma-ray sources from the EGRET third catalogue where observed and the data were analyzed using an excess method. Upper limits were obtained for 10 EGRET sources. Then we investigated the EAS event rates for these 10 sources and obtained a flux for each of them using parameters of our experiment results and simulations. Finally we investigated the gamma-ray spectrum in the UHE range using these fluxes with reported fluxes of the EGRET sources.

Ultra-high-energy (E$>$100 Tev) Extensive Air Showers (EASs) have been monitored for a period of five years (1997 - 2003), using a small array of scintillator detectors in Tehran, Iran. The data have been analyzed to take in to account of the dependence of source counts with respect to the zenith angle. During a calendar year different sources come in the field of view of the detector at varying zenith angles. Because of varying thickness of the overlaying atmosphere, the shower count rate is extremely dependent on zenith angle which have been carefully analyzed over time. High energy gamma-ray sources from EGRET third catalogue where observed and the data were analyzed using an excess method. A number of upper limits for a number of EGRET sources were obtained, including 6 AGNs or probably AGNs and 4 unidentified sources.

In this model a collimated ultra-relativistic ejecta collides with an amorphous dense cloud surrounding the central engine, producing gamma-rays via synchrotron process. The ejecta is taken as a standard candle, while assuming a gaussian distribution in thickness and density of the surrounding cloud. Due to the cloud high density, the synchrotron emission would be an instantaneous phenomenon (fast cooling synchrotron radiation), so a GRB duration corresponds to the time that the ejecta takes to pass through the cloud. Fitting the model with the observed bimodal distribution of GRBs' durations, the ejecta's initial Lorentz factor, and its initial opening angle are obtained as $\Gamma_{0}\lesssim 10^{3}$, and $\zeta_{0} \approx 10^{-2}$, and the mean density and mean thickness of the surrounding cloud as $\bar{n} \sim 3 \times 10^{17} cm^{-3}$ and $\bar{L} \sim 2 \times 10^{13} cm$. The clouds maybe interpreted as the extremely amorphous envelops of Thorne-Zytkow objects. In this model the two classes of long and short duration GRBs are explained in a unique frame.

The effective compton cross section in a non-degenerate plasma($n\ll\{\frac{(kT/c)^{2}+2mkT}{h^{2}}\}^{^{3/2}}$) is investigated in a wide range of temperatures. The results show a decreasing behavior with temperature especially for $kT\gg m_{e}c^{2}$. The $may be important in phenomena like accretion discs or ultra-relativistic blast waves in GRB models, where the emitted radiation has to pass through a medium containing high energy electrons.

We present quantitative statistical evidence for a $\gamma$-ray emission halo surrounding the Galaxy. Maps of the emission are derived. EGRET data were analyzed in a wavelet-based non-parametric hypothesis testing framework, using a model of expected diffuse (Galactic + isotropic) emission as a null hypothesis. The results show a statistically significant large scale halo surrounding the center of the Milky Way as seen from Earth. The halo flux at high latitudes is somewhat smaller than the isotropic gamma-ray flux at the same energy, though of the same order (O(10^(-7)--10^(-6)) ph/cm^2/s/sr above 1 GeV).

The region near the Galactic center may have experienced recurrent episodes of injection of energy in excess of $\sim$ 10$^{55}$ ergs due to repeated starbursts involving more than $\sim$ 10$^4$ supernovae. This hypothesis can be tested by measurements of $\gamma$-ray lines produced by the decay of radioactive isotopes and positron annihilation, or by searches for pulsars produced during starbursts. Recent OSSE observations of 511 keV emission extending above the Galactic center led to the suggestion of a starburst driven fountain from the Galactic center. We present EGRET observations that might support this picture.