import tensorflow as tf

    with self.test_session():
      qdist = distributions.QuantizedDistribution(
          base_dist_cls=distributions.Normal,
          lower_cutoff=0.,
          upper_cutoff=None,
          mu=tf.zeros(
              batch_shape, dtype=tf.float32),
          sigma=tf.ones(
              batch_shape, dtype=tf.float32))

      samps = qdist.sample_n(n=5000, seed=42)
      samps_v = samps.eval()

      # With lower_cutoff = 0, the interval j=0 is (-infty, 0], which holds 1/2
