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