Numerical integration of detector response functions via Monte Carlo simulations