I am trying to write a function using Rcpp to quickly draw samples from a half t distribution. I need to draw n samples, repeated number of times equal to number of Particles. For each n, the arguments controlling the distribution might be different, hence df and sigma being supplied as vectors to the function. The current Rcpp function is below: