From 012140a8f03580d16d34dfdf085a878776dba2a1 Mon Sep 17 00:00:00 2001
From: tobre1 <t.reiter1@live.de>
Date: Mon, 23 Oct 2023 09:35:56 +0200
Subject: [PATCH] Update coned cosine reflection in models

---
 include/Models/PlasmaDamage.hpp | 3 +--
 include/Models/SF6O2Etching.hpp | 2 +-
 2 files changed, 2 insertions(+), 3 deletions(-)

diff --git a/include/Models/PlasmaDamage.hpp b/include/Models/PlasmaDamage.hpp
index 7b051556..930be350 100644
--- a/include/Models/PlasmaDamage.hpp
+++ b/include/Models/PlasmaDamage.hpp
@@ -60,8 +60,7 @@ class DamageIon : public csParticle<DamageIon<T, D>, T> {
     if (NewEnergy > minEnergy) {
       reflect = true;
       auto direction = rayReflectionConedCosine<T, D>(
-          rayInternal::PI / 2. - std::min(incAngle, minAngle), rayDir,
-          geomNormal, Rng);
+          rayDir, geomNormal, Rng, std::min(incAngle, minAngle));
       E = NewEnergy;
       return std::pair<T, rayTriple<T>>{impactEnergy, direction};
     } else {
diff --git a/include/Models/SF6O2Etching.hpp b/include/Models/SF6O2Etching.hpp
index a56ffe78..9affc25c 100644
--- a/include/Models/SF6O2Etching.hpp
+++ b/include/Models/SF6O2Etching.hpp
@@ -228,7 +228,7 @@ class SF6O2Ion : public rayParticle<SF6O2Ion<NumericType, D>, NumericType> {
       E = NewEnergy;
 
       auto direction = rayReflectionConedCosine<NumericType, D>(
-          halfPI - std::min(incAngle, minAngle), rayDir, geomNormal, Rng);
+          rayDir, geomNormal, Rng, std::min(incAngle, minAngle));
 
       return std::pair<NumericType, rayTriple<NumericType>>{0., direction};
     } else {