Turing model that records objective function values

Hi, this might be a bit of a silly question but I’m using Turing to run a Bayesian MCMC algorithm. I’m using a distribution I made up as the distribution of my response variable and I would like the chain object to store the values of the PDF from this distribution. The reason for this is that the PDF takes a long time to compute so I don’t like having to compute it after the algorithm is done, since I’m effectively doing the same computation twice. Any thoughts for a workaround on this? Right now, this is what I have:

@model function turing_model(x)
    α ~ Uniform(-5.0, 5.0)
    β ~ Uniform(-5.0, 5.0)
    γ ~ Uniform(-5.0, 5.0)
    q ~ Uniform(-5.0, 5.0)

    x ~ CustomDist(α, β, γ, q)

Somewhere in there, is there a place I could record the likelihood function values for each value of the parameters? I know how to access the parameter values of each iteration from the output of the sample function, but it doesn’t seem to keep track of the likelihoods. Thoughts?