© 2021 Korean Nuclear SocietyIn this study, 3″ × 3″ NaI(Tl) detector, which is widely used in gamma spectroscopy, was modeled with FLUKA code, and calculations required to determine the detector's energy resolution were reported. Photon beams with isotropic distribution with 59, 81, 302, 356, 511, 662, 835, 1173, 1275, and 1332 keV energy were used as radiation sources. The photon pulse height distribution of the NaI(Tl) without influence of its energy resolution obtained with FLUKA code has been converted into a real NaI(Tl) response function, using the necessary conversion process. The photon pulse height distribution simulated in the conversion process was analyzed using the ROOT data analysis framework. The statistical errors of the simulated data were found in the range of 0.2–1.1%. When the results, obtained with FLUKA and ROOT, are compared with the literature data, it is seen that the results are in good agreement with them. Thus, the applicability of this procedure has been demonstrated for the other energy values mentioned.