Abstract Hexagonal assemblies are used in a high conversion light water reactor (HCLVVR) and a nuclear reactor for cancer therapy. A cross section generation code for hexagonal assemblies is needed to design these cores precisely. For this purpose, a two-dimensional spectrum calculation code HARPY which can treat a hexagonal assembly has been developed. The HARPY code consists of a nodal code for pin-cell coupling in hexagonal geometry and a hexagonal mesh SN transport code and calculates detailed neutron spectra of fuel cells in the assembly. Verification calculations have been performed on a small hexagonal assembly by using SRAC and a Monte Carlo code GMVP to check the calculational accuracy of HARPY. The comparisons of k ∞'s, rod worths and pin power distributions have shown good agreement. This shows the adequacy of HARPY calculation.