Determination of mass composition of high energy cosmic rays is one of the greatest challenge in modern astrophysics. All of previous methods for finding the mass composition of primary cosmic rays in a surface array at least require two independent measurements (e.g. muon and electron components) of extensive air showers (EAS). Here a new statistical parameter is introduced which can be used to determine the mass composition of vertical downward cosmic rays in a simple surface array. The most interesting thing about the new parameter is that it does not need two independent measurements and can be used in a simple surface array which does not have muon detectors.

A large distance between true and reconstructed core locations of an extensive air shower (EAS) may results in great systematic mis-estimation of EAS parameters. The reconstruction of those EASs whose core locations are outside the boundary of a surface array (outside EAS (OEAS)) results in a large distance of the reconstructed core location from the true one, especially when the true core is far outside the array. Although it may not be mentioned, the identification of OEASs is a necessary and important step in the reconstruction procedure of an EAS. In this paper, an existing technique is optimized for the identification of OEASs. The simultaneous use of this technique and a recently developed approach for reconstructing the core location of an EAS can significantly increase the sensitive area of a surface array.

A new accurate method for reconstructing the arrival direction of an extensive air shower (EAS) is described. Compared to existing methods, it is not subject to minimization of a function and, therefore, is fast and stable. This method also does not need to know detailed curvature or thickness structure of an EAS. It can have angular resolution of about 1 degree for a typical surface array in central regions. Also, it has better angular resolution than other methods in the marginal area of arrays.