Rotational transitions of HF that are important as a wavenumber standard have been analyzed by simultaneous fitting of the reported rotational and vibrational-rotational transitions plus rotational measurements of the present study with the non-Born-Oppenheimer effective Hamiltonian expressed with the optimal parameters, i.e., the determinable clusters of the expansion coefficients of the correction functions for the breakdown of the Born-Oppenheimer approximation. Since the spectral fit is made with the analytical solution of the Schrödinger equation, one can easily reproduce the fitting procedure. This is the first example of the non-Born-Oppenheimer analysis of a single isotopologue with physically meaningful parameters based on the traditional concept of the molecular constants. Fitting of a data set for HF has generated 15 independent parameter values that include 4 determinable clusters of the expansion coefficients. These 15 parameters are sufficient to generate 41 Yij coefficients and 72 rotational constants, which provide comprehensive sets of molecular constants of HF.