R2 and R4 in cone shaped optics

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);