How Hamiltonian Monte Carlo handle the discontinuity at zero in hurdle models?

You need to specify a Stan model that marginalizes over the latent binary indicator. Here’s a post with example code for the zero-inflated Poisson case

Here’s @mbjoseph 's super helpful tutorial on marginalization: