We prove the L p boundedness of the circular maximal function on the Heisenberg group H 1 for 2 < p ≤ ∞. The proof is based on the square sum estimate associated with the 2the phase space arising from the vector fields X 1 , X 2 , tX 3 , ∂/∂t on the Heisenberg group, rather than the 2 × 1 cone |(ξ 1 , ξ 2 )| = |ξ 3 | of the frequency space arising from ∂/∂x 1 , ∂/∂x 2 , ∂/∂t on the Euclidean space. Contents 1. Introduction 2. Basic Decompositions 3. Statement of Main Estimates 4. Idea of Proof 5. Torus According to Four Vector Fields ∂ ∂t , tX 3 , X 1 , X 2 6. Square Sum over Radial Decompositions 7. Square Sum over Cubes C k,ℓ,θ 8. Vector-Valued Estimates 9. Littlewood-Paley Inequality Associated with Cubes 10. Appendix References