From d688e787874f35c899bd965109a02b2921f3b300 Mon Sep 17 00:00:00 2001
From: jovoy <50204158+jovoy@users.noreply.github.com>
Date: Thu, 22 Feb 2024 16:51:18 +0100
Subject: [PATCH] Change alpha to beta in hyperbolic fuctions

---
 src/TRestAxionTrueWolterOptics.cxx | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/TRestAxionTrueWolterOptics.cxx b/src/TRestAxionTrueWolterOptics.cxx
index 231d998f..97ba9319 100644
--- a/src/TRestAxionTrueWolterOptics.cxx
+++ b/src/TRestAxionTrueWolterOptics.cxx
@@ -332,10 +332,11 @@ Int_t TRestAxionTrueWolterOptics::SecondMirrorReflection(const TVector3& pos, co
 
     TVector3 vertex(0, 0, fBackVertex[mirror]);
     Double_t focal = fR3[mirror] / TMath::Tan(4 * fAlpha[mirror]);
+    Double_t beta = 3 * fAlpha[mirror];
 
     //// Reflection on second mirror
     fSecondInteractionPosition = REST_Physics::GetHyperbolicVectorIntersection(
-        pos, dir, fAlpha[mirror], fR3[mirror], focal);  // should add this: TVector3(0, 0, -1), vertex,
+        pos, dir, beta, fR3[mirror], focal);  // should add this: TVector3(0, 0, -1), vertex,
 
     if (fSecondInteractionPosition.Z() > GetExitPositionZ() ||
         fSecondInteractionPosition.Z() < (0.5 * fXSep[mirror])) {
@@ -347,7 +348,7 @@ Int_t TRestAxionTrueWolterOptics::SecondMirrorReflection(const TVector3& pos, co
     }
 
     TVector3 hyperNormal =
-        REST_Physics::GetHyperbolicNormal(fSecondInteractionPosition, fAlpha[mirror], fR3[mirror], focal);
+        REST_Physics::GetHyperbolicNormal(fSecondInteractionPosition, beta, fR3[mirror], focal);
 
     fExitDirection = GetVectorReflection(fMiddleDirection, hyperNormal);