# Shared model: latency is log-normal; a trial ends at the cut-off (60 s), so longer latencies are recorded as 60 s
# and marked as censored (event = 0). Control median 25 s; treatment multiplies latency by RATIO; log-SD 0.6.
import numpy as np
CUTOFF = 60.0
def simulate(rng, n, median, sdlog=0.6):
    t = np.exp(np.log(median) + rng.normal(0, sdlog, n))
    return np.minimum(t, CUTOFF), (t < CUTOFF).astype(int)