The  ¹⁶⁰Gd(p,n)¹⁶⁰Tb excitation function was measured between 4-18 MeV using stacked-target activation at Lawrence Berkeley National Laboratory's 88-Inch Cyclotron. Nine copper and eight titanium foils served as proton fluence monitor foils, using the  ^natCu(p,x)⁶⁵Zn,  ^natTi(p,x)⁴⁸V, and  ^natTi(p,x)⁴⁶Sc monitor standards, respectively. Variance minimization using an MCNP v.6.2 model reduced the systematic uncertainties in proton energy and fluence. A priori predictions of the  ¹⁶⁰Gd(p,n) reaction using ALICE, CoH, EMPIRE, and TALYS, as well as the TENDL database, are compared to the experimentally measured values.