Skip to content
Snippets Groups Projects

Update error_bands_pumplin.f

Merged Francesco Giuli requested to merge patch-4 into master
1 file
+ 6
2
Compare changes
  • Side-by-side
  • Inline
+ 6
2
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
Loading