\noindent \textbf{Background:} Detailed spectroscopy of neutron-rich, heavy, deformed nuclei is of broad interest for nuclear astrophysics and nuclear structure. Nuclei in the r-process path and following freeze-out region impact the resulting r-process abundance distribution, and the structure of nuclei midshell in both proton and neutron number helps to understand the evolution of subshell gaps and large deformation in these nuclei. \noindent \textbf{Purpose:} To improve the understanding of the nuclear structure of $^{160}$Gd, specifically the $K^\pi$=$4^+$ bands, as well as study the $\beta$-decay of $^{160}$Eu into $^{160}$Gd. \noindent \textbf{Methods:} High-statistics decay spectroscopy of $^{160}$Gd resulting from the $\beta$-decay of $^{160}$Eu was collected using the GRIFFIN spectrometer at the TRIUMF-ISAC facility. \noindent \textbf{Results:} Two new excited states and ten new transitions were observed in $^{160}$Gd. The $\beta$-decaying half-lives of the low- and high-spin isomer in $^{160}$Eu were determined, and the low-spin state's half-life was measured to be $t_{1/2}=26.0(8)$~s, $\sim$16\% shorter than previous measurements. Lifetimes of the two $K^\pi$=$4^+$ band-heads in $^{160}$Gd were measured for the first time, as well as $\gamma$-$\gamma$ angular correlations and mixing ratios of intense transitions out of those band-heads. \noindent \textbf{Conclusions:} Lifetimes and mixing ratios suggest that the hexadecapole phonon model of the $K^\pi$=$4^+$ band-heads in $^{160}$Gd is preferred over a simple two-state strong mixing scenario, although further theoretical calculations are needed to fully understand these states. Additionally, the 1999.0 keV state in $^{160}$Gd heavily populated in $\beta$-decay is shown to have positive parity, which raises questions regarding the structure of the high-spin $\beta$-decaying state in $^{160}$Eu.
