From 9484fdf75c9e37f62f67e8b874d62b59a68246bd Mon Sep 17 00:00:00 2001 From: Francesco Giuli <francesco.giuli@cern.ch> Date: Fri, 9 Oct 2020 15:56:50 +0000 Subject: [PATCH] Update error_bands_pumplin.f so the bug in plotting the error bands with increased tolerance is fixed --- src/error_bands_pumplin.f | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/src/error_bands_pumplin.f b/src/error_bands_pumplin.f index f2d1fca22..a2e5b7d4e 100644 --- a/src/error_bands_pumplin.f +++ b/src/error_bands_pumplin.f @@ -1,4 +1,3 @@ - subroutine Error_Bands_Pumplin implicit none @@ -56,6 +55,10 @@ C SG: x-dependent fs: double precision GetUmat !function ! double precision DecorVarShift +C for uncertainty bands with increased tolerance + double precision fmin, fedm, errdef + integer npari, nparx, istat + C for theory errors: double precision, allocatable :: TheoVars(:,:,:) @@ -141,7 +144,8 @@ C if(doOffset) then shift = shift_dir * DecorVarShift(iint, j) else - shift = shift_dir * GetUmat(iint,j) + call MNSTAT(fmin, fedm, errdef, npari, nparx, istat) !> MW&FG for scaling with DeltaChi2>1.0 + shift = shift_dir * GetUmat(iint,j)*SQRT(errdef) !> MW&FG for scaling with DeltaChi2>1.0 endif a(i) = a(i) + shift endif -- GitLab