About the computation of finite temperature ensemble averages of hybrid quantum-classical systems with molecular dynamics