This work proposes a parametric equilibrium distribution function to be applied to the gyrokinetic studies of the finite orbit width behavior of guiding centers representing several species encountered in axisymmetric tokamak plasmas, such as fusion products, thermal bulk and energetic particles from ion cyclotron radiation heating and negative neutral beam injections. can be used to fit experimental profiles and it could provide a useful tool for experimental and numerical data analysis. Moreover, it could help one to develop analytical computations for facilitating data interpretation in the light of theoretical models. This distribution function can be easily implemented in gyrokinetic codes, where it could be used to simulate plasma also in the presence of external heating sources.