From c50b855c732855b9c1ea77775f30b66bf355d64e Mon Sep 17 00:00:00 2001 From: David Grote Date: Wed, 8 Jan 2025 14:58:57 -0800 Subject: [PATCH] Fix the energy scaling of the generated photons --- .../BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H b/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H index ac9c213b104..77e2bb42b1e 100644 --- a/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H +++ b/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H @@ -168,10 +168,11 @@ public: auto& uz1 = soa_products_data[0].m_rdata[PIdx::uz][product_index]; // Normalize out the electron velocity and multiply by the photon momentum, E/c + // Also, the photon momentum is normalized by m_e auto u1 = std::sqrt(ux1*ux1 + uy1*uy1 + uz1*uz1); - ux1 *= (p_product_data[i]/PhysConst::c)/u1; - uy1 *= (p_product_data[i]/PhysConst::c)/u1; - uz1 *= (p_product_data[i]/PhysConst::c)/u1; + ux1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e; + uy1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e; + uz1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e; } });