implement fix by M. Bonvini

NNPDF · Dec 8, 2023 · 1f2a5f2 · 1f2a5f2
1 parent ad5c9aa
commit 1f2a5f2
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 4 deletions.
diff --git a/new_grids/run_them_all.sh b/new_grids/run_them_all.sh
@@ -9,19 +9,24 @@ exe=../build/Vrap
 kinfolder=input_kinematics
 output_placeholder=test.pineappl.lz4
 
+pineapploptimize() {
+    pineappl write ${output_placeholder} ${1}
+}
+
+
 ${exe} inputE605nlo.dat input_kinematics/E605.dat
-pineappl optimize ${output_placeholder} E605nlo.pineappl.lz4
+pineapploptimize  E605nlo.pineappl.lz4
 
 ## E886P
 ${exe} inputE866nlo.dat  input_kinematics/E866P.dat
-pineappl optimize ${output_placeholder} E866nlo.pineappl.lz4
+pineapploptimize  E866nlo.pineappl.lz4
 
 ## E886R
 ${exe} inputE866nlo.dat  input_kinematics/E866R.dat
-pineappl optimize ${output_placeholder} E866Rnlo.pineappl.lz4
+pineapploptimize  E866Rnlo.pineappl.lz4
 
 ${exe} inputE866deutnlo.dat  input_kinematics/E866R.dat
-pineappl optimize ${output_placeholder} E866deutRnlo.pineappl.lz4
+pineapploptimize E866deutRnlo.pineappl.lz4
 
 ## E906
 ./run_E906.sh inputE906nlo.dat

diff --git a/src/integration.h b/src/integration.h
@@ -258,6 +258,8 @@ double int_NNLO(double y, double ys, double z){
 
     if ((ys < 1.0e-08) || (ys > 1.0-1.0e-08) || (z > 1.0-1.0e-08)){ return 0.; }
     double tau = Q*Q/E_CM/E_CM;
+    if ((tau * exp(2.*y) >= 1.) || (tau * exp(-2.*y) >= 1.)){ return 0.; } 
+
     if ((tau * exp(2.*y) >= 1.) || (tau * exp(-2.*y) >= 1.) 
                                 || (z <= tau)){ return 0.; }