  While the Quasi-Monte Carlo method of numerical integration achieves smaller
integration error than standard Monte Carlo, its use in particle physics
phenomenology has been hindered by the abscence of a reliable way to estimate
that error. The standard Monte Carlo error estimator relies on the assumption
that the points are generated independently of each other and, therefore, fails
to account for the error improvement advertised by the Quasi-Monte Carlo
method. We advocate the construction of an estimator of stochastic nature,
based on the ensemble of pointsets with a particular discrepancy value. We
investigate the consequences of this choice and give some first empirical
results on the suggested estimators.
