The study of primary mass composition is important in revealing the vital information about the origin, acceleration and propagation mechanisms of high energy cosmic rays as it is difficult to analyse those properties directly. Several parameters may be considered for analysis of mass composition which carry information regarding the identity of the primary particle. The mass composition can be inferred from the study of depth of shower maximum ( X m a x ), shower maximum fluctuations ( Δ X m a x ) and elongation rate ( d X m a x / d E ). A new statistical analysis using higher moments of simulated air shower longitudinal profile using CORSIKA code within energy range 10 15 − 10 19 eV for different primary induced shower is presented here.