diff --git a/CMakeLists.txt b/CMakeLists.txt
index 7448b43..0d67399 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -21,7 +21,7 @@ set(CMAKE_BUILD_TYPE Release)
 
 ######### compiler flags ###########
 set(CMAKE_CXX_EXTENSIONS OFF)
-set(CMAKE_CXX_STANDARD 14)
+set(CMAKE_CXX_STANDARD 17)
 set(CMAKE_CXX_STANDARD_REQUIRED ON)
 if ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "GNU" OR
     "${CMAKE_CXX_COMPILER_ID}" STREQUAL "Clang")
@@ -87,7 +87,7 @@ set_target_properties(catima_static
                         POSITION_INDEPENDENT_CODE ON
                         ARCHIVE_OUTPUT_DIRECTORY ${PROJECT_BINARY_DIR}/lib
                     )
-target_link_libraries(catima ${EXTRA_LIBS} ${GSL_LIBRARIES} nurex::nurex)
+target_link_libraries(catima ${EXTRA_LIBS} ${GSL_LIBRARIES})
 target_link_libraries(catima_static ${EXTRA_LIBS} ${GSL_LIBRARIES})
 
 target_include_directories(catima 
diff --git a/reactions.cpp b/reactions.cpp
index bf1f478..f690af6 100644
--- a/reactions.cpp
+++ b/reactions.cpp
@@ -1,26 +1,38 @@
 #include "catima/reactions.h"
+
 #ifdef NUREX
+#include "nurex/Parametrization.h"
 #include "catima/catima.h"
 #include "catima/abundance_database.h"
+#include "catima/storage.h"
 #include <cmath>
 #include <iostream>
 namespace catima{
     
 double reaction_rate1(Projectile &projectile, const Material &target, const Config &c){
-
-    int num_elements = target.ncomponents();
     int ap = lround(projectile.A);
     int zp = lround(projectile.Z);
-    nurex::Nucleus nurex_projectile = nurex::get_default_nucleus(ap,zp);
-    
     int zt = target.get_element(0).Z;
-    int at = abundance::get_isotope_a(zt,0);
-    nurex::Nucleus nurex_target = nurex::get_default_nucleus(at,zt);
+    int at = abundance::get_isotope_a(zt,0); // most abundand natural isotope mass
 
-    double eout = energy_out(projectile,projectile.T, target,c);
-    std::cout<<eout<<"\n";
-    nurex::GlauberModelOLA_ZeroRange gm(nurex_projectile, nurex_target);
-    double cs = nurex::SigmaR(gm, projectile.T);
+    auto data = _storage.Get(projectile,target,c);
+    Interpolator range_spline(energy_table.values,data.range.data(),energy_table.num);
+
+    auto sigma_r = [&](double th){
+        double stn_sum=0.0, sum=0.0;
+        double e = energy_out(projectile.T, th, range_spline);
+        for(unsigned int i = 0;i<target.ncomponents();i++){
+            stn_sum += target.molar_fraction(i);
+            sum += target.molar_fraction(i)*nurex::SigmaR_Kox(ap,zp,e,at,zt); 
+        }
+        return sum/stn_sum;
+    };
+
+    //nurex::Nucleus nurex_projectile = nurex::get_default_nucleus(ap,zp);
+    //nurex::Nucleus nurex_target = nurex::get_default_nucleus(at,zt);
+    //nurex::GlauberModelOLA_ZeroRange gm(nurex_projectile, nurex_target);
+    //double cs = nurex::SigmaR_Kox(ap,zp,projectile.T,);
+    double cs = catima::integrator.integrate(sigma_r,0,target.thickness());
 
     double rr = reaction_rate(cs,target.number_density_cm2(0));
     std::cout<<rr<<"\n";
diff --git a/reactions.h b/reactions.h
index 4665c49..44bb6fe 100644
--- a/reactions.h
+++ b/reactions.h
@@ -18,7 +18,6 @@
 #define REACTIONS_H
 #include "catima/build_config.h"
 #ifdef NUREX
-#include "nurex/nurex.h"
 #include "catima/structures.h"
 #include "catima/config.h"
 #include "catima/integrator.h"