Collisionless damping of geodesic acoustic mode (GAM) excited in the large safety factor (q) region of a tokamak plasma is investigated taking into account the effects of finite ion Larmor radius and guiding-center drift orbit width as well as parallel electric field contributions. A corresponding analytical expression for the damping rate including higher-order harmonics of ion transit resonances is systematically derived and agrees well with numerical results in its validity regime.