Recently, we reported a new set of Bessel type orbitals, which include the radial part of generalized Bessel functions for LCAO calculations of atomic systems. In this study, to achieve further improvement of the performance of generalized Bessel type basis sets in the Hartree-Fock-Roothaan calculations, different hyperbolic cosine functions inserted into the radial part of those generalized Bessel functions. For this purpose, three different generalized hyperbolic cosine functions have been used to construct the generalized Bessel type hyperbolic cosine basis sets. The accuracies of generalized Bessel type hyperbolic cosine orbitals within the minimal basis sets approach are compared to show their superiority to conventional approaches those in the literature. The performance of the presented basis functions is also compared to the numerical Hartree-Fock results. Our virial ratios are in good agreement to within 8-digits of the -2. It is shown that the results obtained by the new basis sets surpass the quality and accuracy of existing BTOs basis sets