We compute analytical formulae for pricing arithmetic Asian options under jump diffusion CIR processes. To derive the solution, we employ a characteristic function of the underlying asset price process and its integrated process that is not required to take the inversion Fourier or Laplace transform. We conduct numerical tests for validation of proposed formulae to confirm that they provide stable and accurate option prices with much faster computation time than the full Monte Carlo method.