Generating random draws from the Truncated Normal is probably slower than fr9m a vanilla Normal.
I think we are just using scipy trunc_norm. You can see if taking rvs eith similar parameters as in your model is indeed so mich slower.
Can you share a small reproducible snippet that illustrates the slowdown?