model {
   # osservazoni con segnale + background 
   X ~ dpois(lambda)       # nr. eventi osservati
   lambda <- ls + lB       # lambda totale
   ls <- r * T             # lambda solo segnale
                           # intensita' del processo di Poisson segnale
   r ~ dgamma(1, 0.00001)  # prior coniugata gamma, esattamente dexp(0.00001)
                           # (ovvero sostanzialmente piatta)
   lB <- rB * T            # lambda background

   # osservazioni solo background
   lB0 <- rB * TB
   XB ~ dpois(lB0)
   rB ~ dgamma(1, 0.00001) # anche sul backgrond mettiamo prior coniugata vaga
}