Simulating normalising constants with referenced thermodynamic integration: application to COVID-19 model selection