Error calculating Large error uncertainties (read_data.f)
TotalError is computed before the symmetrization, by adding all sources in quadrature.
To fix the issue is to re-set total error by the sum StatError, StatErrorConst, UncorError, UncorConstError in quadrature:
...
Call SetUncorErrors(npoints, StatError,
$ StatErrorConst,UncorError,UncorConstError)
.......
TotalError = StatError^2 + StatErrorConst^2 + UncorError^2 + UncorConstError^2
add to that BETA after symmetrisation:
....
BETA(CompressIdx(i),npoints) =
$ 0.5*( BetaAsym(CompressIdx(i),1,npoints)-
$ BetaAsym(CompressIdx(i),2,npoints))
TotalError = TotalError + BETA(CompressIdx(i),npoints) **2
and re-set the value in the array after the loop over syst. sources:
.....
enddo
JSET(npoints) = NDATASETS ! IndexDataset
E_TOT(npoints) = sqrt(TotalError)
GPlotVarCol(NDATASETS) = PlotVarColumn
GNPlots(NDATASETS) = PlotN
PreviousPlots = 0