To determine this transfer function, I will use the fast analytical circuits techniques or FACTs, as described in my last book on the subject. The principle consists of identifying the resistance \$R\$ driving the energy-storing elements - the \$C\$ and the \$L\$ - to form the time constants of the circuit defined as \$\tau=RC\$ and \$\tau=\frac{L}{R}\$. With these time constants obtained without writing a single line of algebra, it becomes possible to determine the transfer function in a low-entropy form, with well-ordered coefficients.

The difficulty then comes from the factorization of the denominator is one wants to cascade a 2nd-order expression followed by a pole, as in this 3rd-order circuit. You need to compare the times constants and, depending on their values, group certain of them to write the denominator. This document from CoPEC describes the process, it starts with the low-\$Q\$ approximation, page 51.
I used Mathcad for this exercise and I have obtained an expression which is well factored but the leading term does not match the gain at resonance yet:

However, you can see the good matching between the reference brute-force expression and my expression featuring a 2nd-order polynomial followed by an inverted pole.