From a5836d0bb07a24fd945d81dcca85be33ff220753 Mon Sep 17 00:00:00 2001
From: jovoy <50204158+jovoy@users.noreply.github.com>
Date: Wed, 21 Feb 2024 14:40:18 +0100
Subject: [PATCH] R2 and R4 in cone shaped optics

---
 inc/TRestAxionWolterOptics.h | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)

diff --git a/inc/TRestAxionWolterOptics.h b/inc/TRestAxionWolterOptics.h
index eaf00243..e4af9987 100644
--- a/inc/TRestAxionWolterOptics.h
+++ b/inc/TRestAxionWolterOptics.h
@@ -35,9 +35,15 @@ class TRestAxionWolterOptics : public TRestAxionOptics {
     /// Entrance radius R1 in mm. See schematic figure.
     std::vector<Double_t> fR1;  //!
 
+    /// Radius R2 in mm. See schematic figure.
+    std::vector<Double_t> fR2;  //!
+
     /// Radius R3 in mm. See schematic figure.
     std::vector<Double_t> fR3;  //!
 
+    /// Radius R4 in mm. See schematic figure.
+    std::vector<Double_t> fR4;  //!
+
     /// Radius R5 in mm. See schematic figure.
     std::vector<Double_t> fR5;  //!
 
@@ -86,12 +92,24 @@ class TRestAxionWolterOptics : public TRestAxionOptics {
         return r;
     }
 
+    /// It returns a vector with the values of R2
+    std::vector<Double_t> GetR2() {
+        std::vector<Double_t> r = TRestTools::GetColumnFromTable(fOpticsData, 1);
+        return r;
+    }
+
     /// It returns a vector with the values of R3
     std::vector<Double_t> GetR3() {
         std::vector<Double_t> r = TRestTools::GetColumnFromTable(fOpticsData, 2);
         return r;
     }
 
+    /// It returns a vector with the values of R4
+    std::vector<Double_t> GetR4() {
+        std::vector<Double_t> r = TRestTools::GetColumnFromTable(fOpticsData, 3);
+        return r;
+    }
+
     /// It returns a vector with the values of R5
     std::vector<Double_t> GetR5() {
         std::vector<Double_t> r = TRestTools::GetColumnFromTable(fOpticsData, 4);