The SCENIC code package has been developed to integrate self-consistently an anisotropic pressure magnetohydrodynamic equilibrium state with power absorption from ion cyclotron resonance heating and with a guiding center particle distribution function for the energetic particles generated in three-dimensional geometry. The main novelty constitutes the inclusion of the background equilibrium state in the iterative procedure, an approach that has not been previously addressed. Applications to tokamaks and stellarators demonstrate viability of the model considered.