From 1f2a5f203915b90794a6aa54fb3f577152706c20 Mon Sep 17 00:00:00 2001
From: juacrumar <juacrumar@lairen.eu>
Date: Fri, 8 Dec 2023 15:17:35 +0100
Subject: [PATCH] implement fix by M. Bonvini

---
 new_grids/run_them_all.sh | 13 +++++++++----
 src/integration.h         |  2 ++
 2 files changed, 11 insertions(+), 4 deletions(-)

diff --git a/new_grids/run_them_all.sh b/new_grids/run_them_all.sh
index e54640c..dbab0e4 100755
--- 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
index bc222fb..b31453b 100644
--- 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.; }