LC in parallel with R circuit resonance
The equivalent impedance you calculated is only valid if you measure it across \$R\$, disregarding the source, but your output is across \$C\$. To calculate the transfer function first you must account for the impedance (let's consider it a resistance) of the voltage source, which can be modelled as a series resistance. It doesn't make sense otherwise, since a voltage source is considered zero, or very low impedance, thus a resistor across it which is larger that its source will not contribute with anything.
Consider the voltage across \$R\$ to be \$V_1\$. Then, the impedance seen from the source's side is \$Z_1\$:
$$\begin{align*}
Z_{LC}&amp;=sL+\dfrac{1}{sC}\tag{1} \\
Z_1&amp;=\dfrac{1}{\dfrac{1}{Z_{LC}}+\dfrac{1}{R}} \\
&amp;=R\dfrac{s^2+\dfrac{1}{LC}}{s^2+\dfrac{R}{L}s+\dfrac{1}{LC}}\tag{2} \\
&amp;\Rightarrow \\
V_{out}&amp;=V_1\dfrac{Z_1}{Z_1+R_{in}}\tag{3} \\
&amp;\Rightarrow \\
H(s)&amp;=\dfrac{R}{R_{in}+R}\dfrac{s^2+\dfrac{1}{LC}}{s^2+\dfrac{R_{in}}{R_{in}+R}\dfrac{R}{L}s+\dfrac{1}{LC}}\tag{4}
\end{align*}$$
\$(4)\$ shows an attenuated and damped 2nd order lowpass, which should have a response something like this (for some random values of \$L=1\;\text{mH},\;C=1\;\mu\text{F},\;R=100\;\Omega,\;R_{in}=50\;\Omega\$):

But your plots show only a peak resembling a bandpass. So, in the event that your connections are according to your schematic, I would guess that you're only measuring a finite bandwidth around the peak, something like this:

If this is the case, then you need to increase the range, otherwise it's possible you have something different, and it can only be determined by showing the actual values for the elements (including parasitics, for L, at least), source's output impedance, point of measuring, how it's measured, range of frequencies and amplitudes/magnitudes you're measuring.

Given the recent edit, you still haven't said what is the range of frequencies you're mesuring, but here is the reworked response of your circuit with your values:

This further reinforces my guess that you're measuring the output over a limited bandwidth, creating the illusion that it's a bandpass, instead of a very underdamped lowpass.
As for the damping, as your newly found formula shows, it's exactly as shown in my derivation, which means that the usual \$R/L\$ is now \$R_{eq}/L\$, where \$R_{eq}=R_{in}||R\$. In short, both \$R\$ and \$R_{in}\$ contribute. And if \$L\$ had a series resistance, that would have been in there, too. If it makes it easier for you, here's the rewritten transfer function:
$$H(s)=\dfrac{R}{R_{in}+R}\dfrac{s^2+\dfrac{1}{LC}}{s^2+\dfrac{R_{in}||R}{L}s+\dfrac{1}{LC}}$$