Abstract. This work proposes a parametric equilibrium distribution function F eq to be applied to the gyrokinetic studies of the Finite Orbit Width behavior of guiding centers representing several species encountered in axisymmetric tokamak plasma, as fusion products, thermal bulk and energetic particles from Ion Cyclotron Radiation Heating and Negative Neutral Beam Injections.After the analysis of the basic results of orbit theory obtained with a particularly convenient orbit coordinates set, it is shown how the proposed F eq satisfies the two conditions that make it an equilibrium distribution function: (i) it must depend only on the constants of motion and adiabatic invariants, and (ii) the guiding centers must remain confined for suitably long time.Furthermore, the F eq can be modeled, with a proper choice of its parameters, to reproduce the most common distribution functions. A local Maxwellian distribution function is obtained for the thermal plasma in the Zero Orbit Width approximation. For the fusion α particles, F eq can also reproduce the Slowing Down (SD) distribution function. More generally, for supra-thermal particles, when external heatings are present, such as (N)NBI and ICRH, the proposed model distribution function shows similarities with the anisotropic SD and the biMaxwellian distribution functions.F eq can be used to fit experimental profiles and it could provide a useful tool for experimental and numerical data analysis. Moreover, it could help 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 can be used to simulate plasma also in the presence of external heating sources.