I have some data and a model, e.g., 

$H_0: \xi \sim \mathcal N(\mu_1 > \mu_0, \sigma_0)$, 

in other words, this is the similar to the situation with "unequal means, alternative=greater".

How can I calculate the likelihood of the data under this model? Numerical integration could help, but often this integral is divergent (upper limit is infinity). I can put some prior on possible $\mu_1$ values, but it is another problem.