model{
  # for each trigger detected
  for(i in 1:N) {
    part[i] ~ dbern(p.mu)
    part.ind[i] <- part[i] + 1
    T[i] ~ dbern(p.trig[part.ind[i]])
  }
  # where
  p.mu ~ dbeta(beam.r, beam.s)
  next.part ~ dbern(p.mu)
  next.part.ind <- next.part + 1
  next.T ~ dbern(p.trig[next.part.ind])
}